Re: [R] multivariate lme?

From: Douglas Bates <dmbates_at_gmail.com>
Date: Wed 06 Jul 2005 - 22:53:00 EST

Neither lme nor lmer has provision for fitting multivariate models.

On 7/5/05, Constantinos Antoniou <antoniou@central.ntua.gr> wrote:
> 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
>



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 22:56:12 2005

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