Stata Textbook Examples
Experimental Design by Roger Kirk
Chapter 13: Confounded Factorial Designs

Use data file rbcf2-2, page 596.
use http://www.ats.ucla.edu/stat/stata/examples/kirk/rbcf2-2, clear

list

            a         b         g         s         y 
  1.        0         0         0         0         3  
  2.        0         0         0         1         5  
  3.        0         0         0         2         6  
  4.        0         0         0         3         5  
  5.        1         1         0         0        16  
  6.        1         1         0         1        14  
  7.        1         1         0         2        17  
  8.        1         1         0         3        15  
  9.        0         1         1         4        14  
 10.        0         1         1         5        14  
 11.        0         1         1         6        16  
 12.        0         1         1         7        16  
 13.        1         0         1         4         7  
 14.        1         0         1         5         6  
 15.        1         0         1         6         7  
 16.        1         0         1         7        11
Table 13.3-1, page 596.
tabdisp s a b, cellvar(y) by(g) concise

----------+-------------------------
          |         b and a         
          | ---- 0 ---    ---- 1 ---
  g and s |    0     1       0     1
----------+-------------------------
0         |
        0 |    3                  16
        1 |    5                  14
        2 |    6                  17
        3 |    5                  15
----------+-------------------------
1         |
        4 |          7      14      
        5 |          6      14      
        6 |          7      16      
        7 |         11      16      
----------+-------------------------

table a b, cont(sum y) row col

----------+--------------------
          |          b         
        a |     0      1  Total
----------+--------------------
        0 |    19     60     79
        1 |    31     62     93
          | 
    Total |    50    122    172
----------+--------------------
Table 13.3-2, page 598.
anova y g / s|g a b / a*b*s|g

                     Number of obs =      16     R-squared     =  1.0000
                     Root MSE      =    0.00     Adj R-squared =

            Source |  Partial SS    df       MS           F     Prob > F
        -----------+----------------------------------------------------
             Model |      371.00    15  24.7333333   
                   |
                 g |        6.25     1        6.25       1.90     0.2174
               s|g |       19.75     6  3.29166667   
        -----------+----------------------------------------------------
                 a |       12.25     1       12.25       8.40     0.0274
                 b |      324.00     1      324.00     222.17     0.0000
           a*b*s|g |        8.75     6  1.45833333   
        -----------+----------------------------------------------------
                   |
          Residual |        0.00     0
        -----------+----------------------------------------------------
             Total |      371.00    15  24.7333333
Use data file rbcr2-3, page 604.
use http://www.ats.ucla.edu/stat/stata/examples/kirk/rbcf2-3, clear

list

            a         b         c         g         s         y 
  1.        0         0         0         0         0         3  
  2.        0         0         0         0         1         6  
  3.        0         0         0         0         2         3  
  4.        0         0         0         0         3         3  
  5.        0         1         1         0         0         7  
  6.        0         1         1         0         1         8  
  7.        0         1         1         0         2         9  
  8.        0         1         1         0         3         8  
  9.        1         0         1         0         0         2  
 10.        1         0         1         0         1         3  
 11.        1         0         1         0         2         4  
 12.        1         0         1         0         3         3  
 13.        1         1         0         0         0         5  
 14.        1         1         0         0         1         6  
 15.        1         1         0         0         2         5  
 16.        1         1         0         0         3         6  
 17.        0         0         1         1         4         4  
 18.        0         0         1         1         5         5  
 19.        0         0         1         1         6         4  
 20.        0         0         1         1         7         3  
 21.        0         1         0         1         4         7  
 22.        0         1         0         1         5         8  
 23.        0         1         0         1         6         7  
 24.        0         1         0         1         7         6  
 25.        1         0         0         1         4         1  
 26.        1         0         0         1         5         2  
 27.        1         0         0         1         6         2  
 28.        1         0         0         1         7         2  
 29.        1         1         1         1         4        10  
 30.        1         1         1         1         5        10  
 31.        1         1         1         1         6         9  
 32.        1         1         1         1         7        11
Table 13.5-1, page 604.
tabdisp s a b if g==0, cellvar(y) by(c)

----------+-------------------------
          |         b and a         
          | ---- 0 ---    ---- 1 ---
  c and s |    0     1       0     1
