求大神帮忙检查下面程序哪里有问题——(SAS响应面设计试验之图形分析) 20
程序dataGasValue;X2=9;doX1=4to12by0.1;doX3=0.2to5by0.1;Y1=-40.1388+7.3477*X1+9.6052*X2-...
程序
data Gas Value;
X2=9;
do X1=4 to 12 by 0.1;do X3=0.2 to 5 by 0.1;
Y1=-40.1388+7.3477*X1+9.6052*X2-6.4517*X3-0.5777*X1*X2-0.5102*X1*X3+0.8929*X2*X3+0.0016*X1**2-0.4080*X2**2+0.6241*X3**2;
output;end;end;
run;
proc g3grid data=Gas Value out=gridfile;
grid X3*X1=Y1/naxis 1=20 naxis 2=20;
run;
goptions reset=all ftext=swiss htext=1.85 colors=(black r b black r b black r b black r b);
axis 1 label=(f='宋体' c=black h=1.85 c=blue width=1.5;
axis 2 label=(A=90f='宋体' c=black h=1.85 '配料浓度(%)')c=blue width=1.5;
symbol height=1.5 width=2.0;
proc g3d data=gridfile;
plot X3*X1=Y1/caxis=blue ctext=black cbottom=red ctop=black tilt=70 rotate=70 xticknum=7 yticknum=7 zticknum=7 gird;
run;
proc gcontour data=Gas Value;
plot X3*X1=Y1/levels=23 to 37 by 2 autolable nolegend haxis=axis 1 vaxis=axis 2 caxis=b grid;
如下图,一张是书上给出的程序,程序没问题会出来另一张图上的两个图,可现在不知道哪里出了问题了 展开
data Gas Value;
X2=9;
do X1=4 to 12 by 0.1;do X3=0.2 to 5 by 0.1;
Y1=-40.1388+7.3477*X1+9.6052*X2-6.4517*X3-0.5777*X1*X2-0.5102*X1*X3+0.8929*X2*X3+0.0016*X1**2-0.4080*X2**2+0.6241*X3**2;
output;end;end;
run;
proc g3grid data=Gas Value out=gridfile;
grid X3*X1=Y1/naxis 1=20 naxis 2=20;
run;
goptions reset=all ftext=swiss htext=1.85 colors=(black r b black r b black r b black r b);
axis 1 label=(f='宋体' c=black h=1.85 c=blue width=1.5;
axis 2 label=(A=90f='宋体' c=black h=1.85 '配料浓度(%)')c=blue width=1.5;
symbol height=1.5 width=2.0;
proc g3d data=gridfile;
plot X3*X1=Y1/caxis=blue ctext=black cbottom=red ctop=black tilt=70 rotate=70 xticknum=7 yticknum=7 zticknum=7 gird;
run;
proc gcontour data=Gas Value;
plot X3*X1=Y1/levels=23 to 37 by 2 autolable nolegend haxis=axis 1 vaxis=axis 2 caxis=b grid;
如下图,一张是书上给出的程序,程序没问题会出来另一张图上的两个图,可现在不知道哪里出了问题了 展开
2个回答
推荐律师服务:
若未解决您的问题,请您详细描述您的问题,通过百度律临进行免费专业咨询