[R] help with lme function

[Olivier Martin]

Hi all,

I have some difficulties with the lme function and so this is my problem.
Supoose i have the following model
        y_(ijk)=beta_j + e_i + epsilon_(ijk)

where beta_j are  fixed effects, e_i is a random effect and 
epsilon_(ijk) is the error.

If i want to estimate a such model, i execute
 >lme(y~vec.J , random~1 |vec .I )

where y is the vector of my data, vec.J  is a  factor object and vec.I 
is the vector for the i indice.

Now , this is an other model
        y_(ijk)=beta_j + e_(ij) + epsilon_(ijk)

The  parameters estimations can be resolved with
 >lme(y~vec.J , random~1 |vec .IJ )

My problem concern the following model
        y_(ijk)=beta_j + e_i + e_(ij) + epsilon_(ijk)

where e_i and e_(ij) are the random effects.

In this case, what is the command to estimate the parameters?


Regards,
Olivier.


-- 
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
Olivier MARTIN
PhD student                  phone : (33) 04 76 61 53 55
Projet IS2                               06 08 67 93 42		
INRIA Rhone-Alpes            fax   : (33) 04 76 61 54 77
655, Av. de l'Europe
Montbonnot                   e-mail: olivier.martin at inrialpes.fr
38334 Saint Ismier cedex     web   : http://www.inrialpes.fr/is2/people/omartin
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-


-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._