Differences

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

--- repeated_measure_anova [2024/05/08 08:16] – [demo 2] hkimscil
+++ repeated_measure_anova [2024/05/13 08:39] (current) – [Repeated Measure ANOVA] hkimscil
@@ Line 81: / Line 81: @@
 {{:pasted:20240501-082722.png}}
 ----
-{{:pasted:20240501-082738.png}}
+{{:pasted:20240513-083858.png}}
 ----
   * but, $\text{SS}_\text{{within}}$ can be partitioned as
@@ Line 262: / Line 262: @@
 ===== demo 2 =====
 see [[:r:repeated measure anova]]
-====== two way ======
+===== Twoway repeated measure anova=====
-<code>
+see [[:r:twoway repeated measure anova]]
-demo1  <- read.csv("https://stats.idre.ucla.edu/stat/data/demo1.csv")
-demo1
-str(demo1) ## 모든 변인이 int이므로 (숫자) factor로 바꿔야 한다
-## Convert variables to factor
-demo1 <- within(demo1, {
-  group <- factor(group)
-  time <- factor(time)
-  id <- factor(id)
-}) ## 이제 pulse만 제외하고 모두 factor로 변환된 데이터
-str(demo1)
-par(cex = .6)
-with(demo1, interaction.plot(time, group, pulse,
-  ylim = c(5, 20), lty= c(1, 12), lwd = 3,
-  ylab = "mean of pulse", xlab = "time", trace.label = "group"))
-demo1.aov <- aov(pulse ~ group * time + Error(id), data = demo1)
-summary(demo1.aov)
-</code>
-<code>
-> summary(demo1.aov)
-Error: id
-          Df Sum Sq Mean Sq F value  Pr(>F)
-group      1 155.04  155.04    3721 1.3e-09 ***
-Residuals  6   0.25    0.04
----
-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 0.0833 0.04167       1  0.397
-group:time  2 0.0833 0.04167       1  0.397
-Residuals  12 0.5000 0.04167
-</code>
-{{:pasted:20200611-142331.png?350}}
-===== demo2 =====
-<code>
-demo2 <- read.csv("https://stats.idre.ucla.edu/stat/data/demo2.csv")
-## Convert variables to factor
-demo2 <- within(demo2, {
-    group <- factor(group)
-    time <- factor(time)
-    id <- factor(id)
-})
-demo2
-with(demo2, interaction.plot(time, group, pulse,
- ylim = c(10, 40), lty = c(1, 12), lwd = 3,
- ylab = "mean of pulse", xlab = "time", trace.label = "group"))
-demo2.aov <- aov(pulse ~ group * time + Error(id), data = demo2)
-summary(demo2.aov)
-</code>
-{{:pasted:20200611-151520.png?350}}
-<code>
-> demo2 <- read.csv("https://stats.idre.ucla.edu/stat/data/demo2.csv")
-> ## Convert variables to factor
-> demo2 <- within(demo2, {
-+     group <- factor(group)
-+     time <- factor(time)
-+     id <- factor(id)
-+ })
-> demo2
-   id group pulse time
-   1     1    14    1
-   1     1    19    2
-   1     1    29    3
-   2     1    15    1
-   2     1    25    2
-   2     1    26    3
-   3     1    16    1
-   3     1    16    2
-   3     1    31    3
-  4     1    12    1
-  4     1    24    2
-  4     1    32    3
-  5     2    10    1
-  5     2    21    2
-  5     2    24    3
-  6     2    17    1
-  6     2    26    2
-  6     2    35    3
-  7     2    19    1
-  7     2    22    2
-  7     2    32    3
-  8     2    15    1
-  8     2    23    2
-  8     2    34    3
->
-> with(demo2, interaction.plot(time, group, pulse,
-+  ylim = c(10, 40), lty = c(1, 12), lwd = 3,
-+  ylab = "mean of pulse", xlab = "time", trace.label = "group"))
->
-> demo2.aov <- aov(pulse ~ group * time + Error(id), data = demo2)
-> summary(demo2.aov)
-Error: id
-          Df Sum Sq Mean Sq F value Pr(>F)
-group      1  15.04   15.04   0.836  0.396
-Residuals  6 107.92   17.99
-Error: Within
-           Df Sum Sq Mean Sq F value   Pr(>F)
-time        2  978.2   489.1  53.684 1.03e-06 ***
-group:time  2    1.1     0.5   0.059    0.943
-Residuals  12  109.3     9.1
----
-Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
->
-</code>
-===== demo 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)
-})
-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>
-{{:pasted:20200611-151755.png?350}}
-<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)
-+ })
->
-> 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>
 ====== reference ======
   * [[http://wwwstage.valpo.edu/other/dabook/ch12/c12-1.htm|Repeated measures one-way ANOVA]] by Akkelin