Differences

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

--- twoway_repeated_measure_anova [2023/04/30 23:36] – hkimscil
+++ twoway_repeated_measure_anova [2024/05/08 08:33] (current) – removed hkimscil
@@ Line 1: / Line 1: @@
-{{:r:10_rmanova.csv}}
-<code>
-acne <- read.csv("http://commres.net/wiki/_media/r/10_rmanova.csv")
-str(acne)
-acne
-# install.packages("reshape")
-library(reshape)
-acne.re <- reshape(acne, direction="long", varying=3:6, sep="")
-str(acne.re)
-acne.re$group <- factor(acne.re$group)
-acne.re$id <- factor(acne.re$id)
-acne.re$time <- factor(acne.re$time)
-str(acne.re)
-acne.re
-attach(acne.re)
-acne.re.anova <- aov(month~group*time+Error(id/time), data=acne.re)
-summary(acne.re.anova)
-interaction.plot(acne.re$time, acne.re$group, acne.re$month)
-</code>
-{{:pasted:20220510-102845.png}}
-<code>
-> acne <- read.csv("http://commres.net/wiki/_media/r/10_rmanova.csv")
-> str(acne)
-'data.frame':	14 obs. of  6 variables:
- $ group : int  1 1 1 1 1 1 1 2 2 2 ...
- $ id    : int  1 2 3 4 5 6 7 8 9 10 ...
- $ month0: int  60 52 62 58 65 58 53 55 55 60 ...
- $ month1: int  41 38 36 34 34 42 38 42 54 55 ...
- $ month3: int  25 23 22 21 28 26 25 33 46 46 ...
- $ month6: int  16 12 14 13 18 16 21 22 26 23 ...
-> acne
-   group id month0 month1 month3 month6
-      1  1     60     41     25     16
-      1  2     52     38     23     12
-      1  3     62     36     22     14
-      1  4     58     34     21     13
-      1  5     65     34     28     18
-      1  6     58     42     26     16
-      1  7     53     38     25     21
-      2  8     55     42     33     22
-      2  9     55     54     46     26
-     2 10     60     55     46     23
-     2 11     63     45     40     25
-     2 12     52     41     35     22
-     2 13     61     38     32     18
-     2 14     58     43     39     21
->
-> # install.packages("reshape")
-> library(reshape)
-Warning message:
-패키지 ‘reshape’는 R 버전 4.1.3에서 작성되었습니다
->
-> acne.re <- reshape(acne, direction="long", varying=3:6, sep="")
->
-> str(acne.re)
-'data.frame':	56 obs. of  4 variables:
- $ group: int  1 1 1 1 1 1 1 2 2 2 ...
- $ id   : int  1 2 3 4 5 6 7 8 9 10 ...
- $ time : num  0 0 0 0 0 0 0 0 0 0 ...
- $ month: int  60 52 62 58 65 58 53 55 55 60 ...
- - attr(*, "reshapeLong")=List of 4
-  ..$ varying:List of 1
-  .. ..$ month: chr [1:4] "month0" "month1" "month3" "month6"
-  .. ..- attr(*, "v.names")= chr "month"
-  .. ..- attr(*, "times")= num [1:4] 0 1 3 6
-  ..$ v.names: chr "month"
-  ..$ idvar  : chr "id"
-  ..$ timevar: chr "time"
-> acne.re$group <- factor(acne.re$group)
-> acne.re$id <- factor(acne.re$id)
-> acne.re$time <- factor(acne.re$time)
-> str(acne.re)
-'data.frame':	56 obs. of  4 variables:
- $ group: Factor w/ 2 levels "1","2": 1 1 1 1 1 1 1 2 2 2 ...
- $ id   : Factor w/ 14 levels "1","2","3","4",..: 1 2 3 4 5 6 7 8 9 10 ...
- $ time : Factor w/ 4 levels "0","1","3","6": 1 1 1 1 1 1 1 1 1 1 ...
- $ month: int  60 52 62 58 65 58 53 55 55 60 ...
- - attr(*, "reshapeLong")=List of 4
-  ..$ varying:List of 1
-  .. ..$ month: chr [1:4] "month0" "month1" "month3" "month6"
-  .. ..- attr(*, "v.names")= chr "month"
-  .. ..- attr(*, "times")= num [1:4] 0 1 3 6
-  ..$ v.names: chr "month"
-  ..$ idvar  : chr "id"
-  ..$ timevar: chr "time"
->
-> acne.re
-     group id time month
-.0      1  1    0    60
-.0      1  2    0    52
-.0      1  3    0    62
-.0      1  4    0    58
-.0      1  5    0    65
-.0      1  6    0    58
-.0      1  7    0    53
-.0      2  8    0    55
-.0      2  9    0    55
-.0     2 10    0    60
-.0     2 11    0    63
-.0     2 12    0    52
-.0     2 13    0    61
-.0     2 14    0    58
-.1      1  1    1    41
-.1      1  2    1    38
-.1      1  3    1    36
-.1      1  4    1    34
-.1      1  5    1    34
-.1      1  6    1    42
-.1      1  7    1    38
-.1      2  8    1    42
-.1      2  9    1    54
-.1     2 10    1    55
-.1     2 11    1    45
-.1     2 12    1    41
-.1     2 13    1    38
-.1     2 14    1    43
-.3      1  1    3    25
-.3      1  2    3    23
-.3      1  3    3    22
-.3      1  4    3    21
-.3      1  5    3    28
-.3      1  6    3    26
-.3      1  7    3    25
-.3      2  8    3    33
-.3      2  9    3    46
-.3     2 10    3    46
-.3     2 11    3    40
-.3     2 12    3    35
-.3     2 13    3    32
-.3     2 14    3    39
-.6      1  1    6    16
-.6      1  2    6    12
-.6      1  3    6    14
-.6      1  4    6    13
-.6      1  5    6    18
-.6      1  6    6    16
-.6      1  7    6    21
-.6      2  8    6    22
-.6      2  9    6    26
-.6     2 10    6    23
-.6     2 11    6    25
-.6     2 12    6    22
-.6     2 13    6    18
-.6     2 14    6    21
->
-> attach(acne.re)
-The following objects are masked from acne.re (pos = 4):
-    group, id, month, time
-> acne.re.anova <- aov(month~group*time+Error(id/time), data=acne.re)
-> summary(acne.re.anova)
-Error: id
-          Df Sum Sq Mean Sq F value   Pr(>F)
-group      1  707.2   707.2   19.71 0.000808 ***
-Residuals 12  430.6    35.9
----
-Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
-Error: id:time
-           Df Sum Sq Mean Sq F value  Pr(>F)
-time        3  11366    3789  308.78 < 2e-16 ***
-group:time  3    396     132   10.77 3.4e-05 ***
-Residuals  36    442      12
----
-Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
->
-> interaction.plot(acne.re$time, acne.re$group, acne.re$month)
->
-</code>