问题标签 [metafor]

For questions regarding programming in ECMAScript (JavaScript/JS) and its various dialects/implementations (excluding ActionScript). Note JavaScript is NOT the same as Java! Please include all relevant tags on your question; e.g., [node.js], [jquery], [json], [reactjs], [angular], [ember.js], [vue.js], [typescript], [svelte], etc.

0 投票
1 回答
277 浏览

r - 双变量模型和“predict.rma()”函数(“metafor”包)

假设我有以下模型结果:

MAT_e我想使用预测函数和Naddl主持人绘制回归线:

但我得到这种类型的错误:

我想这是因为我没有指定 predict() 函数应该考虑哪个主持人。请注意,上面的函数适用于单变量模型,例如仅使用 MAT_e。

0 投票
1 回答
650 浏览

r - metafor 与 glmulti 相结合:如何排除不良交互

我正在拟合以下内容:

这段代码的目的是只包含这两个交互:“Naddl:MAT_e”和“Naddl:Fert”。因此,我exclude=c()用来从完整模型 ( level=2) 中过滤掉所有其他不需要的成对交互。

这在理论上应该与以下内容相同:

但是,当我添加exclude=c()公式时出现此错误:

我错过了什么exclude=c()吗?是否有更优雅的方式来指定“glmulti”中的交互项?

0 投票
2 回答
384 浏览

r - 避免来自 ifelse() 的 FALSE 部分的警告

我创建了一个用户函数来计算 Hedges 对g的 c() 偏差校正(Hedges,1981)。它直接基于metafor::.cmicalc(). library(metafor)它这样做:

当应用于包含值 <= 1 的向量时,会lgamma()生成警告,因为lgamma(0)(以及任何负值)会生成NaN. 所以,我的解决方案(以及做什么metafor::.cmicalc())包括一个ifelse()声明:

但是,这是我似乎找不到解决方案的问题warnings(),即使所有值都正确呈现为,它仍然会生成NA.

例子:

我知道(例如,从这个答案)即使条件是(并且相反,如果我改变条件和参数的顺序),第三个(FALSE)参数也会被评估......但我根本不知道如何解决这个问题(除了可能隐藏之前和之后的警告......)。ifelse()TRUE

(注意:我也试过dplyr::case_when()了,但这是完全相同的问题。)

0 投票
1 回答
225 浏览

r - 如何在森林图上隐藏 95% CI

我正在尝试一种方法来隐藏我们在森林图右侧看到的信息(ES(95% CI)。我正在使用 R 中的 metafor 包。

代码如下

在手册中找不到此选项。有人成功了吗?

0 投票
0 回答
72 浏览

r - 如何进一步划分亚组以进行荟萃分析中的简单效应分析?

我正在使用metafor在 R 中进行元分析。我学会了如何进行亚组分析并从网站上绘制结果,如下所示:

我的问题是,我可以进一步划分这三个子组吗?例如,以年划分以及分配:70s系统,70s随机,70s交替;80年代系统,80年代随机,80年代交替;90 年代系统,90 年代随机,90 年代交替。

而且,我希望这是在元回归中处理交互的方法。

0 投票
1 回答
847 浏览

r - R:Dplyr 'group_by' 似乎在 Shiny 中不起作用

我刚开始尝试使用Shiny,但让dplyr 'group_by' 函数工作时遇到了问题。

我的一般想法是,我有一个包含许多数值变量的数据集和一个可以对这些变量进行分组的变量(在我的例子中是 19 个国家,所以变量范围从 1 到 19)。现在我想使用metafor包对每个国家的两个变量之间的相关性进行荟萃分析。所以我想分别计算每个国家的两个选定变量之间的相关性,然后将它们放在一起进行荟萃分析,最后显示这些影响的森林图。

我以 mtcars 数据集为例说明我已经走了多远:

用户界面

服务器

在示例中,我按“齿轮”(有 3 个级别)进行分组,并选择了变量“排量”、“马力”、“后轴比”、“重量”和“1/4 英里时间”。

最终输出显示选择小部件以及森林图。但是,所有级别的“齿轮”的相关性都是相同的。因此,我怀疑“group_by”功能无法按预期工作。

我尝试按照网上的建议使用“group_by_”,但变化不大。

关于“group_by”如何与 Shiny 一起工作的任何想法?

0 投票
1 回答
345 浏览

r - 使用 rma.glmm 计算估计值:当 2 个估计值相同时出错

我正在使用 rma.glmm 来计算 2 个不同研究的元比例。例如,在 6 个月时出现 (xi) 和没有 (xii) 不良事件的个体比例:

但是,偶然在特定时间点,两个不同人群中的两个比例相同(11.1%),我收到一个错误:

rma.glmm (measure = "PLO", xi = xi, ni = ni, data = NNPNNP12monthAT) 中的错误:无法拟合 ML 模型。

我知道估计值将等于 11.1(因为这在两个群体中都是如此)......但我想要带有置信区间的输出,关于我可以做些什么来获取这些信息的任何建议?

0 投票
1 回答
994 浏览

r - 为什么文本没有添加到 R metafor-package 中的森林图?

我的数据:

我使用 metafor-package 进行荟萃分析。我想在我的森林图中添加一些文本作为情节图的“头部标签”。我已遵循 metafor-guide 但未添加文本。我不知道为什么。

我的森林剧情:

我的森林图如下所示:

没有文字的森林图

你能帮助我吗?在此先感谢,C.

0 投票
1 回答
794 浏览

r - 如何将 tau^2 添加到 metafor R 包中的 Forest Plot?

我正在使用 metafor 包在 R 中进行元分析。我按照http://www.metafor-project.org/doku.php/plots:forest_plot_with_subgroups上的教程进行操作

如何将 tau^2 添加到“所有研究的 RE 模型”?

在此先感谢,C.

0 投票
1 回答
154 浏览

r - 如何将研究名称添加到 R-metafor 中 L'abbé 图中的点?

是否可以将研究的名称添加到 L'abbé 图中代表研究的“点”?我正在使用 metafor 包。

研究名称可在

谢谢你,C。