[R] converting proc mixed to lme for a random effects meta-analysis

From: Lucia Costanzo <lcostanz_at_uoguelph.ca>
Date: Tue, 19 Jun 2007 08:13:30 -0400


I would like to convert the following SAS code for a Random Effects meta-analysis model for use in R but, I am running into difficulties. The results are not similar, R should be reporting 0.017 for the between-study variance component, 0.478 for the estimated parameter and 0.130 for the standard error of the estimated parameter. I think it is the weighting causing problems. Would anyone have any suggestions or tips?

Thank you,
Lucia

re.teo<-lme(y~1, data=genData2, random =~1, method="ML", weights=varFixed(~w))

data tacrine;

    input study y w;
    cards;

    1 0.284 14.63
    2 0.224 17.02
    3 0.360  9.08
    4 0.785 33.03
    5 0.492  5.63

    ;
run;

*Random Effects using log-odds for tacrine example table 4.29;
DATA remlma;

    SET tacrine;
    var=1/w;
    col = _n_;
    row = _n_;
    value = var;
run;

*random effects for tacrine example;

PROC MIXED data = remlma method=reml order=data;

   CLASS study;
   MODEL y = / solution;
   RANDOM study / gdata = remlma;
   REPEATED diag;
RUN;



R-help_at_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 Tue 19 Jun 2007 - 12:19:17 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Tue 19 Jun 2007 - 13:32:03 GMT.

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