From: souvik banerjee <bansouvik_at_gmail.com>

Date: Thu 10 Aug 2006 - 18:54:12 EST

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 and provide commented, minimal, self-contained, reproducible code. Received on Thu Aug 10 18:58:22 2006

Date: Thu 10 Aug 2006 - 18:54:12 EST

Hi

Is there any way to fit a bivaraite longitudinal mixed model using R. I have
a data set with col names

resp1 (Y_ij1), resp2 (Y_ij2), timepts (t_ij), unit(i)

I want to fit the following two models

Y_ij1, Y_ij2 | U_i = u_i ~ N(alpha + u_i + beta1*t_ij, Sigma) U_i ~ iid N(0, sigu^2)

Sigma = bivariate AR structure

alpha and beta are vectors of order 2.

Model 2

Y_ij, Y_ij2 | U_i = u_i ~ N(alpha + u_i + beta1*t_ij + U_i2*t_ij, Sigma)

U_i=(U_i1,U_i2) ~ iid N(0, Omega)

Here alpha and beta are vectors of order 2.

Moreover can we assume a structure for Omega.

[[alternative HTML version deleted]]

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 and provide commented, minimal, self-contained, reproducible code. Received on Thu Aug 10 18:58:22 2006

Archive maintained by Robert King, hosted by
the discipline of
statistics at the
University of Newcastle,
Australia.

Archive generated by hypermail 2.1.8, at Mon 14 Aug 2006 - 04:19:07 EST.

*
Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help.
Please read the posting
guide before posting to the list.
*