HTML5   发布时间:2022-04-25  发布网站:大佬教程  code.js-code.com
大佬教程收集整理的这篇文章主要介绍了R:如何在包meta的metaprop中指定置信区间大佬教程大佬觉得挺不错的,现在分享给大家,也给大家做个参考。
library(Meta)
event <- c(81,15,1)
n <- c(263,148,20,29)
#
m1 <- Metaprop(event,n,sm="PLOGIT",method.ci="SA")

我有兴趣在荟萃分析中结合比例.在上面的例子中,我有4个研究,每个研究报告一个比例.并且似乎Metaprop计算4项研究中的每一项的CI.但是,由于4项研究已经报告了CI(在原始论文中),是否有办法将实际的,报告的CI纳入Meta分析计算,而不是让Metaprop自己计算它们?

如果他们允许我指定CI,我也愿意探索其他包.

解决方法

首先,请记住,您可能犯了一个错误.我没有看到不从Metaprop函数获取CI的原因.您应该从包中获取计算出的CI,特别是如果您使用的是随机效应模型.平均个体效应和CI的计算方法与固定效应的不同.无论如何,这里是:

您可以在以下位置找到较低的置信区间:

> m1$lower
[1] 0.25219013 0.05273001 0.00000000 0.00000000

以及信心上限:

> m1$upper
[1] 0.3637795 0.1499727 0.0000000 0.1008922

您还可以访问TE中的各个平均效果

> m1$TE
[1] -0.8095575 -2.1822989 -3.7135721 -3.3322045

然后,您可以在这些向量中手动添加自己的计算值,它们将显示在您的林图中.

假设您自己的较低值是myLCL:

myLCL <- seq(-1,-.25,.25) # define lower confidence limits

m1$lower <- myLCL # assign them to the respective location
m1$lower
[1] -1.00 -0.75 -0.50 -0.25

你的上限信心是myUCL:

myUCL <- seq(.25,1,.25) # define upper confidence limits

m1$upper <- myUCL # assign them to the respective location
m1$upper
[1] 0.25 0.50 0.75 1.00

大佬总结

以上是大佬教程为你收集整理的R:如何在包meta的metaprop中指定置信区间全部内容,希望文章能够帮你解决R:如何在包meta的metaprop中指定置信区间所遇到的程序开发问题。

如果觉得大佬教程网站内容还不错,欢迎将大佬教程推荐给程序员好友。

本图文内容来源于网友网络收集整理提供,作为学习参考使用,版权属于原作者。
如您有任何意见或建议可联系处理。小编QQ:384754419,请注明来意。
标签: