Workshop
Workshop copied to clipboard
十五期热图——标记相关系数热图显著性
原代码
p_value <- ifelse(is.na(cortest$P),"",cortest$P)
pheatmap(cortest$r,display_numbers = matrix(ifelse(p_value < 0.05,"*",""),nrow(p_value)))
做出来的结果会导致原矩阵里面的科学记数法数字错误识别,且提前制造的空白值会被判定<0.05。
或许可以把原步骤合并以避免该问题,如下:
p_value <- cortest$P
pheatmap(cortest$r,display_numbers = matrix(ifelse(is.na(p_value),"",ifelse(p_value < 0.05,"*","")),nrow(p_value)))