请各位MATLAB高手,不吝赐教,想用matlab编写logit模型的效用函数,然后用极大似然估计求解效用函数的系数值。其中,mydata都是输入数据,X为想求解的各个变量的系数,我实在是头疼,不懂啊!
mydata=xlsread('D:\wo1.xls');
a=mydata(:,6);
b=mydata(:,7);
c=mydata(:,8);
d=mydata(:,9);
e=mydata(:,10);
f=mydata(:,11);
g=mydata(:,12);
h=mydata(:,13);
X0=[1,1,1,1,1,1,1,1,1,1,1]
options=optimset('largescale','off','display','iter','tolx',1e-5,'tolfun',1e-5);
[X,fval,exitflag,output,grad,hessian]=fminunc('fun1',X0,options);
function fun1=myfun(x)
f= x(1)+x(2)*a+x(3)*b+x(4)*c+x(5)*d+x(6)*e+x(7)*(ln(x(8))+exp(x(9)*f)+exp(x(10))*g+exp(x(11))*h);
fun1=-sum(log(f));


雷达卡


京公网安备 11010802022788号







