发布时间 : 星期一 文章乌龟蛋建模更新完毕开始阅读
Local optimal solution found.
Objective value: 0.3085469E+08 Extended solver steps: 3 Total solver iterations: 176
Variable Value Reduced Cost
X0 0.000000
Y0 0.000000
Z0 0.000000
N -581250.0
Row Slack or Price
1 1.000000
2 0.000000
3 390.6250
4 97656.25
5 0.000000
6 156250.0
7 0.000000
Logit
T=[27.2 27.7 28.3 28.4 29.9]'; Male=[2 17 26 19 27]'; Total=[27 24 30 27 28]'; Proport=Male./Toal;
[b,dev,stats]=glmfit(T,[Male Total],'binomial','logit'); logitFit=glmval(b,T,'logit');
30.00000 19.37500 50.00000 1.000000 0.3085469E+08 0.6250000 0.000000 0.000000 20.62500 0.000000 1.000000 Surplus Dual plot(T,Proport,'o',T,logitFit,'r-');
xlabel('T');ylabel('Proportion of MALE') b,bi=stats.se,dev b =
-61.3183 2.2110 bi =
12.0224 0.4309
dev =
14.8629
10.90.80.7Proportion of MALE0.60.50.40.30.20.102727.52828.5T2929.530[yhat,dylo,dyhi]=glmval(b,T,'logit',stats)
yhat =
0.2354 0.4818
0.7780 0.8138 0.9918
dylo =
0.1065 0.1133 0.1091 0.1094 0.0371
dyhi =
0.1551 0.1152 0.0808 0.0753
0.0068 y=[1.0496 0.8871 1.8721 0.8650 3.2934]'; x1=[27.2 27.7 28.3 28.4 29.9]'; x=[ones(5,1) x1];
[b,bint,r,rint,stats]=regress(y,x,0.05) b =
-23.4877 0.8863
bint =
-49.5757 2.6003 -0.0351 1.8076 r =
0.4310 -0.1746 0.2787 -0.8171 0.2819
rint =
-0.8864 1.7485 -2.0732 1.7240 -1.6769 2.2342 -1.7248 0.0906 -0.4631 1.0269
stats =
0.7575 9.3710 0.0549 0.3470 x 27.2 27.7 28.3 28.4 29.9 y 1.0496 0.8871 1.8721 0.8650 3.2934 27.2 27 2 0.7407 27.7 24 17 0.7083 28.3 30 26 0.8667 28.4 27 19 0.7037 29.9 28 27 0.9642