On Fri, 17 Feb 2006, Barry Zajdlik wrote:
> Hello all,
>
> Thanks for the responses but I am still annoyed by this seemingly simple
> problem; I recorded sessionInfo() as below.
>
> x<-rep(0.02,10)
>> var(x)
> [1] 1.337451e-35
Well, yes, it is accurate only to 35 digits and this is a bit unfortuna
On Thu, 16 Feb 2006, Jacqueline Hall wrote:
> Hi all,
>
> Does anyone know if there is an implementation of Brookmeyer & Crowley's
> confidence interval for the median survival time in R?
> Reference : Brookmeyer & Crowley, "A confidence interval for the median
> survival time" (1982) Biometircs
On Thu, 16 Feb 2006, Arne Jol wrote:
> Dear R-Help,
>
> I wonder if there is a prediction function for a clogit model which can be
> used in the same way as the predict function for the multinom model.
>
> In prediction('multinommodel',testset ...) it is possible to predict the
> class or the cla
de.html
--
??
Deparment of Sociology
Fudan University
Thomas Lumley Assoc. Professor, Biostatistics
[EMAIL PROTECTED] University of Washington, Seattle__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman
On Wed, 15 Feb 2006, Lei Liu wrote:
> Hi,
>
> I am using coxph in simulations and I want to store the "robust se" (or
> "se2" in frailty models) for each replicate. Is there a function to retrieve
> it, like vcov() for the variance estimate? Thanks!
>
There isn't a function, so you need to extra
tal( ~vol, srv1 ) )
>
> I only obtain the total,
>
>> print( svytotal( ~vol, srv1 ) )
>total SE
> vol 2377 34.464
>
> or worse,
>
> print( svytotal( ~vol + strata, srv1 ) )
> total SE
> vol 2377.0 34.464
> strataA 42.0 0.000
>
On Mon, 13 Feb 2006, Stephen Richards wrote:
> Can survreg() handle interval-censored data like the documentation
> says? I ask because the command:
>
> survreg(Surv(start, stop, event) ~ 1, data = heart)
>
> fails with the error message
>
> Invalid survival type
>
> yet the documentation
consider the distributions given
by http://mathworld.wolfram.com/EfronsDice.html)
However, the estimate is not the difference in sample medians. It is the
median pairwise difference.
-thomas
Thomas Lumley Assoc. Professor, Biostatistics
[EMAIL PROTECTED
On Thu, 26 Jan 2006, Douglas Grove wrote:
> Hi,
>
> I'm trying to figure out if there's an automated way to get
> read.table to read in my data and *not* convert the character
> columns into anything, just leave them alone. What I'm referring
> to as 'character columns' are columns in the data th
On Thu, 26 Jan 2006, Dale Steele wrote:
> The data looks like:
> -->
> icsrvepf
> fpevrsci
> ics
> p
>
> f
> ic
> <--
>
> I would like to convert the about to a matrix of the form:
>
> i c s r v e p f
> 1 1 1 1 1 1 1 1
> 1 1 1 1 1 1 1 1
> 1 1 1 0 0 0 0 0
> 0 0 0 0 0 0 1 0
> 0 0 0 0 0 0 0 0
>
On Wed, 25 Jan 2006, singyee ling wrote:
> Dear R-users,
>
> I am sorry if this is obvious. I am testing the proportional hazard
> assumptions using cox.zph. If i am not wrong, a g(t) function must be
> assumed. Four possibilities available in R are "km","identity" and "rank".
> may i know what f
On Tue, 24 Jan 2006, Taka Matzmoto wrote:
> Hi R users
> I have a simple question to ask.
> Why am I getting FALSE on this
>
>> is.integer(10)
> [1] FALSE
>
> 10 is a integer number.
is.integer() asks how a number is stored, not whether the number happens
to be an integer. Explicit numbers type
On Wed, 25 Jan 2006, roger bos wrote:
>
> Does anyone have code that keeps generating random data until the memory is
> full and then tells you how much memory was successfully used? I could try
> writing it, but if someone has already done it, thats all the better!
>
In recent versions of R gc()
>>>> c( 1, 2, 2, 3,3,3,1 )
>>>
>>> One-liner:
>>>
>>> > table(ids)[ids]
>>> ids
>>> ID1 ID2 ID2 ID3 ID3 ID3 ID5
>>> 1 2 2 3 3 3 1
>>>
>>> 'table(ids)' computes t
i) ){
>
>i <- i + 1
>}
>
>spotrep[spot : (spot + i-1)] <- i
>spot <- spot + i
>#cat("spot : ",spot,"\n")
>}
>
>obj$genes$spotrep <- spotrep[order(ord)]
>
>obj
#cat("spot : ",spot,"\n")
>}
>
>obj$genes$spotrep <- spotrep[order(ord)]
>
>obj
>
> }
>
> __
> R-help@stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
>
ly than programming would, and have
generally been correct.
If you want Stata, you know where to find it (and it's a good choice for
many problems).
-thomas
Thomas Lumley Assoc. Professor, Biostatistics
[EMAIL PROTECTED] University of Washington, Seattle
__
On Mon, 23 Jan 2006, Linda Lei wrote:
> Thank you guys.
> But I tried the commands and I still get:
>
>> aml1<-aml[aml$group==1,]
>> aml1
> [1] time status x
> <0 rows> (or 0-length row.names)
>> esf.fit <- survfit(Surv(aml1$weeks,status) ~ 1)
> Error in Surv(aml1$weeks, status) : Time variable
On Fri, 20 Jan 2006, Linda Lei wrote:
> Hi there,
>
>
>
> I have a question about one command sentence when I follow the example
> in the book of "Survival analysis in S":
>
> > aml1<-aml[aml$group==1]
If this is really what the book says you should probably complain to the
author. However, if
alues)
>
>
> Thank you in advance for any hints.
> Best regards,
> P. Olsson
>
> [[alternative HTML version deleted]]
>
> __
> R-help@stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLE
On Wed, 18 Jan 2006, Robert Michael Rausch wrote:
> I would like to generate a contour-plot according to a master plot. The
> problem is that the rainbow-palette included in R does not answer this
> purpose. I need a darker blue, no turquoise, relatively less green, more
> yellow and more red. H
in symbols (...)
> are not behaving as expected.
>
The ... argument is also passed to .Internal, and presumably the code there
gives the warning.
-thomas
Thomas Lumley Assoc. Professor, Biostatistics
[EMAIL PROTECTED] University of Washington, Seattle
__
why lmList matches the linear fits rather than the GLM
> fits would be greatly appreciated. I'm using R2.2.1 with lme version 0.98-1
> in Windows XP.
>
> Daniel Farewell
> Cardiff University
>
> ______
> R-help@stat.math.ethz.ch mailing list
> ht
rs.
> I wonder if
>
> (i) there isn't a way of forcing each variable to be numeric or integer,
> maintaining it's original values (instead of "Stata labels" as "R
> levels"). Or,
Yes. If you read the help page for read.dta() it tells you
and create
SPSS variable and value labels.
--
Chuck Cleland, Ph.D.
NDRI, Inc.
71 West 23rd Street, 8th floor
New York, NY 10010
tel: (212) 845-4495 (Tu, Th)
tel: (732) 452-1424 (M, W, F)
fax: (917) 438-0894
__________
R-help@stat.math.ethz.ch mailing list
https
On Mon, 9 Jan 2006, Michael Dewey wrote:
>
> Further to that I feel that (perhaps because they do not like to blow their
> own trumpet too much) the authors of books on R do not stress how much most
> questioners could gain by buying and reading at least one of the many books
> on R. When I started
>
> Any suggestions?
>
> Thanks.
>
> Harsh,
>
>
>
>
>
> -
>
> [[alternative HTML version deleted]]
>
> __________
> R-help@stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listi
On Wed, 4 Jan 2006, Roger Bivand wrote:
> Could I ask for comments on:
>
> source(url("http://spatial.nhh.no/R/etc/capabilities.R";), echo=TRUE)
>
> as a reproduction of the Stata capabilities session? Both the t test and
> the chi-square from our side point up oddities. I didn't succeed on
> putti
us what that
criterion is).
These functions produce one or more best models of each size, and for
models of the same size all the commonly-used criteria reduce to ranking
by residual sum of squares, which is what leaps() and regsubsets() do.
-thomas
Thomas Lumley Assoc
On Tue, 3 Jan 2006, Verena Hoffmann wrote:
> Hello!
>
> I want to compare two Kaplan-Meier-Curves by using the Logrank-Test:
>
> logrank(Surv(time[b], status[b]) ~ group[b])
>
> This way I only get the value of the test-statistic, but not the p-value.
>
> Does anybody know how I can get the p-valu
hes, the language was substantially reworked in both
versions 8 and 9 to allow the vendor to implement better graphics and
linear mixed models.
On the question of which system really is easier to learn I can only
comment that this isn't the only question where education, as a field,
would
On Tue, 3 Jan 2006, Brian Perron wrote:
> Hello all:
>
> I am interested in computing what the multilevel modeling literature
> calls a multiple membership model. More specifically, I am working with
> a data set involving clients and providers. The clients are the
> lower-level units who ar
On Mon, 2 Jan 2006, Philippe Grosjean wrote:
>
> That said, I think one should interpret Mitchell's paper in a different
> way. Obviously, he is an unconditional and happy Stata user (he even
> wrote a book about graphs programming in Stata). His claim in favor of
> Stata (versus SAS and SPSS, and
file in
the sources.
Description: Functions for reading and writing data stored by
statistical packages such as Minitab, S, SAS, SPSS, Stata,
Systat, ..., and for reading and writing .dbf (dBase) files.
License: GPL version 2 or later
Packaged: Fri Dec 9 14:46:20 2005; ripley
B
On Thu, 15 Dec 2005, Phineas wrote:
> How many distinct values can rnorm return?
2^32-1. This is described in help(Random)
> I assume that rnorm manipulates runif in some way, runif uses the Mersenne
> Twister, which has a period of 2^19937 - 1. Given that runif returns a 64
> bit precision fl
o read the posting guide! http://www.R-project.org/posting-guide.html
Thomas Lumley Assoc. Professor, Biostatistics
[EMAIL PROTECTED] University of Washington, Seattle__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailma
On Tue, 13 Dec 2005, May, Roel wrote:
> Hi all,
>
> I am using a 'tricked' Cox Hazard regression model for discrete choice
> habitat modelling.
> However, I'm having a hard time understanding the meaning of the first
> line the following part of the summary() output:
>
> Rsquare= 0.307 (max possib
-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Thomas Lumley Assoc. Professor, Biostatistics
[EMAIL PROTECTED] University of Washington, Seattle__
R-help@stat.math.ethz.ch mailing list
On Mon, 28 Nov 2005, Hong Ooi wrote:
>
> Hi,
>
> I was just wondering if there was any available R code that could handle
> general constrained optimisation problems. At the moment I'm using
> nlminb and optim, both of which allow box constraints on the parameters,
> but ideally I'd like to be able
___
> R-help@stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
>
Thomas Lumley Assoc. Professor, Biostatistics
[EMAIL PR
On Tue, 22 Nov 2005, giovanni parrinello wrote:
> Dear All,
> I am trying to apply the function 'cox.zph' of the library survival, but I
> receive this error message:
> not found the object 'residuals.coxph'.
> I have re-installed the library 'survival' without any change and also a
> search wit
On Fri, 18 Nov 2005, ronggui wrote:
> x<-1
> f<-function(){
> x<-3
> eval(substitute(x+y,list(y=10)))
> }
> f() #13
>
> x<-1
> f<-function(){
> x<-3
> eval(substitute(x+y,list(y=10)), envir = sys.frame(sys.parent()))
> }
> f() #11
>
> x<-1
> f<-function(){
> x<-3
> eval.parent(substitute(x+y,list(
> 4700 Keele Streethttp://www.math.yorku.ca/SCS/friendly.html
> Toronto, ONT M3J 1P3 CANADA
>
> __
> R-help@stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide! http:/
ames, variables, varnames, extras,
> extranames, :
>invalid variable type
> ===
>
> - Vivien Chen -
> [[alternative HTML version deleted]]
>
> __
> R-help@stat.math.ethz.ch mailing list
> http
On Tue, 15 Nov 2005, Claus Atzenbeck wrote:
> On Tue, 15 Nov 2005, Peter Dalgaard wrote:
>
However, for Wilcoxon tests, the formula for effect sizes is:
r = Z / sqrt(N)
I wonder how I can calculate the Z-score in R for a Wilcoxon test.
>>>
>>> Does anyone know how to calculate
accelerated failure
models to right/left/interval censored data.
-thomas
Thomas Lumley Assoc. Professor, Biostatistics
[EMAIL PROTECTED] University of Washington, Seattle
__
R-help@stat.math.ethz.ch mailing list
not be made
consistent across even the basic vector types.
-thomas
> Ron
>
>>>> Thomas Lumley <[EMAIL PROTECTED]> 11/10/05 11:02 PM >>>
> On Thu, 10 Nov 2005, Ron Ophir wrote:
>
>> Thanks Thomas,
>>
>> "...For logical subscript
On Thu, 10 Nov 2005, Thomas Lumley wrote:
> On Thu, 10 Nov 2005, Ron Ophir wrote:
>
>> Thanks Thomas,
>>
>> "...For logical subscripts you could argue that the
>> ambiguity isn't present and that if the index was NA the element should
>> just b
On Thu, 10 Nov 2005, Ron Ophir wrote:
> Thanks Thomas,
>
> "...For logical subscripts you could argue that the
> ambiguity isn't present and that if the index was NA the element should
> just be set to NA. This change might be worth making."
>
> I see you got my point. NA should return NA no matte
On Thu, 10 Nov 2005, Bickel, David wrote:
> Are there any R functions or packages that can compute distributions,
> expectations, or quantiles of order statistics (or sample quantiles or
> extreme values) for a given distribution such as a normal distribution?
> Both exact and asymptotic calculati
On Thu, 10 Nov 2005, Ron Ophir wrote:
> Hi,
> I am trying to apply two different functions on on a vector as follow:
> a<-c(NA,1,2,3,-3,-4,-6)
> if a>0 I would like to raise it by the power of 2: 2^a and if the a<0 I
> would like to have the inverse value, i.e., -1/2^a.
> so I thought of doing it
On Mon, 7 Nov 2005, Jarrett Byrnes wrote:
> I'm attempting to analyze some survey data comparing multiple docks. I
> surveyed all of the slips within each dock, but as slips are nested
> within docks, getting multiple samples per slip, and don't really
> represent any meaningful gradient, slip is
On Sun, 6 Nov 2005, Kjetil Brinchmann halvorsen wrote:
> Why doesn't paste behave in apply as sum?
>
Because sum maps a vector of inputs to a single output, but paste does
not, unless you use collapse=
-thomas
__
R-help@stat.math.ethz.ch mail
> PLEASE do read the posting guide!
> http://www.R-project.org/posting-guide.html
>
> __
> R-help@stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide! http://www.R-project.org/pos
On Fri, 4 Nov 2005, Uwe Ligges wrote:
> Dr Carbon wrote:
>
>> At the risk of being beaten about the face and body, can somebody explain
>> why the middle example: log2(2^3); floor(log2(2^3)) is different than
>> examples 1 and 3?
>
>
> Because
>
> > log2(2^3) - 3
> [1] -4.440892e-16
>
This is a
> month10
> day 06
> svn rev 35749
> language R
>
> --
> Ken Kelley, Ph.D.
> Inquiry Methodology Program
> Indiana University
> 201 North Rose Avenue, Room 4004
> Bloomington, Indiana 47405
> http://www.indiana.edu/~kenkel
>
> ___
plotmath that starts
## How to combine "math" and numeric variables :
-thomas
Thomas Lumley Assoc. Professor, Biostatistics
[EMAIL PROTECTED] University of Washington, Seattle
__
R-help@stat.m
On Mon, 31 Oct 2005, Jarrett Byrnes wrote:
> Hey, all. Quick question. I'm attempting to use some of the great
> graphs generated in R for an upcoming talk that I'm writing in
> Powerpoint. Copying and pasting (I'm using OSX) yields graphs that
> look great in Powerpoint - until I resize them.
On Mon, 31 Oct 2005, Jeffrey Stratford wrote:
> Hi, I have no experience with R and I'm finding the manuals a bit obtuse
> and written as if I already understood R.
>
> I'm trying to import a csv file from a floppy and it's not working. The
> code I'm using is
>
> read.table("F:\GEORGIA\species_ri
.R-project.org/posting-guide.html
Thomas Lumley Assoc. Professor, Biostatistics
[EMAIL PROTECTED] University of Washington, Seattle__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do rea
; [4,] -0.5898438 -0.2539063 3.592773 2.0351563
> 0.7849121
> [5,] 1.500 -1.1289063 -2.070313 -0.4003906
> 2.1386719
>>
>
> __
> R-help@stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE d
On Fri, 28 Oct 2005, Richard Reiss wrote:
>
> Does anyone know of an R package that I can use to uncensor a normal or
> log-normal dataset? I'm particularly interested in the MLE method of
> Cohen (1959), "Simplified estimators for the normal distribution when
> samples are single censored or trun
mp(alphas=possible.alphas[3], betas=possible.betas[11],
> timep=lifetimes)
>
>
> ### I'd appreciate any kind of advice.
> ### Thanks a lot in advance.
> ### Roland
>
>
> +
> This mail has been sent through the MPI for Demographic Rese...{{dropp
Tel 0543-731836
Tel/Fax 0543-731612
**
[[alternative HTML version deleted]]
Thomas Lumley Assoc. Professor, Biostatistics
[EMAIL PROTECTED] University of Washington, Seattle__
R-help@s
f lines.
Lewis S, Clarke M. Forest plots: trying to see the wood and the trees.
BMJ 2001;322:1479-80."
-thomas
Thomas Lumley Assoc. Professor, Biostatistics
[EMAIL PROTECTED] University of Washington, Seattle
___
in: make.packages.html()
> Error: .onLoad failed in 'loadNamespace' for 'tcltk'
> Error: package 'tcltk' could not be loaded
>
> __
> R-help@stat.math.ethz.ch mailing list
>
On Wed, 19 Oct 2005, Tim Churches wrote:
> MJ Price, Social Medicine wrote:
>> I have been trying to find a function to calculate the Breslow-Day test for
>> homogeneity of the odds ratio in R. I know the test can be preformed in SAS
>> but i was wondering if anyone could help me to perform this i
On Fri, 14 Oct 2005, Ido M. Tamir wrote:
> Hello,
> i am trying to subset a dataframe multiple times:
> something like:
>
> stats <- by(df, list(items), ttestData)
>
> ttestData <- function(df){
>t.test( df[,c(2,3,4), df[,c(5,6,7)]
> }
>
> While this works for small data, it is to slow for my
> R-help@stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
>
Thomas Lumley Assoc. Professor, Biostatistics
[EMAIL PROTECTED] University o
On Tue, 11 Oct 2005, joerg van den hoff wrote:
> many thanks to thomas and gabor for their help. both solutions solve my
> problem perfectly.
>
> but just as an attempt to improve my understanding of the inner workings of R
> (similar problems are sure to come up ...) two more question:
>
> 1.
>
orce evaluation of the variable n such that
> the lapply call works? I suspect it has something to do with eval and
> specifying the correct evaluation frame, but how?
>
>
> many thanks
>
> joerg
>
> ______
> R-help@stat.math.e
284)
and to a ratio estimator
svyratio(~as.numeric(y1*(id2>1)), ~as.numeric(id2>1), design=dmu284)
All three give the same mean estimator and standard error.
-thomas
Thomas Lumley Assoc. Professor, Biostatistics
[EMAIL PROTECTED]
On Fri, 7 Oct 2005, Peter Muhlberger wrote:
>
> The max function won't do the trick because I need the entire matrix. I
> could do one cell at a time, but this is part of a ML routine that needs to
> be evaluated hundreds of thousands of times, so I can't afford to slow it
> down that much.
pmax,
On Fri, 7 Oct 2005, Peter Muhlberger wrote:
> Does anyone know how -log(x) can equal 743 but -log(x+0)=Inf? That's what
> the following stream of calculations suggest:
>
> Browse[2]> -log ( 1e-323+yMat2 - yMat1 * logitShape(matrix(parsList$Xs,
> nrow = numXs, ncol=numOfCurves), matrix(means, nr
On Thu, 6 Oct 2005, Pyy-Martikainen Marjo wrote:
> I specify a model with strata by covariate interactions. I would like to
> conduct a Wald test
> for the null hypothesis "no differences between any covariate effects in the 3
> groups".
The "survey" package has a function regTermTest (which isn
rr(X,Y)= ",(covxy/(varx*vary)^.5)))
>
> print("NOTE -sqrt(varx*vary)<=covxy<=sqrt(varx*vary)")
> #print(A)
> }
> BINORMAL(varx=1,vary=1,covxy=0,meanx=0,meany=0)
>
> thanx
>
> /allan
>
> ___
On Tue, 4 Oct 2005, David L. Van Brunt, Ph.D. wrote:
> Hello, all.
>
> I wanted to use the "survey" package to analyze data from the National
> Ambulatory Medical Care Survey, and am having some difficulty translating
> the analysis keywords from one package (Stata) to the other (R). The data
> we
On Sat, 1 Oct 2005, Jim Hurd wrote:
>
> Which provides data in DTA (STATA), XPT (SAS), and POR (SPSS) formats all
> of which I have tried to read with the foreign package but I am not able to
> load any of them. I have 2 gb of RAM, but R crashes when the memory gets
> just over 1 GB. I am using Wi
On Tue, 4 Oct 2005, Karin Lagesen wrote:
>
> First, how does boxplot determine the size of the box? And is the line
> inside the box the mean or the median (or something completely
> different?) And how does it determine how long out the whiskers should
> go?
Part of the problem is that there are
On Fri, 30 Sep 2005, Erin Hodgess wrote:
> Dear R People:
>
> I have R Version 2.1.1. for Windows in binary form.
>
> I would like to look at the C code for massdist. It is part
> of the density function.
>
The .C() call shows that massdist is in the base package. This means that
it is likely t
On Fri, 30 Sep 2005, Mike Prager wrote:
> Recent R function names seem to be using CaseOfTheLetters to mark words
> rather than dots as was done previously. Is the use of dots in function
> names deprecated, or is that simply a style choice? Will function names
> with dots cause problems in futu
On Thu, 29 Sep 2005, Spencer Graves wrote:
> Questions like this are best directed to the package maintainer(s).
> From help(package="leaps"), I learned that Thomas Lumley is the author and
> maintainer for "leaps"; I'm including him as a '
On Thu, 29 Sep 2005, Karin Lagesen wrote:
>
> I have a file like this:
>
>
> a 0.1
> a 0.2
> a 0.9
> b 0.5
> b 0.9
> b 0.7
> c 0.6
> c 0.99
> c 0.88
>
> Which I would like to get to be the following matrix:
>
> 0.1 0.20.30.4 .
On Wed, 28 Sep 2005, Denis Chabot wrote:
> But what about another analogy, that of polynomials? You may not be sure what
> degree polynomial to use, and you have not decided before analysing your
> data. You fit different polynomials to your data, checking if added degrees
> increase r2 suffici
listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Thomas Lumley Assoc. Professor, Biostatistics
[EMAIL PROTECTED] University of Washington, Seattle__
R-help@stat.math.ethz.ch mailing list
h
On Tue, 27 Sep 2005, Ferran Carrascosa wrote:
> Hi everyone,
>
> I would like some package of R or any help to solve the next problem
> with a weighting fatcors:
>
> Giving a data matrix with dichotomous (2 or more) variables in columns
> and individuals in rows, and also a theorical distribution
On Tue, 27 Sep 2005, Simon.Bond wrote:
>
> Does anyone have any suggestions for a way round this problem for the
> future (I had to resort to using Stata), or maybe more realistically, how
> much work it would take to build an extendible version of the gee
> "algorithm"?
>
I don't think it would
On Tue, 27 Sep 2005, Samuel Bertrand wrote:
> Hello,
>
> I am using the 'regsubsets' function
> (from leaps package)
> to get the best linear models
> to explain 1 variable
> from 1 to 5 explanatory variables
> (exhaustive search).
>
> Is there anyone who can tell me
> on which criterion is based
On Mon, 26 Sep 2005, Denis Chabot wrote:
>
> But the mgcv manual warns that p-level for the smooth can be
> underestimated when df are estimated by the model. Most of the time
> my p-levels are so small that even doubling them would not result in
> a value close to the P=0.05 threshold, but I have
he problem, and show us that one,
> e.g. by making it available via http://... if it is still large,
> otherwise (if it's small), maybe even posting the result of
> dump(..).
>
> Regards,
> Martin
>
> __
> R-help@stat.math.ethz.ch mailing list
> https://s
On Thu, 22 Sep 2005, Anne wrote:
> DeaR list
> I would like to know if there is a direct method to compute the
> expected remaining survival time for a subject having survived up to
> time t. survexp gives me the probabilty for subject S to survive up to
> day D
You can't estimate expected surv
On Tue, 13 Sep 2005, Zhen Zhang wrote:
> Hello everyone,
>
> I tried to use coxph.detail() to get the hazard function. But a warning
> messge always returns to me, even in the example provided by its help
> document:
>
>> ?coxph.detail
>> fit <- coxph(Surv(futime,fustat) ~ age + rx + ecog.ps, o
On Tue, 13 Sep 2005, [EMAIL PROTECTED] wrote:
> A recent exchange on the 'octave' list led to the following
> paper being cited, which I had not met before:
>
> What Every Computer Scientist Should Know About
> Floating-Point Arithmetic, by David Goldberg,
>
> So I'm writing to bring it to the n
On Mon, 12 Sep 2005, Chris Buddenhagen wrote:
>
> Is there a simple means of doing this multiple comparison test in R?
>
p.adjust()
-thomas
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the
2 (0)16 322831
>> <http://www.kuleuven.be/ucs> http://www.kuleuven.be/ucs
>>
>>
>>
>> Disclaimer: http://www.kuleuven.be/cwis/email_disclaimer.htm
>>
>>
>> [[alternative HTML version deleted]]
>>
>> _______
On Sun, 11 Sep 2005, Arkady Sherman wrote:
>
> "c:\Program Files\R\rw2011\bin\R.exe" --no-save < test > out.txt
>
> the file "c:/temp/foo.txt" will contain nothing.
> But I'd like it should contain the warning message "Foo warning".
> Is the behavior a bug of R or there is another way to get it wor
http://faculty.washington.edu/tlumley/survey/NEWS
-thomas
Thomas Lumley Assoc. Professor, Biostatistics
[EMAIL PROTECTED] University of Washington, Seattle
___
R-packages mailing list
[EMAIL PROTECTED]
https://stat.ethz.
On Fri, 9 Sep 2005, Heinz Tuechler wrote:
> Dear Martin,
> Thank you for your answer. As I said, I appreciate this change. The
> documentation does not explain precisely, how variables with labels are
> treated now. It only tells "If SPSS value labels are converted to factors
> the underlying numer
n SPSS (we can't be sure that this is
always true since the file format is undocumented).
> However I could not find a description of this supposed change.
It's in the ChangeLog file in the foreign package.
-thomas
Thomas Lumley Assoc. Prof
X1 or X2...
>
> Many thanks in advance,
>
> Bernard
>
>
> -
>
>
> [[alternative HTML version deleted]]
>
> __
> R-help@stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
>
301 - 400 of 1258 matches
Mail list logo