def myloglike(cube, ndim, nparams):
efac = np.ones(npsr)
equad = np.zeros(npsr)
A = 10 ** cube[0]
gam = cube[1]
loglike = PALLikelihoods.modelIndependentFullPTAPL(psr, F, Diag, f, A, gam, Ared, gred, efac, equad, ORF)
# print efac, rho, loglike
return loglike
def myloglike(cube, ndim, nparams):
efac = np.ones(npsr)
equad = np.zeros(npsr)
A = 10**cube[0]
gam = cube[1]
loglike = PALLikelihoods.modelIndependentFullPTAPL(psr, F, Diag, f, A, gam, \
Ared, gred, efac, equad, ORF)
#print efac, rho, loglike
return loglike
def loglike(x):
A = x[0] * 1e-14
gam = x[1]
efac = x[2:(2 + npsr)]
equad = 10**x[(2 + npsr):(2 + 2 * npsr)]
Ared = 10**x[(2 + 2 * npsr):(2 + 3 * npsr)]
gred = x[(2 + 3 * npsr):(2 + 4 * npsr)]
loglike = PALLikelihoods.modelIndependentFullPTAPL(psr, F, Diag, f, A, gam, \
Ared, gred, efac, equad, ORF)
return loglike
def loglike(x):
A = x[0]*1e-14
gam = x[1]
efac = x[2:(2+npsr)]
equad = 10**x[(2+npsr):(2+2*npsr)]
Ared = 10**x[(2+2*npsr):(2+3*npsr)]
gred = x[(2+3*npsr):(2+4*npsr)]
loglike = PALLikelihoods.modelIndependentFullPTAPL(psr, F, Diag, f, A, gam, \
Ared, gred, efac, equad, ORF)
return loglike
def myloglike(cube, ndim, nparams):
efac = np.zeros(npsr)
equad = np.zeros(npsr)
Ared = np.zeros(npsr)
gred = np.zeros(npsr)
A = cube[0]
gam = cube[1]
for ii in range(npsr):
efac[ii] = cube[ii + 2]
equad[ii] = 10 ** cube[ii + 2 + npsr]
Ared[ii] = 10 ** cube[ii + 2 + 2 * npsr]
gred[ii] = cube[ii + 2 + 3 * npsr]
loglike = PALLikelihoods.modelIndependentFullPTAPL(psr, F, Diag, f, A, gam, Ared, gred, efac, equad, ORF)
# print efac, rho, loglike
return loglike
def myloglike(cube, ndim, nparams):
efac = np.zeros(npsr)
equad = np.zeros(npsr)
Ared = np.zeros(npsr)
gred = np.zeros(npsr)
A = cube[0]
gam = cube[1]
for ii in range(npsr):
efac[ii] = cube[ii + 2]
equad[ii] = 10**cube[ii + 2 + npsr]
Ared[ii] = 10**cube[ii + 2 + 2 * npsr]
gred[ii] = cube[ii + 2 + 3 * npsr]
loglike = PALLikelihoods.modelIndependentFullPTAPL(psr, F, Diag, f, A, gam, \
Ared, gred, efac, equad, ORF)
#print efac, rho, loglike
return loglike