----------+-------------------------
0         |
        0 |    3                   5
        1 |    6                   6
        2 |    3                   5
        3 |    3                   6
----------+-------------------------
1         |
        0 |          2       7      
        1 |          3       8      
        2 |          4       9      
        3 |          3       8      
----------+-------------------------

tabdisp s a b if g==1, cellvar(y) by(c)

----------+-------------------------
          |         b and a         
          | ---- 0 ---    ---- 1 ---
  c and s |    0     1       0     1
----------+-------------------------
0         |
        4 |          1       7      
        5 |          2       8      
        6 |          2       7      
        7 |          2       6      
----------+-------------------------
1         |
        4 |    4                  10
        5 |    5                  10
        6 |    4                   9
        7 |    3                  11
----------+-------------------------

table a b, cont(sum y) row col

----------+--------------------
          |          b         
        a |     0      1  Total
----------+--------------------
        0 |    31     60     91
        1 |    19     62     81
          | 
    Total |    50    122    172
----------+--------------------

table a c, cont(sum y) row

----------+-----------
          |     c     
        a |    0     1
----------+-----------
        0 |   43    48
        1 |   29    52
          | 
    Total |   72   100
----------+-----------

table b c, cont(sum y)

----------+-----------
          |     c     
        b |    0     1
----------+-----------
        0 |   22    28
        1 |   50    72
----------+-----------

table a c b, cont(sum y)

----------+-------------------------
          |         b and c         
          | ---- 0 ---    ---- 1 ---
        a |    0     1       0     1
----------+-------------------------
        0 |   15    16      28    32
        1 |    7    12      22    40
----------+-------------------------
Table 13.5-2, page 607.
anova y g / s|g a b c a*b a*c b*c / a*b*c*s|g

                     Number of obs =      32     R-squared     =  1.0000
                     Root MSE      =    0.00     Adj R-squared =

            Source |  Partial SS    df       MS           F     Prob > F
        -----------+----------------------------------------------------
             Model |      235.50    31  7.59677419   
                   |
                 g |       3.125     1       3.125       2.94     0.1372
               s|g |       6.375     6      1.0625   
        -----------+----------------------------------------------------
                 a |       3.125     1       3.125       4.64     0.0451
                 b |      162.00     1      162.00     240.49     0.0000
                 c |       24.50     1       24.50      36.37     0.0000
               a*b |       6.125     1       6.125       9.09     0.0074
               a*c |      10.125     1      10.125      15.03     0.0011
               b*c |        8.00     1        8.00      11.88     0.0029
         a*b*c*s|g |      12.125    18  .673611111   
        -----------+----------------------------------------------------
                   |
          Residual |        0.00     0
        -----------+----------------------------------------------------
             Total |      235.50    31  7.59677419
Use data file rbpf2-3, page 609.
use http://www.ats.ucla.edu/stat/stata/examples/kirk/rbpf2-3, clear

list

            a         b         c         g         s         y 
  1.        0         0         0         0         0         3  
  2.        0         0         1         0         0         4  
  3.        1         1         0         0         0         5  
  4.        1         1         1         0         0        10  
  5.        0         1         0         0         1         7  
  6.        0         1         1         0         1         7  
  7.        1         0         0         0         1         1  
  8.        1         0         1         0         1         2  
  9.        0         0         0         1         2         6  
 10.        0         1         0         1         2         8  
 11.        1         0         1         1         2         3  
 12.        1         1         1         1         2        10  
 13.        0         0         1         1         3         5  
 14.        0         1         1         1         3         8  
 15.        1         0         0         1         3         2  
 16.        1         1         0         1         3         6  
 17.        0         0         0         2         4         3  
 18.        0         1         1         2         4         9  
 19.        1         0         0         2         4         2  
 20.        1         1         1         2         4         9  
 21.        0         0         1         2         5         4  
 22.        0         1         0         2         5         7  
 23.        1         0         1         2         5         4  
 24.        1         1         0         2         5         5  
 25.        0         0         0         3         6         3  
 26.        0         1         1         3         6         8  
 27.        1         0         1         3         6         3  
 28.        1         1         0         3         6         6  
 29.        0         0         1         3         7         3  
 30.        0         1         0         3         7         6  
 31.        1         0         0         3         7         2  
 32.        1         1         1         3         7        11 
