image.png方法来源于上面这篇文章,不得不说,这篇文章运用了非常多复杂的方法去阐述关注的科学问题,真不愧是出自Broadinstitute实验室的。我这里暂时只讲下文章中一种比较新颖的比较不同条件下亚群丰度变化的方法。首先我们先了解下Dirichlet-multinomialregression。让我们从数学层面开始:假设从正常组织取了samplei,正常组织本身包含了p种celltype,假设各种celltype出现的概率为Cellprobability,samplei中各种celltype出现的数目为Numberofcelltypesamplei中共有N个cell,全部细胞的总和就是