search for: lfit

Displaying 7 results from an estimated 7 matches for "lfit".

Did you mean: fit
2012 Jun 05
1
model.frame and predvars
...nd lm does not. I use a call to model.frame in predict.coxph, which is why I went that route, but never noted the difference till now (preparing for my course in Nashville). Could someone shed light on the rationale for non-preservation? Terry T. Simple example > library(survival) > lfit <- lm(time ~ factor(ph.ecog) + ns(age, 3), data=lung) > ltemp <- model.frame(lfit, data=lung[1:2,]) > ltemp time factor(ph.ecog) ns(age, 3).1 ns(age, 3).2 ns(age, 3).3 1 306 1 -0.1428571 0.4285714 0.7142857 2 455 0 0.0000000 0.0000000...
2010 Apr 06
0
Strange error
...get an error I have never seen before. fit <- lmekin(icam1 ~ factor(center) + age + factor(sex), random= ~1|iid, data=chaidata, varlist=kmat) Error in Y - fitted : non-numeric argument to binary operator Add the recover option, and the offending lines are fitted <- c(X %*% lfit$coef) #fitted, on the original scale residuals <- Y - fitted > options(error=recover) > fit <- lmekin(icam1 ~ factor(center) + age + factor(sex), + random= ~1|iid, data=chaidata, varlist=kmat) Error in Y - fitted : non-numeric argument to binary operator Enter a fra...
2012 Feb 03
1
coxme with frailty--variance of random effect?
Dear all, This probably stems from my lack of understanding of the model, but I do not understand the variance of the random effect reported in coxme. Consider the following toy example: #------------------------------- BEGINNING OF CODE ------------------------------------------------ library(survival) library(coxme) #--- Generate toy data: d <- data.frame(id = c(1:100), #
2009 Jun 29
0
nlsList {nlme} - control arguments problem
...t;port", convergence message: relative convergence (4) ################################## ### nlsList analysis ........................................................................................ ### Boundary condition on dplt is not enforced ............................................. Lfit.nls ? ?<- nlsList(y~SSpowrDplt(trial,min,dplt,dif,p),data=df.gr,control=list(algorithm='port',lower=c(0.0,0.0,0.0,-10),maxiter=100)) ...... example output....... >Call: > ?Model: y ~ SSpowrDplt(trial, min, dplt, dif, p) | id > ? Data: df.gr >Coefficients: > ? min > ?E...
2010 Apr 15
0
nlsList {nlme} - control arguments problem
...ative convergence (4) > ################################## > > ### nlsList analysis > ........................................................................................ > ### Boundary condition on dplt is not enforced > ............................................. > > Lfit.nls &n! bsp; <- > nlsList(y~SSpowrDplt(trial,min,dplt,dif,p),data=df.gr,control=list(algorithm='port',lower=c(0.0,0.0,0.0,-10),maxiter=100) > > ...... example output....... >Call: > Model: y ~ SSpowrDplt(trial, min, dplt, dif, p) | id > Data: df.gr > >Coef...
2013 Jun 05
0
Survival aareg problem
On 06/05/2013 12:33 AM, r-help-request at r-project.org wrote: > Dear friends - I'm on windows 7, R 2.15.2 > > when I run the example for aareg in survival package I see this: > > plot(lfit[4], ylim=c(-4,4)) > error in xy.coords(x, y, xlabel, ylabel, log) : > 'x' is a list, but does not have components 'x' and 'y' > > Is that a matter of an old R? > > Best wishes > > Troels Ring, MD > Aalborg, Denmark It's an error in the sur...
2010 Feb 17
1
Checking the assumptions for a proper GLM model
Hello, Are there any packages/functions available for testing the assumptions underlying assumptions for a good GLM model? Like linktest in STATA and smilar. If not, could somebody please describe their work process when they check the validity of a logit/probit model? Regards, Jay