R语言 如何绘制Weibull分布图
在这篇文章中,我们将看到如何在R编程语言中绘制一个Weibull分布。Weibull分布是一个连续的概率分布,用于分析寿命数据,建立故障时间模型,以及在旧时代没有现代机器的情况下获取产品的可靠性。
要在R语言中绘制Weibull分布,我们需要两个函数,即dweibull和curve()。
- dweibull(): 密度、分布函数、量化函数和随机生成的Weibull分布的参数shape和scale。
dweibull(x, shape, scale= 1) 。
- x – 量子化的向量
- shape – 形状参数。
- scale – 尺度参数。
-
curve(): 在区间[from, to]上绘制与函数相对应的曲线。曲线也可以在变量xname中绘制表达式,默认为x。
curve(expr, from = NULL, to = NULL)
- function – 一个函数的名称,或一个调用,或一个写成x的函数的表达式,它将评估一个与x相同长度的对象。
- from – 函数将被绘制的起始范围。
- to – 绘制函数的结束范围。
例1:
为了绘制Weibull分布的概率密度函数,我们必须在曲线函数中指定形状和比例值以及从和到的范围,如下所示。
curve(dweibull(x, shape=2, scale = 1), from=0, to=5)
输出
Weibull分布
例2:
我们还可以通过为每个分布指定不同的颜色来绘制具有不同形状、尺度和范围的多个Weibull分布,如下图所示。
curve(dweibull(x, shape=2, scale=1),
from=0, to=5, col='blue')
curve(dweibull(x, shape=3, scale=2),
from=0, to=7, col='red', add=TRUE)
curve(dweibull(x, shape=4, scale=3),
from=0, to=10, col='purple', add=TRUE)
输出
多重weibull分布
例3:
我们现在可以用不同的形状、比例和范围为Weibull分布添加一些美感,并包括绘图标题的轴名,如下图所示。
让我们添加一些参数,如main(主标题)、ylab-(y轴标签)、lwd-(线宽)和col(颜色)。此外,我们还可以将形状和比例值传递给图例函数,如下面的代码片断所示。
curve(dweibull(x, shape=3, scale = 1), from=0, to=10,
main = 'Weibull Distribution (shape = 2, scale = 1)', # add title
ylab = 'Density values', # change y-axis label
lwd = 2, # increase line width to 2
col = 'red') # change line color to steelblue
# add legend
legend(2, .7, legend=c("shape=2, scale=1"),
col=c("red"), lty=1, cex=1.5)
输出
带有x轴和y轴名称和图例的Weibull分布