【发布时间】:2019-02-11 07:49:49
【问题描述】:
我试图结合直方图和伽马分布计算的最大似然估计,而直方图表示 x 的观察值,一条线表示伽马分布计算的最大似然估计。示例图如下所示:
但是,直方图在 x 轴上的长度为 13,而 gamma 分布在 x 轴上的长度为 512,它一直抱怨`
错误:美学长度必须为 1 或与数据相同 (13): x, y
den <- density(x) # x's length is 13, however den's length is 512
dat <- data.frame(x = den$x, y = den$y)
dataframe = data.frame(days = days, x = x)
ggplot(data = dat) +
geom_line(aes(x=dat$x, y=dgamma(dat$x,shape, rate)), color="red", size = 1) +
theme_classic()
bplot = ggplot(data = dataframe,aes(x = dataframe$days, y = (dataframe$x)/sum(dataframe$x)));
bplot +
geom_bar(stat="identity", width=0.5)+
geom_line(aes(x=dat$x, y=dgamma(dat$x,shape, rate)), color="blue", size = 1)
【问题讨论】: