Differences

This shows you the differences between two versions of the page.

--- r:factorial_anova [2026/04/28 13:20] – hkimscil
+++ r:factorial_anova [2026/04/28 13:32] (current) – [output interpretation] hkimscil
@@ Line 708: / Line 708: @@
 {{.:pasted:20260428-131558.png}}
 <code>
-> m.aov <- aov(y~a*b, data = dat)
+> m.lm <- (lm(y~a*b, data = dat))
-> summary(m.aov)
+> summary(m.lm)
-            Df Sum Sq Mean Sq F value   Pr(>F)
-a            2  12.57   6.283   6.081 0.007303 **
+Call:
-b            1  12.11  12.107  11.718 0.002227 **
+lm(formula = y ~ a * b, data = dat)
-a:b          2  25.74  12.870  12.456 0.000195 ***
-Residuals   24  24.80   1.033
+Residuals:
+    Min      1Q  Median      3Q     Max
+-1.9469 -0.4908 -0.1742  0.5814  2.0851
+Coefficients:
+             Estimate Std. Error t value Pr(>|t|)
+(Intercept)    9.6868     0.4546  21.309  < 2e-16 ***
+adrg2          0.8050     0.6429   1.252   0.2226
+adrg3          0.6943     0.6429   1.080   0.2909
+bexerc         3.3284     0.6429   5.177 2.65e-05 ***
+adrg2:bexerc  -4.4910     0.9092  -4.940 4.84e-05 ***
+adrg3:bexerc  -1.6826     0.9092  -1.851   0.0766 .
 ---
 Signif. codes:
 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
+Residual standard error: 1.016 on 24 degrees of freedom
+Multiple R-squared:  0.6703,	Adjusted R-squared:  0.6016
+F-statistic: 9.758 on 5 and 24 DF,  p-value: 3.474e-05
+> anova(m.lm)
+Analysis of Variance Table
+Response: y
+          Df Sum Sq Mean Sq F value    Pr(>F)
+a          2 12.566  6.2832  6.0811 0.0073026 **
+b          1 12.107 12.1069 11.7176 0.0022274 **
+a:b        2 25.740 12.8699 12.4561 0.0001948 ***
+Residuals 24 24.797  1.0332
+---
+Signif. codes:
+‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
+>
+> ss.tot <- sum(anova(m.lm)[2])
+> ss.a <- anova(m.lm)[[2]][1]
+> ss.b <- anova(m.lm)[[2]][2]
+> ss.ab <- anova(m.lm)[[2]][3]
+> ss.bet <- ss.a + ss.b + ss.ab
+> ss.bet
+[1] 50.41308
+> ss.tot
+[1] 75.21048
+> ss.bet / ss.tot
+[1] 0.6702933
 >
 </code>