ggplot2のメモ（その3）

正規分布

library(ggplot2)
# 正規分布
# seq(a,b,n):下限a、上限b、分割数nの等差数列を作る。
x <-  seq(-4, 4, length=50)
# mean:平均値、sd：標準偏差(sd = sqrt(sigma^2))
fx <- dnorm(x=x,mean=0,sd = 1)
data <-data.frame(x=x,y=fx)
ggplot(data, aes(x=x,y=fx)) + geom_line()

対数正規分布

library(ggplot2)
# 対数正規分布
# seq(a,b,n):下限a、上限b、分割数nの等差数列を作る。
x <-  seq(0.01, 10, length=1000)
# f(x) = 1/(sqrt(2*pi)sigma x) exp {-(logx-mu)^2 / (2 sigma^2)}
# meanlog:mu、sdlog：sigma
fx <- dlnorm(x=x, meanlog = 0, sdlog=1)
data <-data.frame(x=x,y=fx)
ggplot(data, aes(x=x,y=fx)) + geom_line()
# xを対数にしたもの。上記の正規分布と同じになる。
ggplot(data, aes(x=log(x),y=fx)) + geom_line()

次のグラフは横軸をにして表示したものである。