> dput(can1_scores)
structure(list(Tempo = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 4L, 4L, 4L,
4L, 4L, 5L, 5L, 5L, 5L, 5L, 6L, 6L, 6L, 6L, 6L, 7L, 7L, 7L, 7L,
7L), .Label = c("0", "2", "4", "6", "8", "10", "12"), class = "factor"),
Can1 = c(-3.97903120666332, -2.08094673287043, -3.14865243445412,
-2.86189109535914, -2.43198427698181, -2.6178610707061, -2.43352997837457,
-2.40806466202522, -2.46111031174995, -3.71884011193399,
-2.05034263703473, -2.0052215528943, -3.0429268229614, -1.07030351926172,
-1.59189581930415, -0.488694426925713, -0.382347271019961,
-0.28869621735626, -1.63394987567049, -0.388184676668998,
1.0491443457085, 1.35412326627344, 0.06221094779304, 1.72914975430105,
1.49862908063138, 1.88308178933803, 3.76305713116495, 2.6139946529437,
3.24530431502092, 3.28546336805719, 2.19512933319006, 2.71337547039165,
2.4748149110824, 1.49630700969963, 1.79375734523888, 2.24198697232919,
2.26729170076295, 1.57615670948704, 2.00861522820103, 1.83288136860133
)), .Names = c("Tempo", "Can1"), row.names = c(NA, 40L), class = "data.frame")
str(can1_scores)
'data.frame': 40 obs. of 2 variables:
$ Tempo: Factor w/ 7 levels "0","2","4","6",..: 1 1 1 1 1 1 1 1 1 1 ...
$ Can1 : num -3.98 -2.08 -3.15 -2.86 -2.43 ...
> av1 <- aov(Can1~Tempo, data=can1_scores)
> summary(av1)
Df Sum Sq Mean Sq F value Pr(>F)
Tempo 6 193.01 32.17 88.92 <2e-16 ***
Residuals 33 11.94 0.36
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
> sk1 <- summary(SK(av1, which ='Tempo'))
Levels Means SK(5%)
8 2.9581803 a
10 2.1346768 b
12 1.9853864 b
6 1.1386515 c
4 -0.6363745 d
2 -1.9521381 e
0 -2.8141912 f
There were 14 warnings (use warnings() to see them)