R语言metafor包元分析fsn函数输出结果不匹配问题咨询

1. 「观测显著性水平」与元分析模型显著性不匹配的原因 #

• fsn()函数默认的计算逻辑是基于固定效应模型的逆方差加权合并，采用Z检验计算合并效应的显著性，不会自动继承你用rma()拟合的REML随机效应模型的参数。
• 你调用函数时传入的是单独的效应量a$yi和方差a$vi，fsn()会独立计算固定效应下的合并结果，和你拟合的随机效应模型的效应量、标准误、显著性结果天然存在差异，因此出现不匹配是正常现象。

2. 模型本身不显著但fsn返回大于0的原因 #

• 失安全系数N的返回值完全基于fsn()自身计算的合并显著性结果，与你外部拟合的元分析模型结果无关：
  • 只有当fsn()自己计算的合并效应本身就不显著时，才会返回0，代表不需要额外的零效应研究就已经满足不显著的条件
  • 你提供的示例中，fsn()用固定效应计算得到的合并p值刚好略低于0.05的显著性阈值，因此只需要1个零效应的未发表研究就可以将合并结果拉到不显著，所以返回值为1

补充说明 #

如果需要fsn()的计算与你拟合的随机效应模型保持对齐，可以直接传入rma拟合对象而非单独的yi和vi，即使用fsn(a)调用。但需要注意：经典的失安全系数算法本身是为固定效应元分析设计的，随机效应场景下的失安全系数计算方法目前学界仍存在争议，解读结果时需要谨慎。

内容的提问来源于stack exchange，提问作者user3306720