Parts of Table 13.7-1, page 609.
table a b, cont(sum y) row col

----------+--------------------
          |          b         
        a |     0      1  Total
----------+--------------------
        0 |    31     60     91
        1 |    19     62     81
          | 
    Total |    50    122    172
----------+--------------------

table a c, cont(sum y) row 

----------+-----------
          |     c     
        a |    0     1
----------+-----------
        0 |   43    48
        1 |   29    52
          | 
    Total |   72   100
----------+-----------
Table 13.7-2, page 614.
Note: Ignore F-ratios for g and s|g.
anova y g s|g a b c a*b a*c b*c a*b*c

                     Number of obs =      32     R-squared     =  0.9478
                     Root MSE      = .850317     Adj R-squared =  0.9048

            Source |  Partial SS    df       MS           F     Prob > F
        -----------+----------------------------------------------------
             Model |  223.208333    14  15.9434524      22.05     0.0000
                   |
                 g |        5.25     3        1.75       2.42     0.1016
               s|g |  .958333333     4  .239583333       0.33     0.8530
                 a |       3.125     1       3.125       4.32     0.0531
                 b |      162.00     1      162.00     224.05     0.0000
                 c |       24.50     1       24.50      33.88     0.0000
               a*b |       3.375     1       3.375       4.67     0.0453
               a*c |        6.00     1        6.00       8.30     0.0104
               b*c |  7.04166667     1  7.04166667       9.74     0.0062
             a*b*c |  2.66666667     1  2.66666667       3.69     0.0717
                   |
          Residual |  12.2916667    17  .723039216   
        -----------+----------------------------------------------------
             Total |      235.50    31  7.59677419
Use data file rbcf3-2, page 618.
use http://www.ats.ucla.edu/stat/stata/examples/kirk/rbcf3-2, clear

list

            a         b         g         s         y 
  1.        0         0         0         0         3  
  2.        1         2         0         0         7  
  3.        2         1         0         0         4  
  4.        0         0         0         1         2  
  5.        1         2         0         1         6  
  6.        2         1         0         1         3  
  7.        0         1         1         2         4  
  8.        1         0         1         2         3  
  9.        2         2         1         2         6  
 10.        0         1         1         3         2  
 11.        1         0         1         3         1  
 12.        2         2         1         3         5  
 13.        0         2         2         4         7  
 14.        1         1         2         4         3  
 15.        2         0         2         4         3  
 16.        0         2         2         5         5  
 17.        1         1         2         5         4  
 18.        2         0         2         5         2  
Table 13.8-1, page 618.
tabdisp s a b, cellvar(y) by(g) concise

----------+---------------------------------------------------------
          |                         b and a                         
          | ------- 0 ------    ------- 1 ------    ------- 2 ------
  g and s |    0     1     2       0     1     2       0     1     2
----------+---------------------------------------------------------
0         |
        0 |    3                               4             7      
        1 |    2                               3             6      
----------+---------------------------------------------------------
1         |
        2 |          3             4                               6
        3 |          1             2                               5
----------+---------------------------------------------------------
2         |
        4 |                3             3             7            
        5 |                2             4             5            
----------+---------------------------------------------------------

table a b, cont(sum y) row col

----------+---------------------------
          |             b             
        a |     0      1      2  Total
----------+---------------------------
        0 |     5      6     12     23
        1 |     4      7     13     24
        2 |     5      7     11     23
          | 
    Total |    14     20     36     70
----------+---------------------------
Table 13.8-2, page 620.
anova y g / s|g a b a*b

                     Number of obs =      18     R-squared     =  0.9504
                     Root MSE      = .666667     Adj R-squared =  0.8595

            Source |  Partial SS    df       MS           F     Prob > F
        -----------+----------------------------------------------------
             Model |  51.1111111    11  4.64646465      10.45     0.0046
                   |
                 g |  1.16666667     2  .583333333       0.28     0.7760
               s|g |  6.33333333     3  2.11111111   
        -----------+----------------------------------------------------
                 a |  .111111111     2  .055555556       0.13     0.8847
                 b |  43.1111111     2  21.5555556      48.50     0.0002
               a*b |  .111111111     2  .055555556       0.13     0.8847
                   |
          Residual |  2.66666667     6  .444444444   
        -----------+----------------------------------------------------
             Total |  53.7777778    17  3.16339869 
