На самом деле вырвал мои волосы, потому что я искал, как кодировать это весь день!
Есть ли простой способ использовать вывод из Tukey HSD для создания компактного отображения букв на ggplot? ?
Вот мой код Tukey HSD и вывод для Tukey HSD
model2a<-aov(ChlorophyllSq~Treatment*Depth, data= Chloro_day1)
summary(model2a)
TukeyHSDG3 <- TukeyHSD(model2a)
Вывод для этого:
$`Treatment:Depth`
diff lwr upr p adj
1:0-4:0 -0.22553886 -0.453050360 0.001972637 0.0539570
2:0-4:0 0.08538540 -0.142126102 0.312896895 0.9479148
3:0-4:0 0.13881631 -0.088695190 0.366327807 0.5837634
4:1-4:0 -0.11860442 -0.346115918 0.108907079 0.7605038
1:1-4:0 -0.32883710 -0.556348601 -0.101325604 0.0003326
2:1-4:0 -0.10558438 -0.333095877 0.121927120 0.8530927
3:1-4:0 -0.06751487 -0.295026371 0.159996627 0.9859703
2:0-1:0 0.31092426 0.083412759 0.538435757 0.0009340
3:0-1:0 0.36435517 0.136843672 0.591866669 0.0000363
4:1-1:0 0.10693444 -0.120577057 0.334445940 0.8445680
1:1-1:0 -0.10329824 -0.330809739 0.124213258 0.8669147
2:1-1:0 0.11995448 -0.107557015 0.347465982 0.7497017
3:1-1:0 0.15802399 -0.069487509 0.385535488 0.4094089
3:0-2:0 0.05343091 -0.174080586 0.280942411 0.9966252
4:1-2:0 -0.20398982 -0.431501315 0.023521682 0.1167876
1:1-2:0 -0.41422250 -0.641733997 -0.186711000 0.0000011
2:1-2:0 -0.19096977 -0.418481273 0.036541724 0.1762182
3:1-2:0 -0.15290027 -0.380411767 0.074611230 0.4545645
4:1-3:0 -0.25742073 -0.484932227 -0.029909230 0.0141612
1:1-3:0 -0.46765341 -0.695164909 -0.240141912 0.0000000
2:1-3:0 -0.24440069 -0.471912186 -0.016889189 0.0251239
3:1-3:0 -0.20633118 -0.433842679 0.021180318 0.1079735
1:1-4:1 -0.21023268 -0.437744181 0.017278816 0.0944567
2:1-4:1 0.01302004 -0.214491457 0.240531540 0.9999998
3:1-4:1 0.05108955 -0.176421951 0.278601046 0.9974563
2:1-1:1 0.22325272 -0.004258775 0.450764222 0.0588699
3:1-1:1 0.26132223 0.033810732 0.488833729 0.0118425
3:1-2:1 0.03806951 -0.189441992 0.265581005 0.9996237
Я пытался применить этот пример к мои данные, но я получаю сообщение об ошибке из-за взаимодействия Как создать компактный буквенный дисплей для парного TukeyHSD
glht(model2a, infct=mcp(Treatment*Depth="Tukey")
Любые рекомендации или помощь будет высоко ценится?
Спасибо!