Как сравнить два набора данных с графиком QQ, используя ggplot2?

11

Как новичок в области статистики и R, мне было очень трудно пытаться сгенерировать qqplots с соотношением сторон 1: 1. ggplot2, кажется, предлагает гораздо больший контроль над построением графиков, чем стандартные пакеты графиков R, но я не вижу, как выполнить qqplot в ggplot2 для сравнения двух наборов данных.

Итак, мой вопрос, что такое эквивалент ggplot2 что-то вроде:

qqplot(datset1,dataset2)

r distributions ggplot2 qq-plot

— psychemedia
источник

Документы ggplot2 могут быть полезны: docs.ggplot2.org/current/stat_qq.html

— Чарли,

12

Самое простое, что нужно сделать, это просто посмотреть, как это qqplotработает. Итак, в R тип:

R> qqplot
function (x, y, plot.it = TRUE, xlab = deparse(substitute(x)), 
    ylab = deparse(substitute(y)), ...) 
{
    sx <- sort(x)
    sy <- sort(y)
    lenx <- length(sx)
    leny <- length(sy)
    if (leny < lenx) 
        sx <- approx(1L:lenx, sx, n = leny)$y
    if (leny > lenx) 
        sy <- approx(1L:leny, sy, n = lenx)$y
    if (plot.it) 
        plot(sx, sy, xlab = xlab, ylab = ylab, ...)
    invisible(list(x = sx, y = sy))
}
<environment: namespace:stats>

Так что для генерации сюжета нам просто нужно получить sxи sy, т.е.

x <- rnorm(10);y <- rnorm(20)

sx <- sort(x); sy <- sort(y)
lenx <- length(sx)
leny <- length(sy)
if (leny < lenx)sx <- approx(1L:lenx, sx, n = leny)$y
if (leny > lenx)sy <- approx(1L:leny, sy, n = lenx)$y

require(ggplot2)
g = ggplot() + geom_point(aes(x=sx, y=sy))
g

qqplot используя ggplot2

— csgillespie
источник

2

ggplot2есть stat_qq(), есть ли способ использовать это? Кажется, он предназначен для сравнения одного вектора с теоретическим распределением, я не мог понять, как использовать его для сравнения двух разных векторов.

— Кен Уильямс

7

На самом деле вы можете qqplot()сделать все sort/ length/ approxвычисления для вас: d <- as.data.frame(qqplot(x, y, plot.it=FALSE)); ggplot(d) + geom_point(aes(x=x, y=y))

— Кен Уильямс

9

Я использую это, когда я тоже хочу нормальную линию.

ggplot(data, aes(sample = data$column1)) + stat_qq(color="firebrick2", alpha=1) + geom_abline(intercept = mean(data$column1), slope = sd(data$column1))

— KLDavenport
источник

0

Если ваша первоначальная потребность - просто контролировать соотношение сторон, вот один из способов сделать это:

x <- rnorm(1000)
y <- rnorm(1500, 2)

myqq <- function(x, y, ...) {
  rg <- range(x, y, na.rm=T)
  qqplot(x, y, xlim=rg, ylim=rg, ...)
}

myqq(x, y)

сюжет Myqq

— Кен Уильямс
источник