漏斗 plot，有效样本量为 R

Question

我打算创建一个具有有效样本量的漏斗 plot。 我正在使用 metafor package 的funnel.default() metafor我尝试了以下代码来创建这样一个漏斗 plot

Soil_mineral_nitrogen$inv_n_tilda <-  with(Soil_mineral_nitrogen, (control_mean + treatment_mean) / (control_mean*treatment_mean))


par(mfrow = c(1, 2))
funnel(Soil_mineral_nitrogen$lnrr, Soil_mineral_nitrogen$inv_n_tilda, yaxis="ninv",
       #xlim = c(-3, 3),
       ylab = "Effective sample size (ñ)",
       xlab = "Effect size (RR)", col =Soil_mineral_nitrogen$unique_id, atransf = exp)

但是这段代码给我返回了一个错误，如下所示

Error in funnel.default(Soil_mineral_nitrogen$lnrr, Soil_mineral_nitrogen$inv_n_tilda,  : 
  No sample size information available.

有谁知道如何处理这个错误？

Answer 1

它使用以下代码对我有用

par(mfrow = c(1, 2))
funnel(Soil_mineral_nitrogen$lnrr, Soil_mineral_nitrogen$v, ni = Soil_mineral_nitrogen$inv_n_tilda, 
       yaxis="ni",
       #xlim = c(-3, 3),
       ylab = "Effective sample size (ñ)",
       xlab = "Effect size (RR)", col =Soil_mineral_nitrogen$unique_id, atransf = exp)