reml随机效应模型meta合并出现提示使用不同数据集错误

Question

mvInterceptonly<-mvmeta(coef,vcov,method = "reml") #利用随机效应模型（reml）进行meta合并
summary(mvInterceptonly)#含单个自变量的meta回归模型结果展示
#####PM2.5
PM2.5 <- c(38.14 ,39.74 ,54.12 ,42.18 ,38.15 ,37.32 ,37.96 ,45.16 ,36.35 ,42.07 ,43.99, 
           43.34 ,47.94 )
mvPM2.5<- mvmeta(coef~PM2.5,vcov,method = "reml")
summary(mvPM2.5)#含单个自变量的meta回归模型结果展示
waldtest(mvInterceptonly,mvPM2.5)#Wald检验
Error in lrtest.default(mvInterceptonly, mvPM2.5) : 
  models were not all fitted to the same size of dataset

这个问题怎么解决，我用ml合并以后可以做，但是reml就不行，检查了模型的参数都是一样的

Answer 1

没用过这个包，但是一般不管什么模型里，所有变量都应该在一个数据里的
PM2.5这个为什么不一开始就放到dataset，你这从外部用，可能不规范