我正在进行荟萃分析，并在模型中进行了双弧转换。 我目前正在对主持人变量（实验室方法）进行子组分析。主持人变量有2个组（实验1和实验2）。 lab1中有6个对象，而lab2中有9个对象。我的数据框（一个csv文件）有15行数据。
我正在使用以下代码：
subganal.lab_method=rma(yi, vi, data=ies.da, mods=~lab_method)
pes.da.lab1=rma(yi, vi, data=ies.da, mods=~lab_method=="lab1")
pes.da.lab2=rma(yi, vi, data=ies.da, mods=~lab_method=="lab2")
pes.subg.lab_method=predict(subganal.lab_method,
transf=transf.ipft.hm,
targ=list(ni=dat$samples_total))
dat.samevar=data.frame(estimate=c((pes.da.lab1$b)[1], (pes.da.lab2$b)[1]),
stderror=c((pes.da.lab1$se)[1], (pes.da.lab2$se)[1]),
tau2=subganal.lab_method$tau2)
pes.da.lab_method=rma(estimate, sei=stderror,
method="DL",
data=dat.samevar)
pes.lab_method=predict(pes.da.lab_method, transf=transf.ipft.hm, targ=list(ni=dat$samples_total))
print(pes.subg.lab_method ["lab1"]) #display subgroup 1 summary effect size
print(pes.subg.lab_method["lab2"]) #display subgroup 2 summary effect size
print(subganal.lab_method) #display subgroup analysis results
print(pes.lab_method) #display recomputed summary effect size
我两次收到错误消息“冗余预测变量已从模型中删除”。一次用于lab1，一次用于lab2。 此外，不打印lab1和2的子组摘要效果。它会在应打印值的位置打印“ NULL”。
不会打印总体子组分析。
“ dat”是我的数据框的名称