加载中…
个人资料
  • 博客等级:
  • 博客积分:
  • 博客访问:
  • 关注人气:
  • 获赠金笔:0支
  • 赠出金笔:0支
  • 荣誉徽章:
正文 字体大小:

SAS实现重复测量方差分析

(2013-09-08 15:09:56)
标签:

sas

重复测量

方差分析

教育

分类: SAS与统计学
data ex12_3;
  input t0-t4 g@@;
cards;
120      108      112      120      117      1
118      109      115      126      123      1
119      112      119      124      118      1
121      112      119      126      120      1
127      121      127      133      126      1
121      120      118      131      137      2
122      121      119      129      133      2
128      129      126      135      142      2
117      115      111      123      131      2
118      114      116      123      133      2
131      119      118      135      129      3
129      128      121      148      132      3
123      123      120      143      136      3
123      121      116      145      126      3
125      124      118      142      130      3
;
proc glm;
  class g;
  model t0-t4=g;
  repeated time 5  contrast(1);
run;

**********************************************

                                          The GLM Procedure

                                       Class Level Information

                                    Class         Levels    Values

                                    g                  3    1 2 3


                               Number of Observations Read          15
                               Number of Observations Used          15
                                            The SAS System       15:06 Thursday, September 8, 2013   2

                                          The GLM Procedure

Dependent Variable: t0

                                                 Sum of
         Source                      DF         Squares     Mean Square    F Value    Pr > F

         Model                        2      86.8000000      43.4000000       2.93    0.0919

         Error                       12     177.6000000      14.8000000

         Corrected Total             14     264.4000000


                          R-Square     Coeff Var      Root MSE       t0 Mean

                          0.328290      3.132799      3.847077      122.8000


         Source                      DF       Type I SS     Mean Square    F Value    Pr > F

         g                            2     86.80000000     43.40000000       2.93    0.0919


         Source                      DF     Type III SS     Mean Square    F Value    Pr > F

         g                            2     86.80000000     43.40000000       2.93    0.0919
                                            The SAS System       15:06 Thursday, September 8, 2013   3

                                          The GLM Procedure

Dependent Variable: t1

                                                 Sum of
         Source                      DF         Squares     Mean Square    F Value    Pr > F

         Model                        2     295.6000000     147.8000000       6.03    0.0154

         Error                       12     294.0000000      24.5000000

         Corrected Total             14     589.6000000


                          R-Square     Coeff Var      Root MSE       t1 Mean

                          0.501357      4.180530      4.949747      118.4000


         Source                      DF       Type I SS     Mean Square    F Value    Pr > F

         g                            2     295.6000000     147.8000000       6.03    0.0154


         Source                      DF     Type III SS     Mean Square    F Value    Pr > F

         g                            2     295.6000000     147.8000000       6.03    0.0154
                                            The SAS System       15:06 Thursday, September 8, 2013   4

                                          The GLM Procedure

Dependent Variable: t2

                                                 Sum of
         Source                      DF         Squares     Mean Square    F Value    Pr > F

         Model                        2       0.9333333       0.4666667       0.02    0.9788

         Error                       12     260.4000000      21.7000000

         Corrected Total             14     261.3333333


                          R-Square     Coeff Var      Root MSE       t2 Mean

                          0.003571      3.936613      4.658326      118.3333


         Source                      DF       Type I SS     Mean Square    F Value    Pr > F

         g                            2      0.93333333      0.46666667       0.02    0.9788


         Source                      DF     Type III SS     Mean Square    F Value    Pr > F

         g                            2      0.93333333      0.46666667       0.02    0.9788
                                            The SAS System       15:06 Thursday, September 8, 2013   5

                                          The GLM Procedure

