【发布时间】:2022-11-02 13:29:03
【问题描述】:
我试图计算和绘制相同的jaccard指数
我的虚拟数据和代码
df <- data.frame(tumor_type = c("tumor1", "tumor1", "tumor1", "tumor2", "tumor2", "tumor3", "tumor4", "tumor4"), genes = c("geneA", "geneB", "geneC", "geneA", "geneD", "geneD", "geneA", "geneD"))
D1 <- vegdist(unclass(table(df$tumor_type,df$genes)), method = "jaccard")
pheatmap(1-as.matrix(D1))
对于我的实际数据,我尝试了一个在线工具,我得到了这个
我的目标是获得相似的结果,例如显示每个子类型中存在的元素及其 jacard 分数
【问题讨论】: