Добавление сводной статистики (или даже точки необработанных данных) к коробчатым диаграммам положения, которых избегают,

Хэдли, несомненно, поправит меня, если я ошибаюсь ...

Вот естественный синтаксис:

bplot + geom_point(aes(colour=trt), position=position_dodge(width=.5))

(position = "dodge" будет делать то же самое, без параметра .)

Когда я рисую это, я получаю что-то похожее на position_jitter (), что, по-видимому, тоже есть у вас.

Любопытно, я пошел посмотреть в исходный код, где нашел функцию pos_dodge () . (Введите pos_dodge в приглашении R, чтобы увидеть его ...) Вот его конец:

within(df, {
  xmin <- xmin + width / n * (seq_len(n) - 1) - diff * (n - 1) / (2 * n)
  xmax <- xmin + d_width / n
  x <- (xmin + xmax) / 2
})

n - количество строк в кадре данных. Итак, похоже, что он уклоняется от отдельных точек на дробь, индексированную по строке! Таким образом, первая точка уклоняется шириной / n, вторая уклоняется 2 * ширина / n, а последняя уклоняется n * width / n.

Очевидно, это не то, что вы имели в виду , хотя это то, что вы сказали . Возможно, вы застряли, воссоздавая уклоненную коробчатую диаграмму вручную или используя другую визуализацию, например, фасетирование?

ggplot(to.analyze,aes(inj.site,relief)) + geom_boxplot() + facet_wrap(~ trt)