sdY.est: Estimate trait variance, internal function

In GRealesM/RapidoPGS: A Fast and Light Package to Compute Polygenic Risk Scores

Estimate trait variance, internal function

Description

Estimate trait standard deviation given vectors of variance of coefficients, MAF and sample size

Usage

sdY.est(vbeta, maf, n)

Arguments

vbeta	vector of variance of coefficients
maf	vector of MAF (same length as vbeta)
n	sample size

Details

Estimate is based on var(beta-hat) = var(Y) / (n * var(X)) var(X) = 2* maf*(1-maf) so we can estimate var(Y) by regressing n*var(X) against 1/var(beta) This function is verbatim from its namesake in coloc package (github.com/chr1swallace/coloc/), by Chris Wallace

Value

estimated standard deviation of Y

Author(s)

Chris Wallace

GRealesM/RapidoPGS documentation built on Oct. 15, 2023, 2:43 p.m.