我有2项研究报告了随着时间的推移经历疼痛的个体的轨迹/比例(例如,1个月30%的经历疼痛和70%的患者,在这些分支的2个月内形成2个以上痛苦/无痛等分支。)
我想改变这两项研究中报告的比例(我知道),我有一个问题。说我为第一个时间点制作了数据帧:
study = c("1", "2")
ni = c(539, 89)
xi = c(121, 31)
xii = c(418, 58)
onemonth <-data.frame(study, xi, xii, ni)
然后我想模拟和logit转换我的数据..我不确定是否应该使用随机效应模型来尝试并解释研究之间的差异
res <- rma(measure = "PLO", xi = xi, ni = ni, data = onemonth)
或者如果我应该使用随机效应逻辑回归模型,因为我的结果是二分的。我见过两个报道。
res <- rma.glmm(measure = "PLO", xi = xi, ni = ni, data = onemonth)
然后我会将数据转换回来:
predict(res, transf = transf.ilogit, digits = 2)
非常感谢任何/所有帮助和想法!