专业资料word完美格式matlabmatlabmatlab1.(ptzz=f(p,t)zp0.811.2t609.7387520.7536.598712013.572529.632550.9387518018.9787536.5987580.13875240207512538.2212590.92530022.05544.58104.7725regstatsreglmMfunctionstats=reglm(y,:X,model,varnames)%%%reglm(y,X)%Xnp.yn1%%stats=reglm(y,X)stats.%%stats=reglm(y,X,model)model.model%%'linear'%'interaction'%'quadratic'%'purequadratic'%%stats=reglm(y,X,model,varnames)varnames.varnames%word完美格式%X.X1,X2,.%%:%x=[215250180250180215180215250215215%136.5136.5136.5138.5139.5138.5140.5140.5140.5138.5138.5%y=[6.27.54.85.14.64.62.83.14.34.94.1]，;%reglm(y,x,，quadratic，)%%%Fp%5.000015.02773.00557.61220.0219%5.00001.97420.3948%10.000017.0018%%(RootMSE)0.6284(R-Square)0.8839%(DependentMean)4.7273(AdjR-Sq)0.7678%%%tp%30.94282011.11170.01540.9883%X10.70400.64051.09920.3218%X2-0.848729.1537-0.02910.9779%X1*X2-0.00580.0044-1.31320.2461%X1*X10.00030.00030.83840.4400%X2*X20.00520.10550.04920.9626%%Copyright2009—2010xiezhh.%$Revision:1.0.0.0$$Date:2009/12/2221:41:00$ifnargin<2error(，');endp=size(X,2);%Xifnargin<3||isempty(model)model='linear，;%modelend%varnamesifnargin<4||isempty(varnames)(AdjR-Sq)',...word完美格式varnamel=strcat({'X'},num2str([l:p]'));varnames=makevarnames(varname1,model);%elseifischar(varnames)varname1=cellstr(varnames);elseifiscell(varnames)varname1=varnames(:);elseerror('varnamesendifsize(varname1,1)~=perror('X');elsevarnames=makevarnames(varname1,model);%endendST=regstats(y,X,model);%regstatsSTf=ST.fstat;%Ft=ST.tstat;%t%fprintf('\n');fprintf('');fprintf('\n');邙rintf('%s%7s%15s%15s%15s%12s','','','');fprintf('\n');fmt='%s%13.4f%17.4f%17.4f%16.4f%12.4f';fprintf(fmt,'',f.dfr,f.ssr,f.ssr/f.dfr,f.f,f.pval);fprintf('\n');fmt='%s%13.4f%17.4f%17.4f';fprintf(fmt,'',f.dfe,f.sse,f.sse/f.dfe);fprintf('\n');fmt='%s%13.4f%17.4f';fprintf(fmt,'',f.dfe+f.dfr,f.sse+f.ssr);fprintf('\n');fprintf('\n');%fmt='%22s%15.4f%25s%10.4f';邙rintf(fmt,'(RootMSE)',sqrt(ST.mse),'(R-Square)',ST.rsquare);fprintf('\n');fprintf(fmt,'(DependentMean)',mean(y),'ST.adjrsquare);word完美格式fprintf('\n');fprintf('\n');%tfprintf('')；fprintf('\n');邙rintf('%8s%18s%15s%15s%12s','','','','t','p');fprintf('\n');fori=l:size(t.beta,l)ifi==1fmt='%8s%20.4f%17.4f%17.4f%12.4f\n';fprintf(fmt,'',t.beta(i),t.se(i),t.t(i),t.pval(i));elsefmt='%10s%20.4f%17.4f%17.4f%12.4f\n';fprintf(fmt,varnames{i-l},t.beta(i),t.se(i),t.t(i),t.pval(i));endendifnargout==1stats=ST;%end%functionvarnames=makevarnames(varnamel,model)%p=size(varnamel,l);varname2=[];fori=1:p-1varname2=[varname2;strcat(varnamel(i),'*',varnamel(i+1:end))];endvarname3=strcat(varnamel,'*',varnamel);switchmodelcase'linear'varnames=varnamel;case'interaction'varnames=[varname1;varname2];case'quadratic'varnames=[varname1;varname2;varname3];case'purequadratic'varnames=[varname1;varname3];endreglm>>z=[9.7387520.7536.5987513.572529.632550.9387518.9787536.5987580.1387520.7512538.2212590.92522.05544.58104.7725];word完美格式>>[p,t]=meshgrid([0.811.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)ppttpp.*ttpp."2tt."2]*stats.beta;>>mesh(pnew,tnew,reshape(zhat,[20,20]));>>holdon>>plot3(p(:),t(:),z(:),，k*，)Fp5.000011548.91472309.782993.47390.00009.0000222.394224.710514.000011771.3089(RootMSE)4.9710(R-Squar...
