Пары матрицы с qq-графиками

Question

Полезно для просмотра ваших данных. Когда у вас есть несколько переменных, вы можете сформировать матрицу рассеянного экрана, например, pairs(). Матрица рассеяния дает вам набор двумерных маргинальных прогнозов ваших данных.

set.seed(8092) 
X <- matrix(rnorm(80), ncol=4) 
pairs(X) 

Вы также можете иметь данные для различных групп и хотим сравнить их распределения. Два таких распределения можно сравнить с qq-графиком.

set.seed(4415) 
group1 <- rnorm(20) 
group2 <- rnorm(20) 
qqplot(group1, group2) 
abline(c(0,1)) 

Если у вас есть несколько групп, было бы удобно, если бы там были пары типа участка, который будет отображать матрицу q¯q-участков.

colnames(X) <- c("group1", "group2", "group3", "group4") 
qq.pairs(X) 

Есть ли такая функция? Есть ли простой способ закодировать его с нуля?

Answer 1

Вы можете добавить любую функцию, которая вам нравится, через аргументы *panel. Например:

set.seed(8092) 
X <- matrix(rnorm(80), ncol=4) 

panel.qq <- function(x, y, ...) { 
    usr <- par("usr"); on.exit(par(usr)) 
    par(usr = c(0, 1, 0, 1), new = TRUE) 
    qqplot(x, y, xlab = deparse(substitute(x)), ylab = deparse(substitute(y))) 
    abline(c(0,1), ...) 
} 

pairs(X, lower.panel = panel.qq) 

Конечно, вы можете использовать panel=panel.qq сделать все квадраты ¯qq участки

Answer 2

Я получил сообщение об ошибке с просто положить qqplot в функции панели которой она появляется @rawr решена используя «new = TRUE». Мой подход был вытаскивать х и у значения и использовать points:

panel.qq <- function(x,y, ...) 
{ usr <- par("usr"); on.exit(par(usr)) 
    par(usr = c(usr[1:2], -2,2)) 
    QQ <- qqplot(x,y, plot.it = FALSE) 
    x <- QQ$x; 
    y <- QQ$y; 
    points(x, y, col = "cyan", ...) 
    abline(0,1) 
} 

set.seed(8092) 
X <- matrix(rnorm(80), ncol=4) 
pairs(X, upper.panel=panel.qq)