R - как ограничить вывод для Hmisc Rcorr?

Question

У меня есть два фрейма данных motivation_on с 60 наблюдениями и motivation_off с 146 наблюдениями, каждый из которых содержит 21 переменную и 1 столбец ID, который находится в первом столбце.

Теперь я хочу знать, как переменные соотносятся друг с другом, поэтому я использую:

rcorr(as.matrix(motivation_on[2:ncol(motivation_on)]), type = "spearman")

а также

rcorr(as.matrix(motivation_off[2:ncol(motivation_off)]), type = "spearman")

(Подмножество сделано, чтобы избавиться от столбца ID)

Теперь я хочу вычислить соотношение между онлайн и оффлайн переменными, поэтому я попытался:

rcorr(as.matrix(motivation_off[2:ncol(motivation_off)]), as.matrix(motivation_on[2:ncol(motivation_on)]) , type = "spearman")

Теперь я получаю то, что хочу, но кроме того, он также показывает все корреляции для переменных внутри motivation_on и внутри motivation_off, которые я вычислял ранее. Это делает вывод чрезвычайно длинным. Как получить вывод rcorr исключительно для корреляций on_off?

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

x <- as.matrix(mtcars[1:3])
y <- as.matrix(mtcars[4:6])
rcorr(x,y)

То, что я хотел, это корреляционная таблица с: mpg, cyl, disp в виде строк и hp, drat, wt в виде столбцов, а не с полным выводом. моя текущая работа вокруг:

z <- rcorr(x,y)
q <- as.data.frame(z[1])
q[1:3,4:6]

Answer 1

Кажется, что Hmisc's rcorr не предоставляет возможность извлечь то, что я хочу. Можно извлечь результаты только вручную, как я продемонстрировал в моем примере с mtcars. Тем не менее psych::corr.test(x, y) обеспечивает именно то, что я хотел - спасибо @user20650 за то, что указал на это!