在QQ情节中包括Lambda

Question

我想制作一个QQ情节，我达到了这个目标：

plot(-log(data$QQ, 10), -log(data$UNADJ,10), 
     xlab = "expected –logP values", 
     ylab = "observed –logP values"); 
abline(a = 0, b=1,col="red")

其中：

QQ
7.82E-06
1.82E-05
2.16E-05
2.43E-05

UNADJ
7.82E-06
1.82E-05
2.16E-05
2.43E-05

这有效，但我还需要lambda值，而我无法集成该功能。 有人可以帮我制作这个包含lambda输出的情节吗？

Answer 1

我不熟悉w /以下项目，

library(sos)  
???'inverse chi' 

但您可以在 LaplacesDemon 中尝试dist.Inverse.ChiSquare 或 metrumrg中的rinvchisq
或 geoR 中的InvChisquare。

也许其中一个会产生你需要的lambda。