# HG changeset patch # User davidvanzessen # Date 1491985696 14400 # Node ID cfc9a442e59d72371a2961cce688172e1d0f1a2f # Parent b66e8946ba756853079c247c813385baf8ec116b Uploaded diff -r b66e8946ba75 -r cfc9a442e59d sequence_overview.r --- a/sequence_overview.r Wed Apr 12 03:52:33 2017 -0400 +++ b/sequence_overview.r Wed Apr 12 04:28:16 2017 -0400 @@ -85,7 +85,7 @@ cat("", file=main.html, append=T) cat("SequenceFunctionalityIGA1IGA2IGG1IGG2IGG3IGG4IGMIGEUN", file=main.html, append=T) -cat("total IGAtotal IGGtotal IGEtotal IGMnumber of subclassespresent in both IGA and IGGpresent in IGA, IGG and IGMpresent in IGA, IGG and IGEpresent in IGA, IGG, IGM and IGEIGA1+IGA2", file=main.html, append=T) +cat("total IGAtotal IGGtotal IGMtotal IGEnumber of subclassespresent in both IGA and IGGpresent in IGA, IGG and IGMpresent in IGA, IGG and IGEpresent in IGA, IGG, IGM and IGEIGA1+IGA2", file=main.html, append=T) cat("IGG1+IGG2IGG1+IGG3IGG1+IGG4IGG2+IGG3IGG2+IGG4IGG3+IGG4", file=main.html, append=T) cat("IGG1+IGG2+IGG3IGG2+IGG3+IGG4IGG1+IGG2+IGG4IGG1+IGG3+IGG4IGG1+IGG2+IGG3+IGG4", file=main.html, append=T) cat("", file=main.html, append=T) @@ -147,7 +147,9 @@ next } - in.classes = sum(classes > 0) + classes.no.un = classes[-length(classes)] + + in.classes = sum(classes.no.un > 0) matched = matched + in.classes #count in how many subclasses the sequence occurs. @@ -249,6 +251,8 @@ #rw = c(as.character(dat[i,"seq_conc"]), functionality, ca1.html, ca2.html, cg1.html, cg2.html, cg3.html, cg4.html, cm.html, un.html) rw = c(as.character(dat[i,"seq_conc"]), functionality, ca1.html, ca2.html, cg1.html, cg2.html, cg3.html, cg4.html, cm.html, ce.html, un.html) rw = c(rw, ca.n, cg.n, cm.n, ce.n, in.classes, in.ca.cg, in.ca.cg.cm, in.ca.cg.ce, in.ca.cg.cm.ce, in.ca1.ca2, in.cg1.cg2, in.cg1.cg3, in.cg1.cg4, in.cg2.cg3, in.cg2.cg4, in.cg3.cg4, in.cg1.cg2.cg3, in.cg2.cg3.cg4, in.cg1.cg2.cg4, in.cg1.cg3.cg4, in.cg.all) + + cat(tr(rw), file=main.html, append=T)