glm5_4.sas

黄燕 吴平《SAS统计分析及应用》 书中一些例子的SAS程序代码

data data_prg.anova5_4_grass;                                                                                                                        
input turang $  GMD1   GMD2  @@;                                                                                                                
cards;                                                                                                                                  
t1 2.61 1.04 t1 3.04  1.10  t1 2.77  1.22  t1  2.92  1.04 t1  2.70  0.66
t2 2.15 1.57 t2 2.34  3.10  t2 2.98  1.68  t2  1.91  1.61
t3 0.62 1.36 t3 0.57  0.87  t3 0.43  0.76  t3  0.58  0.80 t3  0.47  0.78 
t4 1.47 1.01 t4 0.96  1.88  t4 1.41  0.73  t4  1.59  1.89 t4  1.48  0.86
t5 3.73 2.17 t5 3.80  1.35  t5 3.19  0.73  t5  3.71  2.02 t5  3.94  2.45
t6 3.65 2.33 t6 3.22  1.91  t6 4.21  1.17  t6  3.00  1.05
;                        
proc glm  data=data_prg.anova5_4_grass;                                                                                                                             
class turang;                                                                                                                                
model GMD1 GMD2=turang  ;                                                                                                                     
means turang / duncan regwq;                                                                                                                         
run;                                                                                                                                    
proc glm  data=data_prg.anova5_4_grass;                                                                                                                             
class turang;                                                                                                                                
model GMD1 GMD2=turang ;  
contrast  't1 vs t5 '  turang 1  0  0  0 -1  0 ;
contrast  't2 vs t5 '  turang 0  1  0  0 -1  0 ;
contrast  't3 vs t5 '  turang 0  0  1  0 -1  0 ;
contrast  't4 vs t5 '  turang 0  0  0  1 -1  0 ;
contrast  't6 vs t5 '  turang 0  0  0  0  1 -1 ;
contrast  't1 vs t6 '  turang 1  0  0  0 -1  0 ;
contrast  't2 vs t6 '  turang 0  1  0  0 -1  0 ;
contrast  't3 vs t6 '  turang 0  0  1  0 -1  0 ;
contrast  't4 vs t6 '  turang 0  0  0  1 -1  0 ;
contrast  't1 vs t2 '  turang 1 -1  0  0  0  0 ;
contrast  't3 vs t2 '  turang 0  1 -1  0  0  0 ;
contrast  't4 vs t2 '  turang 0  1  0 -1  0  0 ;
contrast  't5 vs t2 '  turang 0  1  0  0 -1  0 ;
contrast  't6 vs t2 '  turang 0  1  0  0  0 -1 ;
manova h=turang;
run; quit;