Как построить среднее и стандартную ошибку в Boxplot в R

У меня есть два категориальных фактора («Место обитания» и «Местность») и одна непрерывная переменная (T). «Хабитат» имеет два уровня, а «Местность» - восемь уровней. Я хочу изменить усы по умолчанию, чтобы они представляли SE, а медиану - среднюю для каждого коробочного графика. Есть ли способ сделать это и принять во внимание оба категориальных фактора при построении? Спасибо заранее.

Это то, что я сделал с настройкой по умолчанию boxplot ggplot, показывая первый и третий квартили с медианными интервалами.

ggplot(data,aes(x=Locality,y=T)) + 
  geom_boxplot(aes(fill=interaction(Habitat,Locality), 
                   group=interaction(factor(Habitat),Locality)),
               outlier.shape=1,outlier.size=3) + 
  theme_bw() + 
  theme(
    panel.grid.major=element_blank(),
    panel.grid.minor=element_blank(),
    axis.line=element_line(colour='black'),
    legend.position='none',
    axis.text.x=element_text(angle=90,hjust=1,size=12)) + 
  scale_y_continuous('T') + 
  xlab('Locality')

Ответы на вопрос(2)

Ваш ответ на вопрос