python - 用 pvalues 标记热图
问题描述
我有一个来自我运行的一组回归模型和相关 p 值的 beta 系数矩阵。我想制作一个带有星号或星号的 Beta valeus 热图,这些值具有显着性。最好是两颗星,一颗用于划分bonferonni 阈值,另一颗用于划分FDR。如何用符号或其他方式注释 heatmap.2 中的值?我也很乐意接受 seaborn 解决方案。
#example of betas
> bees[1:3, 1:3]
smri_thick_cdk_banksstslh smri_thick_cdk_cdacatelh smri_thick_cdk_cdmdfrlh
Height 0.61403585 0.9574614 0.8578689
Weight 0.07517904 0.6418907 0.2238402
Waist Circumference 0.07336446 0.4155540 0.1593594
#sample of p-values
> peas[1:3, 1:3]
smri_thick_cdk_banksstslh smri_thick_cdk_cdacatelh smri_thick_cdk_cdmdfrlh
Height 0.61403585 0.9574614 0.8578689
Weight 0.07517904 0.6418907 0.2238402
Waist Circumference 0.07336446 0.4155540 0.1593594
#example heatmap
heatmap.2( x=as.matrix(bees[1:3, 1:3]),
dendrogram="none",
keysize = 1,
key.title = NA,
key.ylab = NA,
density.info = "none",
trace = "none",
Colv = "Rowv",
scale = "none",
main = "betas")