analyses:gleser2009
Differences
This shows you the differences between two versions of the page.
Both sides previous revisionPrevious revisionNext revision | Previous revision | ||
analyses:gleser2009 [2021/09/24 16:34] – Wolfgang Viechtbauer | analyses:gleser2009 [2024/06/09 16:19] (current) – Wolfgang Viechtbauer | ||
---|---|---|---|
Line 79: | Line 79: | ||
Now the model described in equation (19.5) can be fitted with: | Now the model described in equation (19.5) can be fitted with: | ||
<code rsplus> | <code rsplus> | ||
- | res <- rma.mv(yi, V, mods = ~ factor(trt) | + | res <- rma.mv(yi, V, mods = ~ 0 + factor(trt), |
res | res | ||
</ | </ | ||
Line 127: | Line 127: | ||
V <- bldiag(lapply(split(dat, | V <- bldiag(lapply(split(dat, | ||
- | res <- rma.mv(yi, V, mods = ~ factor(trt) | + | res <- rma.mv(yi, V, mods = ~ 0 + factor(trt), |
res | res | ||
</ | </ | ||
Line 165: | Line 165: | ||
V <- bldiag(lapply(split(dat, | V <- bldiag(lapply(split(dat, | ||
- | res <- rma.mv(yi, V, mods = ~ factor(trt) | + | res <- rma.mv(yi, V, mods = ~ 0 + factor(trt), |
res | res | ||
</ | </ | ||
Line 202: | Line 202: | ||
V <- bldiag(lapply(split(dat, | V <- bldiag(lapply(split(dat, | ||
- | res <- rma.mv(yi, V, mods = ~ factor(trt) | + | res <- rma.mv(yi, V, mods = ~ 0 + factor(trt), |
res | res | ||
</ | </ | ||
Line 283: | Line 283: | ||
The analysis can now be carried out with: | The analysis can now be carried out with: | ||
<code rsplus> | <code rsplus> | ||
- | res <- rma.mv(yi, V, mods = ~ factor(trt) | + | res <- rma.mv(yi, V, mods = ~ 0 + factor(trt), |
print(res, digits=3) | print(res, digits=3) | ||
</ | </ | ||
Line 291: | Line 291: | ||
Variance Components: none | Variance Components: none | ||
- | Test for Residual Heterogeneity: | + | Test for Residual Heterogeneity: |
QE(df = 4) = 3.945, p-val = 0.414 | QE(df = 4) = 3.945, p-val = 0.414 | ||
- | Test of Moderators (coefficient(s) 1,2): | + | Test of Moderators (coefficient(s) 1,2): |
QM(df = 2) = 252.165, p-val < .001 | QM(df = 2) = 252.165, p-val < .001 | ||
Model Results: | Model Results: | ||
- | estimate | + | estimate |
factor(trt)1 | factor(trt)1 | ||
factor(trt)2 | factor(trt)2 | ||
Line 315: | Line 315: | ||
factor(trt)2 | factor(trt)2 | ||
</ | </ | ||
- | The same results are given on page 367. | + | The same results are given on page 367. |
==== Multiple-Endpoint Studies ==== | ==== Multiple-Endpoint Studies ==== | ||
Line 336: | Line 336: | ||
dat$yi <- round(with(dat, | dat$yi <- round(with(dat, | ||
dat$vi <- round(with(dat, | dat$vi <- round(with(dat, | ||
- | dat$covi <- round(with(dat, | + | dat$covi <- round(with(dat, |
+ | | ||
+ | | ||
</ | </ | ||
The contents of the resulting dataset are: | The contents of the resulting dataset are: | ||
Line 359: | Line 361: | ||
The variance-covariance matrix for the entire dataset can be constructed with: | The variance-covariance matrix for the entire dataset can be constructed with: | ||
<code rsplus> | <code rsplus> | ||
- | V <- bldiag(lapply(split(dat, | + | V <- bldiag(lapply(split(dat, |
+ | | ||
V | V | ||
</ | </ | ||
Line 383: | Line 386: | ||
Finally, we can fit a model allowing for a different treatment effect depending on the outcome with: | Finally, we can fit a model allowing for a different treatment effect depending on the outcome with: | ||
<code rsplus> | <code rsplus> | ||
- | res <- rma.mv(yi, V, mods = ~ outcome | + | res <- rma.mv(yi, V, mods = ~ 0 + outcome, data=dat) |
print(res, digits=3) | print(res, digits=3) | ||
</ | </ | ||
Line 391: | Line 394: | ||
Variance Components: none | Variance Components: none | ||
- | Test for Residual Heterogeneity: | + | Test for Residual Heterogeneity: |
QE(df = 12) = 19.626, p-val = 0.074 | QE(df = 12) = 19.626, p-val = 0.074 | ||
- | Test of Moderators (coefficient(s) 1,2): | + | Test of Moderators (coefficient(s) 1,2): |
QM(df = 2) = 13.005, p-val = 0.001 | QM(df = 2) = 13.005, p-val = 0.001 | ||
Model Results: | Model Results: | ||
- | estimate | + | estimate |
outcomemath | outcomemath | ||
outcomereading | outcomereading |
analyses/gleser2009.txt · Last modified: 2024/06/09 16:19 by Wolfgang Viechtbauer