R中的相关性

我试图生成与中相同的图形，但使用不同的数据。这是我的密码：

library(SciViews)

args <- commandArgs(TRUE)
pdfname <- args[1]
datafile <- args[2]

pdf(pdfname)
eqdata = read.csv(datafile , header = T,sep=",")
(longley.cor <- correlation(eqdata$feqs))
# Synthetic view of the correlation matrix
summary(longley.cor) 
p <- plot(longley.cor)
print(p)
dev.off()   

我使用这个命令

Rscript ./rscript/correlation.R "/home/co.pdf" "/home/data_correlation.csv"

代码输出

我想这样生成

---

您可以尝试

椭圆

软件包中的

plotcorr

功能。帮助页面提供了以下示例：

你在找什么

编辑： 之后可以添加文本，圆圈放置在1个变量网格上。例如：

data(mtcars)
Corrmat  <- cor(mtcars)
cols <- ifelse(Corrmat>0, rgb(0,0,abs(Corrmat)), rgb(abs(Corrmat),0,0))

library(ellipse)
plotcorr(Corrmat,col=cols)

n <- nrow(Corrmat)
for (i in 1:n)
{
    for (j in 1:n)
    {
        text(j,i,round(Corrmat[n-i+1,j],2),col="white",cex=0.6)     
    }
}

数据（mtcars）
Corrmat不幸的是，我不熟悉
SciViews
，但我仍然非常确信，仅为数据帧的一个变量计算
相关性并不是一个好主意。非常感谢您的建议。我想生成文本的相关性。。。所以现在只考虑每个单词的频率。对文本关联有什么想法吗？如果要关联的变量是
'abc'
、`def'等的频率，那么每个变量都需要一个以上的值。在您的示例中，每个变量只有一个值，您无法计算实际没有变化的事物的协方差。@Marius您是对的。我必须格式化像这样的数据>“abc”、“def”、“ghi”、“jkl”、“mno”、“pqr”>10,20,10,12,13,10>12,13,13,14,30,10>40,12,40,44,12,30，并且更接近我的要求。感谢sacha，我需要在每个椭圆上显示相关性。我可以生成上面的图形，但标签有问题。。。在标签而不是空间中，它被点（.）替换。你有什么想法，如何修复它？我仍然在使用SciViews软件包。
data(mtcars)
Corrmat  <- cor(mtcars)
cols <- ifelse(Corrmat>0, rgb(0,0,abs(Corrmat)), rgb(abs(Corrmat),0,0))

library(ellipse)
plotcorr(Corrmat,col=cols)

n <- nrow(Corrmat)
for (i in 1:n)
{
    for (j in 1:n)
    {
        text(j,i,round(Corrmat[n-i+1,j],2),col="white",cex=0.6)     
    }
}