HelenaLC/muscat

pbDS(pb, design = mm, contrast = contrast) error

Opened this issue · 2 comments

Jadeapple commented

I ran pbDS(pb, design = mm, contrast = contrast) and got the following error:
Error in pbDS(pb, design = mm, contrast = contrast) :
Specified filtering options result in no genes in any clusters being tested. To force testing, consider modifying arguments 'min_cells' and/or 'filter'. See '?pbDS' for details.

then I ran pbDS(pb, design = mm, contrast = contrast,min_cells=1,filter="none"), and the same error appeared. what should I do?

HelenaLC commented

Independent of filter=..., the function will drop samples with fewer than min_cells in a the cluster being tested. Then, the rank of the design matrix is checked and at least two samples need to be present per group. So... my intuition is that your clusters aren't represented by at least two replicates in both groups. Maybe you could screenshot & post something like table(sce$sample_id, sce$group_id) and table(sce$sample_id, sce$cluster_id) to better understand what's causing the issue?

Jadeapple commented

Independent of filter=..., the function will drop samples with fewer than min_cells in a the cluster being tested. Then, the rank of the design matrix is checked and at least two samples need to be present per group. So... my intuition is that your clusters aren't represented by at least two replicates in both groups. Maybe you could screenshot & post something like table(sce$sample_id, sce$group_id) and table(sce$sample_id, sce$cluster_id) to better understand what's causing the issue?

you are right，I have no replicates, only one sample per group. So I can't use this method to do differential analysis?