前几天分析了一批单细胞TCR的数据,需要画这样一个图:
行是不同的样本,列是不同的T细胞类型,每一个位置点的大小表示T细胞的数目多少,并且还需要根据T细胞所属的克隆型类型涂色。这个图用来描述TCR是我自己构思的,之前没有在文章中见过。好在画图工具是现成的,用的是余老师的scatterpie
包。后台回复20210722
即可获取本文的测试数据。
1. 导入测试数据
library(tidyverse)
library(scatterpie)
library(reshape2)
df1=read.table("test.txt",header = T,row.names = 1,sep = "\t",stringsAsFactors = F)
看一下数据格式
> head(df1)
sample celltype clonal_expansion
TTAGTTCGTTCGGCAC tumor5 E 1
ATAACGCGTGAGGCTA tumor5 A morethan3
CTCTAATCAACACCTA tumor1 B 1
AGACGTTAGTGTGGCA tumor5 A 1
AGATCTGTCACGGTTA tumor5 C 1
TCAGCTCTCGGCGCAT tumor3 A 1
> table(df1$samp