Table 8.3-1, page 325.
use https://stats.idre.ucla.edu/stat/stata/examples/kirk/ls4, clear tabdisp s c, cellvar(y) by(b) /* levels of a not shown */ ----------+----------------------- | c b and s | 1 2 3 4 ----------+----------------------- 1 | 1 | 1 2 5 9 2 | 2 3 6 8 ----------+----------------------- 2 | 1 | 3 8 9 2 2 | 4 6 8 3 ----------+----------------------- 3 | 1 | 5 10 3 5 2 | 7 11 2 4 ----------+----------------------- 4 | 1 | 7 6 3 6 2 | 10 3 4 7 ----------+----------------------- table b c, cont(sum y) row col /* levels of a not shown */ ----------+---------------------------------- | c b | 1 2 3 4 Total ----------+---------------------------------- 1 | 3 5 11 17 36 2 | 7 14 17 5 43 3 | 12 21 5 9 47 4 | 17 9 7 13 46 | Total | 39 49 40 44 172 ----------+---------------------------------- table a, cont(sum y) ----------+----------- a | sum(y) ----------+----------- 1 | 22 2 | 28 3 | 50 4 | 72 ----------+-----------
Table 8.3-2, page 327.
Note: The term that Kirk calls residual is called a*b*c in this Stata model.
anova y a b c a*b*c Number of obs = 32 R-squared = 0.9193 Root MSE = 1.08972 Adj R-squared = 0.8437 Source | Partial SS df MS F Prob > F -----------+---------------------------------------------------- Model | 216.50 15 14.4333333 12.15 0.0000 | a | 194.50 3 64.8333333 54.60 0.0000 b | 9.25 3 3.08333333 2.60 0.0884 c | 7.75 3 2.58333333 2.18 0.1308 a*b*c | 5.00 6 .833333333 0.70 0.6525 | Residual | 19.00 16 1.1875 -----------+---------------------------------------------------- Total | 235.50 31 7.59677419
Omega-squared computation, page 329.
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 54.6 3 omega squared = 0.8340 effect size = 2.2417
Table 8.4-1, page 331.
use https://stats.idre.ucla.edu/stat/stata/examples/kirk/ls4a, clear tabdisp b c, cellvar(y) ----------+----------------------- | c b | 1 2 3 4 ----------+----------------------- 1 | 3 5 11 17 2 | 7 14 17 5 3 | 12 21 5 9 4 | 17 9 7 13 ----------+-----------------------
Table 8.4-1, page 332.
nonadd y a b c Tukey's test of nonadditivity for randomized block designs F (1,5) = .0001702 Pr > F: .99009558
Figure 8.4-1, page 333.
anova y a b c, noanova predict yhat /* yhat is the fitted value */ predict e, rstandard /* e is the standardized residual */ graph twoway scatter e yhat, ylabel(-2(.5)2) xlabel(0(2)20)
Table 8.11-1, page 353.
use https://stats.idre.ucla.edu/stat/stata/examples/kirk/ls4a, clear tabdisp s b, cellvar(y) /* levels of a not shown */ ----------+----------------- | b s | 1 2 3 ----------+----------------- 1 | 7 14 12 2 | 3 5 11 3 | 6 7 11 4 | 12 13 9 5 | 7 9 8 6 | 8 13 9 ----------+-----------------
Table 8.11-2, page 355.
anova y a b s Number of obs = 18 R-squared = 0.9021 Root MSE = 1.38944 Adj R-squared = 0.7920 Source | Partial SS df MS F Prob > F -----------+---------------------------------------------------- Model | 142.333333 9 15.8148148 8.19 0.0035 | a | 49.7777778 2 24.8888889 12.89 0.0031 b | 34.1111111 2 17.0555556 8.83 0.0094 s | 58.4444444 5 11.6888889 6.05 0.0131 | Residual | 15.4444444 8 1.93055556 -----------+---------------------------------------------------- Total | 157.777778 17 9.28104575
Omega-squared computation, page 354.
omega2 12.89 2 omega squared = 0.5692 effect size = 1.1494
Table 8.11-3, page 355.
Note: The fhcomp 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).
fhcomp a Fisher-Hayter pairwise comparisons for variable a studentized range critical value(.05, 2, 8) = 3.2611823 mean critical grp vs grp group means dif dif ------------------------------------------------------- 1 vs 2 7.3333 11.3333 4.0000* 1.8499 1 vs 3 7.3333 8.6667 1.3333 1.8499 2 vs 3 11.3333 8.6667 2.6667* 1.8499