r-errorbar-3

とある棒グラフを作成。

errorbar-1

このグラフにエラーバーを付けたい場合、

geom_errorbar()
を追加する。

エラーバーをつける

標準偏差(standard deviation, SD)

SDは、

dplyr::summarise(新規作成列名 = sd(参照データの列))

で求める。

dplyr::summarise(mean = mean(Data), sd = sd(Data))

で、このSDの値をエラーバーにするには、

geom_errorbar(aes(ymax = 平均値 + エラーバー値(上), ymin = 平均値 - エラーバーの値(下)), position = 位置, width = 横の線の長さ, size = 先の太さ)

を追加する。

geom_errorbar(aes(ymax = mean + sd, ymin = mean - sd), position = position_dodge(0.9), width = 0.2, size = 0.75 )
errorbar-2

標準誤差(standard error of the mean, SEM)

SEMは、「標準偏差 ÷ ルート(サンプル数)」なので、

dplyr::mutate()

で列を追加する。

mutate(Data, sem = sd/sqrt(6))

で、新たに作ったSEMの列を、グラフのエラーバーの指示で指定する。

r-errorbar-3

 



 

全体のコードは下記↓

# bar graphs
Data$Group <- factor(Data$Group, levels=c("WTM", "PdM" , "WTF" , "PdF" )) Data$Region <-
factor(Data$Region, levels=c("Ipsi-Hipp", "Ipsi-Ctx" , "Cont-Hipp" , "Cont-Ctx" )) Data %>%
ggplot(aes(x = reorder(x = Sex, X = Group), y = mean, fill = Group)) +
geom_col(position = "dodge", colour = "black", size = 0.75) +
geom_errorbar(aes(ymax = mean + sem, ymin = mean - sem), position = position_dodge(0.9), width = 0.2,
size = 0.75 ) +
labs(title = "Datas", x = "Genotype", y = "% area of AT8+ Datas", fill="Genotype") +
scale_x_discrete(limit = c("M", "F")) +
scale_fill_discrete(limit = c("WTM", "PdM", "WTF", "PdM")) +
scale_fill_manual(values = c("#2980B9", "#154360", "#CB4335", "#641E16")) +
facet_wrap(~Region)
にほんブログ村 子育てブログ ワーキングマザー育児へ