Dear listers,

I experinced a problem prdicting the values using the LME with multilevel
data.
I have NA's in my dependent variable and the model is fitted only on the
completed cases.
I want to estimate the predicted values for the rest of the data (those
cases with missing dep. variable)
I extracted a subset from the original file containing the variables used in
the model as well as the second level indicator.
I used the following command

p<-predict(level2,newdata=d.n.new,level=0:1)

where level2 is my LME model.
But, I get the following error:

Error in eval(expr, envir, enclos) : 1 argument passed to "$" which requires
2.

I tried with omitting the level specification (which is 0 by default) and I
transformed the new data to be groupedData with no luck.

I have tried the example from the Pinheiro,Bates book and it works - mine
doesn't. Does anybody have an idea what could be wrong?

Thanks for all the suggestions.

Andrej

_________
Andrej Kveder, M.A.
researcher
Institute of Medical Sciences SRS SASA; Novi trg 2, SI-1000 Ljubljana,
Slovenia
phone: +386 1 47 06 440   fax: +386 1 42 61 493

	[[alternative HTML version deleted]]


