Open gighuarhguggg45 opened 7 months ago
Hi @gighuarhguggg45 ,
Thank you for using our package. When you perform the CalcPerterbSig() and RunMixscale(), have you included the "split.by"
flag as well? You will need to re-run them if not, otherwise the Run_wmvRegDE() does not support the "split.by"
flag. Sorry for the confusion and we will definitely make this more explicit in the vignette in our next update!
Please let me know if the issue persists even if "split.by" is included, and I will look into this further.
Hi @longmanz thank you for your quick response! I confirm I included "split.by" in CalcPerturbSig and RunxMixscale but still got this error in Run_wmvRegDE.
Hi @gighuarhguggg45 , Thank you for confirming! Then this might be a bug. I will look into this and get back to you!
Hi @gighuarhguggg45 , I found a bug and has fixed it in the weighted DE function. However, this does not seem to be related to the error you reported. Could you please try the updated version and let me know if the issue persists? If so, it would be nice if you could provide a minimal reproducible example data to me so that I can test it from my end. My email is ljiang@nygenome.org.
In addition, could you check if your "gene" and "RNA_snn_res.0.2" contains NAs ?
Hello @longmanz ,
Thank you for the mixscale package. I have been able to run the vignette example as well as my own data without issues. However, when I include the split.by argument in all the functions, I get the following error when I get to the Run_wmvRegDE function:
de_res = Run_wmvRegDE(object = seurat_obj, assay = "RNA", slot = "counts", labels = "gene", nt.class.name = "negative_control", PRTB_list = c("GENE_NAME"), logfc.threshold = 0.2, split.by = 'RNA_snn_res.0.2')
Error message:
sessionInfo:
Is there a way to modify that design matrix manually? Could this be related to issues with dependencies?
Thanks