[R] multivariate lme?

From: Constantinos Antoniou <antoniou_at_central.ntua.gr>
Date: Wed 06 Jul 2005 - 07:34:17 EST


Dear all,

is it possible to estimate a multivariate multilevel model in R (I guess the term in R is mixed-effects model). I can estimate univariate models using either lme{nlme} or lmer{lme4}, but when I use a multivariate response, I get the same output as if I was trying to estimate a univariate model.

[For a reproducible example, please download http://mit.edu/costas/ www/parav.txt and run:

>Parav<-read.table("parav.txt",header=T)
>library(nlme)
>attach(Parav)

]

>myResponse <- cbind(accidents,Killed)

I get the same thing as if I used only accidents as the response.

For example:

>parav7.lme <- lme(accidents~alcontrols,random=~1+alcontrols|
department/prefecture,data=Parav)

yields the same output as:

>parav8.lme <- lme(myResponse~alcontrols,random=~1+alcontrols|
department/prefecture,data=Parav)

Thank you,

Costas

PS. I am working on R-2.1.0a on a Mac 10.4.1.

--
Constantinos Antoniou, Ph.D.
Department of Transportation Planning and Engineering
National Technical University of Athens
5, Iroon Polytechniou str. GR-15773, Athens, Greece

______________________________________________
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Received on Wed Jul 06 07:38:45 2005

This archive was generated by hypermail 2.1.8 : Fri 03 Mar 2006 - 03:33:16 EST