Есть ли в R способ зациклить функции multcompview (например, CLD) через списки ANOVA и Tukeys?

Question

В конце концов я хочу построить каждую переменную с помощью CLD. Но в то же время у меня проблемы с получением CLD. У меня есть данные, которые выглядят так;

      df <- data.frame(Treat = rep(LETTERS[1:4], 100, replace = TRUE),
                 A = rnorm(400),
                 B = rnorm(400),
                 C = rnorm(400),
                 D = rnorm(400),
                 E = rnorm(400),
                 F = rnorm(400),
                 G = rnorm(400),
                 H = rnorm(400))

Благодаря очень полезному ответу на предыдущий вопрос здесь, я затем зациклил ANOVA и Tukey для каждой переменной (и создал кадры данных) следующим образом;

      ## perform anova on dependent variables
aov_res <- apply(df[,2:ncol(df)], 2, function(x) aov(x ~ Treat, data = df))

# Apply Tukey's HSD test to the results of each ANOVA test
tukey_res <- sapply(aov_res, function(x) TukeyHSD(x, "Treat", ordered = TRUE))

# Convert the results of each ANOVA test into a tidy data frame using the broom package
aov_res_df <- do.call(rbind, lapply(aov_res, broom::tidy))

# Combine the results of the Tukey HSD tests into a single data frame
tukey_res_df <- as.data.frame(do.call(rbind, Map(cbind, Name = names(tukey_res), tukey_res)))

Теперь мне просто нужно получить CLD для каждой переменной. я обычно используюmultcompView::multcompLetters4()но я знаю, что есть альтернативы. я тоже не могла понятьlapplyилиmapplyдля работы с двумя отдельными списками. В конце концов, я хотел быgeom_boxplot()результаты, и я буду полностью потерян, но маленькими шагами. Любая помощь приветствуется!