给定不属于r的分布的分位函数

Question

我试图为给定的发行版生成分位数函数(F_q)。这是我的密码：

library(GoFKernel)

f_cdf <- function(x) {
  integrate(fstar, 1e-10, x)$value/I
}
f_q <- inverse(f_cdf, lower=1e-10, upper=5-1e-10)

ggplot(x_plot, aes(sample=x))+
  labs(title="Empirical against theoretical quantiles")+
   stat_qq(distribution=f_q) +
   stat_qq_line(distribution=f_q)

当我这样做时，我会得到以下错误：

stat_qq()中的计算失败:未使用的参数(p =分位数)。

stat_qq_line()中的计算失败:未使用的参数(p =分位数)。

当我的绘图代码和从f得到的样本时没有问题。原因是当我将'distribution=f_q‘改为'distribution=qnorm'：质量图给出了“错误”的分布。时，我可以成功地得到一个to图。

因此，我可以肯定，我的“f_q”，即分位数函数有一些缺陷。因此，如果有人知道stat_qq和stat_qq.line的“分发”参数的具体需求，那将是非常有帮助的。

fstar在0到5之间没有定义，范围在0到5之间。我设计了一个从这个分布来模拟的拒绝方案，我试图用分位数图来比较理论分位数和样本分位数。

Answer 1

分发函数的形式可能需要命名为p。尝试：

f_q_p <- function(p) { f_q(p) }

然后

ggplot(x_plot, aes(sample=x))+
  labs(title="Empirical against theoretical quantiles")+
   stat_qq(distribution=f_q_p) +
   stat_qq_line(distribution=f_q_p)