~~~~~~帮帮俺吧,偶是菜鸟啊~~~~~

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

问题是这样的～～～～求最佳组合，大人近来

a,b ,c ,d ,e在最左边，表示的是５个影响因素，　y1，y2　，y3是变量指标　　，block是三个重复实验

一共４*４＝１６个处理，４８个分区

数据偶在下面了，不知道我说清楚了没？没有原程序，下面是我自己写的，肯定不对的，工参考

//////////////////////////　a=1.5 b=1 c=0.5 d=0 e=0这个是最佳组合也是答案，唉，弄不出来//////////////////////////

下面是大概的图形

                                                                  block1            block  2                      block3

obs　a　　b　　c　d　　e　　　y1　y2　y3   y4    y1    y2   y3   y4        y1    y2    y3   y4   

１　　0     0      0     0      0            16.3  16.9................................

.......................数据区域...................................................................................................

16 .......................................................................................................................................

data a;

 do a=0, 0.5,1,1.5;
   do b=0 ,1,2,3;
      do c=0 ,0.5,1,1.5;
         do d=0 ,0.045,0.9,0.135;
           do e=0,3,6,12;
             do block=1 ,2, 3;
      do z=1,2,3,4;
          input y@@;
          output;
        end;
      end;
   end;
 end;
end;
end;
end;
cards;

16.3 16.9  8.2  1.53   16.9 13.8   7.2  1.24 16.6 16.1 8.3  1.47
16.9 18.0  9.3  1.62   17.2 20.9   9.1  1.75 18.0 19.4 8.9  1.57
17   18.3  9.2  1.62   16.9 16.5   9.8  1.56 18.1 18.5 9.8  1.56
17   20.1  10   1.76   17.4 20.3   9.8  1.73 17.3 19.6 10.9 1.76
17.9 17.2  9.2  1.48   17.1 17.8   9.9  1.62 19.4 17.6 9.9  1.42
18   17.6  9.2  1.49   17.3 17     8.7  1.48 15.6 16.8 9.1  1.66
17.6 19.2  11.6 1.75   16.9 16.9   10.3 1.62 16.2 15.5 10.9 1.63
17   17.6  10.3 1.64   16   17.9   9.4  1.71 15.6 16.7 9.0  1.65
20.4 23    12.1 1.72   19.4 20.3   11.4 1.64 19.7 21.6 10.3 1.62
18.5 19.9  9.6  1.59   17.3 17.6   9.0  1.54 18.1 18.7 9.5  1.56
17.7 17    9.8  1.51   19.4 20.9   10.1 1.6  19.3 21.4 10.6 1.66
22.7 20.9 11.3  1.42   17.4 15.1   10.5 1.47 20.3 18.6 11.2 1.47
21.5 20.8 12.5  1.55   19.6 21.4   10.7 1.63 17.8 20.5 10.9 1.76
21.4 21.8 12.4  1.6    22.4 22.8   12.6 1.58 20.9 21.5 12.7 1.64
18.1 19.6 8.7   1.57   17.1 20.2   10.0 1.77 17.4 18.8 9.1  1.60
;
proc glm  ;
class a b c d e;
model y =a b c d e ;
manova h=a b c d e;
means  a b c /regewq ;
run;

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：MANOVA output block outpu model 菜鸟

循环输出是弄不出来了,下面是我用最笨的方法弄出来的,结果是对了,有高手给用循环弄弄不?

data a;
input g a b c d e y1 y2 y3 y4@@;
cards;
1 0 0 0 0 0 16.3 16.9 8.2 1.53
1 0 1 0.5 0.045 3 16.9 18 9.3 1.62
1 0 2 1 0.9 6 17 18.3 9.2 1.62
1 0 3 1.5 0.135 12 17 20.1 10 1.76
1 0.5 0 0 0.9 0 17.9 17.2 9.2 1.48
1 0.5 1 0.5 0.135 3 18 17.6 9.2 1.49
1 0.5 2 1 0 6 17.6 19.2 11.6 1.75
1 0.5 3 1.5 0.045 12 17 17.6 10.3 1.64
1 1 0 0 0.135 0 19.3 15 8.1 1.19
1 1 1 0.5 0.9 3 20.4 23 12.1 1.72
1 1 2 1 0.045 6 18.5 19.9 9.6 1.59
1 1 3 1.5 0 12 17.7 17 9.8 1.51
1 1.5 0 0 0.045 0 22.7 20.9 11.3 1.42
1 1.5 1 0.5 0 3 21.5 20.8 12.5 1.55
1 1.5 2 1 0.135 6 21.4 21.8 12.4 1.6
1 1.5 3 1.5 0.9 12 18.1 19.6 8.7 1.57
2 0 0 0 0 0 16.9 13.8 7.2 1.24
2 0 1 0.5 0.045 3 17.2 20.9 9.1 1.75
2 0 2 1 0.9 6 16.9 16.5 9.8 1.56
2 0 3 1.5 0.135 12 17.4 20.3 9.8 1.73
2 0.5 0 0 0.9 0 17.1 17.8 9.9 1.62
2 0.5 1 0.5 0.135 3 17.3 17 8.7 1.48
2 0.5 2 1 0 6 16.9 16.9 10.3 1.62
2 0.5 3 1.5 0.045 12 16 17.9 9.4 1.71
2 1 0 0 0.135 0 18.2 14.1 7.7 1.2
2 1 1 0.5 0.9 3 19.4 20.3 11.4 1.64
2 1 2 1 0.045 6 17.3 17.6 9 1.54
2 1 3 1.5 0 12 19.4 20.9 10.1 1.6
2 1.5 0 0 0.045 0 17.4 15.1 10.5 1.47
2 1.5 1 0.5 0 3 19.6 21.4 10.7 1.63
2 1.5 2 1 0.135 6 22.4 22.8 12.6 1.58
2 1.5 3 1.5 0.9 12 17.1 20.2 10 1.77
3 0 0 0 0 0 16.6 16.1 8.3 1.47
3 0 1 0.5 0.045 3 18 19.4 8.9 1.57
3 0 2 1 0.9 6 18.1 18.5 9.8 1.56
3 0 3 1.5 0.135 12 17.3 19.6 10.9 1.76
3 0.5 0 0 0.9 0 19.4 17.6 9.9 1.42
3 0.5 1 0.5 0.135 3 15.6 16.8 9.1 1.66
3 0.5 2 1 0 6 16.2 15.5 10.9 1.63
3 0.5 3 1.5 0.045 12 15.6 16.7 9 1.65
3 1 0 0 0.135 0 16.7 15.9 8.7 1.47
3 1 1 0.5 0.9 3 19.7 21.6 10.3 1.62
3 1 2 1 0.045 6 18.1 18.7 9.5 1.56
3 1 3 1.5 0 12 19.3 21.4 10.6 1.66
3 1.5 0 0 0.045 0 20.3 18.6 11.2 1.47
3 1.5 1 0.5 0 3 17.8 20.5 10.9 1.76
3 1.5 2 1 0.135 6 20.9 21.5 12.7 1.64
3 1.5 3 1.5 0.9 12 17.4 18.8 9.1 1.6
;
run;
proc glm data=a;
class g a b c d e;
model y1-y4=a b c d e;
means a b c d e/ regwq;
run;