[R-sig-ME] CAN'T HANDLE augPred ()

Douglas Bates bates at stat.wisc.edu
Thu Feb 10 17:12:59 CET 2011


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
>




More information about the R-sig-mixed-models mailing list