For the model comparison benchmark, we also included effects on two different cell types. For this, we assumed \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{{{{\boldsymbol{\mu }}}}}}={{{{{\mathrm{log}}}}}}(1000,1000,\ldots ,1000)$$\end{document}μ=log(1000,1000,…,1000) for all control samples, and an increase of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{{{{\boldsymbol{\mu }}}}}}^{\prime} =({\mu }_{1}^{{\prime} },{\mu }_{2}^{{\prime} })$$\end{document}μ′=(μ1′,μ2′) on the first two cell types, leading to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{{{{\boldsymbol{\mu }}}}}}={{{{{\mathrm{log}}}}}}(1000+{\mu }_{1}^{{\prime} },1000+{\mu }_{2}^{{\prime} },\frac{K\cdot 1000\mbox{-}(2000+{\mu }_{1}^{{\prime} }+{\mu }_{2}^{{\prime} })}{K-2})$$\end{document}μ=log(1000+μ1′,1000+μ2′,K⋅1000-(2000+μ1′+μ2′)K−2).