2017-11-14 39 views
0

これは対数正規分布に従う一連のデータを持っています。私はmuとsigmaの値を見つけましたが、ログ正規のPDFとともにデータのヒストグラムをプロットしようとしています。しかし、線をプロットすると0以下になるようですが、意味が分からないようですが、ここで間違っているのは、線が5 yの値より高くなるようです。R - ggplot2ヒストグラムの対数正規曲線をプロットする

おかげ M

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 

enter image description here

+0

本当に 'log = TRUE'ですか?これは、<1が負になるので、 'log(確率)'を取っています。また、 'dlnorm'が'密度 'を与えるので、これは正しくありません。現在のプロットにカウントが表示されます。 –

答えて

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..を使用周波数の代わりにカウント。

+0

ありがとう、私が非常にばかげて、y軸が過去1に近づいていたことに気づかない –

関連する問題