Warning: file_get_contents(/data/phpspider/zhask/data//catemap/4/r/77.json): failed to open stream: No such file or directory in /data/phpspider/zhask/libs/function.php on line 167

Warning: Invalid argument supplied for foreach() in /data/phpspider/zhask/libs/tag.function.php on line 1116

Notice: Undefined index: in /data/phpspider/zhask/libs/function.php on line 180

Warning: array_chunk() expects parameter 1 to be array, null given in /data/phpspider/zhask/libs/function.php on line 181

Warning: file_get_contents(/data/phpspider/zhask/data//catemap/8/perl/9.json): failed to open stream: No such file or directory in /data/phpspider/zhask/libs/function.php on line 167

Warning: Invalid argument supplied for foreach() in /data/phpspider/zhask/libs/tag.function.php on line 1116

Notice: Undefined index: in /data/phpspider/zhask/libs/function.php on line 180

Warning: array_chunk() expects parameter 1 to be array, null given in /data/phpspider/zhask/libs/function.php on line 181
两组数据的组合corrplot?_R_R Corrplot - Fatal编程技术网

两组数据的组合corrplot?

两组数据的组合corrplot?,r,r-corrplot,R,R Corrplot,我想知道我是否可以为两组数据生成一个组合corrplot?也就是说,一组的上三角,另一组的下三角。我举一个例子如下: library(corrplot) mydata <- replicate(5, rnorm(20)) colnames(mydata) <- c('x1','x2','x3','x4','x5') mydata <- transform(mydata, group = c(rep('A',10),rep('B',10))) corrplot(cor(m

我想知道我是否可以为两组数据生成一个组合corrplot?也就是说,一组的上三角,另一组的下三角。我举一个例子如下:

library(corrplot)

mydata <- replicate(5, rnorm(20))

colnames(mydata) <- c('x1','x2','x3','x4','x5')

mydata <- transform(mydata, group = c(rep('A',10),rep('B',10)))

corrplot(cor(mydata[which(mydata$group=='A'),c(1:5)]), method="number", type="upper", title="Group A", mar=c(1,0,1,0))
corrplot(cor(mydata[which(mydata$group=='B'),c(1:5)]), method="number", type="lower", title="Group B", mar=c(1,0,1,0))
库(corrplot)

mydata您可以将一个相关矩阵的值替换为另一个相关矩阵的值

cor_A <- cor(mydata[which(mydata$group=='A'),c(1:5)])
cor_B <- cor(mydata[which(mydata$group=='B'),c(1:5)])

# replace the values
cor_A[lower.tri(cor_A)] <- cor_B[lower.tri(cor_B)]

corrplot(cor_A, method="number", title="Group A & B", mar=c(1,0,1,0))

coru A这里有一种可能性:

library(corrplot)    
set.seed(1)
mydata <- replicate(5, rnorm(20))   
colnames(mydata) <- c('x1','x2','x3','x4','x5')    
mydata <- transform(mydata, group = c(rep('A',10),rep('B',10)))

cor_A = cor(mydata[which(mydata$group=='A'),c(1:5)])
cor_B = cor(mydata[which(mydata$group=='B'),c(1:5)])
my_cor = upper.tri(cor_A)*cor_A + lower.tri(cor_B) * cor_B

corrplot(cor(mydata[which(mydata$group=='A'),c(1:5)]), method="number", type="upper", title="Group A", mar=c(1,0,1,0))
corrplot(cor(mydata[which(mydata$group=='B'),c(1:5)]), method="number", type="lower", title="Group B", mar=c(1,0,1,0))
corrplot(my_cor, method="number", title="Group A (upper) and B (lower)", mar=c(1,0,1,0))
库(corrplot)
种子(1)
我的数据
library(corrplot)    
set.seed(1)
mydata <- replicate(5, rnorm(20))   
colnames(mydata) <- c('x1','x2','x3','x4','x5')    
mydata <- transform(mydata, group = c(rep('A',10),rep('B',10)))

cor_A = cor(mydata[which(mydata$group=='A'),c(1:5)])
cor_B = cor(mydata[which(mydata$group=='B'),c(1:5)])
my_cor = upper.tri(cor_A)*cor_A + lower.tri(cor_B) * cor_B

corrplot(cor(mydata[which(mydata$group=='A'),c(1:5)]), method="number", type="upper", title="Group A", mar=c(1,0,1,0))
corrplot(cor(mydata[which(mydata$group=='B'),c(1:5)]), method="number", type="lower", title="Group B", mar=c(1,0,1,0))
corrplot(my_cor, method="number", title="Group A (upper) and B (lower)", mar=c(1,0,1,0))