Use data file rbpf3-2, page 621.
use http://www.ats.ucla.edu/stat/stata/examples/kirk/rbpf3-2, clear

list

            a         b         g        ab         s         y 
  1.        0         0         0         0         0         3  
  2.        1         2         0         0         0         7  
  3.        2         1         0         0         0         4  
  4.        0         1         0         1         1         4  
  5.        1         0         0         1         1         3  
  6.        2         2         0         1         1         6  
  7.        0         2         0         2         2         7  
  8.        1         1         0         2         2         3  
  9.        2         0         0         2         2         3  
 10.        0         0         1         0         3         2  
 11.        1         1         1         0         3         4  
 12.        2         2         1         0         3         5  
 13.        0         2         1         1         4         5  
 14.        1         0         1         1         4         1  
 15.        2         1         1         1         4         3  
 16.        0         1         1         2         5         2  
 17.        1         2         1         2         5         6  
 18.        2         0         1         2         5         2 
Table 13.8-3, page 621.
table a b, cont(sum y) row col

----------+---------------------------
          |             b             
        a |     0      1      2  Total
----------+---------------------------
        0 |     5      6     12     23
        1 |     4      7     13     24
        2 |     5      7     11     23
          | 
    Total |    14     20     36     70
----------+---------------------------

table a g, cont(sum y)

----------+-----------
          |     g     
        a |    0     1
----------+-----------
        0 |   14     9
        1 |   13    11
        2 |   13    10
----------+-----------
Table 13.8-4, page 624.

Note: Ignore F-ratios for g and s|ab.
anova y g s|ab a b a*b

                     Number of obs =      18     R-squared     =  0.9773
                     Root MSE      = .552771     Adj R-squared =  0.9034

            Source |  Partial SS    df       MS           F     Prob > F
        -----------+----------------------------------------------------
             Model |  52.5555556    13  4.04273504      13.23     0.0116
                   |
                 g |  .711111111     1  .711111111       2.33     0.2018
              s|ab |  2.22222222     4  .555555556       1.82     0.2884
                 a |  .111111111     2  .055555556       0.18     0.8403
                 b |  43.1111111     2  21.5555556      70.55     0.0008
               a*b |  2.88888889     4  .722222222       2.36     0.2126
                   |
          Residual |  1.22222222     4  .305555556   
        -----------+----------------------------------------------------
             Total |  53.7777778    17  3.16339869     
Use data file rbpf32-2, page 643.
use http://www.ats.ucla.edu/stat/stata/examples/kirk/rbpf32-2, clear

list    

            a         b         c         g         s         y 
  1.        0         0         0         0         0         3  
  2.        0         1         1         0         0         7  
  3.        1         0         1         0         0         4  
  4.        1         1         0         0         0         7  
  5.        2         0         1         0         0         3  
  6.        2         1         0         0         0         6  
  7.        0         0         1         0         1         4  
  8.        0         1         0         0         1         7  
  9.        1         0         0         0         1         3  
 10.        1         1         1         0         1         9  
 11.        2         0         0         0         1         3  
 12.        2         1         1         0         1         8  
 13.        0         0         1         1         2         2  
 14.        0         1         0         1         2         5  
 15.        1         0         0         1         2         2  
 16.        1         1         1         1         2        10  
 17.        2         0         1         1         2         3  
 18.        2         1         0         1         2         6  
 19.        0         0         0         1         3         1  
 20.        0         1         1         1         3        10  
 21.        1         0         1         1         3         3  
 22.        1         1         0         1         3         6  
 23.        2         0         0         1         3         2  
 24.        2         1         1         1         3        11  
 25.        0         0         1         2         4         3  
 26.        0         1         0         2         4         6  
 27.        1         0         1         2         4         4  
 28.        1         1         0         2         4         7  
 29.        2         0         0         2         4         2  
 30.        2         1         1         2         4         9  
 31.        0         0         0         2         5         2  
 32.        0         1         1         2         5         8  
 33.        1         0         0         2         5         3  
 34.        1         1         1         2         5         8  
 35.        2         0         1         2         5         3  
 36.        2         1         0         2         5         7 
