Hi,
I'm analyzing a longitudinal data set with 387 cows were
observed in 63 days divided into 6 groups, and every 30 days
was found to produce milk. Does not aim to model the time
using regression. Only compare the groups differ in terms of
milk production. There are many missing observations.
Because the data are correlated I used the SAS program:
proc mixed data=univar method=reml;
class RACA GRUPO APELIDO Dias;
model Prod = GRUPO / solution DDFM=BW;
repeated Dias / type=arh(1) subject=APELIDO r rcorr;
lsmeans GRUPO / pdiff adjust=tukey;
run ;
But, I want use R. What would be the equivalent in R?
Thank you.
--------------------------------------
Silvano Cesar da Costa
Departamento de Estat?stica
Universidade Estadual de Londrina
Fone: 3371-4346
Longitudinal data
3 messages · Uwe Ligges, Silvano
On 27.12.2011 14:43, Silvano wrote:
Hi, I'm analyzing a longitudinal data set with 387 cows were observed in 63 days divided into 6 groups, and every 30 days was found to produce milk. Does not aim to model the time using regression. Only compare the groups differ in terms of milk production. There are many missing observations. Because the data are correlated I used the SAS program: proc mixed data=univar method=reml; class RACA GRUPO APELIDO Dias; model Prod = GRUPO / solution DDFM=BW; repeated Dias / type=arh(1) subject=APELIDO r rcorr; lsmeans GRUPO / pdiff adjust=tukey; run ;
See package SASmixed that includes some nice examples how to move from SAS to R. Uwe Ligges
But, I want use R. What would be the equivalent in R? Thank you. -------------------------------------- Silvano Cesar da Costa Departamento de Estat?stica Universidade Estadual de Londrina Fone: 3371-4346
______________________________________________ R-help at r-project.org 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.
Hi Uwe, was a great suggestion. Thanks, -------------------------------------- Silvano Cesar da Costa Departamento de Estat?stica Universidade Estadual de Londrina Fone: 3371-4346 -------------------------------------- ----- Original Message ----- From: "Uwe Ligges" <ligges at statistik.tu-dortmund.de> To: "Silvano" <silvano at uel.br> Cc: <r-help at r-project.org> Sent: Tuesday, December 27, 2011 1:13 PM Subject: Re: [R] Longitudinal data
On 27.12.2011 14:43, Silvano wrote:
Hi, I'm analyzing a longitudinal data set with 387 cows were observed in 63 days divided into 6 groups, and every 30 days was found to produce milk. Does not aim to model the time using regression. Only compare the groups differ in terms of milk production. There are many missing observations. Because the data are correlated I used the SAS program: proc mixed data=univar method=reml; class RACA GRUPO APELIDO Dias; model Prod = GRUPO / solution DDFM=BW; repeated Dias / type=arh(1) subject=APELIDO r rcorr; lsmeans GRUPO / pdiff adjust=tukey; run ;
See package SASmixed that includes some nice examples how to move from SAS to R. Uwe Ligges
But, I want use R. What would be the equivalent in R? Thank you. -------------------------------------- Silvano Cesar da Costa Departamento de Estat?stica Universidade Estadual de Londrina Fone: 3371-4346
______________________________________________ R-help at r-project.org 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.