Table 12.2-1, page 516.
use https://stats.idre.ucla.edu/stat/stata/examples/kirk/spf2-4, clear
tabdisp s b, cellvar(y) by(a) concise
----------+-----------------------
| b
a and s | 1 2 3 4
----------+-----------------------
1 |
1 | 3 4 7 7
2 | 6 5 8 8
3 | 3 4 7 9
4 | 3 3 6 8
----------+-----------------------
2 |
5 | 1 2 5 10
6 | 2 3 6 10
7 | 2 4 5 9
8 | 2 3 6 11
----------+-----------------------
table a b, cont(sum y) row col
----------+----------------------------------
| b
a | 1 2 3 4 Total
----------+----------------------------------
1 | 15 16 28 32 91
2 | 7 12 22 40 81
|
Total | 22 28 50 72 172
----------+----------------------------------
Table 12.2-2, page 517.
Note: The repeated option is used to compute the Huynh-Feldt values.
anova y a / s|a b a*b /, repeated(b)
Number of obs = 32 R-squared = 0.9613
Root MSE = .712 Adj R-squared = 0.9333
Source | Partial SS df MS F Prob > F
-----------+----------------------------------------------------
Model | 226.375 13 17.4134615 34.35 0.0000
|
a | 3.125 1 3.125 2.00 0.2070
s|a | 9.375 6 1.5625
-----------+----------------------------------------------------
b | 194.50 3 64.8333333 127.89 0.0000
a*b | 19.375 3 6.45833333 12.74 0.0001
|
Residual | 9.125 18 .506944444
-----------+----------------------------------------------------
Total | 235.50 31 7.59677419
Between-subjects error term: s|a
Levels: 8 (6 df)
Lowest b.s.e. variable: s
Covariance pooled over: a (for repeated variable)
Repeated variable: b
Huynh-Feldt epsilon = 0.9432
Greenhouse-Geisser epsilon = 0.5841
Box's conservative epsilon = 0.3333
------------ Prob > F ------------
Source | df F Regular H-F G-G Box
-----------+----------------------------------------------------
b | 3 127.89 0.0000 0.0000 0.0000 0.0000
a*b | 3 12.74 0.0001 0.0002 0.0019 0.0118
Residual | 18
-----------+----------------------------------------------------
Figure 12.2-1, page 518.
egen ybar_a1 = mean(y) if a == 1, by(b) egen ybar_a2 = mean(y) if a == 2, by(b) graph twoway (scatter ybar_a1 b, connect(l) sort mlabel(a) mlabposition(12)) /// (scatter ybar_a2 b, connect(l) sort mlabel(a) mlabposition(12))
egen ybar_b1 = mean(y) if b == 1, by(a) egen ybar_b2 = mean(y) if b == 2, by(a) egen ybar_b3 = mean(y) if b == 3, by(a) egen ybar_b4 = mean(y) if b == 4, by(a) graph twoway (scatter ybar_b1 a, connect(L) sort mlabel(b)mlabposition(12)) /// (scatter ybar_b2 a, connect(L) sort mlabel(b) mlabposition(12)) /// (scatter ybar_b3 a, connect(L) sort mlabel(b) mlabposition(12)) /// (scatter ybar_b4 a, connect(L) sort mlabel(b) mlabposition(12))
Omega-squareds, page 519.
Note: The omega2 command can be downloaded by typing search contrast (see How can I use the search command to search for programs and get additional help? for more information about using search).
omega2 2 1 omega squared = 0.0303 effect size = 0.1768 omega2 127.89 3 omega squared = 0.9225 effect size = 3.4490 omega2 12.74 3 omega squared = 0.5240 effect size = 1.0491
Table 12.4-1, page 528.
anova y a / s|a b a*b /, repeated(b) noanova /* repeat anova with no output */
matrix list e(Srep)
symmetric e(Srep)[4,4]
c1 c2 c3 c4
r1 1.25
r2 .66666667 .66666667
r3 .58333333 .33333333 .5
r4 0 -.16666667 .16666667 .66666667
Table 12.8-1, page 542.
Note: The entry for s1 a1 c1 b4 is incorrect in the book. It should have the value 7. The data file, spf22-4.dta, has the correct value.
use https://stats.idre.ucla.edu/stat/stata/examples/kirk/spf22-4, clear
tabdisp s b, cellvar(y) by(a c) concise /* note corrected value */
----------+-----------------------
a, c and | b
s | 1 2 3 4
----------+-----------------------
1 |
1 |
1 | 3 4 7 7
2 | 6 5 8 8
----------+-----------------------
1 |
2 |
3 | 3 4 7 9
4 | 3 3 6 8
----------+-----------------------
2 |
1 |
5 | 1 2 5 10
6 | 2 3 6 10
----------+-----------------------
2 |
2 |
7 | 2 4 5 9
8 | 2 3 6 11
----------+-----------------------
table a b c, cont(sum y)
----------+-------------------------------------------------
| c and b
| ---------- 1 --------- ---------- 2 ---------
a | 1 2 3 4 1 2 3 4
----------+-------------------------------------------------
1 | 9 9 15 15 6 7 13 17
2 | 3 5 11 20 4 7 11 20
----------+-------------------------------------------------
table c b, cont(sum y)
----------+-----------------------
| b
c | 1 2 3 4
----------+-----------------------
1 | 12 14 26 35
2 | 10 14 24 37
----------+-----------------------
table a b, cont(sum y) row col
----------+----------------------------------
| b
a | 1 2 3 4 Total
----------+----------------------------------
1 | 15 16 28 32 92
2 | 7 12 22 40 81
|
Total | 22 28 50 73 173
----------+----------------------------------
table a c, cont(sum y) row
----------+-----------
| c
a | 1 2
----------+-----------
1 | 48 43
2 | 39 42
|
Total | 88 85
----------+-----------
Table 12.8-2, page 545.
anova y a c a*c / s|a*c b a*b b*c a*b*c /, repeated(b)
Number of obs = 32 R-squared = 0.9798
Root MSE = .629153 Adj R-squared = 0.9479
Source | Partial SS df MS F Prob > F
-----------+----------------------------------------------------
Model | 230.75 19 12.1447368 30.68 0.0000
|
a | 3.125 1 3.125 1.72 0.2594
c | .125 1 .125 0.07 0.8058
a*c | 2.00 1 2.00 1.10 0.3528
s|a*c | 7.25 4 1.8125
-----------+----------------------------------------------------
b | 194.50 3 64.8333333 163.79 0.0000
a*b | 19.375 3 6.45833333 16.32 0.0002
b*c | 1.375 3 .458333333 1.16 0.3660
a*b*c | 3.00 3 1.00 2.53 0.1068
|
Residual | 4.75 12 .395833333
-----------+----------------------------------------------------
Total | 235.50 31 7.59677419
Between-subjects error term: s|a*c
Levels: 8 (4 df)
Lowest b.s.e. variable: s
Covariance pooled over: a*c (for repeated variable)
Repeated variable: b
Huynh-Feldt epsilon = 2.1181
*Huynh-Feldt epsilon reset to 1.0000
Greenhouse-Geisser epsilon = 0.6367
Box's conservative epsilon = 0.3333
------------ Prob > F ------------
Source | df F Regular H-F G-G Box
-----------+----------------------------------------------------
b | 3 163.79 0.0000 0.0000 0.0000 0.0002
a*b | 3 16.32 0.0002 0.0002 0.0018 0.0156
b*c | 3 1.16 0.3660 0.3660 0.3609 0.3425
a*b*c | 3 2.53 0.1068 0.1068 0.1448 0.1872
Residual | 12
-----------+----------------------------------------------------
matrix list e(Srep)
symmetric e(Srep)[4,4]
c1 c2 c3 c4
r1 1.25
r2 .5 .5
r3 .5 .25 .5
r4 .375 0 .5 .75
Figure 12.10-1, page 555.
use https://stats.idre.ucla.edu/stat/stata/examples/kirk/spf2-22, clear
tabdisp s b c, cellvar(y) by(a) concise
----------+-------------------------
| c and b
| ---- 1 --- ---- 2 ---
a and s | 1 2 1 2
----------+-------------------------
1 |
1 | 3 7 4 7
2 | 6 8 5 8
3 | 3 7 4 9
4 | 3 6 3 8
----------+-------------------------
2 |
5 | 1 5 2 10
6 | 2 6 3 10
7 | 2 5 4 9
8 | 2 6 3 11
----------+-------------------------
table a c b, cont(sum y)
----------+-------------------------
| b and c
| ---- 1 --- ---- 2 ---
a | 1 2 1 2
----------+-------------------------
1 | 15 16 28 32
2 | 7 12 22 40
----------+-------------------------
table a c b, cont(sum y)
----------+-------------------------
| b and c
| ---- 1 --- ---- 2 ---
a | 1 2 1 2
----------+-------------------------
1 | 15 16 28 32
2 | 7 12 22 40
----------+-------------------------
Table 12.10-2, page 559.
anova y a / s|a b a*b / b*s|a c a*c / c*s|a b*c a*b*c /, repeated(b c)
Number of obs = 32 R-squared = 0.9920
Root MSE = .559017 Adj R-squared = 0.9589
Source | Partial SS df MS F Prob > F
-----------+----------------------------------------------------
Model | 233.625 25 9.345 29.90 0.0002
|
a | 3.125 1 3.125 2.00 0.2070
s|a | 9.375 6 1.5625
-----------+----------------------------------------------------
b | 162.00 1 162.00 199.38 0.0000
a*b | 6.125 1 6.125 7.54 0.0335
b*s|a | 4.875 6 .8125
-----------+----------------------------------------------------
c | 24.50 1 24.50 61.89 0.0002
a*c | 10.125 1 10.125 25.58 0.0023
c*s|a | 2.375 6 .395833333
-----------+----------------------------------------------------
b*c | 8.00 1 8.00 25.60 0.0023
a*b*c | 3.125 1 3.125 10.00 0.0195
|
Residual | 1.875 6 .3125
-----------+----------------------------------------------------
Total | 235.50 31 7.59677419
Between-subjects error term: s|a
Levels: 8 (6 df)
Lowest b.s.e. variable: s
Covariance pooled over: a (for repeated variables)
Repeated variable: b
Huynh-Feldt epsilon = 1.2000
*Huynh-Feldt epsilon reset to 1.0000
Greenhouse-Geisser epsilon = 1.0000
Box's conservative epsilon = 1.0000
------------ Prob > F ------------
Source | df F Regular H-F G-G Box
-----------+----------------------------------------------------
b | 1 199.38 0.0000 0.0000 0.0000 0.0000
a*b | 1 7.54 0.0335 0.0335 0.0335 0.0335
b*s|a | 6
-----------+----------------------------------------------------
Repeated variable: c
Huynh-Feldt epsilon = 1.2000
*Huynh-Feldt epsilon reset to 1.0000
Greenhouse-Geisser epsilon = 1.0000
Box's conservative epsilon = 1.0000
------------ Prob > F ------------
Source | df F Regular H-F G-G Box
-----------+----------------------------------------------------
c | 1 61.89 0.0002 0.0002 0.0002 0.0002
a*c | 1 25.58 0.0023 0.0023 0.0023 0.0023
c*s|a | 6
-----------+----------------------------------------------------
Repeated variables: b*c
Huynh-Feldt epsilon = 1.2000
*Huynh-Feldt epsilon reset to 1.0000
Greenhouse-Geisser epsilon = 1.0000
Box's conservative epsilon = 1.0000
------------ Prob > F ------------
Source | df F Regular H-F G-G Box
-----------+----------------------------------------------------
b*c | 1 25.60 0.0023 0.0023 0.0023 0.0023
a*b*c | 1 10.00 0.0195 0.0195 0.0195 0.0195
Residual | 6
-----------+----------------------------------------------------
matrix list e(Srep)
symmetric e(Srep)[4,4]
c1 c2 c3 c4
r1 1.25
r2 .66666667 .66666667
r3 .58333333 .33333333 .5
r4 0 -.16666667 .16666667 .66666667
Table 12.14-1, page 574.
use https://stats.idre.ucla.edu/stat/stata/examples/kirk/spf2-3a, clear
tabdisp s b, cellvar(y) by(a) concise
----------+-----------------
| b
a and s | 1 2 3
----------+-----------------
1 |
1 | 3 4 7
2 | 6 5 8
----------+-----------------
2 |
3 | 3 4 7
4 | 3 3 6
5 | 1 2 5
----------+-----------------
Table 12.14-2, page 576.
anova y a / s|a b a*b /, repeated(b)
Number of obs = 15 R-squared = 0.9692
Root MSE = .544331 Adj R-squared = 0.9281
Source | Partial SS df MS F Prob > F
-----------+----------------------------------------------------
Model | 55.9555556 8 6.99444444 23.61 0.0005
|
a | 10.6777778 1 10.6777778 3.08 0.1774
s|a | 10.3888889 3 3.46296296
-----------+----------------------------------------------------
b | 32.3555556 2 16.1777778 54.60 0.0001
a*b | .355555556 2 .177777778 0.60 0.5787
|
Residual | 1.77777778 6 .296296296
-----------+----------------------------------------------------
Total | 57.7333333 14 4.12380952
Between-subjects error term: s|a
Levels: 5 (3 df)
Lowest b.s.e. variable: s
Covariance pooled over: a (for repeated variable)
Repeated variable: b
Huynh-Feldt epsilon = 0.7500
Greenhouse-Geisser epsilon = 0.5000
Box's conservative epsilon = 0.5000
------------ Prob > F ------------
Source | df F Regular H-F G-G Box
-----------+----------------------------------------------------
b | 2 54.60 0.0001 0.0008 0.0051 0.0051
a*b | 2 0.60 0.5787 0.5431 0.4950 0.4950
Residual | 6
-----------+----------------------------------------------------
matrix list e(Srep)
symmetric e(Srep)[3,3]
c1 c2 c3
r1 2.3888889
r2 1.1666667 .83333333
r3 1.1666667 .83333333 .83333333