Hi all:
I am a user of "JM" package.
Here's the problem of "sample size".
The warning is:
Error in jointModel(fitLME, fitSURV_death, timeVar = "time", method =
"piecewise-PH-GH") :
sample sizes in the longitudinal and event processes differ.
According to the suggestion of "missing data",I use the same data set(data_JM)
without any missing value.
fitLME <- lme(CD4 ~ time + time:Group_DR, random = ~time | id,data = data_JM);
fitSURV_death <- coxph(Surv(time_fail, death) ~ (Group_DR), data = data_JM, x =
TRUE);
summary(fitLME);
Linear mixed-effects model fit by REML
Data: data_JM
AIC BIC logLik
9102.502 9134.538 -4544.251
... ...
Number of Observations: 721
summary(fitSURV_death)
Call:
coxph(formula = Surv(time_fail_min, death) ~ (Group_DR), data = data_JM,
x = TRUE)
n= 721
coef exp(coef) se(coef) z Pr(>|z|)
Group_DR -0.16854 0.84490 0.04213 -4 6.33e-05 ***
---
Signif. codes: 0 ¡®***¡¯ 0.001 ¡®**¡¯ 0.01 ¡®*¡¯ 0.05 ¡®.¡¯ 0.1 ¡® ¡¯ 1
The sample size of lme and coxph is the same as marked blue above,and I'm sure
that data_JM contains no missing value.
So, what's reason about the warning "Error in jointModel(fitLME, fitSURV_death,
timeVar = "time", method = "piecewise-PH-GH") :
sample sizes in the longitudinal and event processes differ."
Thanks a lot .
My best.
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
