cellxgene_VIP
cellxgene_VIP copied to clipboard
functions for DEG analysis
Dear cellxgene_vip authors, I am looking into the script used for DEG analysis VIPInterface.py and can't figure out where function "compute_diffexp_ttest" come from. I guess it is corresponding to the cellxgene t-test option in the frontend. scripts relevant: if data['DEmethod']=='default': if sum(mask[0]==True)<10 or sum(mask[1]==True)<10: raise ValueError('Less than 10 cells in a group!') #with app.get_data_adaptor(url_dataroot=data['url_dataroot'],dataset=data['dataset']) as scD: scD=data['data_adapter'] if scD is not None: #res = diffDefault.diffexp_ttest(scD,mask[0].to_numpy(),mask[1].to_numpy(),scD.data.shape[1])# shape[cells as rows, genes as columns] res = scD.compute_diffexp_ttest(mask[0].to_numpy(),mask[1].to_numpy(),scD.data.shape[1]-1,0.01)
Thanks for your input!
should be someting linked to https://github.com/sinnamone/cellxgene/blob/main/backend/czi_hosted/compute/diffexp_generic.py and https://gitlab-ce.gwdg.de/loosolab/container/cellxgene/-/blob/seve/standard-endpoint-naming/server/data_anndata/anndata_adaptor.py. I will read the code more carefully.