[英]Make rectangular matrix of correlation values in R, possibly using corrplot
I want to make a matrix of correlation values in R. However, I am not comparing the parameters to themselves, so it isn't your typical corrplot that is a triangle and mirrored over a diagonal axis.我想在 R 中创建一个相关值矩阵。但是,我没有将参数与其自身进行比较,因此不是典型的 corrplot 是三角形并镜像在对角轴上。 I actually already have the correlation values .我实际上已经有了相关值。
I was just hoping to recruit the corrplot
function (from corrplot package ), which produces images like the following:我只是希望招募corrplot
函数(来自corrplot 包),它生成如下图像:
My data, already containing the correlation values I want to plot, is:我的数据已经包含我想要绘制的相关值,是:
X animal.1 animal.2 animal.3 animal.4 animal.5
1 parameter 1 0.10258087 0.3338782 0.150246554 0.07295133 0.2484152
2 parameter 2 0.04205748 0.4062727 -0.002101464 0.12068818 0.2951127
3 parameter 3 0.11264488 0.4114954 0.067145776 0.13361071 0.3246052
4 parameter 4 -0.02261649 -0.2426341 0.108042167 -0.12820517 -0.2005686
5 parameter 5 -0.01576384 -0.2300852 0.112941655 -0.12391976 -0.1906473
6 parameter 6 -0.09749030 -0.3110920 0.021994297 -0.13570257 -0.2557532
I want the plot to look like a corrplot, but it will be a full rectangle.我希望绘图看起来像一个 corrplot,但它将是一个完整的矩形。 I want the circles of varying diameters, and the two different colors for positive or negative values, and the color gradient along the side of the plot.我想要不同直径的圆圈,正值或负值的两种不同颜色,以及沿着图边的颜色渐变。 In essence I want the corrplot theme of a rectangular table of values.本质上,我想要一个矩形值表的 corrplot 主题。
I have also asked to see if asterisks could be placed over the grid squares with significant (p<0.05) p-values.我还要求查看是否可以在具有显着 (p<0.05) p 值的网格方块上放置星号。 I am now attaching a dput file of the p-value matrix我现在附上 p 值矩阵的 dput 文件
structure(list(X = structure(1:6, .Label = c("parameter 1", "parameter 2",
"parameter 3", "parameter 4", "parameter 5", "parameter 6"), class = "factor"),
Animal.1 = c(0.2454906, 0.63471, 0.2019519, 0.7984066, 0.8587147,
0.2698293), Animal.2 = c(0.000103586, 1.62e-06, 1.15e-06,
0.005412082, 0.008451775, 0.000315107), Animal.3 = c(0.08796553,
0.9810688, 0.447827, 0.2211191, 0.2007585, 0.8038392), Animal.4 = c(0.4094627,
0.1713835, 0.1296492, 0.1460429, 0.1601174, 0.1236984), Animal.5 = c(0.004374306,
0.000653099, 0.00016464, 0.02213469, 0.0298038, 0.003315349
)), .Names = c("X", "animal.1", "animal.2", "animal.3", "animal.4",
"animal.5"), class = "data.frame", row.names = c(NA, -6L))
I think you just need to transform your df into a matrix and then run the corrplot
function.我认为您只需要将 df 转换为矩阵,然后运行corrplot
函数即可。
library(corrplot)
#Transform data to matrix
matrix_cor<-as.matrix(df[,-1])
#Set row names as df first column
row.names(matrix_cor)<-df[,1]
corrplot(matrix_cor,
method = "circle")
As you tag this question with ggplot2
, I'm providing you a way to do that using ggplot2
.当你用ggplot2
标记这个问题时,我为你提供了一种使用ggplot2
来做到这一点的方法。
First, you need to reshape both the dataframe containing correlation values (named df
here) and the one with pvalues (named df_pval
here).首先,您需要重塑包含相关值的数据df_pval
此处命名为df
)和具有 pvalues 的数据df_pval
此处命名为df_pval
)。 For doing that, you can for example use the pivot_longer
function from tidyr
package (also part of the tidyverse
package).对于这样做,比如,你可以使用pivot_longer
功能从tidyr
包(也部分tidyverse
包)。
On the dataframe containing pvalues, we will add a column Label
in order to mark with an asterisk all values below 0.05.在包含 pvalues 的数据框上,我们将添加一个列Label
以便用星号标记所有低于 0.05 的值。
library(tidyr)
library(dplyr)
library(ggplot2)
corr_values <- df %>% pivot_longer(-X, names_to = "Animal", values_to = "value")
# Here an extract of the dataframe with correlation values after reshaping it:
X Animal value
<fct> <chr> <dbl>
1 parameter 1 animal.1 0.103
2 parameter 1 animal.2 0.334
3 parameter 1 animal.3 0.150
4 parameter 1 animal.4 0.0730
5 parameter 1 animal.5 0.248
6 parameter 2 animal.1 0.0421
corr_pval <- df_pval %>%
pivot_longer(-X, names_to = "Animal", values_to = "p.value") %>%
mutate(Label = ifelse(p.value < 0.05, "*",NA))
# Here an extract of the dataframe obtained
X Animal p.value Label
<fct> <chr> <dbl> <chr>
1 parameter 1 animal.1 0.245 NA
2 parameter 1 animal.2 0.000104 *
3 parameter 1 animal.3 0.0880 NA
4 parameter 1 animal.4 0.409 NA
5 parameter 1 animal.5 0.00437 *
6 parameter 2 animal.1 0.635 NA
Then, you can use geom_point
and geom_tile
to get the heatmap.然后,您可以使用geom_point
和geom_tile
来获取热图。 You can also customize, color, size, etc.. of your plot by passing various functions such as scale_size_continuous
and scale_color_gradient
.您还可以通过传递各种函数(例如scale_size_continuous
和scale_color_gradient
定义绘图的颜色、大小等。
Finally, you can add asterisk for pvalues
by calling a new dataframe corr_pval
into geom_text
and specify the appropriate aesthetic
argument:最后,您可以通过将新数据帧corr_pval
调用到geom_text
并指定适当的aesthetic
参数来为pvalues
添加星号:
library(ggplot2)
ggplot(corr_values, aes(x = X, y = Animal))+
geom_tile(color = "black", fill = "white")+
geom_point(aes(color = value, size = abs(value)))+
scale_color_gradient2(low = "green", mid = "white", midpoint = 0, high = "red", name = "")+
scale_size_continuous(range = c(5,15), name = "")+
geom_text(data = corr_pval, aes(label = Label), size = 8, vjust = 0.7, hjust = 0.5)
So, you should get something like that:所以,你应该得到这样的东西:
声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.