построить график хи-квадрат с 1 степенью свободы с помощью ggplot

Question

Я провожу проверку гипотез, и мне нужно построить график следующего хи-квадрата $ \ chi_ {1,0.95} ^ 2 $. Я пытался сделать это с помощью ggplot, и это выглядит не очень хорошо. Надеюсь, ты сможешь помочь.

library(ggplot2)
df<-1
p<-0.05
gg<-data.frame(x=seq(0,15,0.1))
gg$y<-dchisq(gg$x,df)
p<-ggplot(gg)+geom_path(aes(x,y))
p

Answer 1

Возможно, вы хотите построить cdf распределения хи-квадрат, чтобы продемонстрировать, где он пересекает 0,95?

Вот график cdfs хи-квадрат с 1–5 степенями свободы, и точка где первая пересекает p = 0,95:

library(ggplot2)

gg <- data.frame(x = rep(seq(0,15,0.1), 5), df = rep(1:5, each = 151))
gg$y <-pchisq(gg$x, gg$df)
p <- ggplot(gg) +
      geom_path(aes(x, y, color = factor(df))) +
      geom_hline(aes(yintercept = 0.95), linetype = 2) +
      geom_vline(aes(xintercept = qchisq(0.95, 1)), linetype = 2) +
      geom_point(aes(x = qchisq(0.95, 1), y = 0.95))
p

введите описание изображения здесь