Differences

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

--- factorial_anova [2020/10/12 17:34] – [예] hkimscil
+++ factorial_anova [2024/05/01 08:36] (current) – [Materials and links] hkimscil
@@ Line 290: / Line 290: @@
 <code>
-de <- read.csv("http://commres.net/wiki/_media/detergent.csv", sep = "\t", header=T)
+de <- read.csv("http://commres.net/wiki/_media/detergent.csv", sep = ",", header=T)
 de
@@ Line 309: / Line 309: @@
 </code>
-{{:pasted:20200529-165842.png}}
+<code>
+> de <- read.csv("http://commres.net/wiki/_media/detergent.csv", sep = ",", header=T)
+> de
+   type w.temp cleanness
+     1      1         4
+     1      1         5
+     1      1         6
+     1      1         5
+     2      1         6
+     2      1         6
+     2      1         4
+     2      1         4
+     1      2         7
+    1      2         9
+    1      2         8
+    1      2        10
+    2      2        13
+    2      2        15
+    2      2        12
+    2      2        12
+    1      3        10
+    1      3        12
+    1      3        11
+    1      3         9
+    2      3        12
+    2      3        13
+    2      3        10
+    2      3        13
+> de$type <- factor(de$type, level=c(1,2), label=c("super", "best"))
+> de$w.temp <- factor(de$w.temp, level=c(1,2,3), label=c("cold", "warm", "hot"))
+> de
+    type w.temp cleanness
+  super   cold         4
+  super   cold         5
+  super   cold         6
+  super   cold         5
+   best   cold         6
+   best   cold         6
+   best   cold         4
+   best   cold         4
+  super   warm         7
+super   warm         9
+super   warm         8
+super   warm        10
+  best   warm        13
+  best   warm        15
+  best   warm        12
+  best   warm        12
+super    hot        10
+super    hot        12
+super    hot        11
+super    hot         9
+  best    hot        12
+  best    hot        13
+  best    hot        10
+  best    hot        13
+> de.anova <- aov(cleanness ~ type * w.temp, data=de)
+> summary(de.anova)
+            Df Sum Sq Mean Sq F value   Pr(>F)
+type         1     24   24.00   15.43 0.000986 ***
+w.temp       2    193   96.50   62.04 8.41e-09 ***
+type:w.temp  2     21   10.50    6.75 0.006496 **
+Residuals   18     28    1.56
+---
+Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
+>
+> with(de, interaction.plot(x.factor=type,
++                           trace.factor=w.temp, response=cleanness,
++                           fun=mean, type="b", legend=T,
++                           ylab="cleanness", main="Interaction Plot (type by temp)",
++                           pch=c(1,19)))
+>
+>
+</code>
+{{:pasted:20240429-083635.png}}
+만약에 손으로 계산했다면 R에서
+<code>
+pf(15.43, 1, 18, lower.tail=F)
+pf(62.04, 2, 18, lower.tail=F)
+pf(6.75, 2, 18, lower.tail=F)
+</code>
+<code>
+> pf(15.43, 1, 18, lower.tail=F)
+[1] 0.000985713
+> pf(62.04, 2, 18, lower.tail=F)
+[1] 8.40729e-09
+> pf(6.75, 2, 18, lower.tail=F)
+[1] 0.006496173
+>
+> </code>
 ===== 예 2.  cookie experiment =====
@@ Line 395: / Line 487: @@
 | a R Squared = .102 (Adjusted R Squared = .066)  ||||||
+데이터 파일
+{{:r:cookies.csv}}
+손으로 계산하기
+{{:r:cookies.xlsx}}
-{{http://commres.net/wiki/_media/cookies.csv|cookies.csv}} 데이터파일 위치
 <code>
 cookies <- read.csv("http://commres.net/wiki/_media/cookies.csv")
@@ Line 664: / Line 759: @@
 ====== Materials and links ======
   * {{:AnovaData.sav}}
-  * http://www.uwsp.edu/psych/stat/13/anova-2w.htm
-  * http://faculty.vassar.edu/lowry/ch16pt1.html
-  * http://faculty.vassar.edu/lowry/ch16pt2.html
-  * http://faculty.vassar.edu/lowry/ch16pt3.html
-  * http://faculty.vassar.edu/lowry/ch16pt4.html
-  * http://faculty.vassar.edu/lowry/ch16pt5.html
-{{tag>factorial anova, statistics, Two-Factor Analysis of Variance, 팩토리얼 아노바, 상호작용효과, 주효과}}
+{{tag>factorial_anova statistics Two-Factor_Analysis_of_Variance, 팩토리얼 아노바, 상호작용효과, 주효과}}