Я совершенно отчаялся, попробовал glht() и другие вещи, чтобы получить ярлыки от TukeyHSD на моем графике, но это просто не сработало.
Моя модель зависит от нескольких факторов .. ..:
model3 <- lm(barley.yield.t.ha ~ site*SHADE*fertilization*irrigation*samppos, data = newplot)
Мне нужен график ggplot2, показывающий barley.yield.t.ha для различных обработок оттенков (второй фактор в модели, "SHADE"). Существует три варианта обработки (0, 40% и 90%) тени, и я хочу, чтобы они были показаны на оси абсцисс (для всех сайтов и на втором графике для каждого сайта) в виде прямоугольников с значительными буквами над ними, аналогично этому графику здесь : https://www.r-graph-gallery.com/84-tukey-test_files/figure-html/unnamed-chunk-3-1.png
Есть предложения? Я получил всевозможные проблемы и сообщения об ошибках, например:
Error in modelparm.default(model, ...) :
dimensions of coefficients and covariance matrix don't match
In addition: Warning message:
In mcp2matrix(model, linfct = linfct) :
covariate interactions found -- default contrast might be inappropriate
или
LABELS <- generate_label_df(TUKEY , "newplot$site")
Error in variable[, 4] : incorrect number of dimensions
Любой, кто сможет мне помочь или дать хороший совет о том, каким образом go / пытаться? Я очень благодарен ...