r - Weibull 函数作为直线(R 中的 ggplot2 的双对数)

Question

希望有R方面比较懂的人帮忙。

问题:我尝试用 ggplott2 绘制 Weibull 概率图已经有一段时间了。 2-2 因子 Weibull 应在双对数图中显示为一条直线。但它不直或线不见了......一个例子是这样的https://www.mathworks.com/help/stats/wblplot.html

背景:我正在使用 flexsurv 和生存包进行一些实时计算。但我想使用 ggplot2 环境来显示它更多的花里胡哨。最后，我想将故障添加到先前计算中拟合的直线上。

欢迎所有建议。 :)

这里是我的示例代码:

library(ggplot2)

ggplot(data.frame(x=c(0,30)), aes(x=x)) +
  stat_function(fun = dweibull, args = list(shape = 2, scale = 15)) 
  #coord_trans(x = "log10", y = "log10") # without the option the function appears, but not as straight line 

Answer 1

QQ图

您正在寻找的(来自那个 Matlab 示例)是一个 Q-Q 图，ggplot 可以通过 geom_qq 或 stat_qq (此处的文档:https://ggplot2.tidyverse.org/reference/geom_qq.html):

ggplot(data.frame(x = rweibull(n = 100, shape = 1.2, scale = 1.5)), aes(sample = x)) + 
  stat_qq(distribution = stats::qweibull, dparams = list(shape = 1.2, scale = 1.5))

祝你好运!