Saya baru mengenal R, dan ingin menambahkan batang kesalahan standar ke grafik batang saya. Saya memiliki kode berikut:

p1 <- ggplot(plot1, aes(x=factor(drink_type), y=value, fill = manipulation)) + 
  stat_summary(fun.y="mean", geom="bar", position="dodge") +
  theme_classic() +
 labs(x = 'Drink type', y = 'Evidence accumulation')

p1 + stat_summary(fun.data = mean_cl_normal, geom = "errorbar", width = .08, position = position_dodge(0.9)) 

Namun, saya tidak yakin apakah 'geom = "errorbar" ini menghasilkan bilah kesalahan standar (yang saya cari). Ini adalah data yang saya kerjakan. Saya ingin tahu apakah saya perlu membuat kolom baru yang merupakan SE dari kolom 'nilai' untuk menampilkan bilah kesalahan standar?

ppt | manipulation |   drink_type           |  value
1   |       n      |   avg_alcohol_drift    |  1.8854094
1   |       p      |   avg_alcohol_drift    |  1.6257274
1   |       n      |   avg_softdrink_drift  |  1.8519074
1   |       p      |   avg_softdrink_drift  |  1.8477995

Bantuan apa pun akan sangat dihargai!

r
2
abc_95 15 Maret 2020, 20:56

1 menjawab

Jawaban Terbaik

Ya seharusnya, Anda dapat mencoba sesuatu seperti di bawah ini, di mana saya menggunakan contoh dataset:

ggplot(ToothGrowth,aes(x=factor(dose),y=len,fill=supp)) + 
stat_summary(fun.y=mean,geom="bar",position=position_dodge()) + 
stat_summary(fun.data=mean_se,geom="errorbar",
position=position_dodge(0.9),width=0.2)

enter image description here

Kita bisa menghitungnya secara manual:

library(dplyr)
se_data = ToothGrowth %>% 
group_by(dose,supp) %>% 
summarise(mean=mean(len),se=sd(len)/sqrt(n()))

# A tibble: 6 x 4
# Groups:   dose [3]
   dose supp   mean    se
  <dbl> <fct> <dbl> <dbl>
1   0.5 OJ    13.2  1.41 
2   0.5 VC     7.98 0.869
3   1   OJ    22.7  1.24 
4   1   VC    16.8  0.795
5   2   OJ    26.1  0.840
6   2   VC    26.1  1.52 

ggplot(se_data,aes(x=factor(dose),y=mean,fill=supp)) + 
geom_bar(stat="identity",position="dodge") +
geom_errorbar(aes(ymin=mean-se,ymax=mean+se),
position=position_dodge(0.9),width=0.2)

enter image description here

2
StupidWolf 15 Maret 2020, 18:22