用pvalues标记热图

Question

我有一组回归模型的β系数矩阵和相关的p值。我想用星号或星号绘制一张贝塔-瓦莱乌斯的热图，其数值达到显著水平。最好是两颗星，一颗星用于标定邦费罗尼阈值，另一颗星用于标定FDR。如何用符号或其他方式注释热图.2中的值？我也很乐意接受一个海生解决方案

#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")