我目前正在使用 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
关于r - 在 R 中的相关矩阵上应用阈值,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/26487376/