R QQ-график не показывает

Question

делать онлайн курс и я получил следующую задачу:

load("skew.RData") 

Использование QQ-графики, сравнить распределение каждого столбца матрицы к нормальному. То есть, используйте qqnorm() для каждого столбца. Чтобы выполнить это быстро, вы можете использовать следующую строку кода для настройки сетки для 3x3 = 9 графиков. («Mfrow» означает, что мы хотим многофигурную сетку, заполненную в строке за строкой. Другой выбор mfcol.)

par(mfrow = c(3,3)) 

После этого вы можете использовать цикл, цикл по столбцам, а также отображать один qqnorm () за один раз. Вы должны заменить текст между ** своим собственным кодом.

for (i in 1:9) { 
    **put your qqnorm call here** 
} 

У меня есть файл calld dat, что-то с 9 столбцами и около 1000 строк.

Может ли кто-нибудь дать мне подсказку, как может выглядеть команда QQ-plot?

Answer 1

сравнить распределение каждого столбца матрицы к нормальному

Предположим, что ваша матрица dat после прочтения в файле данных, вы хотите

load("skew.RData") ## read in matrix `dat` 
dat <- scale(dat) ## standardization 

par(mfrow = c(3, 3)) 
for (i in 1:9) { 
    qqnorm(dat[, i], main = paste0(i, "-th column")) 
    qqline(dat[, i]) 
    }