Question

我有92套同类数据。

我想为这两种组合提供一个相关矩阵。

i.e.,我希望有一个92x92的矩阵。

此种内容(ci,cj)应当与ci和cj相关。

我如何这样做?

Answer 1

例如:

 d &lt- data.frame(x1=rnorm(10),
                 x2=rnorm(10),
                 x3=rnorm(10))
cor(d) # get correlations (returns matrix)

Answer 2

您可使用校准包。

d <- data.frame(x1=rnorm(10),
                 x2=rnorm(10),
                 x3=rnorm(10))
M <- cor(d) # get correlations

library( corrplot ) #package corrplot
corrplot(M, method = "circle") #plot matrix

“entergraph

详情见:http://cran.r-project.org/web/parages/corrplot/vignettes/corrplot-intro.html>。

Answer 3

The cor function will use the columns of the matrix in the calculation of correlation. So, the number of rows must be the same between your matrix x and matrix y. Ex.:

set.seed(1)
x <- matrix(rnorm(20), nrow=5, ncol=4)
y <- matrix(rnorm(15), nrow=5, ncol=3)
COR <- cor(x,y)
COR
image(x=seq(dim(x)[2]), y=seq(dim(y)[2]), z=COR, xlab="x column", ylab="y column")
text(expand.grid(x=seq(dim(x)[2]), y=seq(dim(y)[2])), labels=round(c(COR),2))

enter image description here

Edit:

此处是按单一矩阵计算的相关矩阵的海关行号及栏目标签的例子:

png("corplot.png", width=5, height=5, units="in", res=200)
op <- par(mar=c(6,6,1,1), ps=10)
COR <- cor(iris[,1:4])
image(x=seq(nrow(COR)), y=seq(ncol(COR)), z=cor(iris[,1:4]), axes=F, xlab="", ylab="")
text(expand.grid(x=seq(dim(COR)[1]), y=seq(dim(COR)[2])), labels=round(c(COR),2))
box()
axis(1, at=seq(nrow(COR)), labels = rownames(COR), las=2)
axis(2, at=seq(ncol(COR)), labels = colnames(COR), las=1)
par(op)
dev.off()

Answer 4

参看。它使你能够创建<>互动相关性矩阵:

library(qtlcharts)
data(iris)
iris$Species <- NULL
iplotCorr(iris, reorder=TRUE)

<

It s more impressive when you correlate more variables, like in the package s vignette:

Answer 5

There are other ways to achieve this here: (Plot correlation matrix into a graph), but I like your version with the correlations in the boxes. Is there a way to add the variable names to the x and y column instead of just those index numbers? For me, that would make this a perfect solution. Thanks!

编辑:我试图通过[在盒子里的马克]就这个职位发表评论,但我显然不知道我做什么。然而,我确实设法为自己回答这一问题。

if d is the matrix (or the original data frame) and the column names are what you want, then the following works:

axis(1, 1:dim(d)[2], colnames(d), las=2)
axis(2, 1:dim(d)[2], colnames(d), las=2)

激光=0将姓名贴回正常位置,地雷很长,因此我用激光=2将其排在轴上。

edit2: to suppress the image() function printing numbers on the grid (otherwise they overlap your variable labels), add xaxt= n , e.g.:

image(x=seq(dim(x)[2]), y=seq(dim(y)[2]), z=COR, col=rev(heat.colors(20)), xlab="x column", ylab="y column", xaxt= n )

Edit:

友情链接