Compute matrix of simulated Z scores

simulated_z_score(N0, N1, snps, W, gamma.W, freq,
GenoProbList = make_GenoProbList(snps = snps, W = W, freq = freq),
nrep = 1)

## Arguments

N0 The number of Y=0 The number of Y=1 The snps at which we wish to compute the expected Z Score The true causal SNPs (these need not be in "snps") The log odds ratios of effect of the true causal SNPs (not including gamma0, the intercept term) Haplotype frequencies as a data.frame, with column Probability indicating relative frequency in controls. An list of objects giving the probability of seeing each X,W genotype vector. This can be calculated within the function if no value supplied, or you can pass a pre-calculated version Number of replicates (simulated vectors of Z scores) under this scenario. Default=1

## Examples

freq=fake_freq(nhaps=100,nsnps=5) # fake haplotype frequency data
Z=simulated_z_score(N0=1000,N1=2000,snps=paste0("s",1:5),
W="s1",gamma.W=log(1.5),freq=freq,nrep=3)
Z # causal variant is SNP 1, with OR 1.5#>          [,1]     [,2]     [,3]       [,4]       [,5]
#> [1,] 6.881099 4.001179 3.929409  1.4730519  0.9434700
#> [2,] 5.757769 3.142459 1.585591 -0.4055338 -0.1870463
#> [3,] 5.756526 3.115909 2.273784  0.3630614 -0.4446704