我正在尝试计算变量(g)与量表中其余变量的相关性。它显示特定变量组的NA。 (s1,s2,s3 ......)
变量中的值如下所示:
g = c(1, NA, 4, NA, 3, NA, 3, NA, NA, NA)
s1 = c(NA, 3, NA, 2, NA, 3, NA, 4, 4, NA)
s2 = c(NA, 3, NA, 2, NA, 1, NA, 1, 4, NA)
s3 = c(NA, 2, NA, 2, NA, 4, NA, 2, 5, NA)
它显示与比例中其他变量的相关性。 这些变量的例子:
as1 = c(4, NA, 4, NA, 3, NA, 3, NA, NA, NA)
ag1 = c(1, NA, 4, NA, 3, NA, 3, NA, NA, NA)
我使用了以下代码:数据集“dn”
x <- dn["g"]
y <- dn[1:20] # all variables
cor(y,x, use="pairwise.complete.obs")
输出:
var_name g
ag1 0.66254379
ag2 -0.76100601
ag3 0.76110091
s1 NA
s2 NA
s3 NA
s4 NA
s5 NA
as1 0.65019200
as2 0.45664999
as3 -0.33246108
我也试过“所有”,“all.obs”,“na.or.complete”和“complete.obs”而不是“pairwise.complete.obs”
我无法弄清楚为什么会出现这些NAs。我一直在检查我的代码中的其他可能的错误,到目前为止没有运气。虽然我继续尝试,有没有其他方法可以让r计算相关性尽管有NA?
注意:所有值均为数字