[R-meta] Weighing factors in rma.mv corresponding to those used in lmer
Martineau, Roger
Roger.Martineau at AGR.GC.CA
Tue Mar 27 22:56:27 CEST 2018
Dear list members,
I want to use the rma.mv function and compare the output with that from the lmer function.
I need to compare 3 weighting schemes: (1) the number of replica (n), (2) the inverse of SE and (3) the inverse of variance or 1/SE^2 in each PubID.
(res1 <- lmer(y ~ x1 + x2 + (1|PubID), data=Dataset, weights = n, na.action = "na.omit"))
(res2 <- lmer(y ~ x1 + x2 + (1|PubID), data=Dataset, weights = 1/SE, na.action = "na.omit"))
(res3 <- lmer(y ~ x1 + x2 + (1|PubID), data=Dataset, weights = 1/SE^2, na.action = "na.omit"))
What are the rma.mv scripts corresponding to res1, res2 and res3 ?
Thanks in advance
Roger ☺
Roger Martineau, mv Ph.D.
Centre de recherche et de développement
sur le bovin laitier et le porc
Agriculture et agroalimentaire Canada/Agriculture and Agri-Food Canada
Téléphone/Telephone: 819-780-7319
Télécopieur/Facsimile: 819-564-5507
2000, Rue Collège / 2000, College Street
Sherbrooke (Québec) J1M 0C8
Canada
roger.martineau at agr.gc.ca<mailto:roger.martineau at agr.gc.ca>
[[alternative HTML version deleted]]
More information about the R-sig-meta-analysis
mailing list