estlambda {GenABEL} | R Documentation |
Estimate the inflation factor for a distribution of P-values or 1df chi-square test. The major use of this procedure is the Genomic Control, but can also be used to visualise the distribution of P-values coming from other tests.
estlambda(data, plot = TRUE, proportion = 1.0, ... )
data |
A vector of reals. If all are <=1, it is assumed that this is a vector of P-values, else it is treated as a vector of chi-squares with 1 d.f. |
plot |
Wether the prot should be presented |
proportion |
The proportion of lowest P (Chi2) to be used when estimating the inflation factor Lambda |
... |
arguments passed to plot function |
A list with elements
estimate |
Estimate of Lambda |
se |
Standard error of the estimate |
Yurii Aulchenko
data(srdta) pex <- summary(srdta@gtdata)[,"Pexact"] estlambda(pex) a <- ccfast("bt",srdta) a$lambda