所以我有一组遵循对数正态分布的数据,我已经找到mu和sigma的值,现在我正在尝试绘制数据的直方图和对数正态的PDF。然而,绘制线似乎将其置于0以下似乎没有意义,我在这里做错了,看起来该线应该是5 y值更高。
由于 中号
h <- read.table("data.csv", header=TRUE, sep=",", row.names="id")
library("ggplot2")
d <- ggplot(data=h, aes(h$time)) +
geom_histogram(binwidth = 0.8) +
stat_function(fun = dlnorm, args = list(meanlog = 2.2, sdlog = 0.44, log = TRUE), colour = "red")
d
答案 0 :(得分:0)
能够进行比较:
h <- data.frame(time = rlnorm(100, 2.2, .44))
library("ggplot2")
d <- ggplot(data=h) +
geom_histogram(binwidth = 0.8, aes(x = time, ..density..)) +
stat_function(fun = dlnorm, args = list(meanlog = 2.2, sdlog = 0.44),
colour = "red")
d
根据建议,请关闭log = TRUE
,然后使用..density..
获取频率而非计数。