Use data file crff24-1, page 596.
use https://stats.idre.ucla.edu/stat/stata/examples/kirk/crff24-1, clear
list
a b c d y
1. 0 0 0 0 3
2. 0 0 0 0 6
3. 0 0 1 1 4
4. 0 0 1 1 3
5. 0 1 0 1 7
6. 0 1 0 1 6
7. 0 1 1 0 7
8. 0 1 1 0 8
9. 1 0 0 1 2
10. 1 0 0 1 2
11. 1 0 1 0 2
12. 1 0 1 0 3
13. 1 1 0 0 5
14. 1 1 0 0 6
15. 1 1 1 1 9
16. 1 1 1 1 11
Parts of Table 14.3-1, page 672.
table a c b, cont(sum y)
----------+-------------------------
| b and c
| ---- 0 --- ---- 1 ---
a | 0 1 0 1
----------+-------------------------
0 | 9 7 13 15
1 | 4 5 11 20
----------+-------------------------
table a b, cont(sum y) row col
----------+--------------------
| b
a | 0 1 Total
----------+--------------------
0 | 16 28 44
1 | 9 31 40
|
Total | 25 59 84
----------+--------------------
table a c, cont(sum y) row
----------+-----------
| c
a | 0 1
----------+-----------
0 | 22 22
1 | 15 25
|
Total | 37 47
----------+-----------
table b c, cont(sum y)
----------+-----------
| c
b | 0 1
----------+-----------
0 | 13 12
1 | 24 35
----------+-----------
Table 14.3-2, page 674.
anova y a b c d a*b a*c b*c
Number of obs = 16 R-squared = 0.9189
Root MSE = 1.06066 Adj R-squared = 0.8480
Source | Partial SS df MS F Prob > F
-----------+----------------------------------------------------
Model | 102.00 7 14.5714286 12.95 0.0009
|
a | 1.00 1 1.00 0.89 0.3734
b | 72.25 1 72.25 64.22 0.0000
c | 6.25 1 6.25 5.56 0.0462
d | 1.00 1 1.00 0.89 0.3734
a*b | 6.25 1 6.25 5.56 0.0462
a*c | 6.25 1 6.25 5.56 0.0462
b*c | 9.00 1 9.00 8.00 0.0222
|
Residual | 9.00 8 1.125
-----------+----------------------------------------------------
Total | 111.00 15 7.40
Use data file lsff33-2, page 688.
use https://stats.idre.ucla.edu/stat/stata/examples/kirk/lsff33-2, clear
list
a b c s y
1. 0 0 0 0 3
2. 0 1 1 0 4
3. 0 2 2 0 6
4. 0 0 0 1 1
5. 0 1 1 1 4
6. 0 2 2 1 6
7. 1 0 2 2 7
8. 1 1 0 2 3
9. 1 2 1 2 3
10. 1 0 2 3 5
11. 1 1 0 3 2
12. 1 2 1 3 3
13. 2 0 1 4 4
14. 2 1 2 4 7
15. 2 2 0 4 3
16. 2 0 1 5 2
17. 2 1 2 5 5
18. 2 2 0 5 2
Table 14.8-1, pages 688 and 689.
tabdisp s b, cellvar(y) by(a) concise
----------+-----------------
| b
a and s | 0 1 2
----------+-----------------
0 |
0 | 3 4 6
1 | 1 4 6
----------+-----------------
1 |
2 | 7 3 3
3 | 5 2 3
----------+-----------------
2 |
4 | 4 7 3
5 | 2 5 2
----------+-----------------
table b c, cont(sum y) by(a)
----------+-----------------
| c
a and b | 0 1 2
----------+-----------------
0 |
0 | 4
1 | 8
2 | 12
----------+-----------------
1 |
0 | 12
1 | 5
2 | 6
----------+-----------------
2 |
0 | 6
1 | 12
2 | 5
----------+-----------------
table c b, cont(sum y) row col
----------+---------------------------
| b
c | 0 1 2 Total
----------+---------------------------
0 | 4 5 5 14
1 | 6 8 6 20
2 | 12 12 12 36
|
Total | 22 25 23 70
----------+---------------------------
Table 14.8-2, page 690.
Note: Kirk calls a*b*c residual and calls the Residual b*c*s|a.
anova y a / s|a b c b*c*a
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
|
a | .111111111 2 .055555556 0.03 0.9742
s|a | 6.33333333 3 2.11111111
-----------+----------------------------------------------------
b | .777777778 2 .388888889 0.88 0.4640
c | 43.1111111 2 21.5555556 48.50 0.0002
b*c*a | .777777778 2 .388888889 0.87 0.4640
|
Residual | 2.66666667 6 .444444444
-----------+----------------------------------------------------
Total | 53.7777778 17 3.16339869
Use data file lsff43-2, page 693.
use https://stats.idre.ucla.edu/stat/stata/examples/kirk/lsff43-2, clear
list
a b c d y
1. 0 0 0 0 2
2. 0 1 3 0 3
3. 0 2 2 0 3
4. 0 3 1 0 6
5. 1 0 1 0 2
6. 1 1 0 0 4
7. 1 2 3 0 5
8. 1 3 2 0 4
9. 2 0 2 0 5
10. 2 1 1 0 8
11. 2 2 0 0 7
12. 2 3 3 0 6
13. 3 0 3 0 9
14. 3 1 2 0 9
15. 3 2 1 0 11
16. 3 3 0 0 10
17. 0 0 0 1 1
18. 0 1 3 1 2
19. 0 2 2 1 2
20. 0 3 1 1 3
21. 1 0 1 1 3
22. 1 1 0 1 3
23. 1 2 3 1 4
24. 1 3 2 1 3
25. 2 0 2 1 6
26. 2 1 1 1 6
27. 2 2 0 1 5
28. 2 3 3 1 7
29. 3 0 3 1 8
30. 3 1 2 1 8
31. 3 2 1 1 10
32. 3 3 0 1 7
Table 14.9-1, page 693.
tabdisp a b, cellvar(y) by(d)
----------+-----------------------
| b
d and a | 0 1 2 3
----------+-----------------------
0 |
0 | 2 3 3 6
1 | 2 4 5 4
2 | 5 8 7 6
3 | 9 9 11 10
----------+-----------------------
1 |
0 | 1 2 2 3
1 | 3 3 4 3
2 | 6 6 5 7
3 | 8 8 10 7
----------+-----------------------
table a d, cont(sum y) col
----------+--------------------
| d
a | 0 1 Total
----------+--------------------
0 | 14 8 22
1 | 15 13 28
2 | 26 24 50
3 | 39 33 72
----------+--------------------
table b d, cont(sum y) col
----------+--------------------
| d
b | 0 1 Total
----------+--------------------
0 | 18 18 36
1 | 24 19 43
2 | 26 21 47
3 | 26 20 46
----------+--------------------
Table 14.9-2, page 695.
anova y a b c d a*d b*d c*d
Number of obs = 32 R-squared = 0.9618
Root MSE = .866025 Adj R-squared = 0.9013
Source | Partial SS df MS F Prob > F
-----------+----------------------------------------------------
Model | 226.50 19 11.9210526 15.89 0.0000
|
a | 194.50 3 64.8333333 86.44 0.0000
b | 9.25 3 3.08333333 4.11 0.0320
c | 7.75 3 2.58333333 3.44 0.0517
d | 8.00 1 8.00 10.67 0.0068
a*d | 2.00 3 .666666667 0.89 0.4747
b*d | 2.75 3 .916666667 1.22 0.3441
c*d | 2.25 3 .75 1.00 0.4262
|
Residual | 9.00 12 .75
-----------+----------------------------------------------------
Total | 235.50 31 7.59677419