Alloy 97-1-1-1 Alloy AuCa
Dentist Method 1500°F 1600°F 1700°F 1500°F 1600°F 1700°F
1 1 813 792 792 907 792 835
2 782 698 665 1115 835 870
3 752 620 835 847 560 585
2 1 715 803 813 858 907 882
2 772 782 743 933 792 824
3 835 715 673 698 734 681
3 1 743 627 752 858 762 724
2 813 743 613 824 847 782
3 743 681 743 715 824 681
4 1 792 743 762 894 792 649
2 690 882 772 813 870 858
3 493 707 289 715 813 312
5 1 707 698 715 772 1048 870
2 803 665 752 824 933 835
3 421 483 405 536 405 312
这是我对上述数据的 sas 代码:
data gold;
do dentist=1, 2, 3, 4, 5;
do method=1, 2, 3;
do alloy= 1,2;
do temp=1500, 1600, 1700;
input y @@; output;
end;
end;
end;
end;
cards;
813 792 792 907 792 835
782 698 665 1115 835 870
752 620 835 847 560 585
715 803 813 858 907 882
772 782 743 933 792 824
835 715 673 698 734 681
743 627 752 858 762 724
813 743 613 824 847 782
743 681 743 715 824 681
792 743 762 894 792 649
690 882 772 813 870 858
493 707 289 715 813 312
707 698 715 772 1048 870
803 665 752 824 933 835
421 483 405 536 405 312
;
ODS graphics on;
proc GLM data=gold;
class dentist method alloy temp;
model y=dentist|method|alloy|temp;
run; quit;
我哪里做错了?
这是输出的一部分:
The GLM Procedure
Dependent Variable: y
Source DF Sum of Squares Mean Square F Value Pr > F
Model 89 1891095.556 21248.265 . .
Error 0 0.000 .
Total 89 1891095.556
R-Square Coeff Var Root MSE y Mean
1.000000 . . 741.7778
错误应该是
Residuals 75772.0 16 4735.7
残差/错误不应该为 0,因为整个代码都是错误的。:(
我还需要知道如何为上述代码创建交互图/图表。对我的代码的任何帮助将不胜感激。