pcadapt (version 1.0)

sdevList: Auxiliary function

Description

sdevList returns a list of estimated standard deviations using the kurtosis-based method.

Usage

sdevList(V, K, maf, minmaf)

Arguments

V

a numeric matrix containing the loadings of the Principal Components Analysis.

K

an integer specifying the number of retained principal components.

Examples

x <- replicate(1000, rnorm(100))
x <- floor(abs(x))
freq <- apply(x,2,sum)/(2*dim(x)[1])
minfreq <- pmin(freq,1-freq)
respc <- prcomp(x)
V <- respc$rotation*sqrt(dim(x)[2])
sdevList(V,K=3,maf=minfreq,minmaf=0.05)