CellPhoneDB的使用

CellPhoneDB主要分成method、plot、query和database4个模块。

安装成功CellPhoneDB后,要先下载CellphoneDB database。下载的地方在cellphoneDB的github里。https://github.com/ventolab/CellphoneDB

练习运行:

在github上下载了test data文件:test_meta.txt,test_counts.txt

在linux命令行输入python3进入python程序,

from cellphonedb.src.core.methods import cpdb_analysis_method

means, deconvoluted = cpdb_analysis_method.call(cpdb_file_path = "/mnt/data/ref/cellphonedb/cellphonedb.zip", meta_file_path = "/mnt/data/ref/cellphonedb/test_meta.txt", counts_file_path = "/mnt/data/ref/cellphonedb/test_counts.txt", counts_data = 'ensembl',output_path = "/mnt/data/ref/a")

注意:

test_counts.txt格式,第一行为:Gene  cell1  cell2...

test_meta.txt:

Cell cell_type

d-pos_AAACCTGAGCAGGTCA  NKcells_1

d-pos_AAACCTGGTACCGAGA  NKcells_0

d-pos_AAACCTGTCGCCATAA  NKcells_1

d-pos_AAACGGGTCAGTTGAC  Tcells

d-pos_AAAGATGCATTGAGCT  NKcells_0

d-pos_AAAGATGTCCAAAGTC  NKcells_0

d-pos_AAAGCAAAGAGGACGG  Myeloid

d-pos_AAAGCAACACATTCGA  NKcells_1

d-pos_AAAGTAGAGAGCCCAA  NKcells_0

d-pos_AAAGTAGCAAGCTGAG  NKcells_0

counts_data选项,要根据test_counts.txt的基因名为哪种形式来选择ensembl,gene_name ,hgnc_symbol

运行成功后,在输出目录有simple_analysis_means_result_xxx.txt和simple_analysis_deconvoluted_result_xxx.txt文件产生


CellPhoneDB可视化参考:

https://zhuanlan.zhihu.com/p/552055550

©著作权归作者所有,转载或内容合作请联系作者
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

推荐阅读更多精彩内容