scripts for heteroskedastic tests
Switch branches/tags
Nothing to show
Clone or download
Fetching latest commit…
Cannot retrieve the latest commit at this time.
Type Name Latest commit message Commit time
Failed to load latest commit information.


working scripts for heteroskedastic tests
last updated: July, 2015

simple example

from simulations import *
from BTHsuite import *

[maf, n, intercept, beta, sigma2, logalpha,T,c] = [0.2, 300,0, 0.1, 1., 0.2, 1, 0 ]
[genotype, phenotype] = sim_xy(maf, n, intercept, beta, sigma2, logalpha,T,c)

!# compute Bayes Factors using BTH for the generated data
logBF = BTH(genotype,phenotype)
!# compute Bayes Factors for the shuffled data
[genotypePermuted, phenotypePermuted] = special_shuffle2(genotype,phenotype)
logBFperm = BTH(genotypePermuted, phenotypePermuted)