我试图在R中绘制一个直方图,并用不同分布的密度曲线叠加在上面。对于常规的直方图,它可以正常工作,但是我无法在ggplot2包中实现这个效果。
a <- dataset$age
现在是我的常规直方图代码:
Histogram_for_age <- hist(a, prob=T, xlim=c(0,80), ylim=c(0,0.055), main="Histogram for age with density lines", xlab="age")
mean <- mean(a)
sd <- sd(a)
现在需要绘制密度的线条/曲线:
lines(density(dataset$age), col="blue", lwd=2, lty=1)
curve(dnorm(x, mean = mean, sd = sd), add = T, col="red", lwd=2, lty=2)
curve(dgamma(x, shape =mean^2/sd^2, scale = sd^2/mean), add = T, col="goldenrod", lwd=2, lty=3)
还有一个传说:
legend("topright",
c("actual distribution of age","gaussian distribution", "gamma distribution"),
lty=c(1,2,3),
lwd=c(2,2,2),col=c("blue","red","goldenrod"), cex=0.65)
这是我迄今为止使用ggplot2尝试的内容:
ggplot(dataset, aes(x=age)) +
geom_histogram(aes(y=..density..),
colour="black", fill="white") +
geom_density(alpha=.2, fill="lightblue") + stat_function(fun = dgamma, shape=shape)
ggplot2中哪个参数等同于我的lines()和curve()参数?
geom="point"
然后定义形状。你能提供一些示例数据吗?这样我就可以更容易地给出适当的答案了。 - Jaapdgamma
函数中的shape
参数。您应该在stat_function()
调用中添加arg=list(shape=...)
- ilirstat_bin: binwidth默认为range/30。使用'binwidth = x'来调整它。
- Cajira