COMMunication
RESearch.NET

User Tools

Site Tools

You are here: index » twoway_repeated_measure_anova
Trace:
twoway_repeated_measure_anova

Differences

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

Link to this comparison view

Both sides previous revisionPrevious revision
Next revision		Previous revision
Next revisionBoth sides next revision
twoway_repeated_measure_anova [2022/10/04 08:56] – hkimsciltwoway_repeated_measure_anova [2023/04/30 23:36] – hkimscil
Line 17: Line 17:
 acne.re$time <- factor(acne.re$time) acne.re$time <- factor(acne.re$time)
 str(acne.re) str(acne.re)
 +
 +acne.re
  
 attach(acne.re) attach(acne.re)
Line 27: Line 29:
  
 {{:pasted:20220510-102845.png}} {{: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  1     60     41     25     16
 +2      1  2     52     38     23     12
 +3      1  3     62     36     22     14
 +4      1  4     58     34     21     13
 +5      1  5     65     34     28     18
 +6      1  6     58     42     26     16
 +7      1  7     53     38     25     21
 +8      2  8     55     42     33     22
 +9      2  9     55     54     46     26
 +10     2 10     60     55     46     23
 +11     2 11     63     45     40     25
 +12     2 12     52     41     35     22
 +13     2 13     61     38     32     18
 +14     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
 +1.0      1  1    0    60
 +2.0      1  2    0    52
 +3.0      1  3    0    62
 +4.0      1  4    0    58
 +5.0      1  5    0    65
 +6.0      1  6    0    58
 +7.0      1  7    0    53
 +8.0      2  8    0    55
 +9.0      2  9    0    55
 +10.0     2 10    0    60
 +11.0     2 11    0    63
 +12.0     2 12    0    52
 +13.0     2 13    0    61
 +14.0     2 14    0    58
 +1.1      1  1    1    41
 +2.1      1  2    1    38
 +3.1      1  3    1    36
 +4.1      1  4    1    34
 +5.1      1  5    1    34
 +6.1      1  6    1    42
 +7.1      1  7    1    38
 +8.1      2  8    1    42
 +9.1      2  9    1    54
 +10.1     2 10    1    55
 +11.1     2 11    1    45
 +12.1     2 12    1    41
 +13.1     2 13    1    38
 +14.1     2 14    1    43
 +1.3      1  1    3    25
 +2.3      1  2    3    23
 +3.3      1  3    3    22
 +4.3      1  4    3    21
 +5.3      1  5    3    28
 +6.3      1  6    3    26
 +7.3      1  7    3    25
 +8.3      2  8    3    33
 +9.3      2  9    3    46
 +10.3     2 10    3    46
 +11.3     2 11    3    40
 +12.3     2 12    3    35
 +13.3     2 13    3    32
 +14.3     2 14    3    39
 +1.6      1  1    6    16
 +2.6      1  2    6    12
 +3.6      1  3    6    14
 +4.6      1  4    6    13
 +5.6      1  5    6    18
 +6.6      1  6    6    16
 +7.6      1  7    6    21
 +8.6      2  8    6    22
 +9.6      2  9    6    26
 +10.6     2 10    6    23
 +11.6     2 11    6    25
 +12.6     2 12    6    22
 +13.6     2 13    6    18
 +14.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>
Donate Powered by PHP Valid HTML5 Valid CSS Driven by DokuWiki