включить лямбду в сюжет QQ

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

это работает, но мне также нужно лямбда-значение, и мне не удается интегрировать функцию. Может ли кто-нибудь помочь мне сделать этот сюжет, включая вывод с помощью лямбды?

Answer 1

Я не знаком со следующими предметами, найденными с

library(sos)  
???'inverse chi' 

, но вы можете попробовать dist.Inverse.ChiSquare в LaplacesDemon или rinvchisq в metrumrg
или InvChisquare в geoR .

Возможно, один из них сгенерирует лямбду, которая вам нужна.