Python 2.7 对logistic回归例子的理解
我不熟悉逻辑回归。下面是来自包mypc的示例 项目我还不清楚它的目的。更具体地说, 变量n是[5,5,5,5],在模式中使用:pymc.Binomial 我想在二项拟合中应该同时有0和1。n代表 “1”案例 你能解释一下这个例子的意思吗?谢谢 示例来自: www.map.ox.ac.uk/media/PDF/Patil_et_al_2010.PDFPython 2.7 对logistic回归例子的理解,python-2.7,Python 2.7,我不熟悉逻辑回归。下面是来自包mypc的示例 项目我还不清楚它的目的。更具体地说, 变量n是[5,5,5,5],在模式中使用:pymc.Binomial 我想在二项拟合中应该同时有0和1。n代表 “1”案例 你能解释一下这个例子的意思吗?谢谢 示例来自: www.map.ox.ac.uk/media/PDF/Patil_et_al_2010.PDF import pymc import numpy as np n = 5*np.ones(4,dtype=int) x = np.arra
import pymc
import numpy as np
n = 5*np.ones(4,dtype=int)
x = np.array([-.86,-.3,-.05,.73])
alpha = pymc.Normal('alpha',mu=0,tau=.01)
beta = pymc.Normal('beta',mu=0,tau=.01)
@pymc.deterministic
def theta(a=alpha, b=beta):
"""theta = logit^{-1}(a+b)"""
return pymc.invlogit(a+b*x)
d = pymc.Binomial('d', n=n, p=theta, value=np.array([0.,1.,3.,5.]),\
observed=True)
.......
import pymc
import pymc.Matplot
import mymodel
S = pymc.MCMC(mymodel, db='pickle')
S.sample(iter=10000, burn=5000, thin=2)
pymc.Matplot.plot(S)
import matplotlib.pyplot as plt
plt.show()