Mc@szddklZlZlZdZdZdZdZdZdZ dZ d Z d Z d Z d S( i(tlogtexptfloorcCsddddddg}t|d}|d}d}d }||d t|8}x0td D]"}|d7}||||7}qcW| td |S( NgSv S@g+WUg8@g뒦~4gM&S?g1!e[־g?g@ig?ig% @(tfloatRtrange(tZtcoftxttemptsertj((sK/woldlab/castor/data00/home/georgi/erange-4.0a/cistematic/cisstat/helper.pytgammaln%s   cCsddddddddg}d }t|d }x0td D]"}|d 7}||||7}qAWt|d t|t|d tt|d}|S(Ng ö)$@g&Ԏg ta@gQfg>@~)@gA~g~\>giE>g ?g?ig*hS@g?g@(RRR(RRR RR tresult((sK/woldlab/castor/data00/home/georgi/erange-4.0a/cistematic/cisstat/helper.pytgammaln25s   <cCsttt|dS(Ni(troundRR (R((sK/woldlab/castor/data00/home/georgi/erange-4.0a/cistematic/cisstat/helper.pytgammaEscCsttt|dS(Ni(RRR (R((sK/woldlab/castor/data00/home/georgi/erange-4.0a/cistematic/cisstat/helper.pytgamma2IscCsd}dddddg}|djo dGHdS|d jot|d Sx7||jo)|}|d 7}|i|||qQW||SdS( Nig?g@g@g8@is5error: factorial only works for non-negative integersii i(Rtappend(tNtntoptarrR ((sK/woldlab/castor/data00/home/georgi/erange-4.0a/cistematic/cisstat/helper.pyt factorialMs    cCs;|djo dGHdS|djodSt|dSdS(Nis5error: factorial only works for non-negative integersiigg?(R (R((sK/woldlab/castor/data00/home/georgi/erange-4.0a/cistematic/cisstat/helper.pytfactln^s   cCs"t|t|t||S(N(R(tntk((sK/woldlab/castor/data00/home/georgi/erange-4.0a/cistematic/cisstat/helper.pytbicolnhscCstdtt||S(Ng?(RRR(RR((sK/woldlab/castor/data00/home/georgi/erange-4.0a/cistematic/cisstat/helper.pytbicolscCs(tt|t|t||S(N(RR (tztw((sK/woldlab/castor/data00/home/georgi/erange-4.0a/cistematic/cisstat/helper.pytbetapscCset|}t|}t|}t|}tt||t||||t||S(s hyperGeometric(N, r, n, x): N is population size r is the number in the population identified as a success n is the sample size x is the number in the sample identified as a success (RRR(RtrRR((sK/woldlab/castor/data00/home/georgi/erange-4.0a/cistematic/cisstat/helper.pythyperGeometricts     N(tmathRRRR R RRRRRRRR(((sK/woldlab/castor/data00/home/georgi/erange-4.0a/cistematic/cisstat/helper.pyt#s