The final log likelihood could be used to compare pooling modes. The lnL is reported for each EM, so once for each cell (individual), celltype (celltype) or all (pseudobulk). This could be easily implemented by returning the lnL with TelescopeLikelihood.em() and summing.
The final log likelihood could be used to compare pooling modes. The lnL is reported for each EM, so once for each cell (individual), celltype (celltype) or all (pseudobulk). This could be easily implemented by returning the lnL with
TelescopeLikelihood.em()
and summing.