R による連続確率分布と累積分布の描画

正規分布、t 分布、カイ2乗分布、F 分布、ベータ分布、ワイブル分布、一様分布、ガンマ分布、コーシー分布、指数分布、対数正規分布、べき分布(パレート分布)を描画する(Distributions)。

正規分布

正規分布と標準正規分布を描く。

#正規分布。平均50、標準偏差5。
curve(dnorm(x,50,5),30,70) #確率分布
curve(qnorm(x,50,5)) #累積分布
#標準正規分布
curve(dnorm, -5,5) #確率分布
curve(qnorm) #累積分布

t 分布

ここでは自由度を1としている。

curve(dt(x,1),-1,1) #確率分布
curve(qt(x,1)) #累積分布

カイ2乗分布

ここでは自由度を3としている。

curve(dchisq(x, 3)) #確率分布
curve(qchisq(x, 3)) #累積分布

F 分布

ここでは自由度を5と10としている。

curve(df(x,5,10)) #確率分布
curve(qf(x,5,10)) #累積分布

ベータ分布

ここではパラメータを0.3と0.5としている。

curve(dbeta(x,0.3,0.5)) #確率分布
curve(qbeta(x,0.3,0.5)) #累積分布

ワイブル分布

ここでは shape = 0.8, scale = 2 としている。

curve(dweibull(x,0.8,2)) #確率分布
curve(qweibull(x,0.8,2)) #累積分布

一様分布

curve(dunif(x,-1,1),-1.1,1.1) #確率分布
curve(qunif(x,-1,1)) #累積分布

ガンマ分布

ここでは shape = 2, scale = 1.5 としている。

curve(dgamma(x,2,1.5),0,5) #確率分布
curve(qgamma(x,2,1.5)) #累積分布

コーシー分布

ここでは location = 0, scale = 1 の標準コーシー分布としている。

curve(dcauchy(x,0,1),-1,1) #確率分布
curve(qcauchy(x,0,1)) #累積分布

指数分布

ここでは rate = 2 としている。

curve(dexp(x,2)) #確率分布
curve(qexp(x,2)) #累積分布

対数正規分布

ここではパラメーターを 1, 2 としている。

curve(dlnorm (x,1,2),1,3, log = "xy") #確率分布
curve(qlnorm (x,1,2)) #累積分布

べき分布(パレート分布)

ここでは a を shape、b を scale として関数を定義し、a = 2, b = 1 としている。

#確率密度関数
Fd = function(x, a, b) a/b * (b/x)^(a+1) 
curve(Fd(x,2,1),1,10, log = "xy")
#累積分布関数
Fq = function(x, a, b) 1-(b/x)^a
curve(Fq(x,2,1))