cellxgene_VIP icon indicating copy to clipboard operation
cellxgene_VIP copied to clipboard

functions for DEG analysis

Open jianglch opened this issue 4 months ago • 1 comments

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!

jianglch avatar Mar 07 '24 03:03 jianglch

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.

jianglch avatar Mar 07 '24 07:03 jianglch