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()

f:id:nishiru3:20180630123737p:plain

対数正規分布

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()

f:id:nishiru3:20180630124633p:plain

次のグラフは横軸xlog(x)にして表示したものである。
f:id:nishiru3:20180630124648p:plain