Mercurial > repos > davidvanzessen > mutation_analysis
changeset 33:ac9a4307861a draft
Uploaded
| author | davidvanzessen |
|---|---|
| date | Thu, 16 Apr 2015 09:36:19 -0400 |
| parents | 2a7343e4be5a |
| children | d436daae9d68 |
| files | mutation_analysis.py mutation_analysis.r |
| diffstat | 2 files changed, 2 insertions(+), 2 deletions(-) [+] |
line wrap: on
line diff
--- a/mutation_analysis.py Thu Apr 09 08:58:59 2015 -0400 +++ b/mutation_analysis.py Thu Apr 16 09:36:19 2015 -0400 @@ -168,4 +168,4 @@ o.write("ID\tRGYWC\tWRCY\tWA\tTW\n") first = False continue - o.write(ID + "\t" + str(int(RGYWCount[ID])) + "\t" + str(int(WRCYCount[ID])) + "\t" + str(int(WACount[ID])) + "\t" + str(int(TWCount[ID])) + "\n") + o.write(ID + "\t" + str(round(RGYWCount[ID], 2)) + "\t" + str(round(WRCYCount[ID], 2)) + "\t" + str(round(WACount[ID], 2)) + "\t" + str(round(TWCount[ID], 2)) + "\n")
--- a/mutation_analysis.r Thu Apr 09 08:58:59 2015 -0400 +++ b/mutation_analysis.r Thu Apr 16 09:36:19 2015 -0400 @@ -284,7 +284,7 @@ result = data.frame(matrx) -row.names(result) = c("Number of Mutations (%)", "Transition (%)", "Transversions (%)", "Transitions at G C (%)", "Targeting of C.G (%)", "FR S/R (ratio)", "CDR S/R (ratio)") +row.names(result) = c("Number of Mutations (%)", "Transition (%)", "Transversions (%)", "Transitions at G C (%)", "Targeting of C.G (%)", "FR R/S (ratio)", "CDR R/S (ratio)") write.table(x=result, file="mutations.txt", sep=",",quote=F,row.names=T,col.names=F)