GeneticsDesign (version 1.40.0)

GeneticPower.Quantitative.Numeric: Power of Genetics Study

Description

Compute power of quantitative genetics studies, when the genotype is handled as a numeric value (0,1,2) GeneticPower.Quantitative.Numeric or as a factor GeneticPower.Quantitative.Factor.

Usage

GeneticPower.Quantitative.Numeric( N=1000, delta=1, freq=0.15, minh=c("additive", "dominant", "recessive"), sigma=1, OtherParms=0, alpha=0.05, numtests=1, moi=NULL, rsquared=NULL)
GeneticPower.Quantitative.Factor(  N=1000, delta=1, freq=0.15, minh=c("additive","dominant","recessive"), sigma=1, OtherParms=0, alpha=0.05, numtests=1, moi=NULL,  rsquared=NULL)

Arguments

N

total samples in the analysis

delta

Treatment effect for an individual homozygote for the disease allele ('b') relative to an individual homozygote for the reference allele ('A')

freq

allele frequency of disease allele 'b'

minh

mode of inheritance: "additive","dominant","recessive", Default is "additive". This parameter is OVER-RIDDEN by moi.

sigma

standard deviation of the response phenotype

OtherParms

number of additional parameters (really, DOF) in the model that will reduce your overall DOF

alpha

desired significance level

numtests

number of tests to be corrected by Bonferroni adjustment beforee achieving 'alpha'

moi

continuous value between 0 and 1 (inclusive) specifying the mode of inheritance: 0 for recessive, 0.5 for additive, 1.0 for dominant. This parameter OVER-RIDES minh.

rsquared

fraction of total sum-of-squares explained by fit. This parameter OVER-RIDES delta AND sigma.

Details

The value of moi overrides any value specified for minh. Specifying a minh="recessive" is equivalent to specifying moi=0, minh="additive" is equivalent to moi=0.5, and minh="dominant" is equivalent to moi=1.0.

Examples

GeneticPower.Quantitative.Numeric(
                                  N=50,
                                  freq=0.1,
                                  minh="recessive",
                                  alpha=0.05
                                  )

GeneticPower.Quantitative.Factor(
                                  N=50,
                                  freq=0.1,
                                  minh="recessive",
                                  alpha=0.05
                                 )


##