sc.tl.rank_genes_groups()问题

今天被问到了一个关于sc.tl.rank_genes_groups()的奇怪的问题

复制代码
import scanpy as sc
import pandas as pd
import numpy as np
import seaborn as sns
import matplotlib.pyplot as plt 
# from CellDART import da_cellfraction
# from CellDART.utils import random_mix
from sklearn.manifold import TSNE

ref_adata = sc.read_h5ad("./scRNA.h5ad")
ref_adata
复制代码
ref_adata.var_names_make_unique()  
sc.pp.normalize_total(ref_adata)
sc.pp.log1p(ref_adata)

#PCA and clustering 
sc.tl.pca(ref_adata, svd_solver='arpack')
sc.pp.neighbors(ref_adata, n_neighbors=6, n_pcs=40)
sc.tl.umap(ref_adata)
sc.tl.leiden(ref_adata, resolution = 0.5)
sc.pl.umap(ref_adata, color=['leiden','cellType'])
复制代码
sc.tl.rank_genes_groups(ref_adata, 'cellType', method='wilcoxon')
sc.pl.rank_genes_groups(ref_adata, n_genes=15, sharey=False)


可以看到这里的检验结果是以数字的形式显示,而不是以基因的形式显示的

解决办法

这个需要注意,有时间可以研究一下sc.tl.rank_genes_groups()内部是怎么实现的

相关推荐
李昊哲小课1 分钟前
pandas销售数据分析
人工智能·python·数据挖掘·数据分析·pandas
C嘎嘎嵌入式开发22 分钟前
python之set详谈
开发语言·python
之歆1 小时前
Python-正则表达式-信息提取-滑动窗口-数据分发-文件加载及分析器-浏览器分析-学习笔记
python·学习·正则表达式
往日情怀酿做酒 V17639296381 小时前
pytorch的介绍以及张量的创建
人工智能·pytorch·python
豌豆花下猫2 小时前
Python 潮流周刊#110:JIT 编译器两年回顾,AI 智能体工具大爆发(摘要)
后端·python·ai
June bug2 小时前
【Python基础】变量、运算与内存管理全解析
开发语言·python·职场和发展·测试
蹦蹦跳跳真可爱5893 小时前
Python----OpenCV(几何变换--图像平移、图像旋转、放射变换、图像缩放、透视变换)
开发语言·人工智能·python·opencv·计算机视觉
蹦蹦跳跳真可爱5893 小时前
Python----循环神经网络(Transformer ----Layer-Normalization(层归一化))
人工智能·python·rnn·transformer
m0_625686553 小时前
Day58
python
夜阳朔3 小时前
Conda环境激活失效问题
人工智能·后端·python