est_statistic: estimate Z score at a single SNP

Description Usage Arguments Details Value Author(s)

View source: R/est_zscore.R

Description

Wrapper function to run est_zscore for all snps in snps

Usage

1
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 April 4, 2018, 12:41 a.m.