aidarripoll

Hi, in line with @zhangguy comment, do you have some updates on the incorporation of mixed effect models? I'd like to know if the latest version of diffxpy supports the...

Hi @davidsebfischer, nice to hear this!! I think it could be a good approach to add this feature to your package, both in terms of time and work! I will...

As I mentioned, when calling `variancePartition::dream()` in a `exprObj` of 21,315 genes (rows) and 981 donors (cols), it stops due to this error. I'm using `try()` to skip the error...

Hi, I attach a **reproducible example** which uses **testing data** host at [Molgenis server](https://molgenis26.gcc.rug.nl/downloads/marta-transcriptomics/limma_pseudobulk_error/): * [geneExpr.rds](https://molgenis26.gcc.rug.nl/downloads/marta-transcriptomics/limma_pseudobulk_error/geneExpr.rds): DGEList object --> 500 genes x 981 samples * [aggregate_metadata.rds](https://molgenis26.gcc.rug.nl/downloads/marta-transcriptomics/limma_pseudobulk_error/aggregate_metadata.rds): data.frame object --> Sample...

Hi again, As we had very few cases showing this type of error, I used the workaround you suggested to catch and skip these cases by passing to `dream()` the...

Thanks for your quick answer! I was not considering that `eBayes()` was gathering info across all genes, so, it makes a lot of sense what you've explained. Then, I guess...

Nice, the variancePartition I was using didn't give this warning, so, this is much more informative :) Then, the options would be: 1) If we want to remove the cases...

Wow, I can't wait to read it! Btw, when using the **pseudobulk mean data** for the [same subset of genes in the big dataset](https://molgenis26.gcc.rug.nl/downloads/marta-transcriptomics/limma_pseudobulk_error/geneExpr.rds) (from previously normalized scRNA-seq data), instead...