R语言metap包绘制Albatross图：负值不显示等问题求助

核心问题根源 #

1. 负值不显示+NaN警告：metap::albatros要求输入的pval是0-1区间内的单侧p值，负值仅用于标记负相关方向，但所有值的绝对值必须在(0,1)之间。你的输入中虽然负值绝对值合法，但axes$xlimit = 0.001是错误参数（需为长度2的向量），导致函数x轴转换逻辑出错，生成NaN并抑制负值点显示。
2. 正值分配到负相关区域：错误的xlimit单值设置破坏了坐标轴映射规则，导致点的区域分配逻辑混乱；同时plotpars$yscale = 1000是无效参数（仅接受"linear"/"sqrt"/"log"），触发的默认缩放警告也会干扰绘图逻辑。

修正后的代码 #

library(metap) 

# 确保pval所有值的绝对值在0-1之间（当前输入已满足）
pval <- c(0.03, 0.7, 0.51, 0.22, 0.02, -0.819, -0.916, 0.22, 0.4, 0.78,
          0.2, 0.46, 0.001, 0.1, 0.002, 0.41, -0.18, 0.001, 0.03, 0.03,
          0.86, 0.84, 0.347, 0.74, 0.08, 0.98, 0.39, 0.43, 0.02, 0.69)

# 修正参数错误，需确保ptrend$`No.Cases`长度与pval一致
fit.m <- albatros(pval, ptrend$`No.Cases`, 
                  axes = list(xlimit = c(0, 1), # xlimit改为长度2的边界向量
                              ylimit = c(1,600000), 
                              lefttext = "Negative correlation", 
                              righttext = "Positive correlation"), 
                  contours = list(type = "corr", 
                                  contvals = c(0.05, 0.009, 0.005), 
                                  ltys = 1:3), 
                  plotpars = list(pchs = 23, cols = 'blue'), # 移除无效的yscale设置
                  main = "Validity")

关键修正说明 #

• 负值显示恢复：修正xlimit参数后，函数能正确解析负值的方向标记，负值点会显示在左侧负相关区域。
• 区域分配修复：合法的xlimit向量会恢复正常的坐标轴映射，正值点将全部显示在右侧正相关区域。
• 警告消除：移除无效的yscale参数后，所有默认缩放警告和NaN警告都会消失。

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