similar to: survreg anova: problem with indirect invocation

Hi. I'm using lmer() from the lme4 package (version 0.8-3) and I can't get anova() to group variables properly. I'm fitting the mixed model Response ~ Weight + Experimenter + (1|SUBJECT.NAME) + (1|Date.StudyDay) where Weight is numeric and Experimenter is a factor, ie, > str(data.df) `data.frame': 4266 obs. of 5 variables: $ SUBJECT.NAME : Factor w/ 2133 levels

Does anyone know if the survreg function in the survival package can fit numerical covariates ? When I fit a survival model of the form survreg( Surv(time,censored) ~ x ) then x is always treated as a factor even if it is numeric (and even if I try to force it to be numeric using as.numeric(x). Thus, in the particular example I am analysing, a simple numerical covariate becomes a factor

I am trying to use R to do loglogistic hazard estimation. My plan is to generate a loglogistic hazard sample data and then use survreg to estimate it. If everything is correct, survreg should return the parameters I have used to generate the sample data. I have written the following code to do a time invariant hazard estimation. The output of summary(modloglog) shows the factor loading of

Dear All, I want to test whether a method exists for given object. For example, whether a function "deviance" is defined for an object of the "lm" class. My imperfect understanding leads me to think something like hasMethod("deviance", object) hasMethod("deviance", "lm") existsMethod("deviance", signature(class="lm"))

Dear All, Is there is a fast way of refitting lm() when the design matrix stays constant but the response is different? For example, y1 ~ X y2 ~ X y3 ~ X ...etc. where y1 is the 1st instance of the response vector. Calling lm() every time seems rather wasteful since the QR-decomposition of X needs to be calculated only once. It would be nice if qr() was called only once and then the same

Dear All, I have been trying to simulate data from a fitted glm using the simulate() function (version details at the bottom). This works for lm() fits and even for lmer() fits (in lme4). However, for glm() fits its output does not make sense to me -- am I missing something or is this a bug? Consider the following count data, modelled as gaussian, poisson and binomial responses: counts

Full_Name: Andrew Manners Version: 2.5.1 OS: windows xp prof 2003 Submission from: (NULL) (130.102.0.177) To whom it may concern, I'm trying to get an ANOVA table within survreg but it always produces NA's in the p-value, regardless of the data set. The data set below comes from Tableman and Kim 2004. I had the same problem on a number of my own data sets. I searched the R site for

Dear R Listers- I am attempting to analyse the survival of seeds in cages (exclosures) that differ in their permeability to rainforest mammals. Because I did not observe the moment of seed disappearance, my data is interval censored. This limits my options for analysis (as I understand it) to survreg, in the survival package. Because I repeated the experiment in 8 sites, I have a random

What objective function is maximized by survreg with the default Weibull model? I'm getting finite parameters in a case that has the likelihood maximzed at Infinite, so it can't be a simple maximum likelihood. Consider the following: ############################# > set.seed(3) > Stress <- rep(1:3, each=3) > ch.life <- exp(9-3*Stress) > simLife <- rexp(9,

Dear All, I wish to model random effects that have known between-group covariance structure using the lme() function from library nlme. However, I have yet to get even a simple example to work. No doubt this is because I am confusing my syntax, but I would appreciate any guidance as to how. I have studied Pinheiro & Bates carefully (though it's always possible I've missed

Hi, I make a weibull survival regression using suvreg function. Bu when I try to get the P values from anova, it give me NAs: I'm using R 2.4.0 and survival 2.29 Look: m <- survreg(Surv(tempo,censor)~grupo*peso) anova(m) Df Deviance Resid. Df -2*LL P(>|Chi|) NULL NA NA 148 966.6416 NA grupo -2 25.6334407 146 941.0081 NA

Hi, I make a weibull survival regression using suvreg function. Bu when I try to get the P values from anova, it give me NAs: I'm using R version 2.4.0 Patched (2006-11-25 r39997) and survival 2.30 library Look: m <- survreg(Surv(tempo,censor)~grupo*peso) anova(m) Df Deviance Resid. Df -2*LL P(>|Chi|) NULL NA NA 148 966.6416 NA grupo

Hi, Is it normal to get intercept in the list of covariates in the output of survreg function with standard error, z, p.value etc? Does it mean that intercept was fitted with the covariates? Does Value column represent coefficients or some thing else? Regards, ------------------------------------------------- tmp = survreg(Surv(futime, fustat) ~ ecog.ps + rx, ovarian,

Dear R-Helpers: I have tried everything I can think of and hope not to appear too foolish when my error is pointed out to me. I have some real data (18 points) that look linear on a log-log plot so I used them for a comparison of lm() and survreg. There are no suspensions. survreg.df <- data.frame(Cycles=c(2009000, 577000, 145000, 376000, 37000, 979000, 17420000, 71065000, 46397000,

Can anybody give me a neat example of interval censored data analysis codes in R? Given that suvreg(Surv(c(1,1,NA,3),c(2,NA,2,3),type="interval2")~1) works why does survreg(Surv(data[,1],data[,2],type="interval2")~1) not work where data is : T.1 T.2 Status 1 0.0000000 0.62873036 1 2 0.0000000 2.07039068 1 3 0.0000000

Dear R-users, I know that a few people have asked whether survreg handles left-truncation data and the reply that i have seen so far is that it does. However, when I try to use survreg on left-truncated data, I got the following error message. > survcs3<-survreg(Surv(start,end,status)~AG, data=DPONEcs3, dist="exponential") Error in survreg(Surv(start, end, status) ~ AG, data =

Hi there. I am receiving an unexpected error message when creating a new distribution for the survreg() function in the survival package. I understand the survival.distributions() function and have been following the Cauchy example provided in the help file. My goal is to use survreg to fit a gamma distribution to interval censored data. Here is a simple example of what I'm trying to do.

Hi, I would like to fit parametric survival models to time-to-event data that are left truncated. I have checked the help page for survreg and looked in the R-help archive, and it appears that the R function survreg from the survival library (version 2.16) should allow me to take account of left truncation. However, when I try the command

---------------------------- Original Message ---------------------------- Subject: Different results from survreg with version 2.6.1 and 2.10.1 From: nathalcs at ulrik.uio.no Date: Fri, March 19, 2010 16:00 To: r-help at r-project.org -------------------------------------------------------------------------- Dear all I'm using survreg command in package survival.

Dear R-list I have been trying to make survreg fit a normal regression model with left truncated data, but unfortunately I am not able to figure out how to do it. The following survreg-call seems to work just fine when the observations are right censored: library(survival) n<-100000 #censored observations x<-rnorm(n) y<-rnorm(n,mean=x) d<-data.frame(x,y) d$ym<-pmin(y,0.5)