广义线性模型matlab,广义线性回归分析的MATLAB函数

>> z = [9.73875 20.75 36.59875

13.5725 29.6325 50.93875

18.97875 36.59875 80.13875

20.75125 38.22125 90.925

22.055 44.58 104.7725];

>> [p,t] = meshgrid([0.8 1 1.2],[60:60:300]);

>> stats = reglm(z(:),[p(:), t(:)],'quadratic',{'p','t'});

>> [pnew,tnew] = meshgrid(linspace(0.8,1.2,20),linspace(60,300,20));

>> pp = pnew(:);

>> tt = tnew(:);

>> zhat = [ones(400,1) pp tt pp.*tt pp.^2 tt.^2]*stats.beta;

>> mesh(pnew,tnew,reshape(zhat,[20,20]));

>> hold on

>> plot3(p(:),t(:),z(:),'k*')

拟合结果：

------------------------------------方差分析表------------------------------------

方差来源    自由度            平方和             均方             F值          p值

回归       5.0000       11548.9147        2309.782