更新时间:2022-03-17 03:49:22
您是否正在寻找类似的东西?我无权访问您的数据集,所以我使用了虹膜数据集
Were you looking for something like this? I don't have access to your dataset so I used the iris dataset
library(dplyr); library(ggplot2)
meanSe <- iris %>%
filter(Species == "setosa") %>%
summarise(means = mean(Sepal.Length), sd=sd(Sepal.Length))
#
meanVe <- iris %>%
filter(Species == "versicolor") %>%
summarise(means = mean(Sepal.Length), sd=sd(Sepal.Length))
#
meanVi <- iris %>%
filter(Species == "virginica") %>%
summarise(means = mean(Sepal.Length), sd=sd(Sepal.Length))
#
ggplot(iris, aes(x=Sepal.Length, color=Species, fill=Species)) +
geom_histogram(aes(y=..density..), position="identity", binwidth=.5) +
stat_function(fun = dnorm, color="red", args=list(mean=meanSe$means, sd=meanSe$sd)) +
stat_function(fun = dnorm, color="green", args=list(mean=meanVe$means, sd=meanVe$sd)) +
stat_function(fun = dnorm, color="blue", args=list(mean=meanVi$means, sd=meanVi$sd)) +
theme_bw()
给这个