我想计算威布尔随机变量X大于1000的概率,其中X ~威布尔(拉姆达,g).
使用CDO,我将其计算为:
g = log(log(.7)/log(.5))/log(40/81)
lambda = -log(.7)/400^g
> exp(-lambda*1000^g)
[1] 0.4294352
但使用如下R包,结果不同
> pweibull(1000, shape = g, scale= lambda, lower.tail = FALSE)
[1] 0
为什么结果不同?我使用R包是否错误?