斯皮尔曼的等级相关性

Question

我正在编写一个脚本，在两个向量中读取两个.txt文件。 之后我想制作一个斯皮尔曼的等级相关并绘制结果。 第一个向量值的长度为12-13个字符（例如7.3445555667或10.3445555667），第二个向量值的长度为一个字符（例如1或2）。

编码：

vector1 <- read.table ("D:...path.../mytext1.txt", header=FALSE)
vector2 <- read.table ("D:...path.../mytext2.txt", header=FALSE)
cor.coeff = cor(vector1 , vector2 , method = "spearman")
cor.test(vector1 , vector2 , method = "spearman")
plot(vector1.var, vector2.var)

.txt文件仅包含数字值。

我得到两个错误，第4行第一个就像“'x'必须是一个数字向量”，第二个错误发生在第5行，就像“无法找到对象向量1.var”

我也试过了

 plot(vector1, vector2)

代替

 plot(vector1.var, vector2.var)

但是有一个错误，如“stripchart.default（x1，...）中的错误：无效的plot-method

实施方向见http://www.gardenersown.co.uk/Education/Lectures/R/correl.htm#correlation

Answer 1

str是一个非常有用的函数（参见?str for more），应该经常使用，特别是验证R对象类型。 快速str(vector1)和str(vector2)将告诉您这些列是否被读取为字符而不是数字。 如果是，则使用as.numeric(vector1)来对每个向量中的数据进行类型转换。

此外， names(vector1)和names(vector2)将告诉您列名称是什么，并可能解决您的绘图问题。

Answer 2

我怀疑vector1和vector2是向量。 阅读?read.table我们在Value部分注意到：

值：

  A data frame ('data.frame') containing a representation of the data in the file. 

....

因此，即使您的两个文本文件只包含一个变量，读入的两个对象也将是每个都包含一个组件的数据框。

其次，您的数据文件不包含标题，因此R将组成变量名称。 我没有测试过这个但是IIRC你的vector1和vector2的变量都被称为X1 。 在vector2 （或names(vector1) ）上执行head(vector1)和相同操作，以查看对象在R中的外观。

我可以看出你为什么认为vector1.var可能有效，但你应该意识到，就R而言，它正在寻找一个名为vector1.var的对象。 的. 是R对象名称中的任何其他字符。 如果你打算用. 作为子集或选择运算符，那么你需要阅读R中的子集运算符。这些是$和[和[[ 。 例如，参见R语言定义手册或R手册 。

我怀疑您可以将代码更改为：

vector1 <- read.table ("D:...path.../mytext1.txt", header=FALSE)[, 1]
vector2 <- read.table ("D:...path.../mytext2.txt", header=FALSE)[, 1]
cor.coeff <- cor(vector1 , vector2 , method = "spearman")
cor.test(vector1 , vector2 , method = "spearman")
plot(vector1, vector2)

但我想你的两个文本文件中有什么...