R でワイブル分布をプロットする方法

R でワイブル分布の確率密度関数をプロットするには、次の関数を使用できます。

• dweibull(x,shape,scale = 1) を使用して確率密度関数を作成します。
• Curve(function, from = NULL, to = NULL)確率密度関数をプロットします。

確率密度関数をプロットするには、 dweibull関数で形状とスケールのパラメーターの値を指定するだけでなく、 curve()関数でfromとto の値を指定する必要があります。

たとえば、次のコードは、パラメーターshape = 2およびscale = 1でワイブル分布の確率密度関数をプロットする方法を示しています。プロットのx軸は0から4です。

 curve(dweibull(x, shape=2, scale = 1), from=0, to=4)

タイトルを追加したり、Y 軸ラベルを変更したり、行幅を広げたり、行の色を変更してプロットをより美しくすることもできます。

 curve(dweibull(x, shape=2, scale = 1), from=0, to=4, 
    main = 'Weibull Distribution (shape = 2, scale = 1)', #add title
    ylab = 'Density', #change y-axis label
    lwd = 2, #increase line width to 2
    col = 'steelblue') #change line color to steelblue 

グラフにいくつかの曲線を追加して、さまざまな形状およびスケール パラメーターを持つワイブル分布を比較することもできます。

 curve(dweibull(x, shape=2, scale = 1), from=0, to=4, col='red')
curve(dweibull(x, shape=1.5, scale = 1), from=0, to=4, col='blue', add=TRUE)

legend()関数を使用してプロットに凡例を追加できます。この関数は次の構文を使用します。

legend(x, y=NULL, 凡例, fill,col, bg, lty, cex)

金：

• x、y:凡例の配置に使用される x 座標と y 座標
• 凡例:凡例に入れるテキスト
• fill:凡例内の塗りつぶしの色
• Col:凡例内の線に使用する色のリスト
• bg:凡例の背景色
• lty:線のスタイル
• cex:凡例内のテキストのサイズ

この例では、次の構文を使用して凡例を作成します。

 #create density plots
curve(dweibull(x, shape=2, scale = 1), from=0, to=4, col='red')
curve(dweibull(x, shape=1.5, scale = 1), from=0, to=4, col='blue', add=TRUE)

#add legend
legend(2, .7, legend=c("shape=2, scale=1", "shape=1.5, scale=1"),
       col=c("red", "blue"), lty=1, cex=1.2)