我有以下代码(来自此处: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
的变量进行计算。现在,我们添加我df5
和df6
数据帧:
df5<-select(hsb2, id, female, race, ses, schtyp, prog, math, science, socst)
df6<-select(hsb2, write, read)
df5
数据框仅包含应作为cor()
函数的第一个参数的变量。 df6
仅包含应作为cor()
函数的第二个参数的变量。如何从数据帧df5
和df6
中获得所有这些变量的结果,作为cor()
函数的第一和第二自变量?