包括LAMBDA在QQ圖

我想作一個QQ情節和我到達這個地步：包括LAMBDA在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輸出？

一個小的可重複的例子呢？ – 2012-03-04 14:28:26

你想回答哪個問題：如何找到積分（或者可能是累積分佈函數？），或者如何在QQplot頂部疊加一個圖形？ – 2012-03-04 20:41:54

@ Carl Witthoft我的意思是第一個。據我所知，Lambda是與卡方統計相關的通貨膨脹因素。 λ是通過取數據的卡方統計量的分佈的中值，並將該中值除以相應的（理想的）卡方分佈的中值來計算的。 – user1236418 2012-03-04 21:43:19

我不熟悉瓦特/下列物品，與

library(sos) 
???'inverse chi'

發現，但你可以在在geoR嘗試dist.Inverse.ChiSquare LaplacesDemon 或rinvchisq在metrumrg
或InvChisquare。

也許其中一個會生成你需要的lambda。

2012-03-05 00:36:42

回答