Dependent Variable: t3

                                                 Sum of
         Source                      DF         Squares     Mean Square    F Value    Pr > F

         Model                        2      825.600000      412.800000      17.03    0.0003

         Error                       12      290.800000       24.233333

         Corrected Total             14     1116.400000


                          R-Square     Coeff Var      Root MSE       t3 Mean

                          0.739520      3.723704      4.922736      132.2000


         Source                      DF       Type I SS     Mean Square    F Value    Pr > F

         g                            2     825.6000000     412.8000000      17.03    0.0003


         Source                      DF     Type III SS     Mean Square    F Value    Pr > F

         g                            2     825.6000000     412.8000000      17.03    0.0003
                                            The SAS System       15:06 Thursday, September 8, 2013   6

                                          The GLM Procedure

Dependent Variable: t4

                                                 Sum of
         Source                      DF         Squares     Mean Square    F Value    Pr > F

         Model                        2     540.9333333     270.4666667      17.37    0.0003

         Error                       12     186.8000000      15.5666667

         Corrected Total             14     727.7333333


                          R-Square     Coeff Var      Root MSE       t4 Mean

                          0.743313      3.061662      3.945462      128.8667


         Source                      DF       Type I SS     Mean Square    F Value    Pr > F

         g                            2     540.9333333     270.4666667      17.37    0.0003


         Source                      DF     Type III SS     Mean Square    F Value    Pr > F

         g                            2     540.9333333     270.4666667      17.37    0.0003
                                            The SAS System       15:06 Thursday, September 8, 2013   7

                                          The GLM Procedure
                                Repeated Measures Analysis of Variance

                                  Repeated Measures Level Information

                  Dependent Variable          t0       t1       t2       t3       t4

                       Level of time           1        2        3        4        5


           MANOVA Test Criteria and Exact F Statistics for the Hypothesis of no time Effect
                                  H = Type III SSCP Matrix for time
                                        E = Error SSCP Matrix

                                         S=1    M=1    N=3.5

           Statistic                        Value    F Value    Num DF    Den DF    Pr > F

           Wilks' Lambda               0.01745414     126.66         4         9    <.0001
           Pillai's Trace              0.98254586     126.66         4         9    <.0001
           Hotelling-Lawley Trace     56.29299264     126.66         4         9    <.0001
           Roy's Greatest Root        56.29299264     126.66         4         9    <.0001


           MANOVA Test Criteria and F Approximations for the Hypothesis of no time*g Effect
                                 H = Type III SSCP Matrix for time*g
                                        E = Error SSCP Matrix

                                        S=2    M=0.5    N=3.5

           Statistic                        Value    F Value    Num DF    Den DF    Pr > F

           Wilks' Lambda               0.00845847      22.21         8        18    <.0001
           Pillai's Trace              1.80884126      23.66         8        20    <.0001
           Hotelling-Lawley Trace     20.59966971      22.13         8     10.78    <.0001
           Roy's Greatest Root        13.37581359      33.44         4        10    <.0001

                    NOTE: F Statistic for Roy's Greatest Root is an upper bound.
                            NOTE: F Statistic for Wilks' Lambda is exact.
                                            The SAS System       15:06 Thursday, September 8, 2013   8

                                          The GLM Procedure
                                Repeated Measures Analysis of Variance
                           Tests of Hypotheses for Between Subjects Effects

         Source                      DF     Type III SS     Mean Square    F Value    Pr > F

         g                            2     912.2400000     456.1200000       5.78    0.0174
         Error                       12     946.4800000      78.8733333
                                            The SAS System       15:06 Thursday, September 8, 2013   9

                                          The GLM Procedure
                                Repeated Measures Analysis of Variance
                      Univariate Tests of Hypotheses for Within Subject Effects

                                                                                       Adj Pr > F
   Source                     DF    Type III SS    Mean Square   F Value   Pr > F    G - G    H - F

   time                        4    2336.453333     584.113333    106.56   <.0001   <.0001   <.0001
   time*g                      8     837.626667     104.703333     19.10   <.0001   <.0001   <.0001
   Error(time)                48     263.120000       5.481667


                                 Greenhouse-Geisser Epsilon    0.6787


0

阅读 收藏 喜欢 打印举报/Report
  

新浪BLOG意见反馈留言板 欢迎批评指正

新浪简介 | About Sina | 广告服务 | 联系我们 | 招聘信息 | 网站律师 | SINA English | 产品答疑

新浪公司 版权所有