正規分布、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))