问题描述:对于定量因变量模型logistic:ln(p/1-p)=c+b(1)*x(1)+b(2)*x(2)+e
自变量是时间序列的,因变量取0或1。
刚我查了很久,用eviews可以直接估计,它接受的因变量是0或者1.
用matlab的话,其估计接受的是好像是概率值p,如果是的话如何给p赋值。
网上的一个例子:
建立m函数文件存为logistic1
function f=logistic1(b)
t=[0,5,10,24,33,48,57,72,96,120,144,168,192,216];y=[0,0.028,0.103,0.336,0.450,0.597,0.716,0.778,0.835,0.849,0.816,0.839,0.811,0.816];
f = y-b(1)./(1+b(2).*exp(-b(3).*t));
b0=[10,2,2];
>> b=leastsq('logistic1',b0)
b =
0.8221 13.9173 0.0818
这里y就是概率值,而我的数据里没有。真心求教高手该如何用matlab解决?