可视化概率分布曲线中的拒绝区域

Question

我希望能够0.05在以下卡方密度曲线中可视化 R 中单尾和双尾 p 值的拒绝区域：

curve(dchisq(x,24), xlim=c(5,50), ylim=c(0,.06))

如果有人也有其他测试统计数据的通用方法，那也会很棒。

提前致谢

托尼奥

Answer 1

使用 ggplot 的方法：

df   <- 24
p    <- 0.05
gg   <- data.frame(x=seq(5,50,0.1))
gg$y <- dchisq(gg$x,df)

library(ggplot2)
ggplot(gg) + 
  geom_path(aes(x,y)) +
  geom_linerange(data=gg[gg$x>qchisq(p,df,lower.tail=F),],
                 aes(x, ymin=0, ymax=y),
                 colour="red")