在R中多个变量的相关系数矩阵如何来检验其显著性呢？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

变量相关系数的检验可以用cor.test()，那么多个变量的相关系数矩阵如何来检验其显著性呢？谢谢！
多个变量的两两相关性检验，加星号表示显著性的那种

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

关键词：相关系数矩阵 相关系数 多个变量 相关性检验 test

去github上搜代码啊
  require(psych)这个包要安装好

1. corstar <- function(x, y = NULL, use = "pairwise", method = "pearson", round = 3, row.labels, col.labels, ...) {
   
2.   
   
3.   require(psych)
   
4. 
   
5.   ct <- corr.test(x, y, use, method)    # calculate correlation
   
6.   r <- ct$r                             # get correlation coefs
   
7.   p <- ct$p                             # get p-values
   
8. 
   
9.   stars <- ifelse(p < .001, "***", ifelse(p < .01, "** ", ifelse(p < .05, "*  ", "   "))) # generate significance stars
   
10.   
    
11.   m <- matrix(NA, nrow = nrow(r) * 2, ncol = ncol(r) + 1) # create empty matrix
    
12. 
    
13.   rlab <- if(missing(row.labels)) rownames(r) else row.labels # add row labels
    
14.   clab <- if(missing(col.labels)) {
    
15.     if(is.null(colnames(r)))
    
16.       deparse(substitute(y))
    
17.     else
    
18.       colnames(r)
    
19.   } else {
    
20.     col.labels # add column labels
    
21.   }
    
22.   
    
23.   rows <- 1:nrow(m)                     # row indices
    
24.   cols <- 2:ncol(m)                     # column indices
    
25.   
    
26.   odd <- rows %% 2 == 1                 # odd rows
    
27.   even <- rows %% 2 == 0                # even rows
    
28.   m[odd, 1] <- rlab                     # add variable names
    
29.   m[even, 1] <- rep("", sum(even))      # add blank
    
30.   
    
31.   m[odd, cols] <- paste(format(round(r, round), nsmall = round, ...), stars, sep = "")     # add r coefs
    
32.   m[even, cols] <- paste("(", format(round(p, round), nsmall = round, ...), ")", sep = "") # add p-values
    
33.   
    
34.   colnames(m) <- c(" ", clab)           # add colnames
    
35.   m                                     # return matrix
    
36. }
    
37. 
    
38. corstar(iris[,-5])
    
39.                   Sepal.Length Sepal.Width Petal.Length Petal.Width
    
40. [1,] "Sepal.Length" " 1.000***"  "-0.118   " " 0.872***"  " 0.818***"
    
41. [2,] ""             "(0.000)"    "(0.152)"   "(0.000)"    "(0.000)"  
    
42. [3,] "Sepal.Width"  "-0.118   "  " 1.000***" "-0.428***"  "-0.366***"
    
43. [4,] ""             "(0.152)"    "(0.000)"   "(0.000)"    "(0.000)"  
    
44. [5,] "Petal.Length" " 0.872***"  "-0.428***" " 1.000***"  " 0.963***"
    
45. [6,] ""             "(0.000)"    "(0.000)"   "(0.000)"    "(0.000)"  
    
46. [7,] "Petal.Width"  " 0.818***"  "-0.366***" " 0.963***"  " 1.000***"
    
47. [8,] ""             "(0.000)"    "(0.000)"   "(0.000)"    "(0.000)"  

复制代码

jiangbeilu 发表于 2017-11-20 14:22
去github上搜代码啊
  require(psych)这个包要安装好

谢谢，我是新手，github是啥啊？

5752 发表于 2017-11-20 14:41
谢谢，我是新手，github是啥啊？

就是好多好多程序员和媛集中的地方

jiangbeilu 发表于 2017-11-20 14:22
去github上搜代码啊
  require(psych)这个包要安装好

您好，那做带显著性星号相关矩阵可视图代码是什么您知道吗？