Differences

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

--- repeated_measure_anova [2023/05/03 08:00] – [demo1] hkimscil
+++ repeated_measure_anova [2024/05/08 08:16] – [demo 2] hkimscil
@@ Line 79: / Line 79: @@
 -- Picture about here --
+{{:pasted:20240501-082722.png}}
+----
+{{:pasted:20240501-082738.png}}
+----
   * but, $\text{SS}_\text{{within}}$ can be partitioned as
     * $\text{SS}_{\text{ subjects}}$ and $\text{SS}_{\text{ error}}$
@@ Line 258: / Line 261: @@
 see {{:r:repeated_measures_anova_eg.xlsx}}
 ===== demo 2 =====
-<code>
+see [[:r:repeated measure anova]]
-# 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>
 ====== two way ======
 <code>