Question

如果x是n*m矩阵，则当我使用cor(x)时，每对列之间都有一个m*m相关矩阵。

我如何在cor.test矩阵上使用n*m函数来拥有m*m p值矩阵？

Answer 1

可能已有功能，但这是我的版本。 p_cor_mat在矩阵cor.test的每一对列上运行x并记录p值。然后将它们放到一个正方形矩阵中并返回。

# Set seed
set.seed(42)

# Matrix of data
x <- matrix(runif(120), ncol = 4)

# Function for creating p value matrix
p_cor_mat <- function(x){
  # All combinations of columns
  colcom <- t(combn(1:ncol(x), 2))

  # Calculate p values
  p_vals <- apply(colcom, MAR = 1, function(i)cor.test(x[,i[1]], x[,i[2]])$p.value)

  # Create matrix for result
  p_mat <- diag(ncol(x))

  # Fill upper & lower triangles
  p_mat[colcom] <- p_mat[colcom[,2:1]] <- p_vals

  # Return result
  p_mat
}

# Test function
p_cor_mat(x)
#>           [,1]      [,2]      [,3]      [,4]
#> [1,] 1.0000000 0.4495713 0.9071164 0.8462530
#> [2,] 0.4495713 1.0000000 0.5960786 0.7093539
#> [3,] 0.9071164 0.5960786 1.0000000 0.7466226
#> [4,] 0.8462530 0.7093539 0.7466226 1.0000000

^{由reprex package（v0.2.1）于2019-03-06创建}

在R中使用cor.test函数

1 个答案: