[R] Plot survival curves after coxph() with frailty() random effects terms

2017-06-23 Thread Andreu Ferrero
the next code I get the 2 groups. surv_group <- survfit(Surv(time, status) ~ sex, rats) lines(surv_group[1:2], lwd=2, conf.int=FALSE) However, I am not sure about these 2 curves are well done, appropriate. If any of you could help… -- Andreu Ferrero Gregori [[alt

[R] Can mice() handle crr()? Fine-Gray model, Object has no vcov() method.

2017-01-22 Thread Andreu Ferrero
ance=TRUE)) summary(pool(models.FG)) #Error in pool(models.FG) : Object has no vcov() method. models.FG -- Andreu Ferrero Gregori [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see http

[R] frailtypack - Model did not converge

2016-06-20 Thread Andreu Ferrero
n= 2473 n recurrent events= 406 n terminal events= 195" Any idea??? Cause PC computing takes hours to "fit/non-fit" this model. Thanks, Andreu Ferrero Gregori [[alternative HTML version deleted]] __ R-help@r-pr

[R] glmer() -> corrected AUC optimism by bootstraping technic bootMer() [internal validation of a mixed-effects-model]

2016-03-24 Thread Andreu Ferrero Gmail
> > > I would like to do an internal validation of a discriminative ability of a > mixed effects models. > > Here is my scrip: > > ### > bootMer-> boot AUC# > ### > > library(lme4) > library(lattice) > data(cbpp) > > #fit a model >