loop_calc_pvalue: loop_calc_pvalue()

In gwrpvr: Genome-Wide Regression P-Value (Gwrpv)

Description

calls calc_pvalue()

Usage

loop_calc_pvalue(lowone, highone, lowtwo, hightwo, lowthree, highthree, n0a,
  n1a, n2a, n0, n1, n2, pa, pb, x, mua, mub, sumsqx, siga, sigb, vary, beta,
  skipiter, pvalue)

Arguments

lowone	lower bound outer loop
highone	upper bound outer loop
lowtwo	lower bound middle loop
hightwo	upper bound middle loop
lowthree	lower bound inner loop
highthree	upper bound inner loop
n0a	outer loop index
n1a	middle loop index
n2a	inner loop index
n0	the major allele homozygotes
n1	the major allele heterozygotes
n2	the minor allele zygotes
pa	parameter of the mixture distribution, a real number between zero and one with pa+pb=1
pb	parameter of the mixture distribution, a real number between zero and one with pa+pb=1
x	a zero mean explanatory variable from the SNP data set
mua	parameter of the mixture distribution, can be any real number
mub	parameter of the mixture distribution, can be any real number
sumsqx	sum of the squares of x
siga	parameter of the mixture distribution, can be any real number
sigb	parameter of the mixture distribution, can be any real number
vary	vary <- pa*(mua^2+siga^2)+pb*(mub^2+sigb^2)-(pa*mua+pb*mub)^2
beta	the beta from the regression being tested
skipiter	flag to determine if we can skip some calculations
pvalue	the input pvalue prior to calculating new improved pvalue

Value

pvalue

gwrpvr documentation built on May 2, 2019, 7:59 a.m.