Рисование qqline на графике QQ хи-квадрат в R - PullRequest
Новая
       14

Рисование qqline на графике QQ хи-квадрат в R

0 голосов
/ 30 марта 2019

Я копирую следующий код r help. Однако функция qqline не работает. В чем причина? 

library(mgcViz)
y <- rchisq(500, df = 3)
## Q-Q plot for Chi^2 data against true theoretical distribution:
qqplot(qchisq(ppoints(500), df = 3), y, main = expression("Q-Q plot for" ~~ {chi^2}[nu == 3]))
qqline(y, distribution = function(p) qchisq(p, df = 3), prob = c(0.1, 0.6), col = 2)

Спасибо!

1 Ответ

1 голос
/ 30 марта 2019

Чтобы добавить qqline с mgcViz, вам нужно использовать ggplot: 

library(mgcViz)
y <- rchisq(500, df = 3)
## Q-Q plot for Chi^2 data against true theoretical distribution:
qqplot(qchisq(ppoints(500), df = 3), y, main = expression("Q-Q plot for" ~~ {chi^2}[nu == 3]))

# Add qq line
library(ggplot2)
ggplot2::last_plot() + qqline(y, distribution = function(p) qchisq(p, df = 3), prob = c(0.1, 0.6), col = 2)

qq plot

Чтобы просмотреть перехват и наклон вашей qqline, вы можете сделать следующее: 

my.qqline = qqline(y, distribution = function(p) qchisq(p, df = 3), prob = c(0.1, 0.6), col = 2)
> my.qqline$data
   intercept     slope
1 0.08606285 0.9290453
...