正交试验方差分析-SAS程序说明书.pdf

正交试验方差分析程序说明书 一、程序代码： data ex; input a b c d@@; do i=1 to 3; input y @@; output; end; cards; 1 1 1 1 8.86 8.81 8.88 1 2 2 2 6.09 6.14 6.04 1 3 3 3 7.71 7.76 7.81 2 1 2 3 9.26 9.21 9.18 2 2 3 1 9.61 9.56 9.58 2 3 1 2 9.63 9.61 9.61 3 1 3 2 10.23 10.20 10.18 3 2 1 3 8.06 8.01 7.98 3 3 2 1 9.23 9.18 9.21 ; proc anova;class a b c d; model y=a b c d; means a b c d/Duncan; run; 二、使用说明 （一）a、b 、c、d 为因素代码，本程序可以适用任何因素和水平的正交试验。常见的为四 因素三水平正交试验（L （34 ）），因为其只需要 9 组试验就能确定最佳工艺或配方，所以采 用较多。 （二）do i=1 to n，n 代表的是重复数，n 一般取 3 及以上。n 越大，估计误差的精度就越高。 （三）y 代表的是每组试验对应的指标。 （四）关于cards 后面的数据排列。 本实验是四因素三水平正交试验，如果是五因素四水平正交试验，则需要将正交试验表 换成 L （45 ）正交表，后面将指标附上就行。 （五）anova 为方差分析语句，class 是分类的意思， （六）model y=a b c d ，模型中 a、b 、c、d 就是因素，有几个列入几个。 （七）means a b c d/Duncan ，这是均值显著性分析，采用的是邓肯氏检验。邓肯氏检验常用 于均值之间的显著性检验，查看各处理组之间的显著性。 三、运行结果分析 以本程序为例，说明如果解读结果。 SAS 运行结果图 1 一、 运行 SAS 程序，之后，首先找到上图，从图 1 可以看出，a、b 、c、d 各因素 对指标 y 是极显著的（因为各因素 PrF 值均小于 0.01 ），其影响程度从小到大 依次为 a、b 、c、d （看F 值的大小） SAS 运行结果图 2 二、 从图 2 可以看出，a 因素各水平间的显著性，通过此可以选出最优的水平。大 写 A 、B 、C 不同，代表各 a 因素的各水平之间存在显著。