Функция pair () создает p-диаграмму рассеяния всех p x p пар переменных.
x<-rnorm(100,0,1)
y<-rnorm(100,0,1)
z<-rnorm(100,1,1)
t<-rnorm(100,2,10)
dd<-cbind(x,y,z,t)
pairs(dd)
Но я хотел бы иметь возможность выбирать для «линий» и столбцов сетки. Например, постройте эти пары (x, y) (x, z) (t, y) (t, z). Есть ли функция, которая принимает, например, формулу типа (x + t) ~ (z + y)?
Вы можете использовать конкретные аргументы pairs или создать настраиваемую функцию, которая принимает формулу в качестве входных данных:
pairs(dd, horInd=c(1, 4), verInd=c(2,3))
Пользовательская функция:
my_pairs <- function(df, formula) {
form <- deparse(formula)
s <- strsplit(form, "~")
lhs <- trimws(unlist(strsplit(s[[1]][1], "\\+")))
rhs <- trimws(unlist(strsplit(s[[1]][2], "\\+")))
lhs.ind <- match(lhs, colnames(df))
rhs.ind <- match(rhs, colnames(df))
all_cmbs <- expand.grid(lhs.ind, rhs.ind)
rows <- all_cmbs[,1]
cols <- all_cmbs[,2]
par(mfrow=c(2,2))
for(i in 1:4) {
eval(substitute(plot(df[,rows[i]], df[,cols[i]],
xlab=colnames(df)[rows[i]],
ylab=colnames(df)[cols[i]])))
}
}
my_pairs(dd, x + t ~ y + z)
