Используя данные внизу этого поста, после запуска ановы и теста Тьюки:
library(multcomp)
library(multcompView)
mod1 <- aov(lm(Var1 ~ Location, data=dat))
tky <- TukeyHSD(mod1)
Я пытаюсь создать такой сюжет (взято из этого поста https://stats.stackexchange.com/questions/31547/how-to-obtain-the-results-of-a-tukey-hsd-post-hoc-test-in-a-table-showing-groupe ): 
Я пробовал несколько примеров, в том числе некоторые с использованием пакетов foreign и agricolae, но я не могу заставить что-либо работать. Может ли кто-то продемонстрировать, как сделать этот сюжет, с данными ниже:
> dput(dat)
structure(list(Location = structure(c(1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 5L, 5L, 5L, 5L, 5L,
5L, 5L, 5L, 5L, 5L, 5L, 5L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L,
4L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 7L, 7L, 7L, 7L, 7L, 7L,
7L, 7L, 7L, 7L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 9L,
9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 14L,
14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 15L, 15L, 15L, 15L, 15L,
15L, 15L, 15L, 15L, 15L, 15L, 15L, 15L, 16L, 16L, 16L, 16L, 16L,
16L, 16L, 16L, 16L, 17L, 17L, 17L, 17L, 17L, 17L, 17L, 17L, 17L,
17L, 17L, 17L, 18L, 18L, 18L, 18L, 18L, 18L, 18L, 18L, 18L, 18L,
18L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 10L, 10L,
10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L,
10L, 10L, 10L, 10L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L,
12L, 12L, 20L, 20L, 20L, 20L, 20L, 20L, 20L, 20L, 20L, 20L, 20L,
20L, 19L, 19L, 19L, 19L, 19L, 19L, 19L, 19L, 19L, 19L, 19L, 19L,
19L, 19L, 25L, 25L, 25L, 25L, 25L, 21L, 21L, 21L, 21L, 21L, 21L,
21L, 21L, 21L, 21L, 22L, 22L, 22L, 22L, 22L, 22L, 22L, 22L, 22L,
13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 23L,
23L, 23L, 23L, 23L, 23L, 23L, 23L, 23L, 23L, 23L, 23L, 23L, 23L,
24L, 24L, 24L, 24L, 24L, 24L, 24L, 24L, 24L, 26L, 26L, 26L, 26L,
26L), .Label = c("Al", "Apg", "Blp", "Can", "Cas", "Cf", "Chag",
"Conp", "Jam", "Jor", "Mic", "Mir", "MLF", "Nat", "Nop", "Noy",
"Opp", "Plr", "Rag", "Rog", "SK", "SQ", "Tay", "Tur", "UPP",
"WGG"), class = "factor"), Var1 = c(1.13, 1.15, 1.59, 1.42, 1.09,
1.11, 1.64, 0.92, 1.53, 1.18, 1.91, 1.35, 1.12, 1.02, 1.48, 1.44,
2.08, 1.31, 1.1, 1.38, 1.97, 1.52, 1.83, 1.08, 1.34, 1.39, 1.6,
1.37, 1.3, 1.24, 1.6, 1.43, 1.43, 1.18, 1.38, 1.46, 1.65, 1.39,
1.72, 1.42, 1.43, 1.2, 1.74, 0.62, 1.32, 1.79, 1.19, 1.09, 1.32,
1.35, 1.3, 1.86, 0.98, 1.41, 1.33, 1.51, 1.54, 1.2, 1.59, 0.88,
1.39, 1.33, 1.48, 1.6, 1.37, 1.42, 1.26, 1.44, 1.04, 0.92, 2.1,
1.35, 1.59, 1.25, 1.56, 1.33, 1.19, 0.86, 1.19, 1.97, 1.41, 1.51,
1.16, 1.09, 1.31, 1.15, 1.11, 1.22, 1.2, 1.25, 1.27, 1, 1.53,
1.6, 1.14, 0.84, 0.97, 1.11, 1.16, 1.16, 1.4, 1.27, 1.06, 1.26,
1.94, 1.24, 1.37, 1.43, 1.5, 1.43, 1.59, 1.05, 1.27, 1.36, 1.43,
1.59, 1.53, 1.56, 1.76, 1.47, 1.14, 1.11, 1.58, 1.07, 1.62, 1.48,
1.16, 0.75, 1.2, 1.21, 0.95, 1.54, 1.39, 0.98, 1.5, 1.36, 1.19,
1.24, 1.38, 1.35, 1.27, 1.19, 1.27, 1.45, 2.12, 1.27, 1.34, 1.54,
1.33, 1.38, 1.15, 1.7, 1.35, 1.28, 1.19, 1.48, 0.93, 1.56, 1.55,
1.57, 1.28, 1.31, 1.6, 1.14, 1.24, 1.06, 1.53, 1.35, 1.84, 1.22,
1.33, 1.42, 1.46, 1.29, 1.35, 1.41, 1.33, 1.43, 1.57, 1.9, 1.05,
1.23, 0.96, 1.21, 0.88, 1.11, 1.2, 1.24, 1.15, 1.82, 1.58, 1.37,
1.38, 1.34, 1.29, 1.31, 1.73, 1.56, 1.86, 1.18, 1.58, 1.53, 1.23,
1.35, 1.28, 1.11, 1.4, 1.5, 1.05, 1, 1.55, 1.05, 0.78, 1.16,
0.92, 0.85, 1.27, 1.19, 1.24, 1.1, 0.89, 1.19, 1.43, 1.41, 1.31,
1.36, 1.11, 1.75, 1.54, 1.34, 1.42, 1.32, 1.29, 2, 1.28, 1.46,
1.05, 1.64, 1.27, 1.26, 1.3, 1.35, 1.45, 1.13, 1.34, 1.49, 1.67,
2.07, 1.62, 1.3, 1.72, 1.17, 1.74, 1.52, 1.66, 1.31, 1.35, 1.94,
1.76, 2.16, 1.51, 2.13, 1.89, 2.15, 1.53, 1.89, 1.1, 1.62, 1.9,
1.16, 1.68, 1.34, 1.56, 1.32, 1.96, 0.97, 1.16, 1.62, 1.5, 1.14,
1.16, 1.05, 1.21, 0.72, 1.1)), class = "data.frame", row.names = c(NA,
-285L))