Differences

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

--- repeated_measure_anova [2023/05/01 00:15] – hkimscil
+++ repeated_measure_anova [2023/05/03 08:04] – [demo 2] hkimscil
@@ Line 257: / Line 257: @@
 see {{:r:repeated_measures_anova_eg.xlsx}}
+===== demo 2 =====
+<code>
+# create data
+df <- data.frame(patient=rep(1:5, each=4),
+                 drug=rep(1:4, times=5),
+                 response=c(30, 28, 16, 34,
+, 18, 10, 22,
+, 20, 18, 30,
+, 34, 20, 44,
+, 28, 14, 30))
+# view data
+df
+# within sujbect anova
+within.aov.mod <- aov(response~drug+Error(patient), data=df)
+</code>
+<code>
+> #create data
+> df <- data.frame(patient=rep(1:5, each=4),
++                  drug=rep(1:4, times=5),
++                  response=c(30, 28, 16, 34,
++                             14, 18, 10, 22,
++                             24, 20, 18, 30,
++                             38, 34, 20, 44,
++                             26, 28, 14, 30))
+>
+> #view data
+> df
+   patient drug response
+        1    1       30
+        1    2       28
+        1    3       16
+        1    4       34
+        2    1       14
+        2    2       18
+        2    3       10
+        2    4       22
+        3    1       24
+       3    2       20
+       3    3       18
+       3    4       30
+       4    1       38
+       4    2       34
+       4    3       20
+       4    4       44
+       5    1       26
+       5    2       28
+       5    3       14
+       5    4       30
+>
+> #within sujbect anova
+> within.aov.mod <- aov(response~drug+Error(patient), data=df)
+>
+>
+> summary(within.aov.mod)
+Error: patient
+          Df Sum Sq Mean Sq F value Pr(>F)
+Residuals  1   67.6    67.6
+Error: Within
+          Df Sum Sq Mean Sq F value Pr(>F)
+drug       1   11.6   11.56   0.139  0.714
+Residuals 17 1412.6   83.10
+</code>
+The above is <fc #ff0000>**WRONG**</fc>.
+<code>
+within.aov.mod <- aov(response~factor(drug)+Error(factor(patient)), data=df)
+summary(within.aov.mod)
+</code>
+<code>
+> within.aov.mod <- aov(response~factor(drug)+Error(factor(patient)), data=df)
+>
+>
+> summary(within.aov.mod)
+Error: factor(patient)
+          Df Sum Sq Mean Sq F value Pr(>F)
+Residuals  4  680.8   170.2
+Error: Within
+             Df Sum Sq Mean Sq F value   Pr(>F)
+factor(drug)  3  698.2   232.7   24.76 1.99e-05 ***
+Residuals    12  112.8     9.4
+---
+Signif. codes:
+‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
+>
+</code>
 ====== two way ======
 <code>