library(readxl)
data_anova <- read_excel("test.xlsx",
sheet = "anova-twoway", range = "D5:F53")
data_anova <- as.data.frame(data_anova)
data_anova$diet <- as.factor(data_anova$diet)
data_anova$time <- as.factor(data_anova$time)
data_anova -> my_data
my_data
diet time stability
1 NT1 10-min 89.17712
2 NT1 10-min 89.76035
3 NT1 10-min 89.56072
4 NT1 5-min 86.81088
5 NT1 5-min 87.28091
6 NT1 5-min 84.50821
7 NT2 10-min 88.51183
8 NT2 10-min 89.83705
9 NT2 10-min 89.02299
10 NT2 5-min 85.94260
11 NT2 5-min 86.53363
12 NT2 5-min 84.58238
13 NT3 10-min 90.92160
14 NT3 10-min 90.31614
15 NT3 10-min 86.30768
16 NT3 5-min 88.22052
17 NT3 5-min 87.87284
18 NT3 5-min 86.48140
19 NT4 10-min 91.13180
20 NT4 10-min 91.06331
21 NT4 10-min 90.97900
22 NT4 5-min 87.41313
23 NT4 5-min 87.79204
24 NT4 5-min 87.83304
25 NT5 10-min 90.33636
26 NT5 10-min 90.45067
27 NT5 10-min 90.97430
28 NT5 5-min 86.72079
29 NT5 5-min 87.59405
30 NT5 5-min 87.51170
31 NT6 10-min 90.73797
32 NT6 10-min 90.21814
33 NT6 10-min 90.40041
34 NT6 5-min 86.89038
35 NT6 5-min 86.45951
36 NT6 5-min 87.25631
37 NT7 10-min 90.87250
38 NT7 10-min 90.58467
39 NT7 10-min 90.56533
40 NT7 5-min 85.88389
41 NT7 5-min 86.83931
42 NT7 5-min 87.13404
43 NT8 10-min 90.83993
44 NT8 10-min 91.04827
45 NT8 10-min 90.10980
46 NT8 5-min 86.33665
47 NT8 5-min 86.20297
48 NT8 5-min 84.82913