2022-04-18

将mixcr转换格式后的文件进行拼接

setwd("C:\\Users\\Administrator.DESKTOP-4UQ3Q0K\\Downloads\\vdjtools-1.2.1\\vdjtools-1.2.1\\弦图Convert")
a1<-read.table("Convert.A1.TRA.txt")[,4]
b1<-read.table("Convert.A1.TRB.txt")[,4]
a2<-read.table("Convert.A2.TRA.txt")[,4]
b2<-read.table("Convert.A2.TRB.txt")[,4]
a3<-read.table("Convert.A3.TRA.txt")[,4]
b3<-read.table("Convert.A3.TRB.txt")[,4]
a4<-read.table("Convert.A4.TRA.txt")[,4]
b4<-read.table("Convert.A4.TRB.txt")[,4]
a5<-read.table("Convert.A5.TRA.txt")[,4]
b5<-read.table("Convert.A5.TRB.txt")[,4]
a6<-read.table("Convert.A6.TRA.txt")[,4]
b6<-read.table("Convert.A6.TRB.txt")[,4]
a7<-read.table("Convert.A7.TRA.txt")[,4]
b7<-read.table("Convert.A7.TRB.txt")[,4]
a8<-read.table("Convert.A8.TRA.txt")[,4]
b8<-read.table("Convert.A8.TRB.txt")[,4]
a9<-read.table("Convert.A9.TRA.txt")[,4]
b9<-read.table("Convert.A9.TRB.txt")[,4]
a10<-read.table("Convert.A10.TRA.txt")[,4]
b10<-read.table("Convert.A10.TRB.txt")[,4]
a11<-read.table("Convert.A11.TRA.txt")[,4]
b11<-read.table("Convert.A11.TRB.txt")[,4]
a12<-read.table("Convert.A12.TRA.txt")[,4]
b12<-read.table("Convert.A12.TRB.txt")[,4]

library(tidyverse)
N.TRA<-c(a1,a2,a3,a4,a5,a6)
N.TRB<-rbind(A1.TRB,A2.TRB,A3.TRB,A4.TRB,A5.TRB,A6.TRB)
T.TRA<-rbind(A7.TRA,A8.TRA,A9.TRA,A10.TRA,A11.TRA,A12.TRA)
T.TRB<-rbind(A7.TRB,A8.TRB,A9.TRB,A10.TRB,A11.TRB,A12.TRB)

setwd("D:\\mixcr-3.0.13\\弦图\\拼接")
write.table (N.TRA, file ="N.A序列.txt", sep ="\t", row.names =FALSE, col.names =TRUE, quote =FALSE)
write.table (N.TRB, file ="N.TRB.txt", sep ="\t", row.names =FALSE, col.names =TRUE, quote =FALSE)
write.table (T.TRA, file ="T.TRA.txt", sep ="\t", row.names =FALSE, col.names =TRUE, quote =FALSE)
write.table (T.TRB, file ="T.TRB.txt", sep ="\t", row.names =FALSE, col.names =TRUE, quote =FALSE)