forked from omerwe/LEAP
-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathleapMain.py
25 lines (17 loc) · 1.17 KB
/
leapMain.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
def eigenDecompose(bed, kinshipFile=None, outFile=None, ignore_neig=False):
from . import eigenDecompose
return eigenDecompose.eigenDecompose(bed, kinshipFile, outFile, ignore_neig)
def findRelated(bed, outFile=None, cutoff=0.05, kinshipFile=None):
from . import findRelated
return findRelated.findRelated(bed, outFile, cutoff, kinshipFile)
def calcH2(pheno, prev, eigen, keepArr=None, covar=None, numRemovePCs=10, lowtail=False):
from . import calc_h2
return calc_h2.calc_h2(pheno, prev, eigen, keepArr, covar, numRemovePCs, lowtail)
def probit(bed, pheno, h2, prev, eigen, outFile=None, keepArr=None, covar=None, thresholds=None, nofail=False,
numSkipTopPCs=0, mineig=1e-3, hess=False, recenter=True, maxFixedIters=100, epsilon=1e-3, treatFixedAsRandom=False):
from . import probit
return probit.probit(bed, pheno, h2, prev, eigen, outFile, keepArr, covar, thresholds, nofail,
numSkipTopPCs, mineig, hess, recenter, maxFixedIters, epsilon, treatFixedAsRandom=treatFixedAsRandom)
def leapGwas(bedSim, bedTest, pheno, h2, outFile=None, eigenFile=None, covar=None):
from . import leap_gwas
return leap_gwas.gwas(bedSim, bedTest, pheno, h2, outFile, eigenFile, covar)