单细胞转录组如何将UCell CytoTRACE2 infercnv等结果添加到总的rds中

单细胞转录组如何将UCell CytoTRACE2 infercnv等结果添加到总的rds中


        选取需要的列,之前有的相同的列信息去除。保留第一列barcode,方便合并;


cut -f 1,27- UCell.metadata.tsv >UCell.metadata1.tsv
cut -f 1,27- infercnv.metadata.tsv >infercnv.metadata1.tsv
cut -f 1,27- CytoTRACE2.metadata.tsv >CytoTRACE2.metadata1.tsv

attachments-2024-08-at6uWP2n66d0086ca2f0e.png




合并这些信息: -b barcode 指定共同的列名字

Rscript $scripts/merge_tsv_files.r -i infercnv.metadata1.tsv UCell1.metadata1.tsv CytoTRACE2.metadata1.tsv -b barcode -p metadata




添加metadata信息到rds文件中:

Rscript $scripts/addMetadata_to_seurat_obj.r  -i GSE160269.rds -m metadata.tsv -p GSE160269.add 
  • 发表于 2024-08-29 13:40
  • 阅读 ( 572 )
  • 分类:转录组

0 条评论

请先 登录 后评论
omicsgene
omicsgene

生物信息

707 篇文章

作家榜 »

  1. omicsgene 707 文章
  2. 安生水 352 文章
  3. Daitoue 167 文章
  4. 生物女学霸 120 文章
  5. xun 82 文章
  6. rzx 79 文章
  7. 红橙子 78 文章
  8. CORNERSTONE 72 文章