w= as.numeric()
for(i in 2:10)){
w[i] = cor(a[,1], a[,i], method = "pearson", use = "complete.obs")
}
是否有办法将以下 for 循环转换为一行代码,如下所示?
w[1:9] = cor(a[,1], a[,i], method = "pearson", use = "complete.obs")
最佳答案
关于:
w = sapply(c(2:10), function(x) cor(a[,1], a[,x], method = "pearson", use = 'pairwise.complete.obs') )
假设您使用 2:10
的原因是因为有 10 列。如果是这种情况并且您希望它更加动态:
w = sapply(c(2:ncol(a)), function(x) cor(a[,1], a[,x], method = "pearson", use = "complete.obs") )
最后,如果您不想使用 apply 系列,您可以生成整个 cor
矩阵并从第一行中提取您需要的内容:
w = cor(a, method = 'pearson', use = 'pairwise.complete.obs')[1,-1]
可重现的示例:
> a = iris[,1:4]
> cor(a, method = 'pearson', use = 'pairwise.complete.obs')[1,-1]
Sepal.Width Petal.Length Petal.Width
-0.1175698 0.8717538 0.8179411
编辑: 如果您想检索不同列的相关性(例如 2):
> cor(a, method = 'pearson', use = 'pairwise.complete.obs')[2,-2]
Sepal.Length Petal.Length Petal.Width
-0.1175698 -0.4284401 -0.3661259
此外,根据 Joris Meys 的建议,请确保使用 'pairwise.complete.obs'
关于r - 计算一列与选择的其他列之间的相关性,无需循环,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/35317936/