Deepseek也沒有解決,還是好好學R代碼。
1、讀入數據data
2、計算相關系數
cor_result<-corr.test(data,
? ? ? ? ? ? ? ? method = "pearson")?
# 提取需要的相關系數和p值矩陣
cor_matrix <- cor_result$r [-12:-21,12:21] # 相關系數矩陣
p_value_matrix <- cor_result$p [-12:-21,12:21] # p值矩陣
mm<-range(cor_matrix)? #系數范圍
3、繪圖,有意義的顯示,無意義的相關系數不顯示
corrplot::corrplot(cor_matrix,
? ? ? ? is.corr = F,? #非對稱矩陣
? ? ? ? method = "color",? ? ? # 使用顏色顯示相關系數
? ? ? ? type = "upper",? ? ? ? # 只顯示上三角矩陣
? ? ? ? tl.col = "black",? ? ? # 設置文本顏色
? ? ? ? tl.srt = 45,? ? ? ? ? ? # 設置文本旋轉角度
? ? ? ? addCoef.col = "black",? # 添加相關系數的數值
? ? ? ? p.mat = p_value_matrix, # 添加p值矩陣
? ? ? ? sig.level = 0.05,? ? ? # 顯著性水平
? ? ? ? insig = "blank",? ? ? ? # 顯示不顯著的區域為空白
? ? ? ? diag = FALSE,? ? ? ? ? # 去掉對角線
? ? ? ? number.cex = 1.2,? ? ? ? # 設置相關系數值的字體大小
? ? ? ? number.font = 2,
????? col = COL2('PRGn'))?