Матричная диаграмма - протокол разведочного анализа
pairs(x, ...)
Параметры:
cex.labels = NULL - размер подписей меток на диагонали, NULL - размер подбирается программой;
labels = c("Ущерб","Расстояни","Следование","Локализация") - текст меток;
font.labels = 3 - шрифт подписей 3 - Arial-курсив
label.pos = 0.4 - смещение подписей по вертикали: 0 - низ, 0.5 - середина, 1 - верх
panel = panel.smooth - сглаживающая кривая;
horInd = 1:5 - отображаемые строки на диаграмме
verInd = 2:5 - отображаемые столбцы на диаграмме
row1attop = F, - изменение расположения диагонали подписей
gap = 0.3 - Размер полей между плашками, по умолчанию 1
lower.panel=NULL - убирает нижнюю половину графиков
Пример:
# Функция для расчета коэффициентов корреляции
panel.cor <- function(x, y, digits = 2, prefix = "", cex.cor, ...)
{
usr <- par("usr"); on.exit(par(usr))
par(usr = c(0, 1, 0, 1))
r <- abs(cor(x, y, method = "pearson"))
txt <- format(c(r, 0.123456789), digits=digits)[1]
txt <- paste(prefix, txt, sep = "")
# эта команда позволяет изменять размер шрифта в соответствии
# со значением коэффициента корреляции:
if(missing(cex.cor)) cex.cor <- 0.8/strwidth(txt)
text(0.5, 0.5, txt, cex = cex.cor * r)
}
pairs(data.frame(n$F36, n$F26, n$Sled, n$Lok, n$F161),
labels = c("Ущерб","Расстояни","Следование","Локализация","Штаб"), # Подписи
cex.labels = NULL, # Размер текста подписей
panel = panel.smooth, # Тренд
lower.panel=panel.cor) # на нижней части выводится значение корреляции Пирсона
Конец примера
PerformanceAnalytics
Очень информативная корреляционная матрица
library(PerformanceAnalytics)
chart.Correlation(mtcars[2:5], histogram=T)
