est_statistic: estimate Z score at a single SNP

In mdfortune/simGWAS: Simulating the expected output from a GWAS with a given causal model

Description

Wrapper function to run est_zscore for all snps in snps

Usage

est_statistic(N0, N1, snps, W, gamma1, freq, GenoProbList)

Arguments

N0	The number of Y=0
N1	The number of Y=1
snps	The snps at which we wish to compute the expected Z Score
W	The true causal SNPs (these need not be in "snps")
gamma1	The log odds ratios of effect of the true causal SNPs (not including gamma0, the intercept term)
freq	Frequencies of SNP appearances (computed using snphap)
GenoProbList	An list of objects giving the probability of seeing each X,W genotype vector

Details

Assumes we have a list, GenoProbList, giving the GenoProb values for each X.

Value

The expected Z Score for SNP X, assuming the causal SNPs are W

Author(s)

Mary Fortune and Chris Wallace

mdfortune/simGWAS documentation built on May 22, 2019, 3:25 p.m.