gpt4 book ai didi

r - 如何在 R 中读取相关矩阵并形成散点图矩阵

转载 作者:行者123 更新时间:2023-12-02 04:10:39 24 4
gpt4 key购买 nike

我在Excel中有一个相关矩阵如下:

dfA <- read.table(text=
"beta1 beta2 beta3 beta4 beta5 beta6 X X2 X3
beta1 1.0000 -0.2515 -0.2157 0.7209 -0.7205 0.4679 0.1025 -0.3606 -0.0356
beta2 -0.2515 1.0000 0.9831 0.1629 -0.1654 -0.5595 -0.0316 0.0946 0.0829
beta3 -0.2157 0.9831 1.0000 0.1529 -0.1559 -0.4976 -0.0266 0.0383 0.0738
beta4 0.7209 0.1629 0.1529 1.0000 -1.0000 -0.2753 0.0837 -0.1445 0.0080
beta5 0.4679 -0.5595 -0.4976 -0.2753 1.0000 0.2757 0.0354 -0.3149 -0.0596
beta6 -0.7205 -0.1654 -0.1559 -1.0000 0.2757 1.0000 -0.0837 0.1451 -0.0081
X 0.1025 -0.0316 -0.0266 0.0837 -0.0837 0.0354 1.0000 0.0278 -0.0875
X2 -0.3606 0.0946 0.0383 -0.1445 0.1451 -0.3149 0.0278 1.0000 0.2047
X3 -0.0356 0.0829 0.0738 0.0080 -0.0081 -0.0596 -0.0875 0.2047 1.0000",
header=TRUE)

我只有相关矩阵,而不是形成矩阵的原始数据,因此,我尝试使用以下代码将此矩阵读入 R 中的矩阵:

 B <- as.matrix(dfA)

但是当我尝试使用以下代码形成散点图矩阵时:

library(corrplot)
corrplot(B, method="circle")

我收到错误

Error in corrplot(B, method = "circle") : The matrix is not in [-1, 1]!

请帮我解决这个问题。

最佳答案

corrplot() 解决方案

根据上面 user20650 的评论使用 ggplot 更新我的第一篇文章。 user20650 显示错误的可能来源是舍入错误,导致某些数字超出允许的 [-1,1] 范围,舍入解决了此问题。我也能够使用 corrplot() 生成绘图。

此时,运行 corrplot() 会产生以下绘图:

corMat<-as.matrix(dfA)

library('corrplot')
corrplot(corMat, method='circle')

enter image description here

ggplot() 解决方案

您还可以在 ggplot2 中执行此操作,只需执行一些额外的步骤。我个人认为看起来好多了。

1)我去掉了矩阵下三角中的冗余信息。

corMat[lower.tri(corMat)]<-NA

> print(corMat)
beta1 beta2 beta3 beta4 beta5 beta6 X X2 X3
beta1 1 -0.2515 -0.2157 0.7209 0.4679 -0.7205 0.1025 -0.3606 -0.0356
beta2 NA 1.0000 0.9831 0.1629 -0.5595 -0.1654 -0.0316 0.0946 0.0829
beta3 NA NA 1.0000 0.1529 -0.4976 -0.1559 -0.0266 0.0383 0.0738
beta4 NA NA NA 1.0000 -0.2753 -1.0000 0.0837 -0.1445 0.0080
beta5 NA NA NA NA 1.0000 0.2757 -0.0837 0.1451 -0.0081
beta6 NA NA NA NA NA 1.0000 0.0354 -0.3149 -0.0596
X NA NA NA NA NA NA 1.0000 0.0278 -0.0875
X2 NA NA NA NA NA NA NA 1.0000 0.2047
X3 NA NA NA NA NA NA NA NA 1.0000

2) 然后,我使用 reshape2::melt() 将矩阵转换为长形式,并创建仅显示小数点后两位的值的格式化版本。这对情节很有用。

library(reshape2)
m<-melt(corMat)
m<-data.frame(m[!is.na(m[,3]),]) # get rid of the NA matrix entries
m$value_lab<-sprintf('%.2f',m$value)

数据如下:

> head(m)
Var1 Var2 value value_lab
1 beta1 beta1 1.0000 1.00
10 beta1 beta2 -0.2515 -0.25
11 beta2 beta2 1.0000 1.00
19 beta1 beta3 -0.2157 -0.22
20 beta2 beta3 0.9831 0.98
21 beta3 beta3 1.0000 1.00

3) 最后,我将这些数据输入 ggplot2 - 主要依靠 geom_tile() 来打印矩阵,并依靠 geom_text() 来打印每个图 block 上的标签。如果你愿意的话,你可以把它打扮得更漂亮。

library(ggplot2)
ggplot(m, aes(Var2, Var1, fill = value, label=value_lab),color='blue') +
geom_tile() +
geom_text() +
xlab('')+
ylab('')+
theme_minimal()

enter image description here

关于r - 如何在 R 中读取相关矩阵并形成散点图矩阵,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/36953464/

24 4 0
Copyright 2021 - 2024 cfsdn All Rights Reserved 蜀ICP备2022000587号
广告合作:1813099741@qq.com 6ren.com