search for: gcorr

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=...