Mercurial > repos > jfb > kinatest_scoring_function
diff scoring function/KinatestFisherScoring.R @ 0:ed9068208721 draft
Uploaded
| author | jfb |
|---|---|
| date | Fri, 26 Apr 2019 16:36:21 -0400 |
| parents | |
| children | 9ade210d19be |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/scoring function/KinatestFisherScoring.R Fri Apr 26 16:36:21 2019 -0400 @@ -0,0 +1,1261 @@ +# system.time({ + +args = commandArgs(trailingOnly=TRUE) +TodaysKinase<-as.character(args[1]) + +CharTable<-read.csv("Fisher-Char.csv", header = FALSE,stringsAsFactors = FALSE) +ThisKinSheet<-read.csv("Fisher-Table.csv", header=FALSE, stringsAsFactors=FALSE) +NormalizationScore<-CharTable[2,1] + +#so here's the question, do I want this file to be able to score using any EPM table as well? I guess I have to. + + + + +M6<-args[2] +M5<-args[3] +M4<-args[4] +M3<-args[5] +M2<-args[6] +M1<-args[7] +D0<-args[8] +P1<-args[9] +P2<-args[10] +P3<-args[11] +P4<-args[12] +P5<-args[13] +P6<-args[14] + + +A=1 +C=2 +D=3 +E=4 +F=5 +G=6 +H=7 +I=8 +K=9 +L=10 +M=11 +N=12 +P=13 +Q=14 +R=15 +S=16 +T=17 +V=18 +W=19 +Y=20 +aa_props <- c("A"=A, "C"=C, "D"=D, "E"=E, "F"=F,"G"=G,"H"=H,"I"=I,"K"=K,"L"=L,"M"=M,"N"=N,"P"=P,"Q"=Q,"R"=R, + "S"=S,"T"=T,"V"=V,"W"=W,"Y"=Y,"xY"=Y,"O"=21) +Positionm6<-gsub("[^a-zA-Z]", "", M6) +Positionm6<-toupper(Positionm6) +Positionm6<-unlist(strsplit(Positionm6,"")) +Positionm6<-sapply(Positionm6, function(x) aa_props[x]) + +Positionm5<-gsub("[^a-zA-Z]", "", M5) +Positionm5<-toupper(Positionm5) +Positionm5<-unlist(strsplit(Positionm5,"")) +Positionm5<-sapply(Positionm5, function(x) aa_props[x]) + +Positionm4<-gsub("[^a-zA-Z]", "", M4) +Positionm4<-toupper(Positionm4) +Positionm4<-unlist(strsplit(Positionm4,"")) +Positionm4<-sapply(Positionm4, function(x) aa_props[x]) + +Positionm3<-gsub("[^a-zA-Z]", "", M3) +Positionm3<-toupper(Positionm3) +Positionm3<-unlist(strsplit(Positionm3,"")) +Positionm3<-sapply(Positionm3, function(x) aa_props[x]) + +Positionm2<-gsub("[^a-zA-Z]", "", M2) +Positionm2<-toupper(Positionm2) +Positionm2<-unlist(strsplit(Positionm2,"")) +Positionm2<-sapply(Positionm2, function(x) aa_props[x]) + +Positionm1<-gsub("[^a-zA-Z]", "", M1) +Positionm1<-toupper(Positionm1) +Positionm1<-unlist(strsplit(Positionm1,"")) +Positionm1<-sapply(Positionm1, function(x) aa_props[x]) + +Positiond0<-gsub("[^a-zA-Z]", "", D0) +Positiond0<-toupper(Positiond0) +Positiond0<-unlist(strsplit(Positiond0,"")) +Positiond0<-sapply(Positiond0, function(x) aa_props[x]) + +Positionp1<-gsub("[^a-zA-Z]", "", P1) +Positionp1<-toupper(Positionp1) +Positionp1<-unlist(strsplit(Positionp1,"")) +Positionp1<-sapply(Positionp1, function(x) aa_props[x]) + +Positionp2<-gsub("[^a-zA-Z]", "", P2) +Positionp2<-toupper(Positionp2) +Positionp2<-unlist(strsplit(Positionp2,"")) +Positionp2<-sapply(Positionp2, function(x) aa_props[x]) + +Positionp3<-gsub("[^a-zA-Z]", "", P3) +Positionp3<-toupper(Positionp3) +Positionp3<-unlist(strsplit(Positionp3,"")) +Positionp3<-sapply(Positionp3, function(x) aa_props[x]) + +Positionp4<-gsub("[^a-zA-Z]", "", P4) +Positionp4<-toupper(Positionp4) +Positionp4<-unlist(strsplit(Positionp4,"")) +Positionp4<-sapply(Positionp4, function(x) aa_props[x]) + +Positionp5<-gsub("[^a-zA-Z]", "", P5) +Positionp5<-toupper(Positionp5) +Positionp5<-unlist(strsplit(Positionp5,"")) +Positionp5<-sapply(Positionp5, function(x) aa_props[x]) + +Positionp6<-gsub("[^a-zA-Z]", "", P6) +Positionp6<-toupper(Positionp6) +Positionp6<-unlist(strsplit(Positionp6,"")) +Positionp6<-sapply(Positionp6, function(x) aa_props[x]) + + + + +#what the below code does is stupid. I need to make a combinatorial library of peptides and then calculate their values +#but each peptide needs a value calculated for every different kinase. So I first make it so that every AA at every position is labelled differently +#IE an A at position negative 6 is called m6a (minus 6 Alanine) +#then further on I replace valued like m6a with whatever value is found in the EPM table for both this kinase and every other kinase +#then I will multiply those values together + +if(2==2){ +aa_props2 <- c("1"="A", "2"="C", "3"="D", "4"="E", "5"="F", "6"="G", "7"="H", "8"="I", "9"="K", "10"="L", "11"="M", "12"="N", + "13"="P", "14"="Q", "15"="R", "16"="S", "17"="T", "18"="V", "19"="W", "20"="Y", "21"="O") + +m6_props2 <- c("1"="m6A", "2"="m6C", "3"="m6D", "4"="m6E", "5"="m6F", "6"="m6G", "7"="m6H", "8"="m6I", "9"="m6K", "10"="m6L", "11"="m6M", "12"="m6N", + "13"="m6P", "14"="m6Q", "15"="m6R", "16"="m6S", "17"="m6T", "18"="m6V", "19"="m6W", "20"="m6Y", "21"="m6O") + +m5_props2 <- c("1"="m5A", "2"="m5C", "3"="m5D", "4"="m5E", "5"="m5F", "6"="m5G", "7"="m5H", "8"="m5I", "9"="m5K", "10"="m5L", "11"="m5M", "12"="m5N", + "13"="m5P", "14"="m5Q", "15"="m5R", "16"="m5S", "17"="m5T", "18"="m5V", "19"="m5W", "20"="m5Y", "21"="m5O") + +m4_props2 <- c("1"="m4A", "2"="m4C", "3"="m4D", "4"="m4E", "5"="m4F", "6"="m4G", "7"="m4H", "8"="m4I", "9"="m4K", "10"="m4L", "11"="m4M", "12"="m4N", + "13"="m4P", "14"="m4Q", "15"="m4R", "16"="m4S", "17"="m4T", "18"="m4V", "19"="m4W", "20"="m4Y", "21"="m4O") + +m3_props2 <- c("1"="m3A", "2"="m3C", "3"="m3D", "4"="m3E", "5"="m3F", "6"="m3G", "7"="m3H", "8"="m3I", "9"="m3K", "10"="m3L", "11"="m3M", "12"="m3N", + "13"="m3P", "14"="m3Q", "15"="m3R", "16"="m3S", "17"="m3T", "18"="m3V", "19"="m3W", "20"="m3Y", "21"="m3O") + +m2_props2 <- c("1"="m2A", "2"="m2C", "3"="m2D", "4"="m2E", "5"="m2F", "6"="m2G", "7"="m2H", "8"="m2I", "9"="m2K", "10"="m2L", "11"="m2M", "12"="m2N", + "13"="m2P", "14"="m2Q", "15"="m2R", "16"="m2S", "17"="m2T", "18"="m2V", "19"="m2W", "20"="m2Y", "21"="m2O") + +m1_props2 <- c("1"="m1A", "2"="m1C", "3"="m1D", "4"="m1E", "5"="m1F", "6"="m1G", "7"="m1H", "8"="m1I", "9"="m1K", "10"="m1L", "11"="m1M", "12"="m1N", + "13"="m1P", "14"="m1Q", "15"="m1R", "16"="m1S", "17"="m1T", "18"="m1V", "19"="m1W", "20"="m1Y", "21"="m1O") + +d0_props2 <- c("1"="d0A", "2"="d0C", "3"="d0D", "4"="d0E", "5"="d0F", "6"="d0G", "7"="d0H", "8"="d0I", "9"="d0K", "10"="d0L", "11"="d0M", "12"="d0N", + "13"="d0P", "14"="d0Q", "15"="d0R", "16"="d0S", "17"="d0T", "18"="d0V", "19"="d0W", "20"="d0Y", "21"="d0O") + +p1_props2 <- c("1"="p1A", "2"="p1C", "3"="p1D", "4"="p1E", "5"="p1F", "6"="p1G", "7"="p1H", "8"="p1I", "9"="p1K", "10"="p1L", "11"="p1M", "12"="p1N", + "13"="p1P", "14"="p1Q", "15"="p1R", "16"="p1S", "17"="p1T", "18"="p1V", "19"="p1W", "20"="p1Y", "21"="p1O") + +p2_props2 <- c("1"="p2A", "2"="p2C", "3"="p2D", "4"="p2E", "5"="p2F", "6"="p2G", "7"="p2H", "8"="p2I", "9"="p2K", "10"="p2L", "11"="p2M", "12"="p2N", + "13"="p2P", "14"="p2Q", "15"="p2R", "16"="p2S", "17"="p2T", "18"="p2V", "19"="p2W", "20"="p2Y", "21"="p2O") + +p3_props2 <- c("1"="p3A", "2"="p3C", "3"="p3D", "4"="p3E", "5"="p3F", "6"="p3G", "7"="p3H", "8"="p3I", "9"="p3K", "10"="p3L", "11"="p3M", "12"="p3N", + "13"="p3P", "14"="p3Q", "15"="p3R", "16"="p3S", "17"="p3T", "18"="p3V", "19"="p3W", "20"="p3Y", "21"="p3O") + +p4_props2 <- c("1"="p4A", "2"="p4C", "3"="p4D", "4"="p4E", "5"="p4F", "6"="p4G", "7"="p4H", "8"="p4I", "9"="p4K", "10"="p4L", "11"="p4M", "12"="p4N", + "13"="p4P", "14"="p4Q", "15"="p4R", "16"="p4S", "17"="p4T", "18"="p4V", "19"="p4W", "20"="p4Y", "21"="p4O") + +p5_props2 <- c("1"="p5A", "2"="p5C", "3"="p5D", "4"="p5E", "5"="p5F", "6"="p5G", "7"="p5H", "8"="p5I", "9"="p5K", "10"="p5L", "11"="p5M", "12"="p5N", + "13"="p5P", "14"="p5Q", "15"="p5R", "16"="p5S", "17"="p5T", "18"="p5V", "19"="p5W", "20"="p5Y", "21"="p5O") + +p6_props2 <- c("1"="p6A", "2"="p6C", "3"="p6D", "4"="p6E", "5"="p6F", "6"="p6G", "7"="p6H", "8"="p6I", "9"="p6K", "10"="p6L", "11"="p6M", "12"="p6N", + "13"="p6P", "14"="p6Q", "15"="p6R", "16"="p6S", "17"="p6T", "18"="p6V", "19"="p6W", "20"="p6Y", "21"="p6O") + +# Positionm7<-sapply(Positionm7, function (x) aa_props2[x]) +pPositionm6<-sapply(Positionm6, function (x) m6_props2[x]) +pPositionm5<-sapply(Positionm5, function (x) m5_props2[x]) +pPositionm4<-sapply(Positionm4, function (x) m4_props2[x]) +pPositionm3<-sapply(Positionm3, function (x) m3_props2[x]) +pPositionm2<-sapply(Positionm2, function (x) m2_props2[x]) +pPositionm1<-sapply(Positionm1, function (x) m1_props2[x]) +pPositiond0<-sapply(Positiond0, function (x) d0_props2[x]) +pPositionp1<-sapply(Positionp1, function (x) p1_props2[x]) +pPositionp2<-sapply(Positionp2, function (x) p2_props2[x]) +pPositionp3<-sapply(Positionp3, function (x) p3_props2[x]) +pPositionp4<-sapply(Positionp4, function (x) p4_props2[x]) +pPositionp5<-sapply(Positionp5, function (x) p5_props2[x]) +pPositionp6<-sapply(Positionp6, function (x) p6_props2[x]) +# Positionp7<-sapply(Positionp7, function (x) aa_props2[x]) + +Positionm6<-sapply(Positionm6, function (x) aa_props2[x]) +Positionm5<-sapply(Positionm5, function (x) aa_props2[x]) +Positionm4<-sapply(Positionm4, function (x) aa_props2[x]) +Positionm3<-sapply(Positionm3, function (x) aa_props2[x]) +Positionm2<-sapply(Positionm2, function (x) aa_props2[x]) +Positionm1<-sapply(Positionm1, function (x) aa_props2[x]) +Positiond0<-sapply(Positiond0, function (x) aa_props2[x]) +Positionp1<-sapply(Positionp1, function (x) aa_props2[x]) +Positionp2<-sapply(Positionp2, function (x) aa_props2[x]) +Positionp3<-sapply(Positionp3, function (x) aa_props2[x]) +Positionp4<-sapply(Positionp4, function (x) aa_props2[x]) +Positionp5<-sapply(Positionp5, function (x) aa_props2[x]) +Positionp6<-sapply(Positionp6, function (x) aa_props2[x]) + + +} +#this turns the amino acids from above into positionally appropriate amino acids, that is an A at negative 6 is called m6A, if it's at -5 it's m5A, +#if it's a B at positive 4 it's p4B + + + +ScreenerFilename<-"screener.csv" +screaner<-read.csv(ScreenerFilename, header = FALSE, stringsAsFactors = FALSE) +#screaner <- read.csv("F:/ParkerLab Dropbox/Dropbox (Parker Lab)/Parker Lab Shared Files/Dropbox (Parker Lab)/Parker Lab Shared Files/CURRENT LAB MEMBERS/John B/Parker Lab/GalaxyP tools/scoring function KTID/screener.csv", header=FALSE, stringsAsFactors=FALSE) + +Abl<-screaner[2:25,] +Arg<-screaner[27:50,] +Btk<-screaner[52:75,] +Csk<-screaner[77:100,] +Fyn<-screaner[102:125,] +Hck<-screaner[127:150,] +JAK2<-screaner[152:175,] +Lck<-screaner[177:200,] +Lyn<-screaner[202:225,] +Pyk2<-screaner[227:250,] +Src<-screaner[252:275,] +Syk<-screaner[277:300,] +Yes<-screaner[302:325,] +FLT3<-screaner[327:350,] +ALK<-screaner[352:375,] + + + +ThisKinTable<-ThisKinSheet[2:22,] +#NormalizationScore<-ThisKinSheet[27,1] + +if(1==1){ + ThisKinTable_props <- c("m6A"=ThisKinTable[1,3], "m6C"=ThisKinTable[2,3], "m6D"=ThisKinTable[3,3], "m6E"=ThisKinTable[4,3], "m6F"=ThisKinTable[5,3],"m6G"=ThisKinTable[6,3],"m6H"=ThisKinTable[7,3],"m6I"=ThisKinTable[8,3],"m6K"=ThisKinTable[9,3],"m6L"=ThisKinTable[10,3],"m6M"=ThisKinTable[11,3],"m6N"=ThisKinTable[12,3],"m6P"=ThisKinTable[13,3],"m6Q"=ThisKinTable[14,3],"m6R"=ThisKinTable[15,3], + "m6S"=ThisKinTable[16,3],"m6T"=ThisKinTable[17,3],"m6V"=ThisKinTable[18,3],"m6W"=ThisKinTable[19,3],"m6Y"=ThisKinTable[20,3],"m6O"=1, + "m5A"=ThisKinTable[1,4], "m5C"=ThisKinTable[2,4], "m5D"=ThisKinTable[3,4], "m5E"=ThisKinTable[4,4], "m5F"=ThisKinTable[5,4],"m5G"=ThisKinTable[6,4],"m5H"=ThisKinTable[7,4],"m5I"=ThisKinTable[8,4],"m5K"=ThisKinTable[9,4],"m5L"=ThisKinTable[10,4],"m5M"=ThisKinTable[11,4],"m5N"=ThisKinTable[12,4],"m5P"=ThisKinTable[13,4],"m5Q"=ThisKinTable[14,4],"m5R"=ThisKinTable[15,4], + "m5S"=ThisKinTable[16,4],"m5T"=ThisKinTable[17,4],"m5V"=ThisKinTable[18,4],"m5W"=ThisKinTable[19,4],"m5Y"=ThisKinTable[20,4],"m5O"=1, + "m4A"=ThisKinTable[1,5], "m4C"=ThisKinTable[2,5], "m4D"=ThisKinTable[3,5], "m4E"=ThisKinTable[4,5], "m4F"=ThisKinTable[5,5],"m4G"=ThisKinTable[6,5],"m4H"=ThisKinTable[7,5],"m4I"=ThisKinTable[8,5],"m4K"=ThisKinTable[9,5],"m4L"=ThisKinTable[10,5],"m4M"=ThisKinTable[11,5],"m4N"=ThisKinTable[12,5],"m4P"=ThisKinTable[13,5],"m4Q"=ThisKinTable[14,5],"m4R"=ThisKinTable[15,5], + "m4S"=ThisKinTable[16,5],"m4T"=ThisKinTable[17,5],"m4V"=ThisKinTable[18,5],"m4W"=ThisKinTable[19,5],"m4Y"=ThisKinTable[20,5],"m4O"=1, + "m3A"=ThisKinTable[1,6], "m3C"=ThisKinTable[2,6], "m3D"=ThisKinTable[3,6], "m3E"=ThisKinTable[4,6], "m3F"=ThisKinTable[5,6],"m3G"=ThisKinTable[6,6],"m3H"=ThisKinTable[7,6],"m3I"=ThisKinTable[8,6],"m3K"=ThisKinTable[9,6],"m3L"=ThisKinTable[10,6],"m3M"=ThisKinTable[11,6],"m3N"=ThisKinTable[12,6],"m3P"=ThisKinTable[13,6],"m3Q"=ThisKinTable[14,6],"m3R"=ThisKinTable[15,6], + "m3S"=ThisKinTable[16,6],"m3T"=ThisKinTable[17,6],"m3V"=ThisKinTable[18,6],"m3W"=ThisKinTable[19,6],"m3Y"=ThisKinTable[20,6],"m3O"=1, + "m2A"=ThisKinTable[1,7], "m2C"=ThisKinTable[2,7], "m2D"=ThisKinTable[3,7], "m2E"=ThisKinTable[4,7], "m2F"=ThisKinTable[5,7],"m2G"=ThisKinTable[6,7],"m2H"=ThisKinTable[7,7],"m2I"=ThisKinTable[8,7],"m2K"=ThisKinTable[9,7],"m2L"=ThisKinTable[10,7],"m2M"=ThisKinTable[11,7],"m2N"=ThisKinTable[12,7],"m2P"=ThisKinTable[13,7],"m2Q"=ThisKinTable[14,7],"m2R"=ThisKinTable[15,7], + "m2S"=ThisKinTable[16,7],"m2T"=ThisKinTable[17,7],"m2V"=ThisKinTable[18,7],"m2W"=ThisKinTable[19,7],"m2Y"=ThisKinTable[20,7],"m2O"=1, + "m1A"=ThisKinTable[1,8], "m1C"=ThisKinTable[2,8], "m1D"=ThisKinTable[3,8], "m1E"=ThisKinTable[4,8], "m1F"=ThisKinTable[5,8],"m1G"=ThisKinTable[6,8],"m1H"=ThisKinTable[7,8],"m1I"=ThisKinTable[8,8],"m1K"=ThisKinTable[9,8],"m1L"=ThisKinTable[10,8],"m1M"=ThisKinTable[11,8],"m1N"=ThisKinTable[12,8],"m1P"=ThisKinTable[13,8],"m1Q"=ThisKinTable[14,8],"m1R"=ThisKinTable[15,8], + "m1S"=ThisKinTable[16,8],"m1T"=ThisKinTable[17,8],"m1V"=ThisKinTable[18,8],"m1W"=ThisKinTable[19,8],"m1Y"=ThisKinTable[20,8],"m1O"=1, + "d0A"=ThisKinTable[1,9], "d0C"=ThisKinTable[2,9], "d0D"=ThisKinTable[3,9], "d0E"=ThisKinTable[4,9], "d0F"=ThisKinTable[5,9],"d0G"=ThisKinTable[6,9],"d0H"=ThisKinTable[7,9],"d0I"=ThisKinTable[8,9],"d0K"=ThisKinTable[9,9],"d0L"=ThisKinTable[10,9],"d0M"=ThisKinTable[11,9],"d0N"=ThisKinTable[12,9],"d0P"=ThisKinTable[13,9],"d0Q"=ThisKinTable[14,9],"d0R"=ThisKinTable[15,9], + "d0S"=ThisKinTable[16,9],"d0T"=ThisKinTable[17,9],"d0V"=ThisKinTable[18,9],"d0W"=ThisKinTable[19,9],"d0Y"=ThisKinTable[20,9],"d0O"=1, + "p1A"=ThisKinTable[1,10], "p1C"=ThisKinTable[2,10], "p1D"=ThisKinTable[3,10], "p1E"=ThisKinTable[4,10], "p1F"=ThisKinTable[5,10],"p1G"=ThisKinTable[6,10],"p1H"=ThisKinTable[7,10],"p1I"=ThisKinTable[8,10],"p1K"=ThisKinTable[9,10],"p1L"=ThisKinTable[10,10],"p1M"=ThisKinTable[11,10],"p1N"=ThisKinTable[12,10],"p1P"=ThisKinTable[13,10],"p1Q"=ThisKinTable[14,10],"p1R"=ThisKinTable[15,10], + "p1S"=ThisKinTable[16,10],"p1T"=ThisKinTable[17,10],"p1V"=ThisKinTable[18,10],"p1W"=ThisKinTable[19,10],"p1Y"=ThisKinTable[20,10],"p1O"=1, + "p2A"=ThisKinTable[1,11], "p2C"=ThisKinTable[2,11], "p2D"=ThisKinTable[3,11], "p2E"=ThisKinTable[4,11], "p2F"=ThisKinTable[5,11],"p2G"=ThisKinTable[6,11],"p2H"=ThisKinTable[7,11],"p2I"=ThisKinTable[8,11],"p2K"=ThisKinTable[9,11],"p2L"=ThisKinTable[10,11],"p2M"=ThisKinTable[11,11],"p2N"=ThisKinTable[12,11],"p2P"=ThisKinTable[13,11],"p2Q"=ThisKinTable[14,11],"p2R"=ThisKinTable[15,11], + "p2S"=ThisKinTable[16,11],"p2T"=ThisKinTable[17,11],"p2V"=ThisKinTable[18,11],"p2W"=ThisKinTable[19,11],"p2Y"=ThisKinTable[20,11],"p2O"=1, + "p3A"=ThisKinTable[1,12], "p3C"=ThisKinTable[2,12], "p3D"=ThisKinTable[3,12], "p3E"=ThisKinTable[4,12], "p3F"=ThisKinTable[5,12],"p3G"=ThisKinTable[6,12],"p3H"=ThisKinTable[7,12],"p3I"=ThisKinTable[8,12],"p3K"=ThisKinTable[9,12],"p3L"=ThisKinTable[10,12],"p3M"=ThisKinTable[11,12],"p3N"=ThisKinTable[12,12],"p3P"=ThisKinTable[13,12],"p3Q"=ThisKinTable[14,12],"p3R"=ThisKinTable[15,12], + "p3S"=ThisKinTable[16,12],"p3T"=ThisKinTable[17,12],"p3V"=ThisKinTable[18,12],"p3W"=ThisKinTable[19,12],"p3Y"=ThisKinTable[20,12],"p3O"=1, + "p4A"=ThisKinTable[1,13], "p4C"=ThisKinTable[2,13], "p4D"=ThisKinTable[3,13], "p4E"=ThisKinTable[4,13], "p4F"=ThisKinTable[5,13],"p4G"=ThisKinTable[6,13],"p4H"=ThisKinTable[7,13],"p4I"=ThisKinTable[8,13],"p4K"=ThisKinTable[9,13],"p4L"=ThisKinTable[10,13],"p4M"=ThisKinTable[11,13],"p4N"=ThisKinTable[12,13],"p4P"=ThisKinTable[13,13],"p4Q"=ThisKinTable[14,13],"p4R"=ThisKinTable[15,13], + "p4S"=ThisKinTable[16,13],"p4T"=ThisKinTable[17,13],"p4V"=ThisKinTable[18,13],"p4W"=ThisKinTable[19,13],"p4Y"=ThisKinTable[20,13],"p4O"=1, + "p5A"=ThisKinTable[1,14], "p5C"=ThisKinTable[2,14], "p5D"=ThisKinTable[3,14], "p5E"=ThisKinTable[4,14], "p5F"=ThisKinTable[5,14],"p5G"=ThisKinTable[6,14],"p5H"=ThisKinTable[7,14],"p5I"=ThisKinTable[8,14],"p5K"=ThisKinTable[9,14],"p5L"=ThisKinTable[10,14],"p5M"=ThisKinTable[11,14],"p5N"=ThisKinTable[12,14],"p5P"=ThisKinTable[13,14],"p5Q"=ThisKinTable[14,14],"p5R"=ThisKinTable[15,14], + "p5S"=ThisKinTable[16,14],"p5T"=ThisKinTable[17,14],"p5V"=ThisKinTable[18,14],"p5W"=ThisKinTable[19,14],"p5Y"=ThisKinTable[20,14],"p5O"=1, + "p6A"=ThisKinTable[1,15], "p6C"=ThisKinTable[2,15], "p6D"=ThisKinTable[3,15], "p6E"=ThisKinTable[4,15], "p6F"=ThisKinTable[5,15],"p6G"=ThisKinTable[6,15],"p6H"=ThisKinTable[7,15],"p6I"=ThisKinTable[8,15],"p6K"=ThisKinTable[9,15],"p6L"=ThisKinTable[10,15],"p6M"=ThisKinTable[11,15],"p6N"=ThisKinTable[12,15],"p6P"=ThisKinTable[13,15],"p6Q"=ThisKinTable[14,15],"p6R"=ThisKinTable[15,15], + "p6S"=ThisKinTable[16,15],"p6T"=ThisKinTable[17,15],"p6V"=ThisKinTable[18,15],"p6W"=ThisKinTable[19,15],"p6Y"=ThisKinTable[20,15],"p6O"=1) + + + Abl_props <- c("m6A"=Abl[1,3], "m6C"=Abl[2,3], "m6D"=Abl[3,3], "m6E"=Abl[4,3], "m6F"=Abl[5,3],"m6G"=Abl[6,3],"m6H"=Abl[7,3],"m6I"=Abl[8,3],"m6K"=Abl[9,3],"m6L"=Abl[10,3],"m6M"=Abl[11,3],"m6N"=Abl[12,3],"m6P"=Abl[13,3],"m6Q"=Abl[14,3],"m6R"=Abl[15,3], + "m6S"=Abl[16,3],"m6T"=Abl[17,3],"m6V"=Abl[18,3],"m6W"=Abl[19,3],"m6Y"=Abl[20,3],"m6O"=1, + "m5A"=Abl[1,4], "m5C"=Abl[2,4], "m5D"=Abl[3,4], "m5E"=Abl[4,4], "m5F"=Abl[5,4],"m5G"=Abl[6,4],"m5H"=Abl[7,4],"m5I"=Abl[8,4],"m5K"=Abl[9,4],"m5L"=Abl[10,4],"m5M"=Abl[11,4],"m5N"=Abl[12,4],"m5P"=Abl[13,4],"m5Q"=Abl[14,4],"m5R"=Abl[15,4], + "m5S"=Abl[16,4],"m5T"=Abl[17,4],"m5V"=Abl[18,4],"m5W"=Abl[19,4],"m5Y"=Abl[20,4],"m5O"=1, + "m4A"=Abl[1,5], "m4C"=Abl[2,5], "m4D"=Abl[3,5], "m4E"=Abl[4,5], "m4F"=Abl[5,5],"m4G"=Abl[6,5],"m4H"=Abl[7,5],"m4I"=Abl[8,5],"m4K"=Abl[9,5],"m4L"=Abl[10,5],"m4M"=Abl[11,5],"m4N"=Abl[12,5],"m4P"=Abl[13,5],"m4Q"=Abl[14,5],"m4R"=Abl[15,5], + "m4S"=Abl[16,5],"m4T"=Abl[17,5],"m4V"=Abl[18,5],"m4W"=Abl[19,5],"m4Y"=Abl[20,5],"m4O"=1, + "m3A"=Abl[1,6], "m3C"=Abl[2,6], "m3D"=Abl[3,6], "m3E"=Abl[4,6], "m3F"=Abl[5,6],"m3G"=Abl[6,6],"m3H"=Abl[7,6],"m3I"=Abl[8,6],"m3K"=Abl[9,6],"m3L"=Abl[10,6],"m3M"=Abl[11,6],"m3N"=Abl[12,6],"m3P"=Abl[13,6],"m3Q"=Abl[14,6],"m3R"=Abl[15,6], + "m3S"=Abl[16,6],"m3T"=Abl[17,6],"m3V"=Abl[18,6],"m3W"=Abl[19,6],"m3Y"=Abl[20,6],"m3O"=1, + "m2A"=Abl[1,7], "m2C"=Abl[2,7], "m2D"=Abl[3,7], "m2E"=Abl[4,7], "m2F"=Abl[5,7],"m2G"=Abl[6,7],"m2H"=Abl[7,7],"m2I"=Abl[8,7],"m2K"=Abl[9,7],"m2L"=Abl[10,7],"m2M"=Abl[11,7],"m2N"=Abl[12,7],"m2P"=Abl[13,7],"m2Q"=Abl[14,7],"m2R"=Abl[15,7], + "m2S"=Abl[16,7],"m2T"=Abl[17,7],"m2V"=Abl[18,7],"m2W"=Abl[19,7],"m2Y"=Abl[20,7],"m2O"=1, + "m1A"=Abl[1,8], "m1C"=Abl[2,8], "m1D"=Abl[3,8], "m1E"=Abl[4,8], "m1F"=Abl[5,8],"m1G"=Abl[6,8],"m1H"=Abl[7,8],"m1I"=Abl[8,8],"m1K"=Abl[9,8],"m1L"=Abl[10,8],"m1M"=Abl[11,8],"m1N"=Abl[12,8],"m1P"=Abl[13,8],"m1Q"=Abl[14,8],"m1R"=Abl[15,8], + "m1S"=Abl[16,8],"m1T"=Abl[17,8],"m1V"=Abl[18,8],"m1W"=Abl[19,8],"m1Y"=Abl[20,8],"m1O"=1, + "d0A"=Abl[1,9], "d0C"=Abl[2,9], "d0D"=Abl[3,9], "d0E"=Abl[4,9], "d0F"=Abl[5,9],"d0G"=Abl[6,9],"d0H"=Abl[7,9],"d0I"=Abl[8,9],"d0K"=Abl[9,9],"d0L"=Abl[10,9],"d0M"=Abl[11,9],"d0N"=Abl[12,9],"d0P"=Abl[13,9],"d0Q"=Abl[14,9],"d0R"=Abl[15,9], + "d0S"=Abl[16,9],"d0T"=Abl[17,9],"d0V"=Abl[18,9],"d0W"=Abl[19,9],"d0Y"=Abl[20,9],"d0O"=1, + "p1A"=Abl[1,10], "p1C"=Abl[2,10], "p1D"=Abl[3,10], "p1E"=Abl[4,10], "p1F"=Abl[5,10],"p1G"=Abl[6,10],"p1H"=Abl[7,10],"p1I"=Abl[8,10],"p1K"=Abl[9,10],"p1L"=Abl[10,10],"p1M"=Abl[11,10],"p1N"=Abl[12,10],"p1P"=Abl[13,10],"p1Q"=Abl[14,10],"p1R"=Abl[15,10], + "p1S"=Abl[16,10],"p1T"=Abl[17,10],"p1V"=Abl[18,10],"p1W"=Abl[19,10],"p1Y"=Abl[20,10],"p1O"=1, + "p2A"=Abl[1,11], "p2C"=Abl[2,11], "p2D"=Abl[3,11], "p2E"=Abl[4,11], "p2F"=Abl[5,11],"p2G"=Abl[6,11],"p2H"=Abl[7,11],"p2I"=Abl[8,11],"p2K"=Abl[9,11],"p2L"=Abl[10,11],"p2M"=Abl[11,11],"p2N"=Abl[12,11],"p2P"=Abl[13,11],"p2Q"=Abl[14,11],"p2R"=Abl[15,11], + "p2S"=Abl[16,11],"p2T"=Abl[17,11],"p2V"=Abl[18,11],"p2W"=Abl[19,11],"p2Y"=Abl[20,11],"p2O"=1, + "p3A"=Abl[1,12], "p3C"=Abl[2,12], "p3D"=Abl[3,12], "p3E"=Abl[4,12], "p3F"=Abl[5,12],"p3G"=Abl[6,12],"p3H"=Abl[7,12],"p3I"=Abl[8,12],"p3K"=Abl[9,12],"p3L"=Abl[10,12],"p3M"=Abl[11,12],"p3N"=Abl[12,12],"p3P"=Abl[13,12],"p3Q"=Abl[14,12],"p3R"=Abl[15,12], + "p3S"=Abl[16,12],"p3T"=Abl[17,12],"p3V"=Abl[18,12],"p3W"=Abl[19,12],"p3Y"=Abl[20,12],"p3O"=1, + "p4A"=Abl[1,13], "p4C"=Abl[2,13], "p4D"=Abl[3,13], "p4E"=Abl[4,13], "p4F"=Abl[5,13],"p4G"=Abl[6,13],"p4H"=Abl[7,13],"p4I"=Abl[8,13],"p4K"=Abl[9,13],"p4L"=Abl[10,13],"p4M"=Abl[11,13],"p4N"=Abl[12,13],"p4P"=Abl[13,13],"p4Q"=Abl[14,13],"p4R"=Abl[15,13], + "p4S"=Abl[16,13],"p4T"=Abl[17,13],"p4V"=Abl[18,13],"p4W"=Abl[19,13],"p4Y"=Abl[20,13],"p4O"=1, + "p5A"=Abl[1,14], "p5C"=Abl[2,14], "p5D"=Abl[3,14], "p5E"=Abl[4,14], "p5F"=Abl[5,14],"p5G"=Abl[6,14],"p5H"=Abl[7,14],"p5I"=Abl[8,14],"p5K"=Abl[9,14],"p5L"=Abl[10,14],"p5M"=Abl[11,14],"p5N"=Abl[12,14],"p5P"=Abl[13,14],"p5Q"=Abl[14,14],"p5R"=Abl[15,14], + "p5S"=Abl[16,14],"p5T"=Abl[17,14],"p5V"=Abl[18,14],"p5W"=Abl[19,14],"p5Y"=Abl[20,14],"p5O"=1, + "p6A"=Abl[1,15], "p6C"=Abl[2,15], "p6D"=Abl[3,15], "p6E"=Abl[4,15], "p6F"=Abl[5,15],"p6G"=Abl[6,15],"p6H"=Abl[7,15],"p6I"=Abl[8,15],"p6K"=Abl[9,15],"p6L"=Abl[10,15],"p6M"=Abl[11,15],"p6N"=Abl[12,15],"p6P"=Abl[13,15],"p6Q"=Abl[14,15],"p6R"=Abl[15,15], + "p6S"=Abl[16,15],"p6T"=Abl[17,15],"p6V"=Abl[18,15],"p6W"=Abl[19,15],"p6Y"=Abl[20,15],"p6O"=1) + + Arg_props <- c("m6A"=Arg[1,3], "m6C"=Arg[2,3], "m6D"=Arg[3,3], "m6E"=Arg[4,3], "m6F"=Arg[5,3],"m6G"=Arg[6,3],"m6H"=Arg[7,3],"m6I"=Arg[8,3],"m6K"=Arg[9,3],"m6L"=Arg[10,3],"m6M"=Arg[11,3],"m6N"=Arg[12,3],"m6P"=Arg[13,3],"m6Q"=Arg[14,3],"m6R"=Arg[15,3], + "m6S"=Arg[16,3],"m6T"=Arg[17,3],"m6V"=Arg[18,3],"m6W"=Arg[19,3],"m6Y"=Arg[20,3],"m6O"=1, + "m5A"=Arg[1,4], "m5C"=Arg[2,4], "m5D"=Arg[3,4], "m5E"=Arg[4,4], "m5F"=Arg[5,4],"m5G"=Arg[6,4],"m5H"=Arg[7,4],"m5I"=Arg[8,4],"m5K"=Arg[9,4],"m5L"=Arg[10,4],"m5M"=Arg[11,4],"m5N"=Arg[12,4],"m5P"=Arg[13,4],"m5Q"=Arg[14,4],"m5R"=Arg[15,4], + "m5S"=Arg[16,4],"m5T"=Arg[17,4],"m5V"=Arg[18,4],"m5W"=Arg[19,4],"m5Y"=Arg[20,4],"m5O"=1, + "m4A"=Arg[1,5], "m4C"=Arg[2,5], "m4D"=Arg[3,5], "m4E"=Arg[4,5], "m4F"=Arg[5,5],"m4G"=Arg[6,5],"m4H"=Arg[7,5],"m4I"=Arg[8,5],"m4K"=Arg[9,5],"m4L"=Arg[10,5],"m4M"=Arg[11,5],"m4N"=Arg[12,5],"m4P"=Arg[13,5],"m4Q"=Arg[14,5],"m4R"=Arg[15,5], + "m4S"=Arg[16,5],"m4T"=Arg[17,5],"m4V"=Arg[18,5],"m4W"=Arg[19,5],"m4Y"=Arg[20,5],"m4O"=1, + "m3A"=Arg[1,6], "m3C"=Arg[2,6], "m3D"=Arg[3,6], "m3E"=Arg[4,6], "m3F"=Arg[5,6],"m3G"=Arg[6,6],"m3H"=Arg[7,6],"m3I"=Arg[8,6],"m3K"=Arg[9,6],"m3L"=Arg[10,6],"m3M"=Arg[11,6],"m3N"=Arg[12,6],"m3P"=Arg[13,6],"m3Q"=Arg[14,6],"m3R"=Arg[15,6], + "m3S"=Arg[16,6],"m3T"=Arg[17,6],"m3V"=Arg[18,6],"m3W"=Arg[19,6],"m3Y"=Arg[20,6],"m3O"=1, + "m2A"=Arg[1,7], "m2C"=Arg[2,7], "m2D"=Arg[3,7], "m2E"=Arg[4,7], "m2F"=Arg[5,7],"m2G"=Arg[6,7],"m2H"=Arg[7,7],"m2I"=Arg[8,7],"m2K"=Arg[9,7],"m2L"=Arg[10,7],"m2M"=Arg[11,7],"m2N"=Arg[12,7],"m2P"=Arg[13,7],"m2Q"=Arg[14,7],"m2R"=Arg[15,7], + "m2S"=Arg[16,7],"m2T"=Arg[17,7],"m2V"=Arg[18,7],"m2W"=Arg[19,7],"m2Y"=Arg[20,7],"m2O"=1, + "m1A"=Arg[1,8], "m1C"=Arg[2,8], "m1D"=Arg[3,8], "m1E"=Arg[4,8], "m1F"=Arg[5,8],"m1G"=Arg[6,8],"m1H"=Arg[7,8],"m1I"=Arg[8,8],"m1K"=Arg[9,8],"m1L"=Arg[10,8],"m1M"=Arg[11,8],"m1N"=Arg[12,8],"m1P"=Arg[13,8],"m1Q"=Arg[14,8],"m1R"=Arg[15,8], + "m1S"=Arg[16,8],"m1T"=Arg[17,8],"m1V"=Arg[18,8],"m1W"=Arg[19,8],"m1Y"=Arg[20,8],"m1O"=1, + "d0A"=Arg[1,9], "d0C"=Arg[2,9], "d0D"=Arg[3,9], "d0E"=Arg[4,9], "d0F"=Arg[5,9],"d0G"=Arg[6,9],"d0H"=Arg[7,9],"d0I"=Arg[8,9],"d0K"=Arg[9,9],"d0L"=Arg[10,9],"d0M"=Arg[11,9],"d0N"=Arg[12,9],"d0P"=Arg[13,9],"d0Q"=Arg[14,9],"d0R"=Arg[15,9], + "d0S"=Arg[16,9],"d0T"=Arg[17,9],"d0V"=Arg[18,9],"d0W"=Arg[19,9],"d0Y"=Arg[20,9],"d0O"=1, + "p1A"=Arg[1,10], "p1C"=Arg[2,10], "p1D"=Arg[3,10], "p1E"=Arg[4,10], "p1F"=Arg[5,10],"p1G"=Arg[6,10],"p1H"=Arg[7,10],"p1I"=Arg[8,10],"p1K"=Arg[9,10],"p1L"=Arg[10,10],"p1M"=Arg[11,10],"p1N"=Arg[12,10],"p1P"=Arg[13,10],"p1Q"=Arg[14,10],"p1R"=Arg[15,10], + "p1S"=Arg[16,10],"p1T"=Arg[17,10],"p1V"=Arg[18,10],"p1W"=Arg[19,10],"p1Y"=Arg[20,10],"p1O"=1, + "p2A"=Arg[1,11], "p2C"=Arg[2,11], "p2D"=Arg[3,11], "p2E"=Arg[4,11], "p2F"=Arg[5,11],"p2G"=Arg[6,11],"p2H"=Arg[7,11],"p2I"=Arg[8,11],"p2K"=Arg[9,11],"p2L"=Arg[10,11],"p2M"=Arg[11,11],"p2N"=Arg[12,11],"p2P"=Arg[13,11],"p2Q"=Arg[14,11],"p2R"=Arg[15,11], + "p2S"=Arg[16,11],"p2T"=Arg[17,11],"p2V"=Arg[18,11],"p2W"=Arg[19,11],"p2Y"=Arg[20,11],"p2O"=1, + "p3A"=Arg[1,12], "p3C"=Arg[2,12], "p3D"=Arg[3,12], "p3E"=Arg[4,12], "p3F"=Arg[5,12],"p3G"=Arg[6,12],"p3H"=Arg[7,12],"p3I"=Arg[8,12],"p3K"=Arg[9,12],"p3L"=Arg[10,12],"p3M"=Arg[11,12],"p3N"=Arg[12,12],"p3P"=Arg[13,12],"p3Q"=Arg[14,12],"p3R"=Arg[15,12], + "p3S"=Arg[16,12],"p3T"=Arg[17,12],"p3V"=Arg[18,12],"p3W"=Arg[19,12],"p3Y"=Arg[20,12],"p3O"=1, + "p4A"=Arg[1,13], "p4C"=Arg[2,13], "p4D"=Arg[3,13], "p4E"=Arg[4,13], "p4F"=Arg[5,13],"p4G"=Arg[6,13],"p4H"=Arg[7,13],"p4I"=Arg[8,13],"p4K"=Arg[9,13],"p4L"=Arg[10,13],"p4M"=Arg[11,13],"p4N"=Arg[12,13],"p4P"=Arg[13,13],"p4Q"=Arg[14,13],"p4R"=Arg[15,13], + "p4S"=Arg[16,13],"p4T"=Arg[17,13],"p4V"=Arg[18,13],"p4W"=Arg[19,13],"p4Y"=Arg[20,13],"p4O"=1, + "p5A"=Arg[1,14], "p5C"=Arg[2,14], "p5D"=Arg[3,14], "p5E"=Arg[4,14], "p5F"=Arg[5,14],"p5G"=Arg[6,14],"p5H"=Arg[7,14],"p5I"=Arg[8,14],"p5K"=Arg[9,14],"p5L"=Arg[10,14],"p5M"=Arg[11,14],"p5N"=Arg[12,14],"p5P"=Arg[13,14],"p5Q"=Arg[14,14],"p5R"=Arg[15,14], + "p5S"=Arg[16,14],"p5T"=Arg[17,14],"p5V"=Arg[18,14],"p5W"=Arg[19,14],"p5Y"=Arg[20,14],"p5O"=1, + "p6A"=Arg[1,15], "p6C"=Arg[2,15], "p6D"=Arg[3,15], "p6E"=Arg[4,15], "p6F"=Arg[5,15],"p6G"=Arg[6,15],"p6H"=Arg[7,15],"p6I"=Arg[8,15],"p6K"=Arg[9,15],"p6L"=Arg[10,15],"p6M"=Arg[11,15],"p6N"=Arg[12,15],"p6P"=Arg[13,15],"p6Q"=Arg[14,15],"p6R"=Arg[15,15], + "p6S"=Arg[16,15],"p6T"=Arg[17,15],"p6V"=Arg[18,15],"p6W"=Arg[19,15],"p6Y"=Arg[20,15],"p6O"=1) + + Btk_props <- c("m6A"=Btk[1,3], "m6C"=Btk[2,3], "m6D"=Btk[3,3], "m6E"=Btk[4,3], "m6F"=Btk[5,3],"m6G"=Btk[6,3],"m6H"=Btk[7,3],"m6I"=Btk[8,3],"m6K"=Btk[9,3],"m6L"=Btk[10,3],"m6M"=Btk[11,3],"m6N"=Btk[12,3],"m6P"=Btk[13,3],"m6Q"=Btk[14,3],"m6R"=Btk[15,3], + "m6S"=Btk[16,3],"m6T"=Btk[17,3],"m6V"=Btk[18,3],"m6W"=Btk[19,3],"m6Y"=Btk[20,3],"m6O"=1, + "m5A"=Btk[1,4], "m5C"=Btk[2,4], "m5D"=Btk[3,4], "m5E"=Btk[4,4], "m5F"=Btk[5,4],"m5G"=Btk[6,4],"m5H"=Btk[7,4],"m5I"=Btk[8,4],"m5K"=Btk[9,4],"m5L"=Btk[10,4],"m5M"=Btk[11,4],"m5N"=Btk[12,4],"m5P"=Btk[13,4],"m5Q"=Btk[14,4],"m5R"=Btk[15,4], + "m5S"=Btk[16,4],"m5T"=Btk[17,4],"m5V"=Btk[18,4],"m5W"=Btk[19,4],"m5Y"=Btk[20,4],"m5O"=1, + "m4A"=Btk[1,5], "m4C"=Btk[2,5], "m4D"=Btk[3,5], "m4E"=Btk[4,5], "m4F"=Btk[5,5],"m4G"=Btk[6,5],"m4H"=Btk[7,5],"m4I"=Btk[8,5],"m4K"=Btk[9,5],"m4L"=Btk[10,5],"m4M"=Btk[11,5],"m4N"=Btk[12,5],"m4P"=Btk[13,5],"m4Q"=Btk[14,5],"m4R"=Btk[15,5], + "m4S"=Btk[16,5],"m4T"=Btk[17,5],"m4V"=Btk[18,5],"m4W"=Btk[19,5],"m4Y"=Btk[20,5],"m4O"=1, + "m3A"=Btk[1,6], "m3C"=Btk[2,6], "m3D"=Btk[3,6], "m3E"=Btk[4,6], "m3F"=Btk[5,6],"m3G"=Btk[6,6],"m3H"=Btk[7,6],"m3I"=Btk[8,6],"m3K"=Btk[9,6],"m3L"=Btk[10,6],"m3M"=Btk[11,6],"m3N"=Btk[12,6],"m3P"=Btk[13,6],"m3Q"=Btk[14,6],"m3R"=Btk[15,6], + "m3S"=Btk[16,6],"m3T"=Btk[17,6],"m3V"=Btk[18,6],"m3W"=Btk[19,6],"m3Y"=Btk[20,6],"m3O"=1, + "m2A"=Btk[1,7], "m2C"=Btk[2,7], "m2D"=Btk[3,7], "m2E"=Btk[4,7], "m2F"=Btk[5,7],"m2G"=Btk[6,7],"m2H"=Btk[7,7],"m2I"=Btk[8,7],"m2K"=Btk[9,7],"m2L"=Btk[10,7],"m2M"=Btk[11,7],"m2N"=Btk[12,7],"m2P"=Btk[13,7],"m2Q"=Btk[14,7],"m2R"=Btk[15,7], + "m2S"=Btk[16,7],"m2T"=Btk[17,7],"m2V"=Btk[18,7],"m2W"=Btk[19,7],"m2Y"=Btk[20,7],"m2O"=1, + "m1A"=Btk[1,8], "m1C"=Btk[2,8], "m1D"=Btk[3,8], "m1E"=Btk[4,8], "m1F"=Btk[5,8],"m1G"=Btk[6,8],"m1H"=Btk[7,8],"m1I"=Btk[8,8],"m1K"=Btk[9,8],"m1L"=Btk[10,8],"m1M"=Btk[11,8],"m1N"=Btk[12,8],"m1P"=Btk[13,8],"m1Q"=Btk[14,8],"m1R"=Btk[15,8], + "m1S"=Btk[16,8],"m1T"=Btk[17,8],"m1V"=Btk[18,8],"m1W"=Btk[19,8],"m1Y"=Btk[20,8],"m1O"=1, + "d0A"=Btk[1,9], "d0C"=Btk[2,9], "d0D"=Btk[3,9], "d0E"=Btk[4,9], "d0F"=Btk[5,9],"d0G"=Btk[6,9],"d0H"=Btk[7,9],"d0I"=Btk[8,9],"d0K"=Btk[9,9],"d0L"=Btk[10,9],"d0M"=Btk[11,9],"d0N"=Btk[12,9],"d0P"=Btk[13,9],"d0Q"=Btk[14,9],"d0R"=Btk[15,9], + "d0S"=Btk[16,9],"d0T"=Btk[17,9],"d0V"=Btk[18,9],"d0W"=Btk[19,9],"d0Y"=Btk[20,9],"d0O"=1, + "p1A"=Btk[1,10], "p1C"=Btk[2,10], "p1D"=Btk[3,10], "p1E"=Btk[4,10], "p1F"=Btk[5,10],"p1G"=Btk[6,10],"p1H"=Btk[7,10],"p1I"=Btk[8,10],"p1K"=Btk[9,10],"p1L"=Btk[10,10],"p1M"=Btk[11,10],"p1N"=Btk[12,10],"p1P"=Btk[13,10],"p1Q"=Btk[14,10],"p1R"=Btk[15,10], + "p1S"=Btk[16,10],"p1T"=Btk[17,10],"p1V"=Btk[18,10],"p1W"=Btk[19,10],"p1Y"=Btk[20,10],"p1O"=1, + "p2A"=Btk[1,11], "p2C"=Btk[2,11], "p2D"=Btk[3,11], "p2E"=Btk[4,11], "p2F"=Btk[5,11],"p2G"=Btk[6,11],"p2H"=Btk[7,11],"p2I"=Btk[8,11],"p2K"=Btk[9,11],"p2L"=Btk[10,11],"p2M"=Btk[11,11],"p2N"=Btk[12,11],"p2P"=Btk[13,11],"p2Q"=Btk[14,11],"p2R"=Btk[15,11], + "p2S"=Btk[16,11],"p2T"=Btk[17,11],"p2V"=Btk[18,11],"p2W"=Btk[19,11],"p2Y"=Btk[20,11],"p2O"=1, + "p3A"=Btk[1,12], "p3C"=Btk[2,12], "p3D"=Btk[3,12], "p3E"=Btk[4,12], "p3F"=Btk[5,12],"p3G"=Btk[6,12],"p3H"=Btk[7,12],"p3I"=Btk[8,12],"p3K"=Btk[9,12],"p3L"=Btk[10,12],"p3M"=Btk[11,12],"p3N"=Btk[12,12],"p3P"=Btk[13,12],"p3Q"=Btk[14,12],"p3R"=Btk[15,12], + "p3S"=Btk[16,12],"p3T"=Btk[17,12],"p3V"=Btk[18,12],"p3W"=Btk[19,12],"p3Y"=Btk[20,12],"p3O"=1, + "p4A"=Btk[1,13], "p4C"=Btk[2,13], "p4D"=Btk[3,13], "p4E"=Btk[4,13], "p4F"=Btk[5,13],"p4G"=Btk[6,13],"p4H"=Btk[7,13],"p4I"=Btk[8,13],"p4K"=Btk[9,13],"p4L"=Btk[10,13],"p4M"=Btk[11,13],"p4N"=Btk[12,13],"p4P"=Btk[13,13],"p4Q"=Btk[14,13],"p4R"=Btk[15,13], + "p4S"=Btk[16,13],"p4T"=Btk[17,13],"p4V"=Btk[18,13],"p4W"=Btk[19,13],"p4Y"=Btk[20,13],"p4O"=1, + "p5A"=Btk[1,14], "p5C"=Btk[2,14], "p5D"=Btk[3,14], "p5E"=Btk[4,14], "p5F"=Btk[5,14],"p5G"=Btk[6,14],"p5H"=Btk[7,14],"p5I"=Btk[8,14],"p5K"=Btk[9,14],"p5L"=Btk[10,14],"p5M"=Btk[11,14],"p5N"=Btk[12,14],"p5P"=Btk[13,14],"p5Q"=Btk[14,14],"p5R"=Btk[15,14], + "p5S"=Btk[16,14],"p5T"=Btk[17,14],"p5V"=Btk[18,14],"p5W"=Btk[19,14],"p5Y"=Btk[20,14],"p5O"=1, + "p6A"=Btk[1,15], "p6C"=Btk[2,15], "p6D"=Btk[3,15], "p6E"=Btk[4,15], "p6F"=Btk[5,15],"p6G"=Btk[6,15],"p6H"=Btk[7,15],"p6I"=Btk[8,15],"p6K"=Btk[9,15],"p6L"=Btk[10,15],"p6M"=Btk[11,15],"p6N"=Btk[12,15],"p6P"=Btk[13,15],"p6Q"=Btk[14,15],"p6R"=Btk[15,15], + "p6S"=Btk[16,15],"p6T"=Btk[17,15],"p6V"=Btk[18,15],"p6W"=Btk[19,15],"p6Y"=Btk[20,15],"p6O"=1) + + Csk_props <- c("m6A"=Csk[1,3], "m6C"=Csk[2,3], "m6D"=Csk[3,3], "m6E"=Csk[4,3], "m6F"=Csk[5,3],"m6G"=Csk[6,3],"m6H"=Csk[7,3],"m6I"=Csk[8,3],"m6K"=Csk[9,3],"m6L"=Csk[10,3],"m6M"=Csk[11,3],"m6N"=Csk[12,3],"m6P"=Csk[13,3],"m6Q"=Csk[14,3],"m6R"=Csk[15,3], + "m6S"=Csk[16,3],"m6T"=Csk[17,3],"m6V"=Csk[18,3],"m6W"=Csk[19,3],"m6Y"=Csk[20,3],"m6O"=1, + "m5A"=Csk[1,4], "m5C"=Csk[2,4], "m5D"=Csk[3,4], "m5E"=Csk[4,4], "m5F"=Csk[5,4],"m5G"=Csk[6,4],"m5H"=Csk[7,4],"m5I"=Csk[8,4],"m5K"=Csk[9,4],"m5L"=Csk[10,4],"m5M"=Csk[11,4],"m5N"=Csk[12,4],"m5P"=Csk[13,4],"m5Q"=Csk[14,4],"m5R"=Csk[15,4], + "m5S"=Csk[16,4],"m5T"=Csk[17,4],"m5V"=Csk[18,4],"m5W"=Csk[19,4],"m5Y"=Csk[20,4],"m5O"=1, + "m4A"=Csk[1,5], "m4C"=Csk[2,5], "m4D"=Csk[3,5], "m4E"=Csk[4,5], "m4F"=Csk[5,5],"m4G"=Csk[6,5],"m4H"=Csk[7,5],"m4I"=Csk[8,5],"m4K"=Csk[9,5],"m4L"=Csk[10,5],"m4M"=Csk[11,5],"m4N"=Csk[12,5],"m4P"=Csk[13,5],"m4Q"=Csk[14,5],"m4R"=Csk[15,5], + "m4S"=Csk[16,5],"m4T"=Csk[17,5],"m4V"=Csk[18,5],"m4W"=Csk[19,5],"m4Y"=Csk[20,5],"m4O"=1, + "m3A"=Csk[1,6], "m3C"=Csk[2,6], "m3D"=Csk[3,6], "m3E"=Csk[4,6], "m3F"=Csk[5,6],"m3G"=Csk[6,6],"m3H"=Csk[7,6],"m3I"=Csk[8,6],"m3K"=Csk[9,6],"m3L"=Csk[10,6],"m3M"=Csk[11,6],"m3N"=Csk[12,6],"m3P"=Csk[13,6],"m3Q"=Csk[14,6],"m3R"=Csk[15,6], + "m3S"=Csk[16,6],"m3T"=Csk[17,6],"m3V"=Csk[18,6],"m3W"=Csk[19,6],"m3Y"=Csk[20,6],"m3O"=1, + "m2A"=Csk[1,7], "m2C"=Csk[2,7], "m2D"=Csk[3,7], "m2E"=Csk[4,7], "m2F"=Csk[5,7],"m2G"=Csk[6,7],"m2H"=Csk[7,7],"m2I"=Csk[8,7],"m2K"=Csk[9,7],"m2L"=Csk[10,7],"m2M"=Csk[11,7],"m2N"=Csk[12,7],"m2P"=Csk[13,7],"m2Q"=Csk[14,7],"m2R"=Csk[15,7], + "m2S"=Csk[16,7],"m2T"=Csk[17,7],"m2V"=Csk[18,7],"m2W"=Csk[19,7],"m2Y"=Csk[20,7],"m2O"=1, + "m1A"=Csk[1,8], "m1C"=Csk[2,8], "m1D"=Csk[3,8], "m1E"=Csk[4,8], "m1F"=Csk[5,8],"m1G"=Csk[6,8],"m1H"=Csk[7,8],"m1I"=Csk[8,8],"m1K"=Csk[9,8],"m1L"=Csk[10,8],"m1M"=Csk[11,8],"m1N"=Csk[12,8],"m1P"=Csk[13,8],"m1Q"=Csk[14,8],"m1R"=Csk[15,8], + "m1S"=Csk[16,8],"m1T"=Csk[17,8],"m1V"=Csk[18,8],"m1W"=Csk[19,8],"m1Y"=Csk[20,8],"m1O"=1, + "d0A"=Csk[1,9], "d0C"=Csk[2,9], "d0D"=Csk[3,9], "d0E"=Csk[4,9], "d0F"=Csk[5,9],"d0G"=Csk[6,9],"d0H"=Csk[7,9],"d0I"=Csk[8,9],"d0K"=Csk[9,9],"d0L"=Csk[10,9],"d0M"=Csk[11,9],"d0N"=Csk[12,9],"d0P"=Csk[13,9],"d0Q"=Csk[14,9],"d0R"=Csk[15,9], + "d0S"=Csk[16,9],"d0T"=Csk[17,9],"d0V"=Csk[18,9],"d0W"=Csk[19,9],"d0Y"=Csk[20,9],"d0O"=1, + "p1A"=Csk[1,10], "p1C"=Csk[2,10], "p1D"=Csk[3,10], "p1E"=Csk[4,10], "p1F"=Csk[5,10],"p1G"=Csk[6,10],"p1H"=Csk[7,10],"p1I"=Csk[8,10],"p1K"=Csk[9,10],"p1L"=Csk[10,10],"p1M"=Csk[11,10],"p1N"=Csk[12,10],"p1P"=Csk[13,10],"p1Q"=Csk[14,10],"p1R"=Csk[15,10], + "p1S"=Csk[16,10],"p1T"=Csk[17,10],"p1V"=Csk[18,10],"p1W"=Csk[19,10],"p1Y"=Csk[20,10],"p1O"=1, + "p2A"=Csk[1,11], "p2C"=Csk[2,11], "p2D"=Csk[3,11], "p2E"=Csk[4,11], "p2F"=Csk[5,11],"p2G"=Csk[6,11],"p2H"=Csk[7,11],"p2I"=Csk[8,11],"p2K"=Csk[9,11],"p2L"=Csk[10,11],"p2M"=Csk[11,11],"p2N"=Csk[12,11],"p2P"=Csk[13,11],"p2Q"=Csk[14,11],"p2R"=Csk[15,11], + "p2S"=Csk[16,11],"p2T"=Csk[17,11],"p2V"=Csk[18,11],"p2W"=Csk[19,11],"p2Y"=Csk[20,11],"p2O"=1, + "p3A"=Csk[1,12], "p3C"=Csk[2,12], "p3D"=Csk[3,12], "p3E"=Csk[4,12], "p3F"=Csk[5,12],"p3G"=Csk[6,12],"p3H"=Csk[7,12],"p3I"=Csk[8,12],"p3K"=Csk[9,12],"p3L"=Csk[10,12],"p3M"=Csk[11,12],"p3N"=Csk[12,12],"p3P"=Csk[13,12],"p3Q"=Csk[14,12],"p3R"=Csk[15,12], + "p3S"=Csk[16,12],"p3T"=Csk[17,12],"p3V"=Csk[18,12],"p3W"=Csk[19,12],"p3Y"=Csk[20,12],"p3O"=1, + "p4A"=Csk[1,13], "p4C"=Csk[2,13], "p4D"=Csk[3,13], "p4E"=Csk[4,13], "p4F"=Csk[5,13],"p4G"=Csk[6,13],"p4H"=Csk[7,13],"p4I"=Csk[8,13],"p4K"=Csk[9,13],"p4L"=Csk[10,13],"p4M"=Csk[11,13],"p4N"=Csk[12,13],"p4P"=Csk[13,13],"p4Q"=Csk[14,13],"p4R"=Csk[15,13], + "p4S"=Csk[16,13],"p4T"=Csk[17,13],"p4V"=Csk[18,13],"p4W"=Csk[19,13],"p4Y"=Csk[20,13],"p4O"=1, + "p5A"=Csk[1,14], "p5C"=Csk[2,14], "p5D"=Csk[3,14], "p5E"=Csk[4,14], "p5F"=Csk[5,14],"p5G"=Csk[6,14],"p5H"=Csk[7,14],"p5I"=Csk[8,14],"p5K"=Csk[9,14],"p5L"=Csk[10,14],"p5M"=Csk[11,14],"p5N"=Csk[12,14],"p5P"=Csk[13,14],"p5Q"=Csk[14,14],"p5R"=Csk[15,14], + "p5S"=Csk[16,14],"p5T"=Csk[17,14],"p5V"=Csk[18,14],"p5W"=Csk[19,14],"p5Y"=Csk[20,14],"p5O"=1, + "p6A"=Csk[1,15], "p6C"=Csk[2,15], "p6D"=Csk[3,15], "p6E"=Csk[4,15], "p6F"=Csk[5,15],"p6G"=Csk[6,15],"p6H"=Csk[7,15],"p6I"=Csk[8,15],"p6K"=Csk[9,15],"p6L"=Csk[10,15],"p6M"=Csk[11,15],"p6N"=Csk[12,15],"p6P"=Csk[13,15],"p6Q"=Csk[14,15],"p6R"=Csk[15,15], + "p6S"=Csk[16,15],"p6T"=Csk[17,15],"p6V"=Csk[18,15],"p6W"=Csk[19,15],"p6Y"=Csk[20,15],"p6O"=1) + + Fyn_props <- c("m6A"=Fyn[1,3], "m6C"=Fyn[2,3], "m6D"=Fyn[3,3], "m6E"=Fyn[4,3], "m6F"=Fyn[5,3],"m6G"=Fyn[6,3],"m6H"=Fyn[7,3],"m6I"=Fyn[8,3],"m6K"=Fyn[9,3],"m6L"=Fyn[10,3],"m6M"=Fyn[11,3],"m6N"=Fyn[12,3],"m6P"=Fyn[13,3],"m6Q"=Fyn[14,3],"m6R"=Fyn[15,3], + "m6S"=Fyn[16,3],"m6T"=Fyn[17,3],"m6V"=Fyn[18,3],"m6W"=Fyn[19,3],"m6Y"=Fyn[20,3],"m6O"=1, + "m5A"=Fyn[1,4], "m5C"=Fyn[2,4], "m5D"=Fyn[3,4], "m5E"=Fyn[4,4], "m5F"=Fyn[5,4],"m5G"=Fyn[6,4],"m5H"=Fyn[7,4],"m5I"=Fyn[8,4],"m5K"=Fyn[9,4],"m5L"=Fyn[10,4],"m5M"=Fyn[11,4],"m5N"=Fyn[12,4],"m5P"=Fyn[13,4],"m5Q"=Fyn[14,4],"m5R"=Fyn[15,4], + "m5S"=Fyn[16,4],"m5T"=Fyn[17,4],"m5V"=Fyn[18,4],"m5W"=Fyn[19,4],"m5Y"=Fyn[20,4],"m5O"=1, + "m4A"=Fyn[1,5], "m4C"=Fyn[2,5], "m4D"=Fyn[3,5], "m4E"=Fyn[4,5], "m4F"=Fyn[5,5],"m4G"=Fyn[6,5],"m4H"=Fyn[7,5],"m4I"=Fyn[8,5],"m4K"=Fyn[9,5],"m4L"=Fyn[10,5],"m4M"=Fyn[11,5],"m4N"=Fyn[12,5],"m4P"=Fyn[13,5],"m4Q"=Fyn[14,5],"m4R"=Fyn[15,5], + "m4S"=Fyn[16,5],"m4T"=Fyn[17,5],"m4V"=Fyn[18,5],"m4W"=Fyn[19,5],"m4Y"=Fyn[20,5],"m4O"=1, + "m3A"=Fyn[1,6], "m3C"=Fyn[2,6], "m3D"=Fyn[3,6], "m3E"=Fyn[4,6], "m3F"=Fyn[5,6],"m3G"=Fyn[6,6],"m3H"=Fyn[7,6],"m3I"=Fyn[8,6],"m3K"=Fyn[9,6],"m3L"=Fyn[10,6],"m3M"=Fyn[11,6],"m3N"=Fyn[12,6],"m3P"=Fyn[13,6],"m3Q"=Fyn[14,6],"m3R"=Fyn[15,6], + "m3S"=Fyn[16,6],"m3T"=Fyn[17,6],"m3V"=Fyn[18,6],"m3W"=Fyn[19,6],"m3Y"=Fyn[20,6],"m3O"=1, + "m2A"=Fyn[1,7], "m2C"=Fyn[2,7], "m2D"=Fyn[3,7], "m2E"=Fyn[4,7], "m2F"=Fyn[5,7],"m2G"=Fyn[6,7],"m2H"=Fyn[7,7],"m2I"=Fyn[8,7],"m2K"=Fyn[9,7],"m2L"=Fyn[10,7],"m2M"=Fyn[11,7],"m2N"=Fyn[12,7],"m2P"=Fyn[13,7],"m2Q"=Fyn[14,7],"m2R"=Fyn[15,7], + "m2S"=Fyn[16,7],"m2T"=Fyn[17,7],"m2V"=Fyn[18,7],"m2W"=Fyn[19,7],"m2Y"=Fyn[20,7],"m2O"=1, + "m1A"=Fyn[1,8], "m1C"=Fyn[2,8], "m1D"=Fyn[3,8], "m1E"=Fyn[4,8], "m1F"=Fyn[5,8],"m1G"=Fyn[6,8],"m1H"=Fyn[7,8],"m1I"=Fyn[8,8],"m1K"=Fyn[9,8],"m1L"=Fyn[10,8],"m1M"=Fyn[11,8],"m1N"=Fyn[12,8],"m1P"=Fyn[13,8],"m1Q"=Fyn[14,8],"m1R"=Fyn[15,8], + "m1S"=Fyn[16,8],"m1T"=Fyn[17,8],"m1V"=Fyn[18,8],"m1W"=Fyn[19,8],"m1Y"=Fyn[20,8],"m1O"=1, + "d0A"=Fyn[1,9], "d0C"=Fyn[2,9], "d0D"=Fyn[3,9], "d0E"=Fyn[4,9], "d0F"=Fyn[5,9],"d0G"=Fyn[6,9],"d0H"=Fyn[7,9],"d0I"=Fyn[8,9],"d0K"=Fyn[9,9],"d0L"=Fyn[10,9],"d0M"=Fyn[11,9],"d0N"=Fyn[12,9],"d0P"=Fyn[13,9],"d0Q"=Fyn[14,9],"d0R"=Fyn[15,9], + "d0S"=Fyn[16,9],"d0T"=Fyn[17,9],"d0V"=Fyn[18,9],"d0W"=Fyn[19,9],"d0Y"=Fyn[20,9],"d0O"=1, + "p1A"=Fyn[1,10], "p1C"=Fyn[2,10], "p1D"=Fyn[3,10], "p1E"=Fyn[4,10], "p1F"=Fyn[5,10],"p1G"=Fyn[6,10],"p1H"=Fyn[7,10],"p1I"=Fyn[8,10],"p1K"=Fyn[9,10],"p1L"=Fyn[10,10],"p1M"=Fyn[11,10],"p1N"=Fyn[12,10],"p1P"=Fyn[13,10],"p1Q"=Fyn[14,10],"p1R"=Fyn[15,10], + "p1S"=Fyn[16,10],"p1T"=Fyn[17,10],"p1V"=Fyn[18,10],"p1W"=Fyn[19,10],"p1Y"=Fyn[20,10],"p1O"=1, + "p2A"=Fyn[1,11], "p2C"=Fyn[2,11], "p2D"=Fyn[3,11], "p2E"=Fyn[4,11], "p2F"=Fyn[5,11],"p2G"=Fyn[6,11],"p2H"=Fyn[7,11],"p2I"=Fyn[8,11],"p2K"=Fyn[9,11],"p2L"=Fyn[10,11],"p2M"=Fyn[11,11],"p2N"=Fyn[12,11],"p2P"=Fyn[13,11],"p2Q"=Fyn[14,11],"p2R"=Fyn[15,11], + "p2S"=Fyn[16,11],"p2T"=Fyn[17,11],"p2V"=Fyn[18,11],"p2W"=Fyn[19,11],"p2Y"=Fyn[20,11],"p2O"=1, + "p3A"=Fyn[1,12], "p3C"=Fyn[2,12], "p3D"=Fyn[3,12], "p3E"=Fyn[4,12], "p3F"=Fyn[5,12],"p3G"=Fyn[6,12],"p3H"=Fyn[7,12],"p3I"=Fyn[8,12],"p3K"=Fyn[9,12],"p3L"=Fyn[10,12],"p3M"=Fyn[11,12],"p3N"=Fyn[12,12],"p3P"=Fyn[13,12],"p3Q"=Fyn[14,12],"p3R"=Fyn[15,12], + "p3S"=Fyn[16,12],"p3T"=Fyn[17,12],"p3V"=Fyn[18,12],"p3W"=Fyn[19,12],"p3Y"=Fyn[20,12],"p3O"=1, + "p4A"=Fyn[1,13], "p4C"=Fyn[2,13], "p4D"=Fyn[3,13], "p4E"=Fyn[4,13], "p4F"=Fyn[5,13],"p4G"=Fyn[6,13],"p4H"=Fyn[7,13],"p4I"=Fyn[8,13],"p4K"=Fyn[9,13],"p4L"=Fyn[10,13],"p4M"=Fyn[11,13],"p4N"=Fyn[12,13],"p4P"=Fyn[13,13],"p4Q"=Fyn[14,13],"p4R"=Fyn[15,13], + "p4S"=Fyn[16,13],"p4T"=Fyn[17,13],"p4V"=Fyn[18,13],"p4W"=Fyn[19,13],"p4Y"=Fyn[20,13],"p4O"=1, + "p5A"=Fyn[1,14], "p5C"=Fyn[2,14], "p5D"=Fyn[3,14], "p5E"=Fyn[4,14], "p5F"=Fyn[5,14],"p5G"=Fyn[6,14],"p5H"=Fyn[7,14],"p5I"=Fyn[8,14],"p5K"=Fyn[9,14],"p5L"=Fyn[10,14],"p5M"=Fyn[11,14],"p5N"=Fyn[12,14],"p5P"=Fyn[13,14],"p5Q"=Fyn[14,14],"p5R"=Fyn[15,14], + "p5S"=Fyn[16,14],"p5T"=Fyn[17,14],"p5V"=Fyn[18,14],"p5W"=Fyn[19,14],"p5Y"=Fyn[20,14],"p5O"=1, + "p6A"=Fyn[1,15], "p6C"=Fyn[2,15], "p6D"=Fyn[3,15], "p6E"=Fyn[4,15], "p6F"=Fyn[5,15],"p6G"=Fyn[6,15],"p6H"=Fyn[7,15],"p6I"=Fyn[8,15],"p6K"=Fyn[9,15],"p6L"=Fyn[10,15],"p6M"=Fyn[11,15],"p6N"=Fyn[12,15],"p6P"=Fyn[13,15],"p6Q"=Fyn[14,15],"p6R"=Fyn[15,15], + "p6S"=Fyn[16,15],"p6T"=Fyn[17,15],"p6V"=Fyn[18,15],"p6W"=Fyn[19,15],"p6Y"=Fyn[20,15],"p6O"=1) + + Hck_props <- c("m6A"=Hck[1,3], "m6C"=Hck[2,3], "m6D"=Hck[3,3], "m6E"=Hck[4,3], "m6F"=Hck[5,3],"m6G"=Hck[6,3],"m6H"=Hck[7,3],"m6I"=Hck[8,3],"m6K"=Hck[9,3],"m6L"=Hck[10,3],"m6M"=Hck[11,3],"m6N"=Hck[12,3],"m6P"=Hck[13,3],"m6Q"=Hck[14,3],"m6R"=Hck[15,3], + "m6S"=Hck[16,3],"m6T"=Hck[17,3],"m6V"=Hck[18,3],"m6W"=Hck[19,3],"m6Y"=Hck[20,3],"m6O"=1, + "m5A"=Hck[1,4], "m5C"=Hck[2,4], "m5D"=Hck[3,4], "m5E"=Hck[4,4], "m5F"=Hck[5,4],"m5G"=Hck[6,4],"m5H"=Hck[7,4],"m5I"=Hck[8,4],"m5K"=Hck[9,4],"m5L"=Hck[10,4],"m5M"=Hck[11,4],"m5N"=Hck[12,4],"m5P"=Hck[13,4],"m5Q"=Hck[14,4],"m5R"=Hck[15,4], + "m5S"=Hck[16,4],"m5T"=Hck[17,4],"m5V"=Hck[18,4],"m5W"=Hck[19,4],"m5Y"=Hck[20,4],"m5O"=1, + "m4A"=Hck[1,5], "m4C"=Hck[2,5], "m4D"=Hck[3,5], "m4E"=Hck[4,5], "m4F"=Hck[5,5],"m4G"=Hck[6,5],"m4H"=Hck[7,5],"m4I"=Hck[8,5],"m4K"=Hck[9,5],"m4L"=Hck[10,5],"m4M"=Hck[11,5],"m4N"=Hck[12,5],"m4P"=Hck[13,5],"m4Q"=Hck[14,5],"m4R"=Hck[15,5], + "m4S"=Hck[16,5],"m4T"=Hck[17,5],"m4V"=Hck[18,5],"m4W"=Hck[19,5],"m4Y"=Hck[20,5],"m4O"=1, + "m3A"=Hck[1,6], "m3C"=Hck[2,6], "m3D"=Hck[3,6], "m3E"=Hck[4,6], "m3F"=Hck[5,6],"m3G"=Hck[6,6],"m3H"=Hck[7,6],"m3I"=Hck[8,6],"m3K"=Hck[9,6],"m3L"=Hck[10,6],"m3M"=Hck[11,6],"m3N"=Hck[12,6],"m3P"=Hck[13,6],"m3Q"=Hck[14,6],"m3R"=Hck[15,6], + "m3S"=Hck[16,6],"m3T"=Hck[17,6],"m3V"=Hck[18,6],"m3W"=Hck[19,6],"m3Y"=Hck[20,6],"m3O"=1, + "m2A"=Hck[1,7], "m2C"=Hck[2,7], "m2D"=Hck[3,7], "m2E"=Hck[4,7], "m2F"=Hck[5,7],"m2G"=Hck[6,7],"m2H"=Hck[7,7],"m2I"=Hck[8,7],"m2K"=Hck[9,7],"m2L"=Hck[10,7],"m2M"=Hck[11,7],"m2N"=Hck[12,7],"m2P"=Hck[13,7],"m2Q"=Hck[14,7],"m2R"=Hck[15,7], + "m2S"=Hck[16,7],"m2T"=Hck[17,7],"m2V"=Hck[18,7],"m2W"=Hck[19,7],"m2Y"=Hck[20,7],"m2O"=1, + "m1A"=Hck[1,8], "m1C"=Hck[2,8], "m1D"=Hck[3,8], "m1E"=Hck[4,8], "m1F"=Hck[5,8],"m1G"=Hck[6,8],"m1H"=Hck[7,8],"m1I"=Hck[8,8],"m1K"=Hck[9,8],"m1L"=Hck[10,8],"m1M"=Hck[11,8],"m1N"=Hck[12,8],"m1P"=Hck[13,8],"m1Q"=Hck[14,8],"m1R"=Hck[15,8], + "m1S"=Hck[16,8],"m1T"=Hck[17,8],"m1V"=Hck[18,8],"m1W"=Hck[19,8],"m1Y"=Hck[20,8],"m1O"=1, + "d0A"=Hck[1,9], "d0C"=Hck[2,9], "d0D"=Hck[3,9], "d0E"=Hck[4,9], "d0F"=Hck[5,9],"d0G"=Hck[6,9],"d0H"=Hck[7,9],"d0I"=Hck[8,9],"d0K"=Hck[9,9],"d0L"=Hck[10,9],"d0M"=Hck[11,9],"d0N"=Hck[12,9],"d0P"=Hck[13,9],"d0Q"=Hck[14,9],"d0R"=Hck[15,9], + "d0S"=Hck[16,9],"d0T"=Hck[17,9],"d0V"=Hck[18,9],"d0W"=Hck[19,9],"d0Y"=Hck[20,9],"d0O"=1, + "p1A"=Hck[1,10], "p1C"=Hck[2,10], "p1D"=Hck[3,10], "p1E"=Hck[4,10], "p1F"=Hck[5,10],"p1G"=Hck[6,10],"p1H"=Hck[7,10],"p1I"=Hck[8,10],"p1K"=Hck[9,10],"p1L"=Hck[10,10],"p1M"=Hck[11,10],"p1N"=Hck[12,10],"p1P"=Hck[13,10],"p1Q"=Hck[14,10],"p1R"=Hck[15,10], + "p1S"=Hck[16,10],"p1T"=Hck[17,10],"p1V"=Hck[18,10],"p1W"=Hck[19,10],"p1Y"=Hck[20,10],"p1O"=1, + "p2A"=Hck[1,11], "p2C"=Hck[2,11], "p2D"=Hck[3,11], "p2E"=Hck[4,11], "p2F"=Hck[5,11],"p2G"=Hck[6,11],"p2H"=Hck[7,11],"p2I"=Hck[8,11],"p2K"=Hck[9,11],"p2L"=Hck[10,11],"p2M"=Hck[11,11],"p2N"=Hck[12,11],"p2P"=Hck[13,11],"p2Q"=Hck[14,11],"p2R"=Hck[15,11], + "p2S"=Hck[16,11],"p2T"=Hck[17,11],"p2V"=Hck[18,11],"p2W"=Hck[19,11],"p2Y"=Hck[20,11],"p2O"=1, + "p3A"=Hck[1,12], "p3C"=Hck[2,12], "p3D"=Hck[3,12], "p3E"=Hck[4,12], "p3F"=Hck[5,12],"p3G"=Hck[6,12],"p3H"=Hck[7,12],"p3I"=Hck[8,12],"p3K"=Hck[9,12],"p3L"=Hck[10,12],"p3M"=Hck[11,12],"p3N"=Hck[12,12],"p3P"=Hck[13,12],"p3Q"=Hck[14,12],"p3R"=Hck[15,12], + "p3S"=Hck[16,12],"p3T"=Hck[17,12],"p3V"=Hck[18,12],"p3W"=Hck[19,12],"p3Y"=Hck[20,12],"p3O"=1, + "p4A"=Hck[1,13], "p4C"=Hck[2,13], "p4D"=Hck[3,13], "p4E"=Hck[4,13], "p4F"=Hck[5,13],"p4G"=Hck[6,13],"p4H"=Hck[7,13],"p4I"=Hck[8,13],"p4K"=Hck[9,13],"p4L"=Hck[10,13],"p4M"=Hck[11,13],"p4N"=Hck[12,13],"p4P"=Hck[13,13],"p4Q"=Hck[14,13],"p4R"=Hck[15,13], + "p4S"=Hck[16,13],"p4T"=Hck[17,13],"p4V"=Hck[18,13],"p4W"=Hck[19,13],"p4Y"=Hck[20,13],"p4O"=1, + "p5A"=Hck[1,14], "p5C"=Hck[2,14], "p5D"=Hck[3,14], "p5E"=Hck[4,14], "p5F"=Hck[5,14],"p5G"=Hck[6,14],"p5H"=Hck[7,14],"p5I"=Hck[8,14],"p5K"=Hck[9,14],"p5L"=Hck[10,14],"p5M"=Hck[11,14],"p5N"=Hck[12,14],"p5P"=Hck[13,14],"p5Q"=Hck[14,14],"p5R"=Hck[15,14], + "p5S"=Hck[16,14],"p5T"=Hck[17,14],"p5V"=Hck[18,14],"p5W"=Hck[19,14],"p5Y"=Hck[20,14],"p5O"=1, + "p6A"=Hck[1,15], "p6C"=Hck[2,15], "p6D"=Hck[3,15], "p6E"=Hck[4,15], "p6F"=Hck[5,15],"p6G"=Hck[6,15],"p6H"=Hck[7,15],"p6I"=Hck[8,15],"p6K"=Hck[9,15],"p6L"=Hck[10,15],"p6M"=Hck[11,15],"p6N"=Hck[12,15],"p6P"=Hck[13,15],"p6Q"=Hck[14,15],"p6R"=Hck[15,15], + "p6S"=Hck[16,15],"p6T"=Hck[17,15],"p6V"=Hck[18,15],"p6W"=Hck[19,15],"p6Y"=Hck[20,15],"p6O"=1) + + JAK2_props <- c("m6A"=JAK2[1,3], "m6C"=JAK2[2,3], "m6D"=JAK2[3,3], "m6E"=JAK2[4,3], "m6F"=JAK2[5,3],"m6G"=JAK2[6,3],"m6H"=JAK2[7,3],"m6I"=JAK2[8,3],"m6K"=JAK2[9,3],"m6L"=JAK2[10,3],"m6M"=JAK2[11,3],"m6N"=JAK2[12,3],"m6P"=JAK2[13,3],"m6Q"=JAK2[14,3],"m6R"=JAK2[15,3], + "m6S"=JAK2[16,3],"m6T"=JAK2[17,3],"m6V"=JAK2[18,3],"m6W"=JAK2[19,3],"m6Y"=JAK2[20,3],"m6O"=1, + "m5A"=JAK2[1,4], "m5C"=JAK2[2,4], "m5D"=JAK2[3,4], "m5E"=JAK2[4,4], "m5F"=JAK2[5,4],"m5G"=JAK2[6,4],"m5H"=JAK2[7,4],"m5I"=JAK2[8,4],"m5K"=JAK2[9,4],"m5L"=JAK2[10,4],"m5M"=JAK2[11,4],"m5N"=JAK2[12,4],"m5P"=JAK2[13,4],"m5Q"=JAK2[14,4],"m5R"=JAK2[15,4], + "m5S"=JAK2[16,4],"m5T"=JAK2[17,4],"m5V"=JAK2[18,4],"m5W"=JAK2[19,4],"m5Y"=JAK2[20,4],"m5O"=1, + "m4A"=JAK2[1,5], "m4C"=JAK2[2,5], "m4D"=JAK2[3,5], "m4E"=JAK2[4,5], "m4F"=JAK2[5,5],"m4G"=JAK2[6,5],"m4H"=JAK2[7,5],"m4I"=JAK2[8,5],"m4K"=JAK2[9,5],"m4L"=JAK2[10,5],"m4M"=JAK2[11,5],"m4N"=JAK2[12,5],"m4P"=JAK2[13,5],"m4Q"=JAK2[14,5],"m4R"=JAK2[15,5], + "m4S"=JAK2[16,5],"m4T"=JAK2[17,5],"m4V"=JAK2[18,5],"m4W"=JAK2[19,5],"m4Y"=JAK2[20,5],"m4O"=1, + "m3A"=JAK2[1,6], "m3C"=JAK2[2,6], "m3D"=JAK2[3,6], "m3E"=JAK2[4,6], "m3F"=JAK2[5,6],"m3G"=JAK2[6,6],"m3H"=JAK2[7,6],"m3I"=JAK2[8,6],"m3K"=JAK2[9,6],"m3L"=JAK2[10,6],"m3M"=JAK2[11,6],"m3N"=JAK2[12,6],"m3P"=JAK2[13,6],"m3Q"=JAK2[14,6],"m3R"=JAK2[15,6], + "m3S"=JAK2[16,6],"m3T"=JAK2[17,6],"m3V"=JAK2[18,6],"m3W"=JAK2[19,6],"m3Y"=JAK2[20,6],"m3O"=1, + "m2A"=JAK2[1,7], "m2C"=JAK2[2,7], "m2D"=JAK2[3,7], "m2E"=JAK2[4,7], "m2F"=JAK2[5,7],"m2G"=JAK2[6,7],"m2H"=JAK2[7,7],"m2I"=JAK2[8,7],"m2K"=JAK2[9,7],"m2L"=JAK2[10,7],"m2M"=JAK2[11,7],"m2N"=JAK2[12,7],"m2P"=JAK2[13,7],"m2Q"=JAK2[14,7],"m2R"=JAK2[15,7], + "m2S"=JAK2[16,7],"m2T"=JAK2[17,7],"m2V"=JAK2[18,7],"m2W"=JAK2[19,7],"m2Y"=JAK2[20,7],"m2O"=1, + "m1A"=JAK2[1,8], "m1C"=JAK2[2,8], "m1D"=JAK2[3,8], "m1E"=JAK2[4,8], "m1F"=JAK2[5,8],"m1G"=JAK2[6,8],"m1H"=JAK2[7,8],"m1I"=JAK2[8,8],"m1K"=JAK2[9,8],"m1L"=JAK2[10,8],"m1M"=JAK2[11,8],"m1N"=JAK2[12,8],"m1P"=JAK2[13,8],"m1Q"=JAK2[14,8],"m1R"=JAK2[15,8], + "m1S"=JAK2[16,8],"m1T"=JAK2[17,8],"m1V"=JAK2[18,8],"m1W"=JAK2[19,8],"m1Y"=JAK2[20,8],"m1O"=1, + "d0A"=JAK2[1,9], "d0C"=JAK2[2,9], "d0D"=JAK2[3,9], "d0E"=JAK2[4,9], "d0F"=JAK2[5,9],"d0G"=JAK2[6,9],"d0H"=JAK2[7,9],"d0I"=JAK2[8,9],"d0K"=JAK2[9,9],"d0L"=JAK2[10,9],"d0M"=JAK2[11,9],"d0N"=JAK2[12,9],"d0P"=JAK2[13,9],"d0Q"=JAK2[14,9],"d0R"=JAK2[15,9], + "d0S"=JAK2[16,9],"d0T"=JAK2[17,9],"d0V"=JAK2[18,9],"d0W"=JAK2[19,9],"d0Y"=JAK2[20,9],"d0O"=1, + "p1A"=JAK2[1,10], "p1C"=JAK2[2,10], "p1D"=JAK2[3,10], "p1E"=JAK2[4,10], "p1F"=JAK2[5,10],"p1G"=JAK2[6,10],"p1H"=JAK2[7,10],"p1I"=JAK2[8,10],"p1K"=JAK2[9,10],"p1L"=JAK2[10,10],"p1M"=JAK2[11,10],"p1N"=JAK2[12,10],"p1P"=JAK2[13,10],"p1Q"=JAK2[14,10],"p1R"=JAK2[15,10], + "p1S"=JAK2[16,10],"p1T"=JAK2[17,10],"p1V"=JAK2[18,10],"p1W"=JAK2[19,10],"p1Y"=JAK2[20,10],"p1O"=1, + "p2A"=JAK2[1,11], "p2C"=JAK2[2,11], "p2D"=JAK2[3,11], "p2E"=JAK2[4,11], "p2F"=JAK2[5,11],"p2G"=JAK2[6,11],"p2H"=JAK2[7,11],"p2I"=JAK2[8,11],"p2K"=JAK2[9,11],"p2L"=JAK2[10,11],"p2M"=JAK2[11,11],"p2N"=JAK2[12,11],"p2P"=JAK2[13,11],"p2Q"=JAK2[14,11],"p2R"=JAK2[15,11], + "p2S"=JAK2[16,11],"p2T"=JAK2[17,11],"p2V"=JAK2[18,11],"p2W"=JAK2[19,11],"p2Y"=JAK2[20,11],"p2O"=1, + "p3A"=JAK2[1,12], "p3C"=JAK2[2,12], "p3D"=JAK2[3,12], "p3E"=JAK2[4,12], "p3F"=JAK2[5,12],"p3G"=JAK2[6,12],"p3H"=JAK2[7,12],"p3I"=JAK2[8,12],"p3K"=JAK2[9,12],"p3L"=JAK2[10,12],"p3M"=JAK2[11,12],"p3N"=JAK2[12,12],"p3P"=JAK2[13,12],"p3Q"=JAK2[14,12],"p3R"=JAK2[15,12], + "p3S"=JAK2[16,12],"p3T"=JAK2[17,12],"p3V"=JAK2[18,12],"p3W"=JAK2[19,12],"p3Y"=JAK2[20,12],"p3O"=1, + "p4A"=JAK2[1,13], "p4C"=JAK2[2,13], "p4D"=JAK2[3,13], "p4E"=JAK2[4,13], "p4F"=JAK2[5,13],"p4G"=JAK2[6,13],"p4H"=JAK2[7,13],"p4I"=JAK2[8,13],"p4K"=JAK2[9,13],"p4L"=JAK2[10,13],"p4M"=JAK2[11,13],"p4N"=JAK2[12,13],"p4P"=JAK2[13,13],"p4Q"=JAK2[14,13],"p4R"=JAK2[15,13], + "p4S"=JAK2[16,13],"p4T"=JAK2[17,13],"p4V"=JAK2[18,13],"p4W"=JAK2[19,13],"p4Y"=JAK2[20,13],"p4O"=1, + "p5A"=JAK2[1,14], "p5C"=JAK2[2,14], "p5D"=JAK2[3,14], "p5E"=JAK2[4,14], "p5F"=JAK2[5,14],"p5G"=JAK2[6,14],"p5H"=JAK2[7,14],"p5I"=JAK2[8,14],"p5K"=JAK2[9,14],"p5L"=JAK2[10,14],"p5M"=JAK2[11,14],"p5N"=JAK2[12,14],"p5P"=JAK2[13,14],"p5Q"=JAK2[14,14],"p5R"=JAK2[15,14], + "p5S"=JAK2[16,14],"p5T"=JAK2[17,14],"p5V"=JAK2[18,14],"p5W"=JAK2[19,14],"p5Y"=JAK2[20,14],"p5O"=1, + "p6A"=JAK2[1,15], "p6C"=JAK2[2,15], "p6D"=JAK2[3,15], "p6E"=JAK2[4,15], "p6F"=JAK2[5,15],"p6G"=JAK2[6,15],"p6H"=JAK2[7,15],"p6I"=JAK2[8,15],"p6K"=JAK2[9,15],"p6L"=JAK2[10,15],"p6M"=JAK2[11,15],"p6N"=JAK2[12,15],"p6P"=JAK2[13,15],"p6Q"=JAK2[14,15],"p6R"=JAK2[15,15], + "p6S"=JAK2[16,15],"p6T"=JAK2[17,15],"p6V"=JAK2[18,15],"p6W"=JAK2[19,15],"p6Y"=JAK2[20,15],"p6O"=1) + + Lck_props <- c("m6A"=Lck[1,3], "m6C"=Lck[2,3], "m6D"=Lck[3,3], "m6E"=Lck[4,3], "m6F"=Lck[5,3],"m6G"=Lck[6,3],"m6H"=Lck[7,3],"m6I"=Lck[8,3],"m6K"=Lck[9,3],"m6L"=Lck[10,3],"m6M"=Lck[11,3],"m6N"=Lck[12,3],"m6P"=Lck[13,3],"m6Q"=Lck[14,3],"m6R"=Lck[15,3], + "m6S"=Lck[16,3],"m6T"=Lck[17,3],"m6V"=Lck[18,3],"m6W"=Lck[19,3],"m6Y"=Lck[20,3],"m6O"=1, + "m5A"=Lck[1,4], "m5C"=Lck[2,4], "m5D"=Lck[3,4], "m5E"=Lck[4,4], "m5F"=Lck[5,4],"m5G"=Lck[6,4],"m5H"=Lck[7,4],"m5I"=Lck[8,4],"m5K"=Lck[9,4],"m5L"=Lck[10,4],"m5M"=Lck[11,4],"m5N"=Lck[12,4],"m5P"=Lck[13,4],"m5Q"=Lck[14,4],"m5R"=Lck[15,4], + "m5S"=Lck[16,4],"m5T"=Lck[17,4],"m5V"=Lck[18,4],"m5W"=Lck[19,4],"m5Y"=Lck[20,4],"m5O"=1, + "m4A"=Lck[1,5], "m4C"=Lck[2,5], "m4D"=Lck[3,5], "m4E"=Lck[4,5], "m4F"=Lck[5,5],"m4G"=Lck[6,5],"m4H"=Lck[7,5],"m4I"=Lck[8,5],"m4K"=Lck[9,5],"m4L"=Lck[10,5],"m4M"=Lck[11,5],"m4N"=Lck[12,5],"m4P"=Lck[13,5],"m4Q"=Lck[14,5],"m4R"=Lck[15,5], + "m4S"=Lck[16,5],"m4T"=Lck[17,5],"m4V"=Lck[18,5],"m4W"=Lck[19,5],"m4Y"=Lck[20,5],"m4O"=1, + "m3A"=Lck[1,6], "m3C"=Lck[2,6], "m3D"=Lck[3,6], "m3E"=Lck[4,6], "m3F"=Lck[5,6],"m3G"=Lck[6,6],"m3H"=Lck[7,6],"m3I"=Lck[8,6],"m3K"=Lck[9,6],"m3L"=Lck[10,6],"m3M"=Lck[11,6],"m3N"=Lck[12,6],"m3P"=Lck[13,6],"m3Q"=Lck[14,6],"m3R"=Lck[15,6], + "m3S"=Lck[16,6],"m3T"=Lck[17,6],"m3V"=Lck[18,6],"m3W"=Lck[19,6],"m3Y"=Lck[20,6],"m3O"=1, + "m2A"=Lck[1,7], "m2C"=Lck[2,7], "m2D"=Lck[3,7], "m2E"=Lck[4,7], "m2F"=Lck[5,7],"m2G"=Lck[6,7],"m2H"=Lck[7,7],"m2I"=Lck[8,7],"m2K"=Lck[9,7],"m2L"=Lck[10,7],"m2M"=Lck[11,7],"m2N"=Lck[12,7],"m2P"=Lck[13,7],"m2Q"=Lck[14,7],"m2R"=Lck[15,7], + "m2S"=Lck[16,7],"m2T"=Lck[17,7],"m2V"=Lck[18,7],"m2W"=Lck[19,7],"m2Y"=Lck[20,7],"m2O"=1, + "m1A"=Lck[1,8], "m1C"=Lck[2,8], "m1D"=Lck[3,8], "m1E"=Lck[4,8], "m1F"=Lck[5,8],"m1G"=Lck[6,8],"m1H"=Lck[7,8],"m1I"=Lck[8,8],"m1K"=Lck[9,8],"m1L"=Lck[10,8],"m1M"=Lck[11,8],"m1N"=Lck[12,8],"m1P"=Lck[13,8],"m1Q"=Lck[14,8],"m1R"=Lck[15,8], + "m1S"=Lck[16,8],"m1T"=Lck[17,8],"m1V"=Lck[18,8],"m1W"=Lck[19,8],"m1Y"=Lck[20,8],"m1O"=1, + "d0A"=Lck[1,9], "d0C"=Lck[2,9], "d0D"=Lck[3,9], "d0E"=Lck[4,9], "d0F"=Lck[5,9],"d0G"=Lck[6,9],"d0H"=Lck[7,9],"d0I"=Lck[8,9],"d0K"=Lck[9,9],"d0L"=Lck[10,9],"d0M"=Lck[11,9],"d0N"=Lck[12,9],"d0P"=Lck[13,9],"d0Q"=Lck[14,9],"d0R"=Lck[15,9], + "d0S"=Lck[16,9],"d0T"=Lck[17,9],"d0V"=Lck[18,9],"d0W"=Lck[19,9],"d0Y"=Lck[20,9],"d0O"=1, + "p1A"=Lck[1,10], "p1C"=Lck[2,10], "p1D"=Lck[3,10], "p1E"=Lck[4,10], "p1F"=Lck[5,10],"p1G"=Lck[6,10],"p1H"=Lck[7,10],"p1I"=Lck[8,10],"p1K"=Lck[9,10],"p1L"=Lck[10,10],"p1M"=Lck[11,10],"p1N"=Lck[12,10],"p1P"=Lck[13,10],"p1Q"=Lck[14,10],"p1R"=Lck[15,10], + "p1S"=Lck[16,10],"p1T"=Lck[17,10],"p1V"=Lck[18,10],"p1W"=Lck[19,10],"p1Y"=Lck[20,10],"p1O"=1, + "p2A"=Lck[1,11], "p2C"=Lck[2,11], "p2D"=Lck[3,11], "p2E"=Lck[4,11], "p2F"=Lck[5,11],"p2G"=Lck[6,11],"p2H"=Lck[7,11],"p2I"=Lck[8,11],"p2K"=Lck[9,11],"p2L"=Lck[10,11],"p2M"=Lck[11,11],"p2N"=Lck[12,11],"p2P"=Lck[13,11],"p2Q"=Lck[14,11],"p2R"=Lck[15,11], + "p2S"=Lck[16,11],"p2T"=Lck[17,11],"p2V"=Lck[18,11],"p2W"=Lck[19,11],"p2Y"=Lck[20,11],"p2O"=1, + "p3A"=Lck[1,12], "p3C"=Lck[2,12], "p3D"=Lck[3,12], "p3E"=Lck[4,12], "p3F"=Lck[5,12],"p3G"=Lck[6,12],"p3H"=Lck[7,12],"p3I"=Lck[8,12],"p3K"=Lck[9,12],"p3L"=Lck[10,12],"p3M"=Lck[11,12],"p3N"=Lck[12,12],"p3P"=Lck[13,12],"p3Q"=Lck[14,12],"p3R"=Lck[15,12], + "p3S"=Lck[16,12],"p3T"=Lck[17,12],"p3V"=Lck[18,12],"p3W"=Lck[19,12],"p3Y"=Lck[20,12],"p3O"=1, + "p4A"=Lck[1,13], "p4C"=Lck[2,13], "p4D"=Lck[3,13], "p4E"=Lck[4,13], "p4F"=Lck[5,13],"p4G"=Lck[6,13],"p4H"=Lck[7,13],"p4I"=Lck[8,13],"p4K"=Lck[9,13],"p4L"=Lck[10,13],"p4M"=Lck[11,13],"p4N"=Lck[12,13],"p4P"=Lck[13,13],"p4Q"=Lck[14,13],"p4R"=Lck[15,13], + "p4S"=Lck[16,13],"p4T"=Lck[17,13],"p4V"=Lck[18,13],"p4W"=Lck[19,13],"p4Y"=Lck[20,13],"p4O"=1, + "p5A"=Lck[1,14], "p5C"=Lck[2,14], "p5D"=Lck[3,14], "p5E"=Lck[4,14], "p5F"=Lck[5,14],"p5G"=Lck[6,14],"p5H"=Lck[7,14],"p5I"=Lck[8,14],"p5K"=Lck[9,14],"p5L"=Lck[10,14],"p5M"=Lck[11,14],"p5N"=Lck[12,14],"p5P"=Lck[13,14],"p5Q"=Lck[14,14],"p5R"=Lck[15,14], + "p5S"=Lck[16,14],"p5T"=Lck[17,14],"p5V"=Lck[18,14],"p5W"=Lck[19,14],"p5Y"=Lck[20,14],"p5O"=1, + "p6A"=Lck[1,15], "p6C"=Lck[2,15], "p6D"=Lck[3,15], "p6E"=Lck[4,15], "p6F"=Lck[5,15],"p6G"=Lck[6,15],"p6H"=Lck[7,15],"p6I"=Lck[8,15],"p6K"=Lck[9,15],"p6L"=Lck[10,15],"p6M"=Lck[11,15],"p6N"=Lck[12,15],"p6P"=Lck[13,15],"p6Q"=Lck[14,15],"p6R"=Lck[15,15], + "p6S"=Lck[16,15],"p6T"=Lck[17,15],"p6V"=Lck[18,15],"p6W"=Lck[19,15],"p6Y"=Lck[20,15],"p6O"=1) + + Lyn_props <- c("m6A"=Lyn[1,3], "m6C"=Lyn[2,3], "m6D"=Lyn[3,3], "m6E"=Lyn[4,3], "m6F"=Lyn[5,3],"m6G"=Lyn[6,3],"m6H"=Lyn[7,3],"m6I"=Lyn[8,3],"m6K"=Lyn[9,3],"m6L"=Lyn[10,3],"m6M"=Lyn[11,3],"m6N"=Lyn[12,3],"m6P"=Lyn[13,3],"m6Q"=Lyn[14,3],"m6R"=Lyn[15,3], + "m6S"=Lyn[16,3],"m6T"=Lyn[17,3],"m6V"=Lyn[18,3],"m6W"=Lyn[19,3],"m6Y"=Lyn[20,3],"m6O"=1, + "m5A"=Lyn[1,4], "m5C"=Lyn[2,4], "m5D"=Lyn[3,4], "m5E"=Lyn[4,4], "m5F"=Lyn[5,4],"m5G"=Lyn[6,4],"m5H"=Lyn[7,4],"m5I"=Lyn[8,4],"m5K"=Lyn[9,4],"m5L"=Lyn[10,4],"m5M"=Lyn[11,4],"m5N"=Lyn[12,4],"m5P"=Lyn[13,4],"m5Q"=Lyn[14,4],"m5R"=Lyn[15,4], + "m5S"=Lyn[16,4],"m5T"=Lyn[17,4],"m5V"=Lyn[18,4],"m5W"=Lyn[19,4],"m5Y"=Lyn[20,4],"m5O"=1, + "m4A"=Lyn[1,5], "m4C"=Lyn[2,5], "m4D"=Lyn[3,5], "m4E"=Lyn[4,5], "m4F"=Lyn[5,5],"m4G"=Lyn[6,5],"m4H"=Lyn[7,5],"m4I"=Lyn[8,5],"m4K"=Lyn[9,5],"m4L"=Lyn[10,5],"m4M"=Lyn[11,5],"m4N"=Lyn[12,5],"m4P"=Lyn[13,5],"m4Q"=Lyn[14,5],"m4R"=Lyn[15,5], + "m4S"=Lyn[16,5],"m4T"=Lyn[17,5],"m4V"=Lyn[18,5],"m4W"=Lyn[19,5],"m4Y"=Lyn[20,5],"m4O"=1, + "m3A"=Lyn[1,6], "m3C"=Lyn[2,6], "m3D"=Lyn[3,6], "m3E"=Lyn[4,6], "m3F"=Lyn[5,6],"m3G"=Lyn[6,6],"m3H"=Lyn[7,6],"m3I"=Lyn[8,6],"m3K"=Lyn[9,6],"m3L"=Lyn[10,6],"m3M"=Lyn[11,6],"m3N"=Lyn[12,6],"m3P"=Lyn[13,6],"m3Q"=Lyn[14,6],"m3R"=Lyn[15,6], + "m3S"=Lyn[16,6],"m3T"=Lyn[17,6],"m3V"=Lyn[18,6],"m3W"=Lyn[19,6],"m3Y"=Lyn[20,6],"m3O"=1, + "m2A"=Lyn[1,7], "m2C"=Lyn[2,7], "m2D"=Lyn[3,7], "m2E"=Lyn[4,7], "m2F"=Lyn[5,7],"m2G"=Lyn[6,7],"m2H"=Lyn[7,7],"m2I"=Lyn[8,7],"m2K"=Lyn[9,7],"m2L"=Lyn[10,7],"m2M"=Lyn[11,7],"m2N"=Lyn[12,7],"m2P"=Lyn[13,7],"m2Q"=Lyn[14,7],"m2R"=Lyn[15,7], + "m2S"=Lyn[16,7],"m2T"=Lyn[17,7],"m2V"=Lyn[18,7],"m2W"=Lyn[19,7],"m2Y"=Lyn[20,7],"m2O"=1, + "m1A"=Lyn[1,8], "m1C"=Lyn[2,8], "m1D"=Lyn[3,8], "m1E"=Lyn[4,8], "m1F"=Lyn[5,8],"m1G"=Lyn[6,8],"m1H"=Lyn[7,8],"m1I"=Lyn[8,8],"m1K"=Lyn[9,8],"m1L"=Lyn[10,8],"m1M"=Lyn[11,8],"m1N"=Lyn[12,8],"m1P"=Lyn[13,8],"m1Q"=Lyn[14,8],"m1R"=Lyn[15,8], + "m1S"=Lyn[16,8],"m1T"=Lyn[17,8],"m1V"=Lyn[18,8],"m1W"=Lyn[19,8],"m1Y"=Lyn[20,8],"m1O"=1, + "d0A"=Lyn[1,9], "d0C"=Lyn[2,9], "d0D"=Lyn[3,9], "d0E"=Lyn[4,9], "d0F"=Lyn[5,9],"d0G"=Lyn[6,9],"d0H"=Lyn[7,9],"d0I"=Lyn[8,9],"d0K"=Lyn[9,9],"d0L"=Lyn[10,9],"d0M"=Lyn[11,9],"d0N"=Lyn[12,9],"d0P"=Lyn[13,9],"d0Q"=Lyn[14,9],"d0R"=Lyn[15,9], + "d0S"=Lyn[16,9],"d0T"=Lyn[17,9],"d0V"=Lyn[18,9],"d0W"=Lyn[19,9],"d0Y"=Lyn[20,9],"d0O"=1, + "p1A"=Lyn[1,10], "p1C"=Lyn[2,10], "p1D"=Lyn[3,10], "p1E"=Lyn[4,10], "p1F"=Lyn[5,10],"p1G"=Lyn[6,10],"p1H"=Lyn[7,10],"p1I"=Lyn[8,10],"p1K"=Lyn[9,10],"p1L"=Lyn[10,10],"p1M"=Lyn[11,10],"p1N"=Lyn[12,10],"p1P"=Lyn[13,10],"p1Q"=Lyn[14,10],"p1R"=Lyn[15,10], + "p1S"=Lyn[16,10],"p1T"=Lyn[17,10],"p1V"=Lyn[18,10],"p1W"=Lyn[19,10],"p1Y"=Lyn[20,10],"p1O"=1, + "p2A"=Lyn[1,11], "p2C"=Lyn[2,11], "p2D"=Lyn[3,11], "p2E"=Lyn[4,11], "p2F"=Lyn[5,11],"p2G"=Lyn[6,11],"p2H"=Lyn[7,11],"p2I"=Lyn[8,11],"p2K"=Lyn[9,11],"p2L"=Lyn[10,11],"p2M"=Lyn[11,11],"p2N"=Lyn[12,11],"p2P"=Lyn[13,11],"p2Q"=Lyn[14,11],"p2R"=Lyn[15,11], + "p2S"=Lyn[16,11],"p2T"=Lyn[17,11],"p2V"=Lyn[18,11],"p2W"=Lyn[19,11],"p2Y"=Lyn[20,11],"p2O"=1, + "p3A"=Lyn[1,12], "p3C"=Lyn[2,12], "p3D"=Lyn[3,12], "p3E"=Lyn[4,12], "p3F"=Lyn[5,12],"p3G"=Lyn[6,12],"p3H"=Lyn[7,12],"p3I"=Lyn[8,12],"p3K"=Lyn[9,12],"p3L"=Lyn[10,12],"p3M"=Lyn[11,12],"p3N"=Lyn[12,12],"p3P"=Lyn[13,12],"p3Q"=Lyn[14,12],"p3R"=Lyn[15,12], + "p3S"=Lyn[16,12],"p3T"=Lyn[17,12],"p3V"=Lyn[18,12],"p3W"=Lyn[19,12],"p3Y"=Lyn[20,12],"p3O"=1, + "p4A"=Lyn[1,13], "p4C"=Lyn[2,13], "p4D"=Lyn[3,13], "p4E"=Lyn[4,13], "p4F"=Lyn[5,13],"p4G"=Lyn[6,13],"p4H"=Lyn[7,13],"p4I"=Lyn[8,13],"p4K"=Lyn[9,13],"p4L"=Lyn[10,13],"p4M"=Lyn[11,13],"p4N"=Lyn[12,13],"p4P"=Lyn[13,13],"p4Q"=Lyn[14,13],"p4R"=Lyn[15,13], + "p4S"=Lyn[16,13],"p4T"=Lyn[17,13],"p4V"=Lyn[18,13],"p4W"=Lyn[19,13],"p4Y"=Lyn[20,13],"p4O"=1, + "p5A"=Lyn[1,14], "p5C"=Lyn[2,14], "p5D"=Lyn[3,14], "p5E"=Lyn[4,14], "p5F"=Lyn[5,14],"p5G"=Lyn[6,14],"p5H"=Lyn[7,14],"p5I"=Lyn[8,14],"p5K"=Lyn[9,14],"p5L"=Lyn[10,14],"p5M"=Lyn[11,14],"p5N"=Lyn[12,14],"p5P"=Lyn[13,14],"p5Q"=Lyn[14,14],"p5R"=Lyn[15,14], + "p5S"=Lyn[16,14],"p5T"=Lyn[17,14],"p5V"=Lyn[18,14],"p5W"=Lyn[19,14],"p5Y"=Lyn[20,14],"p5O"=1, + "p6A"=Lyn[1,15], "p6C"=Lyn[2,15], "p6D"=Lyn[3,15], "p6E"=Lyn[4,15], "p6F"=Lyn[5,15],"p6G"=Lyn[6,15],"p6H"=Lyn[7,15],"p6I"=Lyn[8,15],"p6K"=Lyn[9,15],"p6L"=Lyn[10,15],"p6M"=Lyn[11,15],"p6N"=Lyn[12,15],"p6P"=Lyn[13,15],"p6Q"=Lyn[14,15],"p6R"=Lyn[15,15], + "p6S"=Lyn[16,15],"p6T"=Lyn[17,15],"p6V"=Lyn[18,15],"p6W"=Lyn[19,15],"p6Y"=Lyn[20,15],"p6O"=1) + + Pyk2_props <- c("m6A"=Pyk2[1,3], "m6C"=Pyk2[2,3], "m6D"=Pyk2[3,3], "m6E"=Pyk2[4,3], "m6F"=Pyk2[5,3],"m6G"=Pyk2[6,3],"m6H"=Pyk2[7,3],"m6I"=Pyk2[8,3],"m6K"=Pyk2[9,3],"m6L"=Pyk2[10,3],"m6M"=Pyk2[11,3],"m6N"=Pyk2[12,3],"m6P"=Pyk2[13,3],"m6Q"=Pyk2[14,3],"m6R"=Pyk2[15,3], + "m6S"=Pyk2[16,3],"m6T"=Pyk2[17,3],"m6V"=Pyk2[18,3],"m6W"=Pyk2[19,3],"m6Y"=Pyk2[20,3],"m6O"=1, + "m5A"=Pyk2[1,4], "m5C"=Pyk2[2,4], "m5D"=Pyk2[3,4], "m5E"=Pyk2[4,4], "m5F"=Pyk2[5,4],"m5G"=Pyk2[6,4],"m5H"=Pyk2[7,4],"m5I"=Pyk2[8,4],"m5K"=Pyk2[9,4],"m5L"=Pyk2[10,4],"m5M"=Pyk2[11,4],"m5N"=Pyk2[12,4],"m5P"=Pyk2[13,4],"m5Q"=Pyk2[14,4],"m5R"=Pyk2[15,4], + "m5S"=Pyk2[16,4],"m5T"=Pyk2[17,4],"m5V"=Pyk2[18,4],"m5W"=Pyk2[19,4],"m5Y"=Pyk2[20,4],"m5O"=1, + "m4A"=Pyk2[1,5], "m4C"=Pyk2[2,5], "m4D"=Pyk2[3,5], "m4E"=Pyk2[4,5], "m4F"=Pyk2[5,5],"m4G"=Pyk2[6,5],"m4H"=Pyk2[7,5],"m4I"=Pyk2[8,5],"m4K"=Pyk2[9,5],"m4L"=Pyk2[10,5],"m4M"=Pyk2[11,5],"m4N"=Pyk2[12,5],"m4P"=Pyk2[13,5],"m4Q"=Pyk2[14,5],"m4R"=Pyk2[15,5], + "m4S"=Pyk2[16,5],"m4T"=Pyk2[17,5],"m4V"=Pyk2[18,5],"m4W"=Pyk2[19,5],"m4Y"=Pyk2[20,5],"m4O"=1, + "m3A"=Pyk2[1,6], "m3C"=Pyk2[2,6], "m3D"=Pyk2[3,6], "m3E"=Pyk2[4,6], "m3F"=Pyk2[5,6],"m3G"=Pyk2[6,6],"m3H"=Pyk2[7,6],"m3I"=Pyk2[8,6],"m3K"=Pyk2[9,6],"m3L"=Pyk2[10,6],"m3M"=Pyk2[11,6],"m3N"=Pyk2[12,6],"m3P"=Pyk2[13,6],"m3Q"=Pyk2[14,6],"m3R"=Pyk2[15,6], + "m3S"=Pyk2[16,6],"m3T"=Pyk2[17,6],"m3V"=Pyk2[18,6],"m3W"=Pyk2[19,6],"m3Y"=Pyk2[20,6],"m3O"=1, + "m2A"=Pyk2[1,7], "m2C"=Pyk2[2,7], "m2D"=Pyk2[3,7], "m2E"=Pyk2[4,7], "m2F"=Pyk2[5,7],"m2G"=Pyk2[6,7],"m2H"=Pyk2[7,7],"m2I"=Pyk2[8,7],"m2K"=Pyk2[9,7],"m2L"=Pyk2[10,7],"m2M"=Pyk2[11,7],"m2N"=Pyk2[12,7],"m2P"=Pyk2[13,7],"m2Q"=Pyk2[14,7],"m2R"=Pyk2[15,7], + "m2S"=Pyk2[16,7],"m2T"=Pyk2[17,7],"m2V"=Pyk2[18,7],"m2W"=Pyk2[19,7],"m2Y"=Pyk2[20,7],"m2O"=1, + "m1A"=Pyk2[1,8], "m1C"=Pyk2[2,8], "m1D"=Pyk2[3,8], "m1E"=Pyk2[4,8], "m1F"=Pyk2[5,8],"m1G"=Pyk2[6,8],"m1H"=Pyk2[7,8],"m1I"=Pyk2[8,8],"m1K"=Pyk2[9,8],"m1L"=Pyk2[10,8],"m1M"=Pyk2[11,8],"m1N"=Pyk2[12,8],"m1P"=Pyk2[13,8],"m1Q"=Pyk2[14,8],"m1R"=Pyk2[15,8], + "m1S"=Pyk2[16,8],"m1T"=Pyk2[17,8],"m1V"=Pyk2[18,8],"m1W"=Pyk2[19,8],"m1Y"=Pyk2[20,8],"m1O"=1, + "d0A"=Pyk2[1,9], "d0C"=Pyk2[2,9], "d0D"=Pyk2[3,9], "d0E"=Pyk2[4,9], "d0F"=Pyk2[5,9],"d0G"=Pyk2[6,9],"d0H"=Pyk2[7,9],"d0I"=Pyk2[8,9],"d0K"=Pyk2[9,9],"d0L"=Pyk2[10,9],"d0M"=Pyk2[11,9],"d0N"=Pyk2[12,9],"d0P"=Pyk2[13,9],"d0Q"=Pyk2[14,9],"d0R"=Pyk2[15,9], + "d0S"=Pyk2[16,9],"d0T"=Pyk2[17,9],"d0V"=Pyk2[18,9],"d0W"=Pyk2[19,9],"d0Y"=Pyk2[20,9],"d0O"=1, + "p1A"=Pyk2[1,10], "p1C"=Pyk2[2,10], "p1D"=Pyk2[3,10], "p1E"=Pyk2[4,10], "p1F"=Pyk2[5,10],"p1G"=Pyk2[6,10],"p1H"=Pyk2[7,10],"p1I"=Pyk2[8,10],"p1K"=Pyk2[9,10],"p1L"=Pyk2[10,10],"p1M"=Pyk2[11,10],"p1N"=Pyk2[12,10],"p1P"=Pyk2[13,10],"p1Q"=Pyk2[14,10],"p1R"=Pyk2[15,10], + "p1S"=Pyk2[16,10],"p1T"=Pyk2[17,10],"p1V"=Pyk2[18,10],"p1W"=Pyk2[19,10],"p1Y"=Pyk2[20,10],"p1O"=1, + "p2A"=Pyk2[1,11], "p2C"=Pyk2[2,11], "p2D"=Pyk2[3,11], "p2E"=Pyk2[4,11], "p2F"=Pyk2[5,11],"p2G"=Pyk2[6,11],"p2H"=Pyk2[7,11],"p2I"=Pyk2[8,11],"p2K"=Pyk2[9,11],"p2L"=Pyk2[10,11],"p2M"=Pyk2[11,11],"p2N"=Pyk2[12,11],"p2P"=Pyk2[13,11],"p2Q"=Pyk2[14,11],"p2R"=Pyk2[15,11], + "p2S"=Pyk2[16,11],"p2T"=Pyk2[17,11],"p2V"=Pyk2[18,11],"p2W"=Pyk2[19,11],"p2Y"=Pyk2[20,11],"p2O"=1, + "p3A"=Pyk2[1,12], "p3C"=Pyk2[2,12], "p3D"=Pyk2[3,12], "p3E"=Pyk2[4,12], "p3F"=Pyk2[5,12],"p3G"=Pyk2[6,12],"p3H"=Pyk2[7,12],"p3I"=Pyk2[8,12],"p3K"=Pyk2[9,12],"p3L"=Pyk2[10,12],"p3M"=Pyk2[11,12],"p3N"=Pyk2[12,12],"p3P"=Pyk2[13,12],"p3Q"=Pyk2[14,12],"p3R"=Pyk2[15,12], + "p3S"=Pyk2[16,12],"p3T"=Pyk2[17,12],"p3V"=Pyk2[18,12],"p3W"=Pyk2[19,12],"p3Y"=Pyk2[20,12],"p3O"=1, + "p4A"=Pyk2[1,13], "p4C"=Pyk2[2,13], "p4D"=Pyk2[3,13], "p4E"=Pyk2[4,13], "p4F"=Pyk2[5,13],"p4G"=Pyk2[6,13],"p4H"=Pyk2[7,13],"p4I"=Pyk2[8,13],"p4K"=Pyk2[9,13],"p4L"=Pyk2[10,13],"p4M"=Pyk2[11,13],"p4N"=Pyk2[12,13],"p4P"=Pyk2[13,13],"p4Q"=Pyk2[14,13],"p4R"=Pyk2[15,13], + "p4S"=Pyk2[16,13],"p4T"=Pyk2[17,13],"p4V"=Pyk2[18,13],"p4W"=Pyk2[19,13],"p4Y"=Pyk2[20,13],"p4O"=1, + "p5A"=Pyk2[1,14], "p5C"=Pyk2[2,14], "p5D"=Pyk2[3,14], "p5E"=Pyk2[4,14], "p5F"=Pyk2[5,14],"p5G"=Pyk2[6,14],"p5H"=Pyk2[7,14],"p5I"=Pyk2[8,14],"p5K"=Pyk2[9,14],"p5L"=Pyk2[10,14],"p5M"=Pyk2[11,14],"p5N"=Pyk2[12,14],"p5P"=Pyk2[13,14],"p5Q"=Pyk2[14,14],"p5R"=Pyk2[15,14], + "p5S"=Pyk2[16,14],"p5T"=Pyk2[17,14],"p5V"=Pyk2[18,14],"p5W"=Pyk2[19,14],"p5Y"=Pyk2[20,14],"p5O"=1, + "p6A"=Pyk2[1,15], "p6C"=Pyk2[2,15], "p6D"=Pyk2[3,15], "p6E"=Pyk2[4,15], "p6F"=Pyk2[5,15],"p6G"=Pyk2[6,15],"p6H"=Pyk2[7,15],"p6I"=Pyk2[8,15],"p6K"=Pyk2[9,15],"p6L"=Pyk2[10,15],"p6M"=Pyk2[11,15],"p6N"=Pyk2[12,15],"p6P"=Pyk2[13,15],"p6Q"=Pyk2[14,15],"p6R"=Pyk2[15,15], + "p6S"=Pyk2[16,15],"p6T"=Pyk2[17,15],"p6V"=Pyk2[18,15],"p6W"=Pyk2[19,15],"p6Y"=Pyk2[20,15],"p6O"=1) + + Src_props <- c("m6A"=Src[1,3], "m6C"=Src[2,3], "m6D"=Src[3,3], "m6E"=Src[4,3], "m6F"=Src[5,3],"m6G"=Src[6,3],"m6H"=Src[7,3],"m6I"=Src[8,3],"m6K"=Src[9,3],"m6L"=Src[10,3],"m6M"=Src[11,3],"m6N"=Src[12,3],"m6P"=Src[13,3],"m6Q"=Src[14,3],"m6R"=Src[15,3], + "m6S"=Src[16,3],"m6T"=Src[17,3],"m6V"=Src[18,3],"m6W"=Src[19,3],"m6Y"=Src[20,3],"m6O"=1, + "m5A"=Src[1,4], "m5C"=Src[2,4], "m5D"=Src[3,4], "m5E"=Src[4,4], "m5F"=Src[5,4],"m5G"=Src[6,4],"m5H"=Src[7,4],"m5I"=Src[8,4],"m5K"=Src[9,4],"m5L"=Src[10,4],"m5M"=Src[11,4],"m5N"=Src[12,4],"m5P"=Src[13,4],"m5Q"=Src[14,4],"m5R"=Src[15,4], + "m5S"=Src[16,4],"m5T"=Src[17,4],"m5V"=Src[18,4],"m5W"=Src[19,4],"m5Y"=Src[20,4],"m5O"=1, + "m4A"=Src[1,5], "m4C"=Src[2,5], "m4D"=Src[3,5], "m4E"=Src[4,5], "m4F"=Src[5,5],"m4G"=Src[6,5],"m4H"=Src[7,5],"m4I"=Src[8,5],"m4K"=Src[9,5],"m4L"=Src[10,5],"m4M"=Src[11,5],"m4N"=Src[12,5],"m4P"=Src[13,5],"m4Q"=Src[14,5],"m4R"=Src[15,5], + "m4S"=Src[16,5],"m4T"=Src[17,5],"m4V"=Src[18,5],"m4W"=Src[19,5],"m4Y"=Src[20,5],"m4O"=1, + "m3A"=Src[1,6], "m3C"=Src[2,6], "m3D"=Src[3,6], "m3E"=Src[4,6], "m3F"=Src[5,6],"m3G"=Src[6,6],"m3H"=Src[7,6],"m3I"=Src[8,6],"m3K"=Src[9,6],"m3L"=Src[10,6],"m3M"=Src[11,6],"m3N"=Src[12,6],"m3P"=Src[13,6],"m3Q"=Src[14,6],"m3R"=Src[15,6], + "m3S"=Src[16,6],"m3T"=Src[17,6],"m3V"=Src[18,6],"m3W"=Src[19,6],"m3Y"=Src[20,6],"m3O"=1, + "m2A"=Src[1,7], "m2C"=Src[2,7], "m2D"=Src[3,7], "m2E"=Src[4,7], "m2F"=Src[5,7],"m2G"=Src[6,7],"m2H"=Src[7,7],"m2I"=Src[8,7],"m2K"=Src[9,7],"m2L"=Src[10,7],"m2M"=Src[11,7],"m2N"=Src[12,7],"m2P"=Src[13,7],"m2Q"=Src[14,7],"m2R"=Src[15,7], + "m2S"=Src[16,7],"m2T"=Src[17,7],"m2V"=Src[18,7],"m2W"=Src[19,7],"m2Y"=Src[20,7],"m2O"=1, + "m1A"=Src[1,8], "m1C"=Src[2,8], "m1D"=Src[3,8], "m1E"=Src[4,8], "m1F"=Src[5,8],"m1G"=Src[6,8],"m1H"=Src[7,8],"m1I"=Src[8,8],"m1K"=Src[9,8],"m1L"=Src[10,8],"m1M"=Src[11,8],"m1N"=Src[12,8],"m1P"=Src[13,8],"m1Q"=Src[14,8],"m1R"=Src[15,8], + "m1S"=Src[16,8],"m1T"=Src[17,8],"m1V"=Src[18,8],"m1W"=Src[19,8],"m1Y"=Src[20,8],"m1O"=1, + "d0A"=Src[1,9], "d0C"=Src[2,9], "d0D"=Src[3,9], "d0E"=Src[4,9], "d0F"=Src[5,9],"d0G"=Src[6,9],"d0H"=Src[7,9],"d0I"=Src[8,9],"d0K"=Src[9,9],"d0L"=Src[10,9],"d0M"=Src[11,9],"d0N"=Src[12,9],"d0P"=Src[13,9],"d0Q"=Src[14,9],"d0R"=Src[15,9], + "d0S"=Src[16,9],"d0T"=Src[17,9],"d0V"=Src[18,9],"d0W"=Src[19,9],"d0Y"=Src[20,9],"d0O"=1, + "p1A"=Src[1,10], "p1C"=Src[2,10], "p1D"=Src[3,10], "p1E"=Src[4,10], "p1F"=Src[5,10],"p1G"=Src[6,10],"p1H"=Src[7,10],"p1I"=Src[8,10],"p1K"=Src[9,10],"p1L"=Src[10,10],"p1M"=Src[11,10],"p1N"=Src[12,10],"p1P"=Src[13,10],"p1Q"=Src[14,10],"p1R"=Src[15,10], + "p1S"=Src[16,10],"p1T"=Src[17,10],"p1V"=Src[18,10],"p1W"=Src[19,10],"p1Y"=Src[20,10],"p1O"=1, + "p2A"=Src[1,11], "p2C"=Src[2,11], "p2D"=Src[3,11], "p2E"=Src[4,11], "p2F"=Src[5,11],"p2G"=Src[6,11],"p2H"=Src[7,11],"p2I"=Src[8,11],"p2K"=Src[9,11],"p2L"=Src[10,11],"p2M"=Src[11,11],"p2N"=Src[12,11],"p2P"=Src[13,11],"p2Q"=Src[14,11],"p2R"=Src[15,11], + "p2S"=Src[16,11],"p2T"=Src[17,11],"p2V"=Src[18,11],"p2W"=Src[19,11],"p2Y"=Src[20,11],"p2O"=1, + "p3A"=Src[1,12], "p3C"=Src[2,12], "p3D"=Src[3,12], "p3E"=Src[4,12], "p3F"=Src[5,12],"p3G"=Src[6,12],"p3H"=Src[7,12],"p3I"=Src[8,12],"p3K"=Src[9,12],"p3L"=Src[10,12],"p3M"=Src[11,12],"p3N"=Src[12,12],"p3P"=Src[13,12],"p3Q"=Src[14,12],"p3R"=Src[15,12], + "p3S"=Src[16,12],"p3T"=Src[17,12],"p3V"=Src[18,12],"p3W"=Src[19,12],"p3Y"=Src[20,12],"p3O"=1, + "p4A"=Src[1,13], "p4C"=Src[2,13], "p4D"=Src[3,13], "p4E"=Src[4,13], "p4F"=Src[5,13],"p4G"=Src[6,13],"p4H"=Src[7,13],"p4I"=Src[8,13],"p4K"=Src[9,13],"p4L"=Src[10,13],"p4M"=Src[11,13],"p4N"=Src[12,13],"p4P"=Src[13,13],"p4Q"=Src[14,13],"p4R"=Src[15,13], + "p4S"=Src[16,13],"p4T"=Src[17,13],"p4V"=Src[18,13],"p4W"=Src[19,13],"p4Y"=Src[20,13],"p4O"=1, + "p5A"=Src[1,14], "p5C"=Src[2,14], "p5D"=Src[3,14], "p5E"=Src[4,14], "p5F"=Src[5,14],"p5G"=Src[6,14],"p5H"=Src[7,14],"p5I"=Src[8,14],"p5K"=Src[9,14],"p5L"=Src[10,14],"p5M"=Src[11,14],"p5N"=Src[12,14],"p5P"=Src[13,14],"p5Q"=Src[14,14],"p5R"=Src[15,14], + "p5S"=Src[16,14],"p5T"=Src[17,14],"p5V"=Src[18,14],"p5W"=Src[19,14],"p5Y"=Src[20,14],"p5O"=1, + "p6A"=Src[1,15], "p6C"=Src[2,15], "p6D"=Src[3,15], "p6E"=Src[4,15], "p6F"=Src[5,15],"p6G"=Src[6,15],"p6H"=Src[7,15],"p6I"=Src[8,15],"p6K"=Src[9,15],"p6L"=Src[10,15],"p6M"=Src[11,15],"p6N"=Src[12,15],"p6P"=Src[13,15],"p6Q"=Src[14,15],"p6R"=Src[15,15], + "p6S"=Src[16,15],"p6T"=Src[17,15],"p6V"=Src[18,15],"p6W"=Src[19,15],"p6Y"=Src[20,15],"p6O"=1) + + Syk_props <- c("m6A"=Syk[1,3], "m6C"=Syk[2,3], "m6D"=Syk[3,3], "m6E"=Syk[4,3], "m6F"=Syk[5,3],"m6G"=Syk[6,3],"m6H"=Syk[7,3],"m6I"=Syk[8,3],"m6K"=Syk[9,3],"m6L"=Syk[10,3],"m6M"=Syk[11,3],"m6N"=Syk[12,3],"m6P"=Syk[13,3],"m6Q"=Syk[14,3],"m6R"=Syk[15,3], + "m6S"=Syk[16,3],"m6T"=Syk[17,3],"m6V"=Syk[18,3],"m6W"=Syk[19,3],"m6Y"=Syk[20,3],"m6O"=1, + "m5A"=Syk[1,4], "m5C"=Syk[2,4], "m5D"=Syk[3,4], "m5E"=Syk[4,4], "m5F"=Syk[5,4],"m5G"=Syk[6,4],"m5H"=Syk[7,4],"m5I"=Syk[8,4],"m5K"=Syk[9,4],"m5L"=Syk[10,4],"m5M"=Syk[11,4],"m5N"=Syk[12,4],"m5P"=Syk[13,4],"m5Q"=Syk[14,4],"m5R"=Syk[15,4], + "m5S"=Syk[16,4],"m5T"=Syk[17,4],"m5V"=Syk[18,4],"m5W"=Syk[19,4],"m5Y"=Syk[20,4],"m5O"=1, + "m4A"=Syk[1,5], "m4C"=Syk[2,5], "m4D"=Syk[3,5], "m4E"=Syk[4,5], "m4F"=Syk[5,5],"m4G"=Syk[6,5],"m4H"=Syk[7,5],"m4I"=Syk[8,5],"m4K"=Syk[9,5],"m4L"=Syk[10,5],"m4M"=Syk[11,5],"m4N"=Syk[12,5],"m4P"=Syk[13,5],"m4Q"=Syk[14,5],"m4R"=Syk[15,5], + "m4S"=Syk[16,5],"m4T"=Syk[17,5],"m4V"=Syk[18,5],"m4W"=Syk[19,5],"m4Y"=Syk[20,5],"m4O"=1, + "m3A"=Syk[1,6], "m3C"=Syk[2,6], "m3D"=Syk[3,6], "m3E"=Syk[4,6], "m3F"=Syk[5,6],"m3G"=Syk[6,6],"m3H"=Syk[7,6],"m3I"=Syk[8,6],"m3K"=Syk[9,6],"m3L"=Syk[10,6],"m3M"=Syk[11,6],"m3N"=Syk[12,6],"m3P"=Syk[13,6],"m3Q"=Syk[14,6],"m3R"=Syk[15,6], + "m3S"=Syk[16,6],"m3T"=Syk[17,6],"m3V"=Syk[18,6],"m3W"=Syk[19,6],"m3Y"=Syk[20,6],"m3O"=1, + "m2A"=Syk[1,7], "m2C"=Syk[2,7], "m2D"=Syk[3,7], "m2E"=Syk[4,7], "m2F"=Syk[5,7],"m2G"=Syk[6,7],"m2H"=Syk[7,7],"m2I"=Syk[8,7],"m2K"=Syk[9,7],"m2L"=Syk[10,7],"m2M"=Syk[11,7],"m2N"=Syk[12,7],"m2P"=Syk[13,7],"m2Q"=Syk[14,7],"m2R"=Syk[15,7], + "m2S"=Syk[16,7],"m2T"=Syk[17,7],"m2V"=Syk[18,7],"m2W"=Syk[19,7],"m2Y"=Syk[20,7],"m2O"=1, + "m1A"=Syk[1,8], "m1C"=Syk[2,8], "m1D"=Syk[3,8], "m1E"=Syk[4,8], "m1F"=Syk[5,8],"m1G"=Syk[6,8],"m1H"=Syk[7,8],"m1I"=Syk[8,8],"m1K"=Syk[9,8],"m1L"=Syk[10,8],"m1M"=Syk[11,8],"m1N"=Syk[12,8],"m1P"=Syk[13,8],"m1Q"=Syk[14,8],"m1R"=Syk[15,8], + "m1S"=Syk[16,8],"m1T"=Syk[17,8],"m1V"=Syk[18,8],"m1W"=Syk[19,8],"m1Y"=Syk[20,8],"m1O"=1, + "d0A"=Syk[1,9], "d0C"=Syk[2,9], "d0D"=Syk[3,9], "d0E"=Syk[4,9], "d0F"=Syk[5,9],"d0G"=Syk[6,9],"d0H"=Syk[7,9],"d0I"=Syk[8,9],"d0K"=Syk[9,9],"d0L"=Syk[10,9],"d0M"=Syk[11,9],"d0N"=Syk[12,9],"d0P"=Syk[13,9],"d0Q"=Syk[14,9],"d0R"=Syk[15,9], + "d0S"=Syk[16,9],"d0T"=Syk[17,9],"d0V"=Syk[18,9],"d0W"=Syk[19,9],"d0Y"=Syk[20,9],"d0O"=1, + "p1A"=Syk[1,10], "p1C"=Syk[2,10], "p1D"=Syk[3,10], "p1E"=Syk[4,10], "p1F"=Syk[5,10],"p1G"=Syk[6,10],"p1H"=Syk[7,10],"p1I"=Syk[8,10],"p1K"=Syk[9,10],"p1L"=Syk[10,10],"p1M"=Syk[11,10],"p1N"=Syk[12,10],"p1P"=Syk[13,10],"p1Q"=Syk[14,10],"p1R"=Syk[15,10], + "p1S"=Syk[16,10],"p1T"=Syk[17,10],"p1V"=Syk[18,10],"p1W"=Syk[19,10],"p1Y"=Syk[20,10],"p1O"=1, + "p2A"=Syk[1,11], "p2C"=Syk[2,11], "p2D"=Syk[3,11], "p2E"=Syk[4,11], "p2F"=Syk[5,11],"p2G"=Syk[6,11],"p2H"=Syk[7,11],"p2I"=Syk[8,11],"p2K"=Syk[9,11],"p2L"=Syk[10,11],"p2M"=Syk[11,11],"p2N"=Syk[12,11],"p2P"=Syk[13,11],"p2Q"=Syk[14,11],"p2R"=Syk[15,11], + "p2S"=Syk[16,11],"p2T"=Syk[17,11],"p2V"=Syk[18,11],"p2W"=Syk[19,11],"p2Y"=Syk[20,11],"p2O"=1, + "p3A"=Syk[1,12], "p3C"=Syk[2,12], "p3D"=Syk[3,12], "p3E"=Syk[4,12], "p3F"=Syk[5,12],"p3G"=Syk[6,12],"p3H"=Syk[7,12],"p3I"=Syk[8,12],"p3K"=Syk[9,12],"p3L"=Syk[10,12],"p3M"=Syk[11,12],"p3N"=Syk[12,12],"p3P"=Syk[13,12],"p3Q"=Syk[14,12],"p3R"=Syk[15,12], + "p3S"=Syk[16,12],"p3T"=Syk[17,12],"p3V"=Syk[18,12],"p3W"=Syk[19,12],"p3Y"=Syk[20,12],"p3O"=1, + "p4A"=Syk[1,13], "p4C"=Syk[2,13], "p4D"=Syk[3,13], "p4E"=Syk[4,13], "p4F"=Syk[5,13],"p4G"=Syk[6,13],"p4H"=Syk[7,13],"p4I"=Syk[8,13],"p4K"=Syk[9,13],"p4L"=Syk[10,13],"p4M"=Syk[11,13],"p4N"=Syk[12,13],"p4P"=Syk[13,13],"p4Q"=Syk[14,13],"p4R"=Syk[15,13], + "p4S"=Syk[16,13],"p4T"=Syk[17,13],"p4V"=Syk[18,13],"p4W"=Syk[19,13],"p4Y"=Syk[20,13],"p4O"=1, + "p5A"=Syk[1,14], "p5C"=Syk[2,14], "p5D"=Syk[3,14], "p5E"=Syk[4,14], "p5F"=Syk[5,14],"p5G"=Syk[6,14],"p5H"=Syk[7,14],"p5I"=Syk[8,14],"p5K"=Syk[9,14],"p5L"=Syk[10,14],"p5M"=Syk[11,14],"p5N"=Syk[12,14],"p5P"=Syk[13,14],"p5Q"=Syk[14,14],"p5R"=Syk[15,14], + "p5S"=Syk[16,14],"p5T"=Syk[17,14],"p5V"=Syk[18,14],"p5W"=Syk[19,14],"p5Y"=Syk[20,14],"p5O"=1, + "p6A"=Syk[1,15], "p6C"=Syk[2,15], "p6D"=Syk[3,15], "p6E"=Syk[4,15], "p6F"=Syk[5,15],"p6G"=Syk[6,15],"p6H"=Syk[7,15],"p6I"=Syk[8,15],"p6K"=Syk[9,15],"p6L"=Syk[10,15],"p6M"=Syk[11,15],"p6N"=Syk[12,15],"p6P"=Syk[13,15],"p6Q"=Syk[14,15],"p6R"=Syk[15,15], + "p6S"=Syk[16,15],"p6T"=Syk[17,15],"p6V"=Syk[18,15],"p6W"=Syk[19,15],"p6Y"=Syk[20,15],"p6O"=1) + + Yes_props <- c("m6A"=Yes[1,3], "m6C"=Yes[2,3], "m6D"=Yes[3,3], "m6E"=Yes[4,3], "m6F"=Yes[5,3],"m6G"=Yes[6,3],"m6H"=Yes[7,3],"m6I"=Yes[8,3],"m6K"=Yes[9,3],"m6L"=Yes[10,3],"m6M"=Yes[11,3],"m6N"=Yes[12,3],"m6P"=Yes[13,3],"m6Q"=Yes[14,3],"m6R"=Yes[15,3], + "m6S"=Yes[16,3],"m6T"=Yes[17,3],"m6V"=Yes[18,3],"m6W"=Yes[19,3],"m6Y"=Yes[20,3],"m6O"=1, + "m5A"=Yes[1,4], "m5C"=Yes[2,4], "m5D"=Yes[3,4], "m5E"=Yes[4,4], "m5F"=Yes[5,4],"m5G"=Yes[6,4],"m5H"=Yes[7,4],"m5I"=Yes[8,4],"m5K"=Yes[9,4],"m5L"=Yes[10,4],"m5M"=Yes[11,4],"m5N"=Yes[12,4],"m5P"=Yes[13,4],"m5Q"=Yes[14,4],"m5R"=Yes[15,4], + "m5S"=Yes[16,4],"m5T"=Yes[17,4],"m5V"=Yes[18,4],"m5W"=Yes[19,4],"m5Y"=Yes[20,4],"m5O"=1, + "m4A"=Yes[1,5], "m4C"=Yes[2,5], "m4D"=Yes[3,5], "m4E"=Yes[4,5], "m4F"=Yes[5,5],"m4G"=Yes[6,5],"m4H"=Yes[7,5],"m4I"=Yes[8,5],"m4K"=Yes[9,5],"m4L"=Yes[10,5],"m4M"=Yes[11,5],"m4N"=Yes[12,5],"m4P"=Yes[13,5],"m4Q"=Yes[14,5],"m4R"=Yes[15,5], + "m4S"=Yes[16,5],"m4T"=Yes[17,5],"m4V"=Yes[18,5],"m4W"=Yes[19,5],"m4Y"=Yes[20,5],"m4O"=1, + "m3A"=Yes[1,6], "m3C"=Yes[2,6], "m3D"=Yes[3,6], "m3E"=Yes[4,6], "m3F"=Yes[5,6],"m3G"=Yes[6,6],"m3H"=Yes[7,6],"m3I"=Yes[8,6],"m3K"=Yes[9,6],"m3L"=Yes[10,6],"m3M"=Yes[11,6],"m3N"=Yes[12,6],"m3P"=Yes[13,6],"m3Q"=Yes[14,6],"m3R"=Yes[15,6], + "m3S"=Yes[16,6],"m3T"=Yes[17,6],"m3V"=Yes[18,6],"m3W"=Yes[19,6],"m3Y"=Yes[20,6],"m3O"=1, + "m2A"=Yes[1,7], "m2C"=Yes[2,7], "m2D"=Yes[3,7], "m2E"=Yes[4,7], "m2F"=Yes[5,7],"m2G"=Yes[6,7],"m2H"=Yes[7,7],"m2I"=Yes[8,7],"m2K"=Yes[9,7],"m2L"=Yes[10,7],"m2M"=Yes[11,7],"m2N"=Yes[12,7],"m2P"=Yes[13,7],"m2Q"=Yes[14,7],"m2R"=Yes[15,7], + "m2S"=Yes[16,7],"m2T"=Yes[17,7],"m2V"=Yes[18,7],"m2W"=Yes[19,7],"m2Y"=Yes[20,7],"m2O"=1, + "m1A"=Yes[1,8], "m1C"=Yes[2,8], "m1D"=Yes[3,8], "m1E"=Yes[4,8], "m1F"=Yes[5,8],"m1G"=Yes[6,8],"m1H"=Yes[7,8],"m1I"=Yes[8,8],"m1K"=Yes[9,8],"m1L"=Yes[10,8],"m1M"=Yes[11,8],"m1N"=Yes[12,8],"m1P"=Yes[13,8],"m1Q"=Yes[14,8],"m1R"=Yes[15,8], + "m1S"=Yes[16,8],"m1T"=Yes[17,8],"m1V"=Yes[18,8],"m1W"=Yes[19,8],"m1Y"=Yes[20,8],"m1O"=1, + "d0A"=Yes[1,9], "d0C"=Yes[2,9], "d0D"=Yes[3,9], "d0E"=Yes[4,9], "d0F"=Yes[5,9],"d0G"=Yes[6,9],"d0H"=Yes[7,9],"d0I"=Yes[8,9],"d0K"=Yes[9,9],"d0L"=Yes[10,9],"d0M"=Yes[11,9],"d0N"=Yes[12,9],"d0P"=Yes[13,9],"d0Q"=Yes[14,9],"d0R"=Yes[15,9], + "d0S"=Yes[16,9],"d0T"=Yes[17,9],"d0V"=Yes[18,9],"d0W"=Yes[19,9],"d0Y"=Yes[20,9],"d0O"=1, + "p1A"=Yes[1,10], "p1C"=Yes[2,10], "p1D"=Yes[3,10], "p1E"=Yes[4,10], "p1F"=Yes[5,10],"p1G"=Yes[6,10],"p1H"=Yes[7,10],"p1I"=Yes[8,10],"p1K"=Yes[9,10],"p1L"=Yes[10,10],"p1M"=Yes[11,10],"p1N"=Yes[12,10],"p1P"=Yes[13,10],"p1Q"=Yes[14,10],"p1R"=Yes[15,10], + "p1S"=Yes[16,10],"p1T"=Yes[17,10],"p1V"=Yes[18,10],"p1W"=Yes[19,10],"p1Y"=Yes[20,10],"p1O"=1, + "p2A"=Yes[1,11], "p2C"=Yes[2,11], "p2D"=Yes[3,11], "p2E"=Yes[4,11], "p2F"=Yes[5,11],"p2G"=Yes[6,11],"p2H"=Yes[7,11],"p2I"=Yes[8,11],"p2K"=Yes[9,11],"p2L"=Yes[10,11],"p2M"=Yes[11,11],"p2N"=Yes[12,11],"p2P"=Yes[13,11],"p2Q"=Yes[14,11],"p2R"=Yes[15,11], + "p2S"=Yes[16,11],"p2T"=Yes[17,11],"p2V"=Yes[18,11],"p2W"=Yes[19,11],"p2Y"=Yes[20,11],"p2O"=1, + "p3A"=Yes[1,12], "p3C"=Yes[2,12], "p3D"=Yes[3,12], "p3E"=Yes[4,12], "p3F"=Yes[5,12],"p3G"=Yes[6,12],"p3H"=Yes[7,12],"p3I"=Yes[8,12],"p3K"=Yes[9,12],"p3L"=Yes[10,12],"p3M"=Yes[11,12],"p3N"=Yes[12,12],"p3P"=Yes[13,12],"p3Q"=Yes[14,12],"p3R"=Yes[15,12], + "p3S"=Yes[16,12],"p3T"=Yes[17,12],"p3V"=Yes[18,12],"p3W"=Yes[19,12],"p3Y"=Yes[20,12],"p3O"=1, + "p4A"=Yes[1,13], "p4C"=Yes[2,13], "p4D"=Yes[3,13], "p4E"=Yes[4,13], "p4F"=Yes[5,13],"p4G"=Yes[6,13],"p4H"=Yes[7,13],"p4I"=Yes[8,13],"p4K"=Yes[9,13],"p4L"=Yes[10,13],"p4M"=Yes[11,13],"p4N"=Yes[12,13],"p4P"=Yes[13,13],"p4Q"=Yes[14,13],"p4R"=Yes[15,13], + "p4S"=Yes[16,13],"p4T"=Yes[17,13],"p4V"=Yes[18,13],"p4W"=Yes[19,13],"p4Y"=Yes[20,13],"p4O"=1, + "p5A"=Yes[1,14], "p5C"=Yes[2,14], "p5D"=Yes[3,14], "p5E"=Yes[4,14], "p5F"=Yes[5,14],"p5G"=Yes[6,14],"p5H"=Yes[7,14],"p5I"=Yes[8,14],"p5K"=Yes[9,14],"p5L"=Yes[10,14],"p5M"=Yes[11,14],"p5N"=Yes[12,14],"p5P"=Yes[13,14],"p5Q"=Yes[14,14],"p5R"=Yes[15,14], + "p5S"=Yes[16,14],"p5T"=Yes[17,14],"p5V"=Yes[18,14],"p5W"=Yes[19,14],"p5Y"=Yes[20,14],"p5O"=1, + "p6A"=Yes[1,15], "p6C"=Yes[2,15], "p6D"=Yes[3,15], "p6E"=Yes[4,15], "p6F"=Yes[5,15],"p6G"=Yes[6,15],"p6H"=Yes[7,15],"p6I"=Yes[8,15],"p6K"=Yes[9,15],"p6L"=Yes[10,15],"p6M"=Yes[11,15],"p6N"=Yes[12,15],"p6P"=Yes[13,15],"p6Q"=Yes[14,15],"p6R"=Yes[15,15], + "p6S"=Yes[16,15],"p6T"=Yes[17,15],"p6V"=Yes[18,15],"p6W"=Yes[19,15],"p6Y"=Yes[20,15],"p6O"=1) + + FLT3_props <- c("m6A"=FLT3[1,3], "m6C"=FLT3[2,3], "m6D"=FLT3[3,3], "m6E"=FLT3[4,3], "m6F"=FLT3[5,3],"m6G"=FLT3[6,3],"m6H"=FLT3[7,3],"m6I"=FLT3[8,3],"m6K"=FLT3[9,3],"m6L"=FLT3[10,3],"m6M"=FLT3[11,3],"m6N"=FLT3[12,3],"m6P"=FLT3[13,3],"m6Q"=FLT3[14,3],"m6R"=FLT3[15,3], + "m6S"=FLT3[16,3],"m6T"=FLT3[17,3],"m6V"=FLT3[18,3],"m6W"=FLT3[19,3],"m6Y"=FLT3[20,3],"m6O"=1, + "m5A"=FLT3[1,4], "m5C"=FLT3[2,4], "m5D"=FLT3[3,4], "m5E"=FLT3[4,4], "m5F"=FLT3[5,4],"m5G"=FLT3[6,4],"m5H"=FLT3[7,4],"m5I"=FLT3[8,4],"m5K"=FLT3[9,4],"m5L"=FLT3[10,4],"m5M"=FLT3[11,4],"m5N"=FLT3[12,4],"m5P"=FLT3[13,4],"m5Q"=FLT3[14,4],"m5R"=FLT3[15,4], + "m5S"=FLT3[16,4],"m5T"=FLT3[17,4],"m5V"=FLT3[18,4],"m5W"=FLT3[19,4],"m5Y"=FLT3[20,4],"m5O"=1, + "m4A"=FLT3[1,5], "m4C"=FLT3[2,5], "m4D"=FLT3[3,5], "m4E"=FLT3[4,5], "m4F"=FLT3[5,5],"m4G"=FLT3[6,5],"m4H"=FLT3[7,5],"m4I"=FLT3[8,5],"m4K"=FLT3[9,5],"m4L"=FLT3[10,5],"m4M"=FLT3[11,5],"m4N"=FLT3[12,5],"m4P"=FLT3[13,5],"m4Q"=FLT3[14,5],"m4R"=FLT3[15,5], + "m4S"=FLT3[16,5],"m4T"=FLT3[17,5],"m4V"=FLT3[18,5],"m4W"=FLT3[19,5],"m4Y"=FLT3[20,5],"m4O"=1, + "m3A"=FLT3[1,6], "m3C"=FLT3[2,6], "m3D"=FLT3[3,6], "m3E"=FLT3[4,6], "m3F"=FLT3[5,6],"m3G"=FLT3[6,6],"m3H"=FLT3[7,6],"m3I"=FLT3[8,6],"m3K"=FLT3[9,6],"m3L"=FLT3[10,6],"m3M"=FLT3[11,6],"m3N"=FLT3[12,6],"m3P"=FLT3[13,6],"m3Q"=FLT3[14,6],"m3R"=FLT3[15,6], + "m3S"=FLT3[16,6],"m3T"=FLT3[17,6],"m3V"=FLT3[18,6],"m3W"=FLT3[19,6],"m3Y"=FLT3[20,6],"m3O"=1, + "m2A"=FLT3[1,7], "m2C"=FLT3[2,7], "m2D"=FLT3[3,7], "m2E"=FLT3[4,7], "m2F"=FLT3[5,7],"m2G"=FLT3[6,7],"m2H"=FLT3[7,7],"m2I"=FLT3[8,7],"m2K"=FLT3[9,7],"m2L"=FLT3[10,7],"m2M"=FLT3[11,7],"m2N"=FLT3[12,7],"m2P"=FLT3[13,7],"m2Q"=FLT3[14,7],"m2R"=FLT3[15,7], + "m2S"=FLT3[16,7],"m2T"=FLT3[17,7],"m2V"=FLT3[18,7],"m2W"=FLT3[19,7],"m2Y"=FLT3[20,7],"m2O"=1, + "m1A"=FLT3[1,8], "m1C"=FLT3[2,8], "m1D"=FLT3[3,8], "m1E"=FLT3[4,8], "m1F"=FLT3[5,8],"m1G"=FLT3[6,8],"m1H"=FLT3[7,8],"m1I"=FLT3[8,8],"m1K"=FLT3[9,8],"m1L"=FLT3[10,8],"m1M"=FLT3[11,8],"m1N"=FLT3[12,8],"m1P"=FLT3[13,8],"m1Q"=FLT3[14,8],"m1R"=FLT3[15,8], + "m1S"=FLT3[16,8],"m1T"=FLT3[17,8],"m1V"=FLT3[18,8],"m1W"=FLT3[19,8],"m1Y"=FLT3[20,8],"m1O"=1, + "d0A"=FLT3[1,9], "d0C"=FLT3[2,9], "d0D"=FLT3[3,9], "d0E"=FLT3[4,9], "d0F"=FLT3[5,9],"d0G"=FLT3[6,9],"d0H"=FLT3[7,9],"d0I"=FLT3[8,9],"d0K"=FLT3[9,9],"d0L"=FLT3[10,9],"d0M"=FLT3[11,9],"d0N"=FLT3[12,9],"d0P"=FLT3[13,9],"d0Q"=FLT3[14,9],"d0R"=FLT3[15,9], + "d0S"=FLT3[16,9],"d0T"=FLT3[17,9],"d0V"=FLT3[18,9],"d0W"=FLT3[19,9],"d0Y"=FLT3[20,9],"d0O"=1, + "p1A"=FLT3[1,10], "p1C"=FLT3[2,10], "p1D"=FLT3[3,10], "p1E"=FLT3[4,10], "p1F"=FLT3[5,10],"p1G"=FLT3[6,10],"p1H"=FLT3[7,10],"p1I"=FLT3[8,10],"p1K"=FLT3[9,10],"p1L"=FLT3[10,10],"p1M"=FLT3[11,10],"p1N"=FLT3[12,10],"p1P"=FLT3[13,10],"p1Q"=FLT3[14,10],"p1R"=FLT3[15,10], + "p1S"=FLT3[16,10],"p1T"=FLT3[17,10],"p1V"=FLT3[18,10],"p1W"=FLT3[19,10],"p1Y"=FLT3[20,10],"p1O"=1, + "p2A"=FLT3[1,11], "p2C"=FLT3[2,11], "p2D"=FLT3[3,11], "p2E"=FLT3[4,11], "p2F"=FLT3[5,11],"p2G"=FLT3[6,11],"p2H"=FLT3[7,11],"p2I"=FLT3[8,11],"p2K"=FLT3[9,11],"p2L"=FLT3[10,11],"p2M"=FLT3[11,11],"p2N"=FLT3[12,11],"p2P"=FLT3[13,11],"p2Q"=FLT3[14,11],"p2R"=FLT3[15,11], + "p2S"=FLT3[16,11],"p2T"=FLT3[17,11],"p2V"=FLT3[18,11],"p2W"=FLT3[19,11],"p2Y"=FLT3[20,11],"p2O"=1, + "p3A"=FLT3[1,12], "p3C"=FLT3[2,12], "p3D"=FLT3[3,12], "p3E"=FLT3[4,12], "p3F"=FLT3[5,12],"p3G"=FLT3[6,12],"p3H"=FLT3[7,12],"p3I"=FLT3[8,12],"p3K"=FLT3[9,12],"p3L"=FLT3[10,12],"p3M"=FLT3[11,12],"p3N"=FLT3[12,12],"p3P"=FLT3[13,12],"p3Q"=FLT3[14,12],"p3R"=FLT3[15,12], + "p3S"=FLT3[16,12],"p3T"=FLT3[17,12],"p3V"=FLT3[18,12],"p3W"=FLT3[19,12],"p3Y"=FLT3[20,12],"p3O"=1, + "p4A"=FLT3[1,13], "p4C"=FLT3[2,13], "p4D"=FLT3[3,13], "p4E"=FLT3[4,13], "p4F"=FLT3[5,13],"p4G"=FLT3[6,13],"p4H"=FLT3[7,13],"p4I"=FLT3[8,13],"p4K"=FLT3[9,13],"p4L"=FLT3[10,13],"p4M"=FLT3[11,13],"p4N"=FLT3[12,13],"p4P"=FLT3[13,13],"p4Q"=FLT3[14,13],"p4R"=FLT3[15,13], + "p4S"=FLT3[16,13],"p4T"=FLT3[17,13],"p4V"=FLT3[18,13],"p4W"=FLT3[19,13],"p4Y"=FLT3[20,13],"p4O"=1, + "p5A"=FLT3[1,14], "p5C"=FLT3[2,14], "p5D"=FLT3[3,14], "p5E"=FLT3[4,14], "p5F"=FLT3[5,14],"p5G"=FLT3[6,14],"p5H"=FLT3[7,14],"p5I"=FLT3[8,14],"p5K"=FLT3[9,14],"p5L"=FLT3[10,14],"p5M"=FLT3[11,14],"p5N"=FLT3[12,14],"p5P"=FLT3[13,14],"p5Q"=FLT3[14,14],"p5R"=FLT3[15,14], + "p5S"=FLT3[16,14],"p5T"=FLT3[17,14],"p5V"=FLT3[18,14],"p5W"=FLT3[19,14],"p5Y"=FLT3[20,14],"p5O"=1, + "p6A"=FLT3[1,15], "p6C"=FLT3[2,15], "p6D"=FLT3[3,15], "p6E"=FLT3[4,15], "p6F"=FLT3[5,15],"p6G"=FLT3[6,15],"p6H"=FLT3[7,15],"p6I"=FLT3[8,15],"p6K"=FLT3[9,15],"p6L"=FLT3[10,15],"p6M"=FLT3[11,15],"p6N"=FLT3[12,15],"p6P"=FLT3[13,15],"p6Q"=FLT3[14,15],"p6R"=FLT3[15,15], + "p6S"=FLT3[16,15],"p6T"=FLT3[17,15],"p6V"=FLT3[18,15],"p6W"=FLT3[19,15],"p6Y"=FLT3[20,15],"p6O"=1) + + ALK_props <- c("m6A"=ALK[1,3], "m6C"=ALK[2,3], "m6D"=ALK[3,3], "m6E"=ALK[4,3], "m6F"=ALK[5,3],"m6G"=ALK[6,3],"m6H"=ALK[7,3],"m6I"=ALK[8,3],"m6K"=ALK[9,3],"m6L"=ALK[10,3],"m6M"=ALK[11,3],"m6N"=ALK[12,3],"m6P"=ALK[13,3],"m6Q"=ALK[14,3],"m6R"=ALK[15,3], + "m6S"=ALK[16,3],"m6T"=ALK[17,3],"m6V"=ALK[18,3],"m6W"=ALK[19,3],"m6Y"=ALK[20,3],"m6O"=1, + "m5A"=ALK[1,4], "m5C"=ALK[2,4], "m5D"=ALK[3,4], "m5E"=ALK[4,4], "m5F"=ALK[5,4],"m5G"=ALK[6,4],"m5H"=ALK[7,4],"m5I"=ALK[8,4],"m5K"=ALK[9,4],"m5L"=ALK[10,4],"m5M"=ALK[11,4],"m5N"=ALK[12,4],"m5P"=ALK[13,4],"m5Q"=ALK[14,4],"m5R"=ALK[15,4], + "m5S"=ALK[16,4],"m5T"=ALK[17,4],"m5V"=ALK[18,4],"m5W"=ALK[19,4],"m5Y"=ALK[20,4],"m5O"=1, + "m4A"=ALK[1,5], "m4C"=ALK[2,5], "m4D"=ALK[3,5], "m4E"=ALK[4,5], "m4F"=ALK[5,5],"m4G"=ALK[6,5],"m4H"=ALK[7,5],"m4I"=ALK[8,5],"m4K"=ALK[9,5],"m4L"=ALK[10,5],"m4M"=ALK[11,5],"m4N"=ALK[12,5],"m4P"=ALK[13,5],"m4Q"=ALK[14,5],"m4R"=ALK[15,5], + "m4S"=ALK[16,5],"m4T"=ALK[17,5],"m4V"=ALK[18,5],"m4W"=ALK[19,5],"m4Y"=ALK[20,5],"m4O"=1, + "m3A"=ALK[1,6], "m3C"=ALK[2,6], "m3D"=ALK[3,6], "m3E"=ALK[4,6], "m3F"=ALK[5,6],"m3G"=ALK[6,6],"m3H"=ALK[7,6],"m3I"=ALK[8,6],"m3K"=ALK[9,6],"m3L"=ALK[10,6],"m3M"=ALK[11,6],"m3N"=ALK[12,6],"m3P"=ALK[13,6],"m3Q"=ALK[14,6],"m3R"=ALK[15,6], + "m3S"=ALK[16,6],"m3T"=ALK[17,6],"m3V"=ALK[18,6],"m3W"=ALK[19,6],"m3Y"=ALK[20,6],"m3O"=1, + "m2A"=ALK[1,7], "m2C"=ALK[2,7], "m2D"=ALK[3,7], "m2E"=ALK[4,7], "m2F"=ALK[5,7],"m2G"=ALK[6,7],"m2H"=ALK[7,7],"m2I"=ALK[8,7],"m2K"=ALK[9,7],"m2L"=ALK[10,7],"m2M"=ALK[11,7],"m2N"=ALK[12,7],"m2P"=ALK[13,7],"m2Q"=ALK[14,7],"m2R"=ALK[15,7], + "m2S"=ALK[16,7],"m2T"=ALK[17,7],"m2V"=ALK[18,7],"m2W"=ALK[19,7],"m2Y"=ALK[20,7],"m2O"=1, + "m1A"=ALK[1,8], "m1C"=ALK[2,8], "m1D"=ALK[3,8], "m1E"=ALK[4,8], "m1F"=ALK[5,8],"m1G"=ALK[6,8],"m1H"=ALK[7,8],"m1I"=ALK[8,8],"m1K"=ALK[9,8],"m1L"=ALK[10,8],"m1M"=ALK[11,8],"m1N"=ALK[12,8],"m1P"=ALK[13,8],"m1Q"=ALK[14,8],"m1R"=ALK[15,8], + "m1S"=ALK[16,8],"m1T"=ALK[17,8],"m1V"=ALK[18,8],"m1W"=ALK[19,8],"m1Y"=ALK[20,8],"m1O"=1, + "d0A"=ALK[1,9], "d0C"=ALK[2,9], "d0D"=ALK[3,9], "d0E"=ALK[4,9], "d0F"=ALK[5,9],"d0G"=ALK[6,9],"d0H"=ALK[7,9],"d0I"=ALK[8,9],"d0K"=ALK[9,9],"d0L"=ALK[10,9],"d0M"=ALK[11,9],"d0N"=ALK[12,9],"d0P"=ALK[13,9],"d0Q"=ALK[14,9],"d0R"=ALK[15,9], + "d0S"=ALK[16,9],"d0T"=ALK[17,9],"d0V"=ALK[18,9],"d0W"=ALK[19,9],"d0Y"=ALK[20,9],"d0O"=1, + "p1A"=ALK[1,10], "p1C"=ALK[2,10], "p1D"=ALK[3,10], "p1E"=ALK[4,10], "p1F"=ALK[5,10],"p1G"=ALK[6,10],"p1H"=ALK[7,10],"p1I"=ALK[8,10],"p1K"=ALK[9,10],"p1L"=ALK[10,10],"p1M"=ALK[11,10],"p1N"=ALK[12,10],"p1P"=ALK[13,10],"p1Q"=ALK[14,10],"p1R"=ALK[15,10], + "p1S"=ALK[16,10],"p1T"=ALK[17,10],"p1V"=ALK[18,10],"p1W"=ALK[19,10],"p1Y"=ALK[20,10],"p1O"=1, + "p2A"=ALK[1,11], "p2C"=ALK[2,11], "p2D"=ALK[3,11], "p2E"=ALK[4,11], "p2F"=ALK[5,11],"p2G"=ALK[6,11],"p2H"=ALK[7,11],"p2I"=ALK[8,11],"p2K"=ALK[9,11],"p2L"=ALK[10,11],"p2M"=ALK[11,11],"p2N"=ALK[12,11],"p2P"=ALK[13,11],"p2Q"=ALK[14,11],"p2R"=ALK[15,11], + "p2S"=ALK[16,11],"p2T"=ALK[17,11],"p2V"=ALK[18,11],"p2W"=ALK[19,11],"p2Y"=ALK[20,11],"p2O"=1, + "p3A"=ALK[1,12], "p3C"=ALK[2,12], "p3D"=ALK[3,12], "p3E"=ALK[4,12], "p3F"=ALK[5,12],"p3G"=ALK[6,12],"p3H"=ALK[7,12],"p3I"=ALK[8,12],"p3K"=ALK[9,12],"p3L"=ALK[10,12],"p3M"=ALK[11,12],"p3N"=ALK[12,12],"p3P"=ALK[13,12],"p3Q"=ALK[14,12],"p3R"=ALK[15,12], + "p3S"=ALK[16,12],"p3T"=ALK[17,12],"p3V"=ALK[18,12],"p3W"=ALK[19,12],"p3Y"=ALK[20,12],"p3O"=1, + "p4A"=ALK[1,13], "p4C"=ALK[2,13], "p4D"=ALK[3,13], "p4E"=ALK[4,13], "p4F"=ALK[5,13],"p4G"=ALK[6,13],"p4H"=ALK[7,13],"p4I"=ALK[8,13],"p4K"=ALK[9,13],"p4L"=ALK[10,13],"p4M"=ALK[11,13],"p4N"=ALK[12,13],"p4P"=ALK[13,13],"p4Q"=ALK[14,13],"p4R"=ALK[15,13], + "p4S"=ALK[16,13],"p4T"=ALK[17,13],"p4V"=ALK[18,13],"p4W"=ALK[19,13],"p4Y"=ALK[20,13],"p4O"=1, + "p5A"=ALK[1,14], "p5C"=ALK[2,14], "p5D"=ALK[3,14], "p5E"=ALK[4,14], "p5F"=ALK[5,14],"p5G"=ALK[6,14],"p5H"=ALK[7,14],"p5I"=ALK[8,14],"p5K"=ALK[9,14],"p5L"=ALK[10,14],"p5M"=ALK[11,14],"p5N"=ALK[12,14],"p5P"=ALK[13,14],"p5Q"=ALK[14,14],"p5R"=ALK[15,14], + "p5S"=ALK[16,14],"p5T"=ALK[17,14],"p5V"=ALK[18,14],"p5W"=ALK[19,14],"p5Y"=ALK[20,14],"p5O"=1, + "p6A"=ALK[1,15], "p6C"=ALK[2,15], "p6D"=ALK[3,15], "p6E"=ALK[4,15], "p6F"=ALK[5,15],"p6G"=ALK[6,15],"p6H"=ALK[7,15],"p6I"=ALK[8,15],"p6K"=ALK[9,15],"p6L"=ALK[10,15],"p6M"=ALK[11,15],"p6N"=ALK[12,15],"p6P"=ALK[13,15],"p6Q"=ALK[14,15],"p6R"=ALK[15,15], + "p6S"=ALK[16,15],"p6T"=ALK[17,15],"p6V"=ALK[18,15],"p6W"=ALK[19,15],"p6Y"=ALK[20,15],"p6O"=1) + +} +#this creates the apply functions that are used below. Each amino acid at each position gets assigned a value from an EPM table + + +#this part turns amino acids at positions into values from the EPM tables using the above apply function +if(1==1){ + ThisKinTablem6<-sapply(pPositionm6, function(x) ThisKinTable_props[x]) + ThisKinTablem5<-sapply(pPositionm5, function(x) ThisKinTable_props[x]) + ThisKinTablem4<-sapply(pPositionm4, function(x) ThisKinTable_props[x]) + ThisKinTablem3<-sapply(pPositionm3, function(x) ThisKinTable_props[x]) + ThisKinTablem2<-sapply(pPositionm2, function(x) ThisKinTable_props[x]) + ThisKinTablem1<-sapply(pPositionm1, function(x) ThisKinTable_props[x]) + ThisKinTabled0<-sapply(pPositiond0, function(x) ThisKinTable_props[x]) + ThisKinTablep1<-sapply(pPositionp1, function(x) ThisKinTable_props[x]) + ThisKinTablep2<-sapply(pPositionp2, function(x) ThisKinTable_props[x]) + ThisKinTablep3<-sapply(pPositionp3, function(x) ThisKinTable_props[x]) + ThisKinTablep4<-sapply(pPositionp4, function(x) ThisKinTable_props[x]) + ThisKinTablep5<-sapply(pPositionp5, function(x) ThisKinTable_props[x]) + ThisKinTablep6<-sapply(pPositionp6, function(x) ThisKinTable_props[x]) + + Ablm6<-sapply(pPositionm6, function(x) Abl_props[x]) + Ablm5<-sapply(pPositionm5, function(x) Abl_props[x]) + Ablm4<-sapply(pPositionm4, function(x) Abl_props[x]) + Ablm3<-sapply(pPositionm3, function(x) Abl_props[x]) + Ablm2<-sapply(pPositionm2, function(x) Abl_props[x]) + Ablm1<-sapply(pPositionm1, function(x) Abl_props[x]) + Abld0<-sapply(pPositiond0, function(x) Abl_props[x]) + Ablp1<-sapply(pPositionp1, function(x) Abl_props[x]) + Ablp2<-sapply(pPositionp2, function(x) Abl_props[x]) + Ablp3<-sapply(pPositionp3, function(x) Abl_props[x]) + Ablp4<-sapply(pPositionp4, function(x) Abl_props[x]) + Ablp5<-sapply(pPositionp5, function(x) Abl_props[x]) + Ablp6<-sapply(pPositionp6, function(x) Abl_props[x]) + + Argm6<-sapply(pPositionm6, function(x) Arg_props[x]) + Argm5<-sapply(pPositionm5, function(x) Arg_props[x]) + Argm4<-sapply(pPositionm4, function(x) Arg_props[x]) + Argm3<-sapply(pPositionm3, function(x) Arg_props[x]) + Argm2<-sapply(pPositionm2, function(x) Arg_props[x]) + Argm1<-sapply(pPositionm1, function(x) Arg_props[x]) + Argd0<-sapply(pPositiond0, function(x) Arg_props[x]) + Argp1<-sapply(pPositionp1, function(x) Arg_props[x]) + Argp2<-sapply(pPositionp2, function(x) Arg_props[x]) + Argp3<-sapply(pPositionp3, function(x) Arg_props[x]) + Argp4<-sapply(pPositionp4, function(x) Arg_props[x]) + Argp5<-sapply(pPositionp5, function(x) Arg_props[x]) + Argp6<-sapply(pPositionp6, function(x) Arg_props[x]) + + Btkm6<-sapply(pPositionm6, function(x) Btk_props[x]) + Btkm5<-sapply(pPositionm5, function(x) Btk_props[x]) + Btkm4<-sapply(pPositionm4, function(x) Btk_props[x]) + Btkm3<-sapply(pPositionm3, function(x) Btk_props[x]) + Btkm2<-sapply(pPositionm2, function(x) Btk_props[x]) + Btkm1<-sapply(pPositionm1, function(x) Btk_props[x]) + Btkd0<-sapply(pPositiond0, function(x) Btk_props[x]) + Btkp1<-sapply(pPositionp1, function(x) Btk_props[x]) + Btkp2<-sapply(pPositionp2, function(x) Btk_props[x]) + Btkp3<-sapply(pPositionp3, function(x) Btk_props[x]) + Btkp4<-sapply(pPositionp4, function(x) Btk_props[x]) + Btkp5<-sapply(pPositionp5, function(x) Btk_props[x]) + Btkp6<-sapply(pPositionp6, function(x) Btk_props[x]) + + Cskm6<-sapply(pPositionm6, function(x) Csk_props[x]) + Cskm5<-sapply(pPositionm5, function(x) Csk_props[x]) + Cskm4<-sapply(pPositionm4, function(x) Csk_props[x]) + Cskm3<-sapply(pPositionm3, function(x) Csk_props[x]) + Cskm2<-sapply(pPositionm2, function(x) Csk_props[x]) + Cskm1<-sapply(pPositionm1, function(x) Csk_props[x]) + Cskd0<-sapply(pPositiond0, function(x) Csk_props[x]) + Cskp1<-sapply(pPositionp1, function(x) Csk_props[x]) + Cskp2<-sapply(pPositionp2, function(x) Csk_props[x]) + Cskp3<-sapply(pPositionp3, function(x) Csk_props[x]) + Cskp4<-sapply(pPositionp4, function(x) Csk_props[x]) + Cskp5<-sapply(pPositionp5, function(x) Csk_props[x]) + Cskp6<-sapply(pPositionp6, function(x) Csk_props[x]) + + Fynm6<-sapply(pPositionm6, function(x) Fyn_props[x]) + Fynm5<-sapply(pPositionm5, function(x) Fyn_props[x]) + Fynm4<-sapply(pPositionm4, function(x) Fyn_props[x]) + Fynm3<-sapply(pPositionm3, function(x) Fyn_props[x]) + Fynm2<-sapply(pPositionm2, function(x) Fyn_props[x]) + Fynm1<-sapply(pPositionm1, function(x) Fyn_props[x]) + Fynd0<-sapply(pPositiond0, function(x) Fyn_props[x]) + Fynp1<-sapply(pPositionp1, function(x) Fyn_props[x]) + Fynp2<-sapply(pPositionp2, function(x) Fyn_props[x]) + Fynp3<-sapply(pPositionp3, function(x) Fyn_props[x]) + Fynp4<-sapply(pPositionp4, function(x) Fyn_props[x]) + Fynp5<-sapply(pPositionp5, function(x) Fyn_props[x]) + Fynp6<-sapply(pPositionp6, function(x) Fyn_props[x]) + + Hckm6<-sapply(pPositionm6, function(x) Hck_props[x]) + Hckm5<-sapply(pPositionm5, function(x) Hck_props[x]) + Hckm4<-sapply(pPositionm4, function(x) Hck_props[x]) + Hckm3<-sapply(pPositionm3, function(x) Hck_props[x]) + Hckm2<-sapply(pPositionm2, function(x) Hck_props[x]) + Hckm1<-sapply(pPositionm1, function(x) Hck_props[x]) + Hckd0<-sapply(pPositiond0, function(x) Hck_props[x]) + Hckp1<-sapply(pPositionp1, function(x) Hck_props[x]) + Hckp2<-sapply(pPositionp2, function(x) Hck_props[x]) + Hckp3<-sapply(pPositionp3, function(x) Hck_props[x]) + Hckp4<-sapply(pPositionp4, function(x) Hck_props[x]) + Hckp5<-sapply(pPositionp5, function(x) Hck_props[x]) + Hckp6<-sapply(pPositionp6, function(x) Hck_props[x]) + + JAK2m6<-sapply(pPositionm6, function(x) JAK2_props[x]) + JAK2m5<-sapply(pPositionm5, function(x) JAK2_props[x]) + JAK2m4<-sapply(pPositionm4, function(x) JAK2_props[x]) + JAK2m3<-sapply(pPositionm3, function(x) JAK2_props[x]) + JAK2m2<-sapply(pPositionm2, function(x) JAK2_props[x]) + JAK2m1<-sapply(pPositionm1, function(x) JAK2_props[x]) + JAK2d0<-sapply(pPositiond0, function(x) JAK2_props[x]) + JAK2p1<-sapply(pPositionp1, function(x) JAK2_props[x]) + JAK2p2<-sapply(pPositionp2, function(x) JAK2_props[x]) + JAK2p3<-sapply(pPositionp3, function(x) JAK2_props[x]) + JAK2p4<-sapply(pPositionp4, function(x) JAK2_props[x]) + JAK2p5<-sapply(pPositionp5, function(x) JAK2_props[x]) + JAK2p6<-sapply(pPositionp6, function(x) JAK2_props[x]) + + Lckm6<-sapply(pPositionm6, function(x) Lck_props[x]) + Lckm5<-sapply(pPositionm5, function(x) Lck_props[x]) + Lckm4<-sapply(pPositionm4, function(x) Lck_props[x]) + Lckm3<-sapply(pPositionm3, function(x) Lck_props[x]) + Lckm2<-sapply(pPositionm2, function(x) Lck_props[x]) + Lckm1<-sapply(pPositionm1, function(x) Lck_props[x]) + Lckd0<-sapply(pPositiond0, function(x) Lck_props[x]) + Lckp1<-sapply(pPositionp1, function(x) Lck_props[x]) + Lckp2<-sapply(pPositionp2, function(x) Lck_props[x]) + Lckp3<-sapply(pPositionp3, function(x) Lck_props[x]) + Lckp4<-sapply(pPositionp4, function(x) Lck_props[x]) + Lckp5<-sapply(pPositionp5, function(x) Lck_props[x]) + Lckp6<-sapply(pPositionp6, function(x) Lck_props[x]) + + Lynm6<-sapply(pPositionm6, function(x) Lyn_props[x]) + Lynm5<-sapply(pPositionm5, function(x) Lyn_props[x]) + Lynm4<-sapply(pPositionm4, function(x) Lyn_props[x]) + Lynm3<-sapply(pPositionm3, function(x) Lyn_props[x]) + Lynm2<-sapply(pPositionm2, function(x) Lyn_props[x]) + Lynm1<-sapply(pPositionm1, function(x) Lyn_props[x]) + Lynd0<-sapply(pPositiond0, function(x) Lyn_props[x]) + Lynp1<-sapply(pPositionp1, function(x) Lyn_props[x]) + Lynp2<-sapply(pPositionp2, function(x) Lyn_props[x]) + Lynp3<-sapply(pPositionp3, function(x) Lyn_props[x]) + Lynp4<-sapply(pPositionp4, function(x) Lyn_props[x]) + Lynp5<-sapply(pPositionp5, function(x) Lyn_props[x]) + Lynp6<-sapply(pPositionp6, function(x) Lyn_props[x]) + + Pyk2m6<-sapply(pPositionm6, function(x) Pyk2_props[x]) + Pyk2m5<-sapply(pPositionm5, function(x) Pyk2_props[x]) + Pyk2m4<-sapply(pPositionm4, function(x) Pyk2_props[x]) + Pyk2m3<-sapply(pPositionm3, function(x) Pyk2_props[x]) + Pyk2m2<-sapply(pPositionm2, function(x) Pyk2_props[x]) + Pyk2m1<-sapply(pPositionm1, function(x) Pyk2_props[x]) + Pyk2d0<-sapply(pPositiond0, function(x) Pyk2_props[x]) + Pyk2p1<-sapply(pPositionp1, function(x) Pyk2_props[x]) + Pyk2p2<-sapply(pPositionp2, function(x) Pyk2_props[x]) + Pyk2p3<-sapply(pPositionp3, function(x) Pyk2_props[x]) + Pyk2p4<-sapply(pPositionp4, function(x) Pyk2_props[x]) + Pyk2p5<-sapply(pPositionp5, function(x) Pyk2_props[x]) + Pyk2p6<-sapply(pPositionp6, function(x) Pyk2_props[x]) + + Srcm6<-sapply(pPositionm6, function(x) Src_props[x]) + Srcm5<-sapply(pPositionm5, function(x) Src_props[x]) + Srcm4<-sapply(pPositionm4, function(x) Src_props[x]) + Srcm3<-sapply(pPositionm3, function(x) Src_props[x]) + Srcm2<-sapply(pPositionm2, function(x) Src_props[x]) + Srcm1<-sapply(pPositionm1, function(x) Src_props[x]) + Srcd0<-sapply(pPositiond0, function(x) Src_props[x]) + Srcp1<-sapply(pPositionp1, function(x) Src_props[x]) + Srcp2<-sapply(pPositionp2, function(x) Src_props[x]) + Srcp3<-sapply(pPositionp3, function(x) Src_props[x]) + Srcp4<-sapply(pPositionp4, function(x) Src_props[x]) + Srcp5<-sapply(pPositionp5, function(x) Src_props[x]) + Srcp6<-sapply(pPositionp6, function(x) Src_props[x]) + + Sykm6<-sapply(pPositionm6, function(x) Syk_props[x]) + Sykm5<-sapply(pPositionm5, function(x) Syk_props[x]) + Sykm4<-sapply(pPositionm4, function(x) Syk_props[x]) + Sykm3<-sapply(pPositionm3, function(x) Syk_props[x]) + Sykm2<-sapply(pPositionm2, function(x) Syk_props[x]) + Sykm1<-sapply(pPositionm1, function(x) Syk_props[x]) + Sykd0<-sapply(pPositiond0, function(x) Syk_props[x]) + Sykp1<-sapply(pPositionp1, function(x) Syk_props[x]) + Sykp2<-sapply(pPositionp2, function(x) Syk_props[x]) + Sykp3<-sapply(pPositionp3, function(x) Syk_props[x]) + Sykp4<-sapply(pPositionp4, function(x) Syk_props[x]) + Sykp5<-sapply(pPositionp5, function(x) Syk_props[x]) + Sykp6<-sapply(pPositionp6, function(x) Syk_props[x]) + + Yesm6<-sapply(pPositionm6, function(x) Yes_props[x]) + Yesm5<-sapply(pPositionm5, function(x) Yes_props[x]) + Yesm4<-sapply(pPositionm4, function(x) Yes_props[x]) + Yesm3<-sapply(pPositionm3, function(x) Yes_props[x]) + Yesm2<-sapply(pPositionm2, function(x) Yes_props[x]) + Yesm1<-sapply(pPositionm1, function(x) Yes_props[x]) + Yesd0<-sapply(pPositiond0, function(x) Yes_props[x]) + Yesp1<-sapply(pPositionp1, function(x) Yes_props[x]) + Yesp2<-sapply(pPositionp2, function(x) Yes_props[x]) + Yesp3<-sapply(pPositionp3, function(x) Yes_props[x]) + Yesp4<-sapply(pPositionp4, function(x) Yes_props[x]) + Yesp5<-sapply(pPositionp5, function(x) Yes_props[x]) + Yesp6<-sapply(pPositionp6, function(x) Yes_props[x]) + + FLT3m6<-sapply(pPositionm6, function(x) FLT3_props[x]) + FLT3m5<-sapply(pPositionm5, function(x) FLT3_props[x]) + FLT3m4<-sapply(pPositionm4, function(x) FLT3_props[x]) + FLT3m3<-sapply(pPositionm3, function(x) FLT3_props[x]) + FLT3m2<-sapply(pPositionm2, function(x) FLT3_props[x]) + FLT3m1<-sapply(pPositionm1, function(x) FLT3_props[x]) + FLT3d0<-sapply(pPositiond0, function(x) FLT3_props[x]) + FLT3p1<-sapply(pPositionp1, function(x) FLT3_props[x]) + FLT3p2<-sapply(pPositionp2, function(x) FLT3_props[x]) + FLT3p3<-sapply(pPositionp3, function(x) FLT3_props[x]) + FLT3p4<-sapply(pPositionp4, function(x) FLT3_props[x]) + FLT3p5<-sapply(pPositionp5, function(x) FLT3_props[x]) + FLT3p6<-sapply(pPositionp6, function(x) FLT3_props[x]) + + ALKm6<-sapply(pPositionm6, function(x) ALK_props[x]) + ALKm5<-sapply(pPositionm5, function(x) ALK_props[x]) + ALKm4<-sapply(pPositionm4, function(x) ALK_props[x]) + ALKm3<-sapply(pPositionm3, function(x) ALK_props[x]) + ALKm2<-sapply(pPositionm2, function(x) ALK_props[x]) + ALKm1<-sapply(pPositionm1, function(x) ALK_props[x]) + ALKd0<-sapply(pPositiond0, function(x) ALK_props[x]) + ALKp1<-sapply(pPositionp1, function(x) ALK_props[x]) + ALKp2<-sapply(pPositionp2, function(x) ALK_props[x]) + ALKp3<-sapply(pPositionp3, function(x) ALK_props[x]) + ALKp4<-sapply(pPositionp4, function(x) ALK_props[x]) + ALKp5<-sapply(pPositionp5, function(x) ALK_props[x]) + ALKp6<-sapply(pPositionp6, function(x) ALK_props[x]) +} +############################# +#here's my plan: EITHER +#make each peptide with different values for each kinase's EPM +#OR +#make the peptides as I'm doing now, then use apply functions. I think the first will be faster because I'm already spending time in the for loop +#so let's draw this out, the which function up above figures out which positions the kinase likes. Then I need an apply function to turn those into AAs +#I think I need 13 different apply functions, one for each plus or minus position of the AAs, sofor positions m7 every AA because AAm7 +#then I need 13 new apply functions, each with 13*20 values in them, which replace those values with individual values in the EPM table +#then I make the combinatorial library +#then I use an apply function to multiply all values together, which is how fisher does the scoring +######################################## + + +total=length(Positionp6)*length(Positionp5)*length(Positionp4)*length(Positionp3)*(length(Positionp2))*length(Positionp1)* + length(Positiond0)*length(Positionm1)*length(Positionm2)*length(Positionm3)*length(Positionm4)*length(Positionm5)*length(Positionm6) +#this is just a way to doublecheck that the length of the generated peptides vector is correct + +GeneratedPeptides<-rep(NA, times=total*13) +GeneratedPeptides<-matrix(data = GeneratedPeptides,ncol = 13) + + +ThisKinPeptides<-GeneratedPeptides + +AblPeptides<-GeneratedPeptides +ArgPeptides<-GeneratedPeptides +BtkPeptides<-GeneratedPeptides +CskPeptides<-GeneratedPeptides +FynPeptides<-GeneratedPeptides +HckPeptides<-GeneratedPeptides +JAK2Peptides<-GeneratedPeptides +LckPeptides<-GeneratedPeptides +LynPeptides<-GeneratedPeptides +Pyk2Peptides<-GeneratedPeptides +SrcPeptides<-GeneratedPeptides +SykPeptides<-GeneratedPeptides +YesPeptides<-GeneratedPeptides +FLT3Peptides<-GeneratedPeptides +ALKPeptides<-GeneratedPeptides + +#create an empty vector of correct length by finding the number of each AAs per position and multiplying them +count<-0 + +# for (t in 1:length(Positionm7)) { +for (s in 1:length(Positionm6)) { + for (r in 1:length(Positionm5)) { + for (i in 1:length(Positionm4)) { + for (j in 1:length(Positionm3)) { + for (k in 1:length(Positionm2)) { + for (l in 1:length(Positionm1)) { + for (m in 1:length(Positiond0)) { + for (n in 1:length(Positionp1)) { + for (o in 1:length(Positionp2)) { + for (p in 1:length(Positionp3)) { + for (q in 1:length(Positionp4)) { + for (u in 1:length(Positionp5)) { + for (v in 1:length(Positionp6)) { + # for (w in 1:length(Positionp7)) { + # i=1 + # j=1 + # k=1 + # l=1 + # m=1 + # n=1 + # o=1 + # p=1 + # q=1 + # + #for every single position, increment the count number, create a peptide using the AAs at that position + #then put them together into the generated peptides sequencex + count<-count+1 + tabulation<-c(Positionm6[s],Positionm5[r],Positionm4[i],Positionm3[j],Positionm2[k],Positionm1[l],Positiond0[m],Positionp1[n], + Positionp2[o],Positionp3[p],Positionp4[q],Positionp5[u],Positionp6[v]) + # numeration<-c(number14[s],number13[r],number1[i],number2[j],number3[k],number4[l],number5[m],number6[n],number7[o],number8[p],number9[q],number10[u],number11[v]) + #tabulation<-paste(tabulation, sep="", collapse="") + GeneratedPeptides[count,1:13]<-tabulation + # NumeratedPeptides[count,1:13]<-numeration + + ThisKinTableNumeration<-c(ThisKinTablem6[s],ThisKinTablem5[r],ThisKinTablem4[i],ThisKinTablem3[j],ThisKinTablem2[k],ThisKinTablem1[l],ThisKinTabled0[m],ThisKinTablep1[n], + ThisKinTablep2[o],ThisKinTablep3[p],ThisKinTablep4[q],ThisKinTablep5[u],ThisKinTablep6[v]) + + AblNumeration<-c(Ablm6[s],Ablm5[r],Ablm4[i],Ablm3[j],Ablm2[k],Ablm1[l],Abld0[m],Ablp1[n], + Ablp2[o],Ablp3[p],Ablp4[q],Ablp5[u],Ablp6[v]) + ArgNumeration<-c(Argm6[s],Argm5[r],Argm4[i],Argm3[j],Argm2[k],Argm1[l],Argd0[m],Argp1[n], + Argp2[o],Argp3[p],Argp4[q],Argp5[u],Argp6[v]) + BtkNumeration<-c(Btkm6[s],Btkm5[r],Btkm4[i],Btkm3[j],Btkm2[k],Btkm1[l],Btkd0[m],Btkp1[n], + Btkp2[o],Btkp3[p],Btkp4[q],Btkp5[u],Btkp6[v]) + CskNumeration<-c(Cskm6[s],Cskm5[r],Cskm4[i],Cskm3[j],Cskm2[k],Cskm1[l],Cskd0[m],Cskp1[n], + Cskp2[o],Cskp3[p],Cskp4[q],Cskp5[u],Cskp6[v]) + FynNumeration<-c(Fynm6[s],Fynm5[r],Fynm4[i],Fynm3[j],Fynm2[k],Fynm1[l],Fynd0[m],Fynp1[n], + Fynp2[o],Fynp3[p],Fynp4[q],Fynp5[u],Fynp6[v]) + HckNumeration<-c(Hckm6[s],Hckm5[r],Hckm4[i],Hckm3[j],Hckm2[k],Hckm1[l],Hckd0[m],Hckp1[n], + Hckp2[o],Hckp3[p],Hckp4[q],Hckp5[u],Hckp6[v]) + JAK2Numeration<-c(JAK2m6[s],JAK2m5[r],JAK2m4[i],JAK2m3[j],JAK2m2[k],JAK2m1[l],JAK2d0[m],JAK2p1[n], + JAK2p2[o],JAK2p3[p],JAK2p4[q],JAK2p5[u],JAK2p6[v]) + LckNumeration<-c(Lckm6[s],Lckm5[r],Lckm4[i],Lckm3[j],Lckm2[k],Lckm1[l],Lckd0[m],Lckp1[n], + Lckp2[o],Lckp3[p],Lckp4[q],Lckp5[u],Lckp6[v]) + LynNumeration<-c(Lynm6[s],Lynm5[r],Lynm4[i],Lynm3[j],Lynm2[k],Lynm1[l],Lynd0[m],Lynp1[n], + Lynp2[o],Lynp3[p],Lynp4[q],Lynp5[u],Lynp6[v]) + Pyk2Numeration<-c(Pyk2m6[s],Pyk2m5[r],Pyk2m4[i],Pyk2m3[j],Pyk2m2[k],Pyk2m1[l],Pyk2d0[m],Pyk2p1[n], + Pyk2p2[o],Pyk2p3[p],Pyk2p4[q],Pyk2p5[u],Pyk2p6[v]) + SrcNumeration<-c(Srcm6[s],Srcm5[r],Srcm4[i],Srcm3[j],Srcm2[k],Srcm1[l],Srcd0[m],Srcp1[n], + Srcp2[o],Srcp3[p],Srcp4[q],Srcp5[u],Srcp6[v]) + SykNumeration<-c(Sykm6[s],Sykm5[r],Sykm4[i],Sykm3[j],Sykm2[k],Sykm1[l],Sykd0[m],Sykp1[n], + Sykp2[o],Sykp3[p],Sykp4[q],Sykp5[u],Sykp6[v]) + YesNumeration<-c(Yesm6[s],Yesm5[r],Yesm4[i],Yesm3[j],Yesm2[k],Yesm1[l],Yesd0[m],Yesp1[n], + Yesp2[o],Yesp3[p],Yesp4[q],Yesp5[u],Yesp6[v]) + FLT3Numeration<-c(FLT3m6[s],FLT3m5[r],FLT3m4[i],FLT3m3[j],FLT3m2[k],FLT3m1[l],FLT3d0[m],FLT3p1[n], + FLT3p2[o],FLT3p3[p],FLT3p4[q],FLT3p5[u],FLT3p6[v]) + ALKNumeration<-c(ALKm6[s],ALKm5[r],ALKm4[i],ALKm3[j],ALKm2[k],ALKm1[l],ALKd0[m],ALKp1[n], + ALKp2[o],ALKp3[p],ALKp4[q],ALKp5[u],ALKp6[v]) + + ThisKinPeptides[count,1:13]<-ThisKinTableNumeration + + AblPeptides[count,1:13]<-AblNumeration + ArgPeptides[count,1:13]<-ArgNumeration + BtkPeptides[count,1:13]<-BtkNumeration + CskPeptides[count,1:13]<-CskNumeration + FynPeptides[count,1:13]<-FynNumeration + HckPeptides[count,1:13]<-HckNumeration + JAK2Peptides[count,1:13]<-JAK2Numeration + LckPeptides[count,1:13]<-LckNumeration + LynPeptides[count,1:13]<-LynNumeration + Pyk2Peptides[count,1:13]<-Pyk2Numeration + SrcPeptides[count,1:13]<-SrcNumeration + SykPeptides[count,1:13]<-SykNumeration + YesPeptides[count,1:13]<-YesNumeration + FLT3Peptides[count,1:13]<-FLT3Numeration + ALKPeptides[count,1:13]<-ALKNumeration + + + + # } + } + } + } + } + } + } + } + } + } + } + } + } +} +# } + +ThisKinPeptides<-as.data.frame(ThisKinPeptides) +ThisKinPeptides$V1<-as.numeric(ThisKinPeptides$V1) +ThisKinPeptides$V2<-as.numeric(ThisKinPeptides$V2) +ThisKinPeptides$V3<-as.numeric(ThisKinPeptides$V3) +ThisKinPeptides$V4<-as.numeric(ThisKinPeptides$V4) +ThisKinPeptides$V5<-as.numeric(ThisKinPeptides$V5) +ThisKinPeptides$V6<-as.numeric(ThisKinPeptides$V6) +ThisKinPeptides$V7<-as.numeric(ThisKinPeptides$V7) +ThisKinPeptides$V8<-as.numeric(ThisKinPeptides$V8) +ThisKinPeptides$V9<-as.numeric(ThisKinPeptides$V9) +ThisKinPeptides$V10<-as.numeric(ThisKinPeptides$V10) +ThisKinPeptides$V11<-as.numeric(ThisKinPeptides$V11) +ThisKinPeptides$V12<-as.numeric(ThisKinPeptides$V12) +ThisKinPeptides$V13<-as.numeric(ThisKinPeptides$V13) +ThisKinGeneratedScores<-with(ThisKinPeptides, V1*V2*V3*V4*V5*V6*V7*V8*V9*V10*V11*V12*V13) + +AblPeptides<-as.data.frame(AblPeptides) +AblGeneratedScores<-with(AblPeptides, V1*V2*V3*V4*V5*V6*V7*V8*V9*V10*V11*V12*V13) +ArgPeptides<-as.data.frame(ArgPeptides) +ArgGeneratedScores<-with(ArgPeptides, V1*V2*V3*V4*V5*V6*V7*V8*V9*V10*V11*V12*V13) +BtkPeptides<-as.data.frame(BtkPeptides) +BtkGeneratedScores<-with(BtkPeptides, V1*V2*V3*V4*V5*V6*V7*V8*V9*V10*V11*V12*V13) +CskPeptides<-as.data.frame(CskPeptides) +CskGeneratedScores<-with(CskPeptides, V1*V2*V3*V4*V5*V6*V7*V8*V9*V10*V11*V12*V13) +FynPeptides<-as.data.frame(FynPeptides) +FynGeneratedScores<-with(FynPeptides, V1*V2*V3*V4*V5*V6*V7*V8*V9*V10*V11*V12*V13) +HckPeptides<-as.data.frame(HckPeptides) +HckGeneratedScores<-with(HckPeptides, V1*V2*V3*V4*V5*V6*V7*V8*V9*V10*V11*V12*V13) +JAK2Peptides<-as.data.frame(JAK2Peptides) +JAK2GeneratedScores<-with(JAK2Peptides, V1*V2*V3*V4*V5*V6*V7*V8*V9*V10*V11*V12*V13) +LckPeptides<-as.data.frame(LckPeptides) +LckGeneratedScores<-with(LckPeptides, V1*V2*V3*V4*V5*V6*V7*V8*V9*V10*V11*V12*V13) +LynPeptides<-as.data.frame(LynPeptides) +LynGeneratedScores<-with(LynPeptides, V1*V2*V3*V4*V5*V6*V7*V8*V9*V10*V11*V12*V13) +Pyk2Peptides<-as.data.frame(Pyk2Peptides) +Pyk2GeneratedScores<-with(Pyk2Peptides, V1*V2*V3*V4*V5*V6*V7*V8*V9*V10*V11*V12*V13) +SrcPeptides<-as.data.frame(SrcPeptides) +SrcGeneratedScores<-with(SrcPeptides, V1*V2*V3*V4*V5*V6*V7*V8*V9*V10*V11*V12*V13) +SykPeptides<-as.data.frame(SykPeptides) +SykGeneratedScores<-with(SykPeptides, V1*V2*V3*V4*V5*V6*V7*V8*V9*V10*V11*V12*V13) +YesPeptides<-as.data.frame(YesPeptides) +YesGeneratedScores<-with(YesPeptides, V1*V2*V3*V4*V5*V6*V7*V8*V9*V10*V11*V12*V13) +FLT3Peptides<-as.data.frame(FLT3Peptides) +FLT3GeneratedScores<-with(FLT3Peptides, V1*V2*V3*V4*V5*V6*V7*V8*V9*V10*V11*V12*V13) +ALKPeptides<-as.data.frame(ALKPeptides) +ALKGeneratedScores<-with(ALKPeptides, V1*V2*V3*V4*V5*V6*V7*V8*V9*V10*V11*V12*V13) + +AblGeneratedScores<-as.matrix(AblGeneratedScores, ncol=1) +AblNorm<-AblNorm<-1/as.numeric(Abl[22,1]) +AblNormColumn<-as.matrix(rep(AblNorm,times=nrow(AblGeneratedScores)),ncol=1) +AblGeneratedScores<-cbind(AblGeneratedScores,AblNormColumn) +AblGeneratedScores<-as.data.frame(AblGeneratedScores) +AblGeneratedScores<-with(AblGeneratedScores, V1/(V1+V2)) + +ArgGeneratedScores<-as.matrix(ArgGeneratedScores, ncol=1) +ArgNorm<-ArgNorm<-1/as.numeric(Arg[22,1]) +ArgNormColumn<-as.matrix(rep(ArgNorm,times=nrow(ArgGeneratedScores)),ncol=1) +ArgGeneratedScores<-cbind(ArgGeneratedScores,ArgNormColumn) +ArgGeneratedScores<-as.data.frame(ArgGeneratedScores) +ArgGeneratedScores<-with(ArgGeneratedScores, V1/(V1+V2)) + +BtkGeneratedScores<-as.matrix(BtkGeneratedScores, ncol=1) +BtkNorm<-BtkNorm<-1/as.numeric(Btk[22,1]) +BtkNormColumn<-as.matrix(rep(BtkNorm,times=nrow(BtkGeneratedScores)),ncol=1) +BtkGeneratedScores<-cbind(BtkGeneratedScores,BtkNormColumn) +BtkGeneratedScores<-as.data.frame(BtkGeneratedScores) +BtkGeneratedScores<-with(BtkGeneratedScores, V1/(V1+V2)) + +CskGeneratedScores<-as.matrix(CskGeneratedScores, ncol=1) +CskNorm<-CskNorm<-1/as.numeric(Csk[22,1]) +CskNormColumn<-as.matrix(rep(CskNorm,times=nrow(CskGeneratedScores)),ncol=1) +CskGeneratedScores<-cbind(CskGeneratedScores,CskNormColumn) +CskGeneratedScores<-as.data.frame(CskGeneratedScores) +CskGeneratedScores<-with(CskGeneratedScores, V1/(V1+V2)) + +FynGeneratedScores<-as.matrix(FynGeneratedScores, ncol=1) +FynNorm<-FynNorm<-1/as.numeric(Fyn[22,1]) +FynNormColumn<-as.matrix(rep(FynNorm,times=nrow(FynGeneratedScores)),ncol=1) +FynGeneratedScores<-cbind(FynGeneratedScores,FynNormColumn) +FynGeneratedScores<-as.data.frame(FynGeneratedScores) +FynGeneratedScores<-with(FynGeneratedScores, V1/(V1+V2)) + +HckGeneratedScores<-as.matrix(HckGeneratedScores, ncol=1) +HckNorm<-HckNorm<-1/as.numeric(Hck[22,1]) +HckNormColumn<-as.matrix(rep(HckNorm,times=nrow(HckGeneratedScores)),ncol=1) +HckGeneratedScores<-cbind(HckGeneratedScores,HckNormColumn) +HckGeneratedScores<-as.data.frame(HckGeneratedScores) +HckGeneratedScores<-with(HckGeneratedScores, 100*V1/(V1+V2)) + +JAK2GeneratedScores<-as.matrix(JAK2GeneratedScores, ncol=1) +JAK2Norm<-JAK2Norm<-1/as.numeric(JAK2[22,1]) +JAK2NormColumn<-as.matrix(rep(JAK2Norm,times=nrow(JAK2GeneratedScores)),ncol=1) +JAK2GeneratedScores<-cbind(JAK2GeneratedScores,JAK2NormColumn) +JAK2GeneratedScores<-as.data.frame(JAK2GeneratedScores) +JAK2GeneratedScores<-with(JAK2GeneratedScores, 100*V1/(V1+V2)) + +LckGeneratedScores<-as.matrix(LckGeneratedScores, ncol=1) +LckNorm<-LckNorm<-1/as.numeric(Lck[22,1]) +LckNormColumn<-as.matrix(rep(LckNorm,times=nrow(LckGeneratedScores)),ncol=1) +LckGeneratedScores<-cbind(LckGeneratedScores,LckNormColumn) +LckGeneratedScores<-as.data.frame(LckGeneratedScores) +LckGeneratedScores<-with(LckGeneratedScores, 100*V1/(V1+V2)) + +LynGeneratedScores<-as.matrix(LynGeneratedScores, ncol=1) +LynNorm<-LynNorm<-1/as.numeric(Lyn[22,1]) +LynNormColumn<-as.matrix(rep(LynNorm,times=nrow(LynGeneratedScores)),ncol=1) +LynGeneratedScores<-cbind(LynGeneratedScores,LynNormColumn) +LynGeneratedScores<-as.data.frame(LynGeneratedScores) +LynGeneratedScores<-with(LynGeneratedScores, 100*V1/(V1+V2)) + +Pyk2GeneratedScores<-as.matrix(Pyk2GeneratedScores, ncol=1) +Pyk2Norm<-Pyk2Norm<-1/as.numeric(Pyk2[22,1]) +Pyk2NormColumn<-as.matrix(rep(Pyk2Norm,times=nrow(Pyk2GeneratedScores)),ncol=1) +Pyk2GeneratedScores<-cbind(Pyk2GeneratedScores,Pyk2NormColumn) +Pyk2GeneratedScores<-as.data.frame(Pyk2GeneratedScores) +Pyk2GeneratedScores<-with(Pyk2GeneratedScores, 100*V1/(V1+V2)) + +SrcGeneratedScores<-as.matrix(SrcGeneratedScores, ncol=1) +SrcNorm<-SrcNorm<-1/as.numeric(Src[22,1]) +SrcNormColumn<-as.matrix(rep(SrcNorm,times=nrow(SrcGeneratedScores)),ncol=1) +SrcGeneratedScores<-cbind(SrcGeneratedScores,SrcNormColumn) +SrcGeneratedScores<-as.data.frame(SrcGeneratedScores) +SrcGeneratedScores<-with(SrcGeneratedScores, 100*V1/(V1+V2)) + +SykGeneratedScores<-as.matrix(SykGeneratedScores, ncol=1) +SykNorm<-SykNorm<-1/as.numeric(Syk[22,1]) +SykNormColumn<-as.matrix(rep(SykNorm,times=nrow(SykGeneratedScores)),ncol=1) +SykGeneratedScores<-cbind(SykGeneratedScores,SykNormColumn) +SykGeneratedScores<-as.data.frame(SykGeneratedScores) +SykGeneratedScores<-with(SykGeneratedScores, 100*V1/(V1+V2)) + +YesGeneratedScores<-as.matrix(YesGeneratedScores, ncol=1) +YesNorm<-YesNorm<-1/as.numeric(Yes[22,1]) +YesNormColumn<-as.matrix(rep(YesNorm,times=nrow(YesGeneratedScores)),ncol=1) +YesGeneratedScores<-cbind(YesGeneratedScores,YesNormColumn) +YesGeneratedScores<-as.data.frame(YesGeneratedScores) +YesGeneratedScores<-with(YesGeneratedScores, 100*V1/(V1+V2)) + +FLT3GeneratedScores<-as.matrix(FLT3GeneratedScores, ncol=1) +FLT3Norm<-FLT3Norm<-1/as.numeric(FLT3[22,1]) +FLT3NormColumn<-as.matrix(rep(FLT3Norm,times=nrow(FLT3GeneratedScores)),ncol=1) +FLT3GeneratedScores<-cbind(FLT3GeneratedScores,FLT3NormColumn) +FLT3GeneratedScores<-as.data.frame(FLT3GeneratedScores) +FLT3GeneratedScores<-with(FLT3GeneratedScores, 100*V1/(V1+V2)) + +ALKGeneratedScores<-as.matrix(ALKGeneratedScores, ncol=1) +ALKNorm<-ALKNorm<-1/as.numeric(ALK[22,1]) +ALKNormColumn<-as.matrix(rep(ALKNorm,times=nrow(ALKGeneratedScores)),ncol=1) +ALKGeneratedScores<-cbind(ALKGeneratedScores,ALKNormColumn) +ALKGeneratedScores<-as.data.frame(ALKGeneratedScores) +ALKGeneratedScores<-with(ALKGeneratedScores, 100*V1/(V1+V2)) + +ThisKinGeneratedScores1<-as.matrix(ThisKinGeneratedScores, ncol=1) +ThisKinNorm<-ThisKinNorm<-1/as.numeric(NormalizationScore) +ThisKinNormColumn<-as.matrix(rep(ThisKinNorm,times=nrow(ThisKinGeneratedScores1)),ncol=1) +ThisKinGeneratedScores1<-cbind(ThisKinGeneratedScores1,ThisKinNormColumn) +ThisKinGeneratedScores1<-as.data.frame(ThisKinGeneratedScores1) +ThisKinGeneratedScores1<-with(ThisKinGeneratedScores1, V1/(V1+V2)) + +AblThresh<-as.numeric(Abl[24,1]) +AblActive<-unlist(AblGeneratedScores)>AblThresh/100 +ArgThresh<-as.numeric(Arg[24,1]) +ArgActive<-unlist(ArgGeneratedScores)>ArgThresh/100 +BtkThresh<-as.numeric(Btk[24,1]) +BtkActive<-unlist(BtkGeneratedScores)>BtkThresh/100 +CskThresh<-as.numeric(Csk[24,1]) +CskActive<-(CskGeneratedScores)>CskThresh/100 +FynThresh<-as.numeric(Fyn[24,1]) +FynActive<-unlist(FynGeneratedScores)>FynThresh/100 +HckThresh<-as.numeric(Hck[24,1]) +HckActive<-unlist(HckGeneratedScores)>HckThresh/100 +JAK2Thresh<-as.numeric(JAK2[24,1]) +JAk2Active<-unlist(JAK2GeneratedScores)>JAK2Thresh/100 +LckThresh<-as.numeric(Lck[24,1]) +LckActive<-unlist(LckGeneratedScores)>LckThresh/100 +LynThresh<-as.numeric(Lyn[24,1]) +LynActive<-unlist(LynGeneratedScores)>LynThresh/100 +Pyk2Thresh<-as.numeric(Pyk2[24,1]) +Pyk2Active<-unlist(Pyk2GeneratedScores)>Pyk2Thresh/100 +SrcThresh<-as.numeric(Src[24,1]) +SrcActive<-unlist(SrcGeneratedScores)>SrcThresh/100 +SykThresh<-as.numeric(Syk[24,1]) +SykActive<-unlist(SykGeneratedScores)>SykThresh/100 +YesThresh<-as.numeric(Yes[24,1]) +YesActive<-unlist(YesGeneratedScores)>YesThresh/100 +ALKThresh<-as.numeric(ALK[24,1]) +ALKActive<-unlist(ALKGeneratedScores)>ALKThresh/100 +FLT3Thresh<-as.numeric(FLT3[24,1]) +FLT3Active<-unlist(FLT3GeneratedScores)>FLT3Thresh/100 + + +AllActive<-AblActive+ArgActive+BtkActive+CskActive+FynActive+HckActive+JAk2Active+LckActive+LynActive+Pyk2Active+SrcActive+SykActive+YesActive+ALKActive+FLT3Active +Scores<-ThisKinGeneratedScores +ThresholdValues<-ThisKinGeneratedScores1 + +FullMotifs<-rep("Z",times=nrow(GeneratedPeptides)) +for (i in 1:nrow(GeneratedPeptides)) { + motif<-GeneratedPeptides[i,1:13] + motif<-paste(motif,sep = "", collapse = "") + FullMotifs[i]<-motif +} + +Scores<-1/Scores + +PeptidesWithRanks<-cbind.data.frame(FullMotifs,GeneratedPeptides,Scores,ThresholdValues) + +PeptidesWithRanks<-cbind.data.frame(PeptidesWithRanks,AllActive,AblActive,ArgActive,BtkActive,CskActive,FynActive,HckActive,JAk2Active,LckActive,LynActive,Pyk2Active,SrcActive,SykActive,YesActive) +RanksPeptides<-PeptidesWithRanks[order(PeptidesWithRanks$AllActive, + PeptidesWithRanks$Scores, + decreasing = FALSE),] + +RanksPeptides$Scores<-1/RanksPeptides$Scores + +#RanksPeptides$Scores<-1/RanksPeptides$Scores + +#I want to see how much of the sequence space we still have available. I should do so + +write.table(RanksPeptides,file = "output.csv",append = FALSE,row.names = FALSE,col.names = TRUE,sep = ",") + +# }) \ No newline at end of file