UCLA Academic Technology Services HomeServicesClassesContactJobs

Stata Textbook Examples
Experimental Design by Roger Kirk
Chapter 9: Completely Randomized Factorial Design with Two Treatments

Table 9.3-1, page 369.
use http://www.ats.ucla.edu/stat/stata/examples/kirk/crf33, clear

tabdisp order b a, cellvar(y)

----------+---------------------------------------------------------
          |                         a and b                         
          | ------- 1 ------    ------- 2 ------    ------- 3 ------
    order |    1     2     3       1     2     3       1     2     3
----------+---------------------------------------------------------
        1 |   24    44    38      30    35    26      21    41    42
        2 |   33    36    29      21    40    27      18    39    52
        3 |   37    25    28      39    27    36      10    50    53
        4 |   29    27    47      26    31    46      31    36    49
        5 |   42    43    48      34    22    45      20    34    64
----------+---------------------------------------------------------

table a b, cont(sum y)

----------+-----------------
          |        b        
        a |    1     2     3
----------+-----------------
        1 |  165   175   190
        2 |  150   155   180
        3 |  100   200   260
----------+-----------------
Table 9.3-2, page 370.
anova y a b a*b

                     Number of obs =      45     R-squared     =  0.5690
                     Root MSE      = 7.90569     Adj R-squared =  0.4732

            Source |  Partial SS    df       MS           F     Prob > F
        -----------+----------------------------------------------------
             Model |     2970.00     8      371.25       5.94     0.0001
                   |
                 a |      190.00     2       95.00       1.52     0.2324
                 b |  1543.33333     2  771.666667      12.35     0.0001
               a*b |  1236.66667     4  309.166667       4.95     0.0028
                   |
          Residual |     2250.00    36       62.50   
        -----------+----------------------------------------------------
             Total |     5220.00    44  118.636364
Figure 9.3-1, page 371.
Note:  You can download anovaplot from within Stata by typing findit anovaplot (see How can I use the findit command to search for programs and get additional help? for more information about using findit).
anovaplot, scatter(msymbol(i))  /* Figure 9.3-1 (a) */
quietly anova y a b a*b
anovaplot, scatter(msymbol(i))  /* Figure 9.3-1 (b) */ 
Table 9.6-2, page 380.

Note: The sme command can be downloaded by typing findit contrast (see How can I use the findit command to search for programs and get additional help? for more information about using findit).
sme a b
 
Test of a at b(1): F(2/36)  = 3.7066667
Test of a at b(2): F(2/36)  = 1.6266667
Test of a at b(3): F(2/36)  = 6.08


Critical value of F for alpha = .05 using ...
--------------------------------------------------
Dunn's procedure              = 4.0941238
Marascuilo & Levin            = 4.5974255
per family error rate         = 4.5974255
simultaneous test procedure   = 6.5295994

sme b a
 
Test of b at a(1): F(2/36)  = .50666667
Test of b at a(2): F(2/36)  = .82666667
Test of b at a(3): F(2/36)  = 20.906667


Critical value of F for alpha = .05 using ...
--------------------------------------------------
Dunn's procedure              = 4.0941238
Marascuilo & Levin            = 4.5974255
per family error rate         = 4.5974255
simultaneous test procedure   = 6.5295994
Omega-squared computations, page 397.

Note: The omega2 command can be downloaded by typing findit contrast (see How can I use the findit command to search for programs and get additional help? for more information about using findit).
omega2 1.52 2  /* see Kirk, page 399, for formula */

 omega squared = 0.0226
 effect size   = 0.1520

omega2 12.35 2

 omega squared = 0.3353
 effect size   = 0.7102

omega2 4.95 4

 omega squared = 0.2599
 effect size   = 0.5925
Table 8-14-1, page 423.
use http://www.ats.ucla.edu/stat/stata/examples/kirk/crf33a, clear

tabdisp order b a, cellvar(y)

----------+---------------------------------------------------------
          |                         a and b                         
          | ------- 1 ------    ------- 2 ------    ------- 3 ------
    order |    1     2     3       1     2     3       1     2     3
----------+---------------------------------------------------------
        1 |   24    44    38      30    35    26      21    41    42
        2 |   33    36    29      21    40    27      18    39    52
        3 |   37    25    28      39    27    36      10    50    53
        4 |   29    27    47      26    31    46      31    36    49
        5 |         43    48      34    22    45      20    34    64
----------+---------------------------------------------------------

table a b, cont(freq mean y)

----------+--------------------
          |          b         
        a |     1      2      3
----------+--------------------
        1 |     5      5      5
          | 30.75     35     38
          | 
        2 |     5      5      5
          |    30     31     36
          | 
        3 |     5      5      5
          |    20     40     52
----------+--------------------
SSA, page 424.
anova y a b a*b, sequential

                     Number of obs =      44     R-squared     =  0.5844
                     Root MSE      = 7.83536     Adj R-squared =  0.4894

            Source |    Seq. SS     df       MS           F     Prob > F
        -----------+----------------------------------------------------
             Model |  3021.13636     8  377.642045       6.15     0.0001
                   |
                 a |  187.505411     2  93.7527056       1.53     0.2313
                 b |  1716.21429     2  858.107143      13.98     0.0000
               a*b |  1117.41667     4  279.354167       4.55     0.0046
                   |
          Residual |     2148.75    35  61.3928571   
        -----------+----------------------------------------------------
             Total |  5169.88636    43  120.229915 
SSB and SSAB, page 425.
anova y b a a*b, sequential

                     Number of obs =      44     R-squared     =  0.5844
                     Root MSE      = 7.83536     Adj R-squared =  0.4894

            Source |    Seq. SS     df       MS           F     Prob > F
        -----------+----------------------------------------------------
             Model |  3021.13636     8  377.642045       6.15     0.0001
                   |
                 b |  1713.33874     2  856.669372      13.95     0.0000
                 a |  190.380952     2  95.1904762       1.55     0.2264
               a*b |  1117.41667     4  279.354167       4.55     0.0046
                   |
          Residual |     2148.75    35  61.3928571   
        -----------+----------------------------------------------------
             Total |  5169.88636    43  120.229915

How to cite this page

Report an error on this page

UCLA Researchers are invited to our Statistical Consulting Services
We recommend others to our list of Other Resources for Statistical Computing Help
These pages are Copyrighted (c) by UCLA Academic Technology Services


The content of this web site should not be construed as an endorsement of any particular web site, book, or software product by the University of California.