Question

我是Logistic回归的新手。以下是来自包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.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()

了解逻辑回归示例

0 个答案: