CAN'T HANDLE augPred ()
On Sun, Feb 6, 2011 at 3:44 PM, Federico Bonofiglio <bonoricus at gmail.com> wrote:
Hi Dears, the model below has been ran standardly, but while using augPred I receive an error message
cd4CORR$call
*lme.formula(fixed = cd4 ~ rt + age + pharmac + haartd, random = ~rt | ? ?id, correlation = corGaus(form = ~rt), subset = which(cd4 < ? ?2200), method = "REML", na.action = na.omit) *
Use the data argument in the call to lme - otherwise it is difficult to determine the original data frame in order to do the augmentation.
cd4CORR$dims
$N [1] *7023* $Q [1] 1 $qvec [1] 2 0 0 $ngrps ?id ? X ? y *868* ? 1 ? 1 $ncol [1] 2 6 1
plot(augPred(cd4CORR))
*Errore in augPred.lme(cd4CORR) : ?Data in cd4CORR call must evaluate to a data frame Errore in plot(augPred(cd4CORR)) : ?error in evaluating the argument 'x' in selecting a method for function 'plot' * ..any hint is well appreciated... and I thank in advance.. Federico -- *Little u can do against ignorance,....it will always disarm u: is the 2nd principle of thermodinamics made manifest, ...entropy in expansion.**....But setting order is the real quest 4 truth, ......and the mission of a (temporally) wise dude. * ? ? ? ?[[alternative HTML version deleted]]
_______________________________________________ R-sig-mixed-models at r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models