獲得資料幀中一個變數與所有其他變數之間的引導相關區間-有解無憂

我想獲得一個變數（v）和所有其他變數（x，y，z）之間的引導相關引數（間隔和 rho）。如果我只對相關性本身感興趣（而不是 boostrap），我會使用以下公式。

df[,-1] %>% map(.,~cor.test(df$v, .))

我們如何在 bootstrap 中使用類似的想法？

v<-rnorm(20,50,2)
x<-rnorm(20,2)
y<-rnorm(20,10,2)
z<-rnorm(20,10,2)

df<-data.frame(v=v,x=x,y=y,z=z)

uj5u.com熱心網友回復：

撰寫一個函式來運行測驗，就像在一個運行案例中一樣。然后提取 rho 的估計值和置信區間。該函式回傳這些值。

使用 packageboot運行以上R次數，在下面的代碼中運行 5 次。

set.seed(2022)
v<-rnorm(20,50,2)
x<-rnorm(20,2)
y<-rnorm(20,10,2)
z<-rnorm(20,10,2)

df<-data.frame(v=v,x=x,y=y,z=z)

library(magrittr)
library(boot)

boot_cor <- function(data, i) {
  d <- data[i, ]
  cor_list <- d[-1] %>% purrr::map(~cor.test(v, .))
  pval <- sapply(cor_list, `[[`, 'estimate')
  conf.int <- sapply(cor_list, `[[`, 'conf.int')
  out <- rbind(pval, conf.int)
  out
}

b <- boot(df, boot_cor, R = 5)
b$t
#>            [,1]        [,2]      [,3]        [,4]        [,5]       [,6]        [,7]       [,8]      [,9]
#> [1,] 0.10901729 -0.35040802 0.5261551  0.47729667  0.04408784 0.75941797 -0.28443711 -0.6456858 0.1808467
#> [2,] 0.17520115 -0.28978685 0.5732758 -0.01138502 -0.45163044 0.43331887  0.20238775 -0.2637551 0.5918977
#> [3,] 0.45951068  0.02132649 0.7496046 -0.42373326 -0.72947047 0.02312342 -0.09795514 -0.5180212 0.3601784
#> [4,] 0.08274794 -0.37344763 0.5067140 -0.26406702 -0.63265841 0.20206624 -0.32949464 -0.6737737 0.1323194
#> [5,] 0.18779791 -0.27781014 0.5819556  0.11815477 -0.34226148 0.53281673  0.16289648 -0.3013469 0.5647101

# maybe more readable
array(b$t, dim = c(5, 3, 3),
      dimnames = list(NULL, c("rho", "lower", "upper"), names(df[-1])))
#> , , x
#> 
#>             rho       lower     upper
#> [1,] 0.10901729 -0.35040802 0.5261551
#> [2,] 0.17520115 -0.28978685 0.5732758
#> [3,] 0.45951068  0.02132649 0.7496046
#> [4,] 0.08274794 -0.37344763 0.5067140
#> [5,] 0.18779791 -0.27781014 0.5819556
#> 
#> , , y
#> 
#>              rho       lower      upper
#> [1,]  0.47729667  0.04408784 0.75941797
#> [2,] -0.01138502 -0.45163044 0.43331887
#> [3,] -0.42373326 -0.72947047 0.02312342
#> [4,] -0.26406702 -0.63265841 0.20206624
#> [5,]  0.11815477 -0.34226148 0.53281673
#> 
#> , , z
#> 
#>              rho      lower     upper
#> [1,] -0.28443711 -0.6456858 0.1808467
#> [2,]  0.20238775 -0.2637551 0.5918977
#> [3,] -0.09795514 -0.5180212 0.3601784
#> [4,] -0.32949464 -0.6737737 0.1323194
#> [5,]  0.16289648 -0.3013469 0.5647101

^{由reprex 包（v2.0.1）于 2022-06-04 創建}

轉載請註明出處，本文鏈接：https://www.uj5u.com/qiye/488135.html

標籤：r 字典应用开机

上一篇：這個“jQuery110...”得到回應是什么，我該如何處理它？

下一篇：只保留嵌套字典的第一個元素