def unit_test_basic(): print spbetaln(10, 20) print spbetainc(0.5, 10, 3) print spbetainc(0.5, 10, 3) print spbetainc(10, 3, 0.5) print spbtdtri(2,3,0.3) print fq(0.3, 2, 3) print tq(0.3, 10)
def fq(p, v1, v2): """ Fisher's F distribution Quantiles\n Returns q """ #q = special.spbtdtri(v1/2.0, v2/2.0, p) #q = betaq(p,v1/2,v2/2) q = spbtdtri(v1/2.0,v2/2.0,p) #if(isarray(q) == False): # print 'array' # for i in range(len(q)): # q[i] = v2/v1 * q[i]/(1-q[i]) #else: # print 'not array' q = float(v2)/float(v1) * float(q / (1-q)) #http://docs.scipy.org/doc/scipy/reference/generated/scipy.special.spbtdtri.html #pth quantile of the beta distribution #q = v2/v1 * q/1-q return q