def nu2phr(nu): phr = np.sqrt(2.0 / nu) * spesh.gamma((nu + 1.0) / 2.0) / \ spesh.gamma(nu / 2.0) phr = sps.t.pdf(0.0, nu) / sps.norm.pdf(0.0) return phr # read in data