Рассчитать результаты для всех переменных из двух фреймов данных

Question

У меня есть этот код (отсюда: https://stats.idre.ucla.edu/r/faq/how-can-i-generate-bootstrap-statistics-in-r/):

library(boot)
library(tidyverse)


hsb2 <- read.table("https://stats.idre.ucla.edu/stat/data/hsb2.csv", sep=",", header=T)


fc <- function(d, i){
  d2 <- d[i,]
  return(cor(d2$write, d2$math))
}


#turn off set.seed() if you want the results to vary
set.seed(626)
fc
bootcorr <- boot(hsb2, fc, R=500)
temp1<-boot.ci(boot.out = bootcorr, type = "perc")
result<-temp1$percent[4:5]
result

Этот код рассчитывает только для первого аргумента cor() function, который равен d2$write, и второго аргумента, который d2$math переменных. Теперь давайте добавим I df5 и df6 фреймы данных:

df5<-select(hsb2, id, female, race, ses, schtyp, prog, math, science, socst)
df6<-select(hsb2, write, read)

df5 dataframe содержит только переменные, которые должны быть первым аргументом для функции cor(). df6 содержит только переменные, которые должны быть вторым аргументом для функции cor(). Как получить результаты для всех этих переменных из фреймов данных df5 и df6 в качестве первого и второго аргументов для соответствующей функции cor()?