我目前正在使用R处理数据集。我为变量创建了相关性martix(Pearson)。但是现在我想为矩阵中显示的值设置一个阈值。我正在尝试以下代码:
cor_relation = cor(mydata_frame, use="all.obs", method="pearson")
我得到以下输出:
200605_s_at 202592_at 202958_at
200605_s_at 1.000000000 0.295065389 0.169772244
202592_at 0.695065389 1.000000000 -0.534394180
202958_at 0.869772244 -0.534394180 1.000000000
我想找到以下输出(当我将阈值设置为0.6时):
200605_s_at 202592_at 202958_at
200605_s_at 1.000000000 NA NA
202592_at 0.695065389 1.000000000 NA
202958_at 0.869772244 NA 1.000000000
在此先感谢您的帮助!
备用:
cor_relation[abs(cor_relation) < 0.6] <- NA
本文收集自互联网,转载请注明来源。
如有侵权,请联系[email protected] 删除。
我来说两句