|
47
|
1 args <- commandArgs(trailingOnly = TRUE)
|
|
|
2
|
|
|
3 naive.file = args[1]
|
|
|
4 shm.file = args[2]
|
|
69
|
5 output.file.ca = args[3]
|
|
|
6 output.file.cg = args[4]
|
|
|
7 output.file.cm = args[5]
|
|
47
|
8
|
|
|
9 naive = read.table(naive.file, sep="\t", header=T, quote="", fill=T)
|
|
|
10 shm.merge = read.table(shm.file, sep="\t", header=T, quote="", fill=T)
|
|
|
11
|
|
|
12
|
|
|
13 final = merge(naive, shm.merge[,c("Sequence.ID", "best_match")], by.x="ID", by.y="Sequence.ID")
|
|
|
14 print(paste("nrow final:", nrow(final)))
|
|
|
15 names(final)[names(final) == "best_match"] = "Sample"
|
|
|
16 final.numeric = final[,sapply(final, is.numeric)]
|
|
|
17 final.numeric[is.na(final.numeric)] = 0
|
|
|
18 final[,sapply(final, is.numeric)] = final.numeric
|
|
69
|
19
|
|
|
20 final.ca = final[grepl("ca", final$Sample),]
|
|
|
21 final.cg = final[grepl("cg", final$Sample),]
|
|
|
22 final.cm = final[grepl("cm", final$Sample),]
|
|
|
23
|
|
|
24
|
|
47
|
25
|
|
69
|
26 #print(paste("nrow final:", nrow(final)))
|
|
|
27 #final2 = final
|
|
|
28 #final2$Sample = gsub("[0-9]", "", final2$Sample)
|
|
|
29 #final = rbind(final, final2)
|
|
|
30 #final$Replicate = 1
|
|
|
31
|
|
|
32 write.table(final.ca, output.file.ca, quote=F, sep="\t", row.names=F, col.names=T)
|
|
|
33 write.table(final.cg, output.file.cg, quote=F, sep="\t", row.names=F, col.names=T)
|
|
|
34 write.table(final.cm, output.file.cm, quote=F, sep="\t", row.names=F, col.names=T)
|
