我想进行将个人数据和汇总数据与分层相关回归相结合的研究。因此，我想使用ecoreg包。
以下是我插入的数据的代码和摘要。
pDV n pIV
0.02974330 17920 0.02779018
0.02843174 13295 0.02564874
0.02761653 11370 0.02515391
0.01849075 10005 0.01589205
0.02738854 7850 0.02267516
0.02488095 8400 0.01904762
0.02459336 7685 0.02121015
0.02352941 6375 0.01866667
0.02079295 5675 0.01797357
eco(cbind(pDV, n) ~ pIV, data = agg)
但是我一直收到同样的错误消息。
checkdata（X，Y，补充，ndim）出错： X和Y的值必须介于0和1之间。
我的所有值都在0到1之间，n除外。但是，在ecoreg文件中，它们也使用值大于1的n，并且能够获得正确的输出。如果我删除了n。我能够得到输出，但是我没有得到我想要的输出，例如：
Call:
eco(formula = cbind(pIV) ~ pDV, data = agg)
Aggregate In-sample Estimates:
W1 W2
posterior mean 0.948122 0.008483
相反，我希望在ecoreg文件中拥有与Jackson相同的输出。即，利用比值比和对数似然：
eco(cbind(disease, n) ~ p.smoke, data = agg)
Call:
eco(formula = cbind(disease, n) ~ p.smoke, data = agg)
Aggregate-level odds ratios:
OR l95 u95
(Intercept) 0.31 0.3 0.32
p.smoke 3.61 2.7 4.74
No individual-level covariates
-2 x log-likelihood: 2141
所以我希望得到与上面类似的输出，我不明白为什么n会在我的控制台中出现问题，但它似乎适用于ecoreg指南。
有人可以告诉我我做错了什么。这将是非常有帮助的！你可以在这里找到ecoreg的指南：https://cran.r-project.org/web/packages/ecoreg/vignettes/ecoreg-guide.pdf