Differences

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

--- r:repeated_measures_anova [2020/06/04 00:49] – hkimscil
+++ r:repeated_measures_anova [2020/06/04 00:53] – [e.g. 3] hkimscil
@@ Line 167: / Line 167: @@
 >
 </code>
+====== e.g. 3 ======
+<code>
+demo3 <- read.csv("https://stats.idre.ucla.edu/stat/data/demo3.csv")
+## Convert variables to factor
+demo3 <- within(demo3, {
+  group <- factor(group)
+  time <- factor(time)
+  id <- factor(id)
+})
+demo3
+par(cex = .6)
+with(demo3, interaction.plot(time, group, pulse,
+  ylim = c(10, 60), lty = c(1, 12), lwd = 3,
+  ylab = "mean of pulse", xlab = "time", trace.label = "group"))
+demo3.aov <- aov(pulse ~ group * time + Error(id), data = demo3)
+summary(demo3.aov)
+</code>
+{{:r:pasted:20200604-005114.png}}
+<code>
+> demo3 <- read.csv("https://stats.idre.ucla.edu/stat/data/demo3.csv")
+> ## Convert variables to factor
+> demo3 <- within(demo3, {
++     group <- factor(group)
++     time <- factor(time)
++     id <- factor(id)
++ })
+>
+> demo3
+   id group pulse time
+   1     1    35    1
+   1     1    25    2
+   1     1    16    3
+   2     1    32    1
+   2     1    23    2
+   2     1    12    3
+   3     1    36    1
+   3     1    22    2
+   3     1    14    3
+  4     1    34    1
+  4     1    21    2
+  4     1    13    3
+  5     2    57    1
+  5     2    43    2
+  5     2    22    3
+  6     2    54    1
+  6     2    46    2
+  6     2    26    3
+  7     2    55    1
+  7     2    46    2
+  7     2    23    3
+  8     2    60    1
+  8     2    47    2
+  8     2    25    3
+>
+> par(cex = .6)
+>
+> with(demo3, interaction.plot(time, group, pulse,
++                              ylim = c(10, 60), lty = c(1, 12), lwd = 3,
++                              ylab = "mean of pulse", xlab = "time", trace.label = "group"))
+>
+> demo3.aov <- aov(pulse ~ group * time + Error(id), data = demo3)
+> summary(demo3.aov)
+Error: id
+          Df Sum Sq Mean Sq F value  Pr(>F)
+group      1 2035.0  2035.0   343.1 1.6e-06 ***
+Residuals  6   35.6     5.9
+---
+Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
+Error: Within
+           Df Sum Sq Mean Sq F value   Pr(>F)
+time        2 2830.3  1415.2   553.8 1.52e-12 ***
+group:time  2  200.3   100.2    39.2 5.47e-06 ***
+Residuals  12   30.7     2.6
+---
+Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
+>
+</code>
+====== e.g. 4 ======
+<code>
+demo4 <- read.csv("https://stats.idre.ucla.edu/stat/data/demo4.csv")
+## Convert variables to factor
+demo4 <- within(demo4, {
+  group <- factor(group)
+  time <- factor(time)
+  id <- factor(id)
+})
+demo4
+par(cex = .6)
+with(demo4, interaction.plot(time, group, pulse,
+  ylim = c(10, 60), lty = c(1, 12), lwd = 3,
+  ylab = "mean of pulse", xlab = "time", trace.label = "group"))
+demo4.aov <- aov(pulse ~ group * time + Error(id), data = demo4)
+summary(demo4.aov)
+</code>
+<code>
+>
+> demo4 <- read.csv("https://stats.idre.ucla.edu/stat/data/demo4.csv")
+> ## Convert variables to factor
+> demo4 <- within(demo4, {
++     group <- factor(group)
++     time <- factor(time)
++     id <- factor(id)
++ })
+>
+> demo4
+   id group pulse time
+   1     1    35    1
+   1     1    25    2
+   1     1    12    3
+   2     1    34    1
+   2     1    22    2
+   2     1    13    3
+   3     1    36    1
+   3     1    21    2
+   3     1    18    3
+  4     1    35    1
+  4     1    23    2
+  4     1    15    3
+  5     2    31    1
+  5     2    43    2
+  5     2    57    3
+  6     2    35    1
+  6     2    46    2
+  6     2    58    3
+  7     2    37    1
+  7     2    48    2
+  7     2    51    3
+  8     2    32    1
+  8     2    45    2
+  8     2    53    3
+>
+> par(cex = .6)
+>
+> with(demo4, interaction.plot(time, group, pulse,
++                              ylim = c(10, 60), lty = c(1, 12), lwd = 3,
++                              ylab = "mean of pulse", xlab = "time", trace.label = "group"))
+>
+> demo4.aov <- aov(pulse ~ group * time + Error(id), data = demo4)
+> summary(demo4.aov)
+Error: id
+          Df Sum Sq Mean Sq F value   Pr(>F)
+group      1 2542.0    2542     629 2.65e-07 ***
+Residuals  6   24.3       4
+---
+Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
+Error: Within
+           Df Sum Sq Mean Sq F value   Pr(>F)
+time        2      1     0.5   0.079    0.925
+group:time  2   1736   868.2 137.079 5.44e-09 ***
+Residuals  12     76     6.3
+---
+Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
+>
+</code>