Parts of Table 13.11-2, page 643.
table a b, cont(sum y) row col

----------+--------------------
          |          b         
        a |     0      1  Total
----------+--------------------
        0 |    15     43     58
        1 |    19     47     66
        2 |    16     47     63
          | 
    Total |    50    137    187
----------+--------------------

table a c, cont(sum y) row

----------+-----------
          |     c     
        a |    0     1
----------+-----------
        0 |   24    34
        1 |   28    38
        2 |   26    37
          | 
    Total |   78   109
----------+-----------
Table 13.11-3, page 647.

Note: Ignore F-ratios for g and s|g.
anova y g s|g a b c a*b a*c b*c a*b*c

                     Number of obs =      36     R-squared     =  0.9287
                     Root MSE      = 1.00204     Adj R-squared =  0.8687

            Source |  Partial SS    df       MS           F     Prob > F
        -----------+----------------------------------------------------
             Model |  248.561111    16  15.5350694      15.47     0.0000
                   |
                 g |  .388888889     2  .194444444       0.19     0.8256
               s|g |        1.20     3         .40       0.40     0.7557
                 a |  2.72222222     2  1.36111111       1.36     0.2816
                 b |      210.25     1      210.25     209.39     0.0000
                 c |  26.6944444     1  26.6944444      26.59     0.0001
               a*b |         .50     2         .25       0.25     0.7821
               a*c |  .055555556     2  .027777778       0.03     0.9728
               b*c |        4.50     1        4.50       4.48     0.0477
             a*b*c |  .033333333     2  .016666667       0.02     0.9836
                   |
          Residual |  19.0777778    19  1.00409357   
        -----------+----------------------------------------------------
             Total |  267.638889    35   7.6468254 
Use data file lscf3-2, page 653.
use http://www.ats.ucla.edu/stat/stata/examples/kirk/lscf3-2, clear

list  

            a         b         g         s         y 
  1.        0         0         0         0         3  
  2.        1         1         0         0         3  
  3.        2         2         0         0         6  
  4.        0         0         0         1         2  
  5.        1         1         0         1         4  
  6.        2         2         0         1         5  
  7.        1         0         1         2         3  
  8.        2         1         1         2         4  
  9.        0         2         1         2         7  
 10.        1         0         1         3         1  
 11.        2         1         1         3         3  
 12.        0         2         1         3         5  
 13.        2         0         2         4         3  
 14.        0         1         2         4         4  
 15.        1         2         2         4         7  
 16.        2         0         2         5         2  
 17.        0         1         2         5         2  
 18.        1         2         2         5         6  
Table 13.12-1, page 653.
tabdisp s b, cellvar(y) by(g) concise

----------+-----------------
          |        b        
  g and s |    0     1     2
----------+-----------------
0         |
        0 |    3     3     6
        1 |    2     4     5
----------+-----------------
1         |
        2 |    3     4     7
        3 |    1     3     5
----------+-----------------
2         |
        4 |    3     4     7
        5 |    2     2     6
----------+-----------------   

table a b, cont(sum y) row col

----------+---------------------------
          |             b             
        a |     0      1      2  Total
----------+---------------------------
        0 |     5      6     12     23
        1 |     4      7     13     24
        2 |     5      7     11     23
          | 
    Total |    14     20     36     70
----------+---------------------------
Table 13.12-2, page 655.
anova y g / s|g a b a*b

                     Number of obs =      18     R-squared     =  0.9628
                     Root MSE      =  .57735     Adj R-squared =  0.8946

            Source |  Partial SS    df       MS           F     Prob > F
        -----------+----------------------------------------------------
             Model |  51.7777778    11  4.70707071      14.12     0.0020
                   |
                 g |  1.16666667     2  .583333333       0.25     0.7936
               s|g |        7.00     3  2.33333333   
        -----------+----------------------------------------------------
                 a |  .111111111     2  .055555556       0.17     0.8503
                 b |  43.1111111     2  21.5555556      64.67     0.0001
               a*b |  1.44444444     2  .722222222       2.17     0.1958
                   |
          Residual |        2.00     6  .333333333   
        -----------+----------------------------------------------------
             Total |  53.7777778    17  3.16339869

How to cite this page

Report an error on this page or leave a comment

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.