Displaying 4 results from an estimated 4 matches for "gcorr".
Did you mean:
xcorr
2012 Sep 21
1
translating SAS proc mixed into R lme()
...everal times a year). I need to allow slope
and intercept vary.
SAS codes are:
proc mixed data = survey method=reml;
class subject var1 var3 var2 time;
model score = var2 score_base var4 var5 var3 var6 var7 var1 time/ noint
solution;
random intercept timecontinious / subject=subject type=un g gcorr v vcorr;
run;
Thank you a lot!
[[alternative HTML version deleted]]
2006 Jun 28
1
lme - Random Effects Struture
Thanks for the help Dimitris,
However I still have a question, this time I'll be more specific,
the following is my SAS code
proc mixed data=Reg;
class ID;
model y=Time Time*x1 Time*x2 Time*x3 /S;
random intercept Time /S type=UN subject=ID G GCORR V;
repeated /subject = ID R RCORR;
run; **
(Type =UN for random effects)
The eqivalent lme statement I am using is :
reglme <- lme(y ~ Time+Time*x1+Time*x2+Time*x3, data=Reg, random = ~ Time |
ID)
When I compare the results, the values differ by considerable margin; I
suppose...
2008 May 07
0
Help with Mixed effect modeling in R
...tement by the GROUP=grp option;
title 'RANDOM COEFFICIENT MODEL WITH DIAGONAL WITHIN-PATIENT';
title1 'DIFFERENT D MATRIX FOR BOTH GENDERS';
proc mixed method=ml data=dent1;
class pt grp;
model y = grp grp*t / noint solution;
random intercept t / type=un group=grp subject=pt g gcorr v vcorr;
run;
The key to specifying different covariance structure for the random
effects seems to be the highlighted portion in the code. What would be
it's equivalent in R?
In R, I tried the following
Model1 <- lme(y ~ g1+Tg1+g2+Tg2-1,random =
pdBlocked(list(pdSymm(~g1+Tg1-1),...
2007 Apr 06
0
translating sas proc mixed to lme()
...of the model, but a few things are still different.
It is a 2-level bivariate model (some call it a pseudo-3-level model).
PROC MIXED DATA=psdata.bivar COVTEST METHOD = ml;
CLASS cluster_ID individual_id variable_id ;
MODEL y = Dp Dq / SOLUTION NOINT;
RANDOM Dp Dq / SUBJECT = cluster_ID TYPE=UN G GCORR;
REPEATED variable_id / SUBJECT = individual_ID(cluster_ID) TYPE=UN R RCORR;
RUN;
Here is my try:
dta = sqlQuery(odbcConnect("sasodbc", believeNRows=FALSE), "select * from
psdata.bivar")
v = lme(Y~Dp+Dq-1, data=dta, random=~Dp+Dq-1|cluster_ID, method="ML",
weights=...