Estimate the inflation factor for a distribution of P-values

estlambda {GenABEL}R Documentation

Estimate the inflation factor for a distribution of P-values

Description

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.

Usage

estlambda(data, plot = TRUE, proportion = 1.0, ... )

Arguments

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

Value

A list with elements

estimate Estimate of Lambda
se Standard error of the estimate

Author(s)

Yurii Aulchenko

See Also

ccfast, qtscore

Examples

data(srdta)
pex <- summary(srdta@gtdata)[,"Pexact"]
estlambda(pex)
a <- ccfast("bt",srdta)
lambda(a)

[Package GenABEL version 1.6-7 Index]