日韩性视频-久久久蜜桃-www中文字幕-在线中文字幕av-亚洲欧美一区二区三区四区-撸久久-香蕉视频一区-久久无码精品丰满人妻-国产高潮av-激情福利社-日韩av网址大全-国产精品久久999-日本五十路在线-性欧美在线-久久99精品波多结衣一区-男女午夜免费视频-黑人极品ⅴideos精品欧美棵-人人妻人人澡人人爽精品欧美一区-日韩一区在线看-欧美a级在线免费观看

歡迎訪問 生活随笔!

生活随笔

當前位置: 首頁 > 编程资源 > 编程问答 >内容正文

编程问答

scanpy和umap版本不匹配问题总结

發布時間:2023/12/20 编程问答 37 豆豆
生活随笔 收集整理的這篇文章主要介紹了 scanpy和umap版本不匹配问题总结 小編覺得挺不錯的,現在分享給大家,幫大家做個參考.

我目前python所有的環境基本上都是python3.6.10, 這里安裝的scanpy默認是1.7.2,但是安裝完,你會發現umap圖出現很奇怪的現象,下面是測試的結果

%load_ext rpy2.ipython %%R -o counts -o metasuppressMessages(library(splatter))params <- newSplatParams()params <- setParam(params, "nGenes", 5000) params <- setParam(params, "batchCells", c(500,500,500)) params <- setParam(params, "batch.facLoc", 0.5) params <- setParam(params, "batch.facScale", 0.5) params <- setParam(params, "group.prob", c(1/3,1/3,1/3))sim <- splatSimulate(params, method="groups", verbose=FALSE)counts = data.frame(counts(sim)) meta = data.frame(colData(sim)) import anndata import scanpy as sc adata = anndata.AnnData(X=counts.values.T, obs=meta)sc.pp.normalize_per_cell(adata,counts_per_cell_after=10000) sc.pp.log1p(adata) sc.tl.pca(adata,svd_solver='arpack') sc.pp.neighbors(adata, n_neighbors=9) sc.tl.umap(adata) sc.pl.umap(adata,color=['Batch','Group'])

結果如下
這個圖很明顯是不對的,而且louvain也沒有安裝

解決方法如下

conda install -c vtraag louvain # 可選,因為louvain 也要安裝的 pip uninstall umap-learn pip install umap-learn==0.4.6 pip uninstall numba pip install numba==0.51.2

重新測試代碼

%load_ext rpy2.ipython%%R -o counts -o metasuppressMessages(library(splatter))params <- newSplatParams()params <- setParam(params, "nGenes", 5000) params <- setParam(params, "batchCells", c(500,500,500)) params <- setParam(params, "batch.facLoc", 0.5) params <- setParam(params, "batch.facScale", 0.5) params <- setParam(params, "group.prob", c(1/3,1/3,1/3))sim <- splatSimulate(params, method="groups", verbose=FALSE)counts = data.frame(counts(sim)) meta = data.frame(colData(sim))import anndata import scanpy as sc adata = anndata.AnnData(X=counts.values.T, obs=meta)sc.pp.normalize_per_cell(adata,counts_per_cell_after=10000) sc.pp.log1p(adata) sc.tl.pca(adata,svd_solver='arpack') sc.pp.neighbors(adata, n_neighbors=9) sc.tl.umap(adata) sc.pl.umap(adata,color=['Batch','Group'])

結果如下
# 然后進行louvain

sc.tl.louvain(adata) sc.pl.umap(adata,color=["louvain"])

結果如下
這個結果才是對的。

總結

以上是生活随笔為你收集整理的scanpy和umap版本不匹配问题总结的全部內容,希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯,歡迎將生活随笔推薦給好友。