getsigma: Get the best estimate of sigma conditional on s

In danjlawson/robustarchitecture: Estimate Genomic Architecture parameters using robust statistics

Description

Estimate sigma using the provided SNPs and a given s. This is done by ensuring that the whole distribution is scaled to obtain the theoretical proportion of SNPs above the second threshold, given the first, thereby ensuring that if s is correct then there should be no relationship between MAF and whether beta is over each threshold.

This uses a binary search of the possible thresholds that impact the choice, and is therefore logarithmic in the number of SNPs.

Usage

getsigma(ratio, thresholds)

Arguments

ratio	The ratio of beta/betahat(f,s,sigma=1)
thresholds	The architecture thresholds to be used in s.d. units.

Value

The estimate of sigma

danjlawson/robustarchitecture documentation built on Dec. 19, 2021, 8:09 p.m.