0
|
1 # system.time({
|
|
2
|
6
|
3 options(warn=-1)
|
|
4
|
5
|
5 args = commandArgs(trailingOnly=TRUE)
|
3
|
6 #args = c("ABL","ABC","ABC","ABC","ABC","ABC","ABC","ABC","ABC","ABC","ABC","ABC","ABC","ABC")
|
0
|
7 TodaysKinase<-as.character(args[1])
|
|
8
|
|
9 CharTable<-read.csv("Fisher-Char.csv", header = FALSE,stringsAsFactors = FALSE)
|
|
10 ThisKinSheet<-read.csv("Fisher-Table.csv", header=FALSE, stringsAsFactors=FALSE)
|
|
11 NormalizationScore<-CharTable[2,1]
|
|
12
|
|
13 #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.
|
|
14
|
|
15
|
|
16
|
|
17
|
3
|
18 M6<-as.character(args[2])
|
|
19 M5<-as.character(args[3])
|
|
20 M4<-as.character(args[4])
|
|
21 M3<-as.character(args[5])
|
|
22 M2<-as.character(args[6])
|
|
23 M1<-as.character(args[7])
|
|
24 D0<-as.character(args[8])
|
|
25 P1<-as.character(args[9])
|
|
26 P2<-as.character(args[10])
|
|
27 P3<-as.character(args[11])
|
|
28 P4<-as.character(args[12])
|
|
29 P5<-as.character(args[13])
|
|
30 P6<-as.character(args[14])
|
0
|
31
|
|
32
|
|
33 A=1
|
|
34 C=2
|
|
35 D=3
|
|
36 E=4
|
|
37 F=5
|
|
38 G=6
|
|
39 H=7
|
|
40 I=8
|
|
41 K=9
|
|
42 L=10
|
|
43 M=11
|
|
44 N=12
|
|
45 P=13
|
|
46 Q=14
|
|
47 R=15
|
|
48 S=16
|
|
49 T=17
|
|
50 V=18
|
|
51 W=19
|
|
52 Y=20
|
|
53 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,
|
|
54 "S"=S,"T"=T,"V"=V,"W"=W,"Y"=Y,"xY"=Y,"O"=21)
|
|
55 Positionm6<-gsub("[^a-zA-Z]", "", M6)
|
|
56 Positionm6<-toupper(Positionm6)
|
|
57 Positionm6<-unlist(strsplit(Positionm6,""))
|
|
58 Positionm6<-sapply(Positionm6, function(x) aa_props[x])
|
|
59
|
|
60 Positionm5<-gsub("[^a-zA-Z]", "", M5)
|
|
61 Positionm5<-toupper(Positionm5)
|
|
62 Positionm5<-unlist(strsplit(Positionm5,""))
|
|
63 Positionm5<-sapply(Positionm5, function(x) aa_props[x])
|
|
64
|
|
65 Positionm4<-gsub("[^a-zA-Z]", "", M4)
|
|
66 Positionm4<-toupper(Positionm4)
|
|
67 Positionm4<-unlist(strsplit(Positionm4,""))
|
|
68 Positionm4<-sapply(Positionm4, function(x) aa_props[x])
|
|
69
|
|
70 Positionm3<-gsub("[^a-zA-Z]", "", M3)
|
|
71 Positionm3<-toupper(Positionm3)
|
|
72 Positionm3<-unlist(strsplit(Positionm3,""))
|
|
73 Positionm3<-sapply(Positionm3, function(x) aa_props[x])
|
|
74
|
|
75 Positionm2<-gsub("[^a-zA-Z]", "", M2)
|
|
76 Positionm2<-toupper(Positionm2)
|
|
77 Positionm2<-unlist(strsplit(Positionm2,""))
|
|
78 Positionm2<-sapply(Positionm2, function(x) aa_props[x])
|
|
79
|
|
80 Positionm1<-gsub("[^a-zA-Z]", "", M1)
|
|
81 Positionm1<-toupper(Positionm1)
|
|
82 Positionm1<-unlist(strsplit(Positionm1,""))
|
|
83 Positionm1<-sapply(Positionm1, function(x) aa_props[x])
|
|
84
|
|
85 Positiond0<-gsub("[^a-zA-Z]", "", D0)
|
|
86 Positiond0<-toupper(Positiond0)
|
|
87 Positiond0<-unlist(strsplit(Positiond0,""))
|
|
88 Positiond0<-sapply(Positiond0, function(x) aa_props[x])
|
|
89
|
|
90 Positionp1<-gsub("[^a-zA-Z]", "", P1)
|
|
91 Positionp1<-toupper(Positionp1)
|
|
92 Positionp1<-unlist(strsplit(Positionp1,""))
|
|
93 Positionp1<-sapply(Positionp1, function(x) aa_props[x])
|
|
94
|
|
95 Positionp2<-gsub("[^a-zA-Z]", "", P2)
|
|
96 Positionp2<-toupper(Positionp2)
|
|
97 Positionp2<-unlist(strsplit(Positionp2,""))
|
|
98 Positionp2<-sapply(Positionp2, function(x) aa_props[x])
|
|
99
|
|
100 Positionp3<-gsub("[^a-zA-Z]", "", P3)
|
|
101 Positionp3<-toupper(Positionp3)
|
|
102 Positionp3<-unlist(strsplit(Positionp3,""))
|
|
103 Positionp3<-sapply(Positionp3, function(x) aa_props[x])
|
|
104
|
|
105 Positionp4<-gsub("[^a-zA-Z]", "", P4)
|
|
106 Positionp4<-toupper(Positionp4)
|
|
107 Positionp4<-unlist(strsplit(Positionp4,""))
|
|
108 Positionp4<-sapply(Positionp4, function(x) aa_props[x])
|
|
109
|
|
110 Positionp5<-gsub("[^a-zA-Z]", "", P5)
|
|
111 Positionp5<-toupper(Positionp5)
|
|
112 Positionp5<-unlist(strsplit(Positionp5,""))
|
|
113 Positionp5<-sapply(Positionp5, function(x) aa_props[x])
|
|
114
|
|
115 Positionp6<-gsub("[^a-zA-Z]", "", P6)
|
|
116 Positionp6<-toupper(Positionp6)
|
|
117 Positionp6<-unlist(strsplit(Positionp6,""))
|
|
118 Positionp6<-sapply(Positionp6, function(x) aa_props[x])
|
|
119
|
|
120
|
|
121
|
|
122
|
|
123 #what the below code does is stupid. I need to make a combinatorial library of peptides and then calculate their values
|
|
124 #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
|
|
125 #IE an A at position negative 6 is called m6a (minus 6 Alanine)
|
|
126 #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
|
|
127 #then I will multiply those values together
|
|
128
|
|
129 if(2==2){
|
|
130 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",
|
|
131 "13"="P", "14"="Q", "15"="R", "16"="S", "17"="T", "18"="V", "19"="W", "20"="Y", "21"="O")
|
|
132
|
|
133 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",
|
|
134 "13"="m6P", "14"="m6Q", "15"="m6R", "16"="m6S", "17"="m6T", "18"="m6V", "19"="m6W", "20"="m6Y", "21"="m6O")
|
|
135
|
|
136 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",
|
|
137 "13"="m5P", "14"="m5Q", "15"="m5R", "16"="m5S", "17"="m5T", "18"="m5V", "19"="m5W", "20"="m5Y", "21"="m5O")
|
|
138
|
|
139 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",
|
|
140 "13"="m4P", "14"="m4Q", "15"="m4R", "16"="m4S", "17"="m4T", "18"="m4V", "19"="m4W", "20"="m4Y", "21"="m4O")
|
|
141
|
|
142 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",
|
|
143 "13"="m3P", "14"="m3Q", "15"="m3R", "16"="m3S", "17"="m3T", "18"="m3V", "19"="m3W", "20"="m3Y", "21"="m3O")
|
|
144
|
|
145 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",
|
|
146 "13"="m2P", "14"="m2Q", "15"="m2R", "16"="m2S", "17"="m2T", "18"="m2V", "19"="m2W", "20"="m2Y", "21"="m2O")
|
|
147
|
|
148 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",
|
|
149 "13"="m1P", "14"="m1Q", "15"="m1R", "16"="m1S", "17"="m1T", "18"="m1V", "19"="m1W", "20"="m1Y", "21"="m1O")
|
|
150
|
|
151 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",
|
|
152 "13"="d0P", "14"="d0Q", "15"="d0R", "16"="d0S", "17"="d0T", "18"="d0V", "19"="d0W", "20"="d0Y", "21"="d0O")
|
|
153
|
|
154 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",
|
|
155 "13"="p1P", "14"="p1Q", "15"="p1R", "16"="p1S", "17"="p1T", "18"="p1V", "19"="p1W", "20"="p1Y", "21"="p1O")
|
|
156
|
|
157 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",
|
|
158 "13"="p2P", "14"="p2Q", "15"="p2R", "16"="p2S", "17"="p2T", "18"="p2V", "19"="p2W", "20"="p2Y", "21"="p2O")
|
|
159
|
|
160 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",
|
|
161 "13"="p3P", "14"="p3Q", "15"="p3R", "16"="p3S", "17"="p3T", "18"="p3V", "19"="p3W", "20"="p3Y", "21"="p3O")
|
|
162
|
|
163 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",
|
|
164 "13"="p4P", "14"="p4Q", "15"="p4R", "16"="p4S", "17"="p4T", "18"="p4V", "19"="p4W", "20"="p4Y", "21"="p4O")
|
|
165
|
|
166 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",
|
|
167 "13"="p5P", "14"="p5Q", "15"="p5R", "16"="p5S", "17"="p5T", "18"="p5V", "19"="p5W", "20"="p5Y", "21"="p5O")
|
|
168
|
|
169 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",
|
|
170 "13"="p6P", "14"="p6Q", "15"="p6R", "16"="p6S", "17"="p6T", "18"="p6V", "19"="p6W", "20"="p6Y", "21"="p6O")
|
|
171
|
|
172 # Positionm7<-sapply(Positionm7, function (x) aa_props2[x])
|
|
173 pPositionm6<-sapply(Positionm6, function (x) m6_props2[x])
|
|
174 pPositionm5<-sapply(Positionm5, function (x) m5_props2[x])
|
|
175 pPositionm4<-sapply(Positionm4, function (x) m4_props2[x])
|
|
176 pPositionm3<-sapply(Positionm3, function (x) m3_props2[x])
|
|
177 pPositionm2<-sapply(Positionm2, function (x) m2_props2[x])
|
|
178 pPositionm1<-sapply(Positionm1, function (x) m1_props2[x])
|
|
179 pPositiond0<-sapply(Positiond0, function (x) d0_props2[x])
|
|
180 pPositionp1<-sapply(Positionp1, function (x) p1_props2[x])
|
|
181 pPositionp2<-sapply(Positionp2, function (x) p2_props2[x])
|
|
182 pPositionp3<-sapply(Positionp3, function (x) p3_props2[x])
|
|
183 pPositionp4<-sapply(Positionp4, function (x) p4_props2[x])
|
|
184 pPositionp5<-sapply(Positionp5, function (x) p5_props2[x])
|
|
185 pPositionp6<-sapply(Positionp6, function (x) p6_props2[x])
|
|
186 # Positionp7<-sapply(Positionp7, function (x) aa_props2[x])
|
|
187
|
|
188 Positionm6<-sapply(Positionm6, function (x) aa_props2[x])
|
|
189 Positionm5<-sapply(Positionm5, function (x) aa_props2[x])
|
|
190 Positionm4<-sapply(Positionm4, function (x) aa_props2[x])
|
|
191 Positionm3<-sapply(Positionm3, function (x) aa_props2[x])
|
|
192 Positionm2<-sapply(Positionm2, function (x) aa_props2[x])
|
|
193 Positionm1<-sapply(Positionm1, function (x) aa_props2[x])
|
|
194 Positiond0<-sapply(Positiond0, function (x) aa_props2[x])
|
|
195 Positionp1<-sapply(Positionp1, function (x) aa_props2[x])
|
|
196 Positionp2<-sapply(Positionp2, function (x) aa_props2[x])
|
|
197 Positionp3<-sapply(Positionp3, function (x) aa_props2[x])
|
|
198 Positionp4<-sapply(Positionp4, function (x) aa_props2[x])
|
|
199 Positionp5<-sapply(Positionp5, function (x) aa_props2[x])
|
|
200 Positionp6<-sapply(Positionp6, function (x) aa_props2[x])
|
|
201
|
|
202
|
|
203 }
|
|
204 #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,
|
|
205 #if it's a B at positive 4 it's p4B
|
|
206
|
|
207
|
|
208
|
|
209 ScreenerFilename<-"screener.csv"
|
|
210 screaner<-read.csv(ScreenerFilename, header = FALSE, stringsAsFactors = FALSE)
|
|
211 #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)
|
|
212
|
|
213 Abl<-screaner[2:25,]
|
|
214 Arg<-screaner[27:50,]
|
|
215 Btk<-screaner[52:75,]
|
|
216 Csk<-screaner[77:100,]
|
|
217 Fyn<-screaner[102:125,]
|
|
218 Hck<-screaner[127:150,]
|
|
219 JAK2<-screaner[152:175,]
|
|
220 Lck<-screaner[177:200,]
|
|
221 Lyn<-screaner[202:225,]
|
|
222 Pyk2<-screaner[227:250,]
|
|
223 Src<-screaner[252:275,]
|
|
224 Syk<-screaner[277:300,]
|
|
225 Yes<-screaner[302:325,]
|
|
226 FLT3<-screaner[327:350,]
|
|
227 ALK<-screaner[352:375,]
|
|
228
|
|
229
|
|
230
|
|
231 ThisKinTable<-ThisKinSheet[2:22,]
|
|
232 #NormalizationScore<-ThisKinSheet[27,1]
|
|
233
|
|
234 if(1==1){
|
|
235 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],
|
|
236 "m6S"=ThisKinTable[16,3],"m6T"=ThisKinTable[17,3],"m6V"=ThisKinTable[18,3],"m6W"=ThisKinTable[19,3],"m6Y"=ThisKinTable[20,3],"m6O"=1,
|
|
237 "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],
|
|
238 "m5S"=ThisKinTable[16,4],"m5T"=ThisKinTable[17,4],"m5V"=ThisKinTable[18,4],"m5W"=ThisKinTable[19,4],"m5Y"=ThisKinTable[20,4],"m5O"=1,
|
|
239 "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],
|
|
240 "m4S"=ThisKinTable[16,5],"m4T"=ThisKinTable[17,5],"m4V"=ThisKinTable[18,5],"m4W"=ThisKinTable[19,5],"m4Y"=ThisKinTable[20,5],"m4O"=1,
|
|
241 "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],
|
|
242 "m3S"=ThisKinTable[16,6],"m3T"=ThisKinTable[17,6],"m3V"=ThisKinTable[18,6],"m3W"=ThisKinTable[19,6],"m3Y"=ThisKinTable[20,6],"m3O"=1,
|
|
243 "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],
|
|
244 "m2S"=ThisKinTable[16,7],"m2T"=ThisKinTable[17,7],"m2V"=ThisKinTable[18,7],"m2W"=ThisKinTable[19,7],"m2Y"=ThisKinTable[20,7],"m2O"=1,
|
|
245 "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],
|
|
246 "m1S"=ThisKinTable[16,8],"m1T"=ThisKinTable[17,8],"m1V"=ThisKinTable[18,8],"m1W"=ThisKinTable[19,8],"m1Y"=ThisKinTable[20,8],"m1O"=1,
|
|
247 "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],
|
|
248 "d0S"=ThisKinTable[16,9],"d0T"=ThisKinTable[17,9],"d0V"=ThisKinTable[18,9],"d0W"=ThisKinTable[19,9],"d0Y"=ThisKinTable[20,9],"d0O"=1,
|
|
249 "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],
|
|
250 "p1S"=ThisKinTable[16,10],"p1T"=ThisKinTable[17,10],"p1V"=ThisKinTable[18,10],"p1W"=ThisKinTable[19,10],"p1Y"=ThisKinTable[20,10],"p1O"=1,
|
|
251 "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],
|
|
252 "p2S"=ThisKinTable[16,11],"p2T"=ThisKinTable[17,11],"p2V"=ThisKinTable[18,11],"p2W"=ThisKinTable[19,11],"p2Y"=ThisKinTable[20,11],"p2O"=1,
|
|
253 "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],
|
|
254 "p3S"=ThisKinTable[16,12],"p3T"=ThisKinTable[17,12],"p3V"=ThisKinTable[18,12],"p3W"=ThisKinTable[19,12],"p3Y"=ThisKinTable[20,12],"p3O"=1,
|
|
255 "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],
|
|
256 "p4S"=ThisKinTable[16,13],"p4T"=ThisKinTable[17,13],"p4V"=ThisKinTable[18,13],"p4W"=ThisKinTable[19,13],"p4Y"=ThisKinTable[20,13],"p4O"=1,
|
|
257 "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],
|
|
258 "p5S"=ThisKinTable[16,14],"p5T"=ThisKinTable[17,14],"p5V"=ThisKinTable[18,14],"p5W"=ThisKinTable[19,14],"p5Y"=ThisKinTable[20,14],"p5O"=1,
|
|
259 "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],
|
|
260 "p6S"=ThisKinTable[16,15],"p6T"=ThisKinTable[17,15],"p6V"=ThisKinTable[18,15],"p6W"=ThisKinTable[19,15],"p6Y"=ThisKinTable[20,15],"p6O"=1)
|
|
261
|
|
262
|
|
263 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],
|
|
264 "m6S"=Abl[16,3],"m6T"=Abl[17,3],"m6V"=Abl[18,3],"m6W"=Abl[19,3],"m6Y"=Abl[20,3],"m6O"=1,
|
|
265 "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],
|
|
266 "m5S"=Abl[16,4],"m5T"=Abl[17,4],"m5V"=Abl[18,4],"m5W"=Abl[19,4],"m5Y"=Abl[20,4],"m5O"=1,
|
|
267 "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],
|
|
268 "m4S"=Abl[16,5],"m4T"=Abl[17,5],"m4V"=Abl[18,5],"m4W"=Abl[19,5],"m4Y"=Abl[20,5],"m4O"=1,
|
|
269 "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],
|
|
270 "m3S"=Abl[16,6],"m3T"=Abl[17,6],"m3V"=Abl[18,6],"m3W"=Abl[19,6],"m3Y"=Abl[20,6],"m3O"=1,
|
|
271 "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],
|
|
272 "m2S"=Abl[16,7],"m2T"=Abl[17,7],"m2V"=Abl[18,7],"m2W"=Abl[19,7],"m2Y"=Abl[20,7],"m2O"=1,
|
|
273 "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],
|
|
274 "m1S"=Abl[16,8],"m1T"=Abl[17,8],"m1V"=Abl[18,8],"m1W"=Abl[19,8],"m1Y"=Abl[20,8],"m1O"=1,
|
|
275 "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],
|
|
276 "d0S"=Abl[16,9],"d0T"=Abl[17,9],"d0V"=Abl[18,9],"d0W"=Abl[19,9],"d0Y"=Abl[20,9],"d0O"=1,
|
|
277 "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],
|
|
278 "p1S"=Abl[16,10],"p1T"=Abl[17,10],"p1V"=Abl[18,10],"p1W"=Abl[19,10],"p1Y"=Abl[20,10],"p1O"=1,
|
|
279 "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],
|
|
280 "p2S"=Abl[16,11],"p2T"=Abl[17,11],"p2V"=Abl[18,11],"p2W"=Abl[19,11],"p2Y"=Abl[20,11],"p2O"=1,
|
|
281 "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],
|
|
282 "p3S"=Abl[16,12],"p3T"=Abl[17,12],"p3V"=Abl[18,12],"p3W"=Abl[19,12],"p3Y"=Abl[20,12],"p3O"=1,
|
|
283 "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],
|
|
284 "p4S"=Abl[16,13],"p4T"=Abl[17,13],"p4V"=Abl[18,13],"p4W"=Abl[19,13],"p4Y"=Abl[20,13],"p4O"=1,
|
|
285 "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],
|
|
286 "p5S"=Abl[16,14],"p5T"=Abl[17,14],"p5V"=Abl[18,14],"p5W"=Abl[19,14],"p5Y"=Abl[20,14],"p5O"=1,
|
|
287 "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],
|
|
288 "p6S"=Abl[16,15],"p6T"=Abl[17,15],"p6V"=Abl[18,15],"p6W"=Abl[19,15],"p6Y"=Abl[20,15],"p6O"=1)
|
|
289
|
|
290 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],
|
|
291 "m6S"=Arg[16,3],"m6T"=Arg[17,3],"m6V"=Arg[18,3],"m6W"=Arg[19,3],"m6Y"=Arg[20,3],"m6O"=1,
|
|
292 "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],
|
|
293 "m5S"=Arg[16,4],"m5T"=Arg[17,4],"m5V"=Arg[18,4],"m5W"=Arg[19,4],"m5Y"=Arg[20,4],"m5O"=1,
|
|
294 "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],
|
|
295 "m4S"=Arg[16,5],"m4T"=Arg[17,5],"m4V"=Arg[18,5],"m4W"=Arg[19,5],"m4Y"=Arg[20,5],"m4O"=1,
|
|
296 "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],
|
|
297 "m3S"=Arg[16,6],"m3T"=Arg[17,6],"m3V"=Arg[18,6],"m3W"=Arg[19,6],"m3Y"=Arg[20,6],"m3O"=1,
|
|
298 "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],
|
|
299 "m2S"=Arg[16,7],"m2T"=Arg[17,7],"m2V"=Arg[18,7],"m2W"=Arg[19,7],"m2Y"=Arg[20,7],"m2O"=1,
|
|
300 "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],
|
|
301 "m1S"=Arg[16,8],"m1T"=Arg[17,8],"m1V"=Arg[18,8],"m1W"=Arg[19,8],"m1Y"=Arg[20,8],"m1O"=1,
|
|
302 "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],
|
|
303 "d0S"=Arg[16,9],"d0T"=Arg[17,9],"d0V"=Arg[18,9],"d0W"=Arg[19,9],"d0Y"=Arg[20,9],"d0O"=1,
|
|
304 "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],
|
|
305 "p1S"=Arg[16,10],"p1T"=Arg[17,10],"p1V"=Arg[18,10],"p1W"=Arg[19,10],"p1Y"=Arg[20,10],"p1O"=1,
|
|
306 "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],
|
|
307 "p2S"=Arg[16,11],"p2T"=Arg[17,11],"p2V"=Arg[18,11],"p2W"=Arg[19,11],"p2Y"=Arg[20,11],"p2O"=1,
|
|
308 "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],
|
|
309 "p3S"=Arg[16,12],"p3T"=Arg[17,12],"p3V"=Arg[18,12],"p3W"=Arg[19,12],"p3Y"=Arg[20,12],"p3O"=1,
|
|
310 "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],
|
|
311 "p4S"=Arg[16,13],"p4T"=Arg[17,13],"p4V"=Arg[18,13],"p4W"=Arg[19,13],"p4Y"=Arg[20,13],"p4O"=1,
|
|
312 "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],
|
|
313 "p5S"=Arg[16,14],"p5T"=Arg[17,14],"p5V"=Arg[18,14],"p5W"=Arg[19,14],"p5Y"=Arg[20,14],"p5O"=1,
|
|
314 "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],
|
|
315 "p6S"=Arg[16,15],"p6T"=Arg[17,15],"p6V"=Arg[18,15],"p6W"=Arg[19,15],"p6Y"=Arg[20,15],"p6O"=1)
|
|
316
|
|
317 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],
|
|
318 "m6S"=Btk[16,3],"m6T"=Btk[17,3],"m6V"=Btk[18,3],"m6W"=Btk[19,3],"m6Y"=Btk[20,3],"m6O"=1,
|
|
319 "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],
|
|
320 "m5S"=Btk[16,4],"m5T"=Btk[17,4],"m5V"=Btk[18,4],"m5W"=Btk[19,4],"m5Y"=Btk[20,4],"m5O"=1,
|
|
321 "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],
|
|
322 "m4S"=Btk[16,5],"m4T"=Btk[17,5],"m4V"=Btk[18,5],"m4W"=Btk[19,5],"m4Y"=Btk[20,5],"m4O"=1,
|
|
323 "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],
|
|
324 "m3S"=Btk[16,6],"m3T"=Btk[17,6],"m3V"=Btk[18,6],"m3W"=Btk[19,6],"m3Y"=Btk[20,6],"m3O"=1,
|
|
325 "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],
|
|
326 "m2S"=Btk[16,7],"m2T"=Btk[17,7],"m2V"=Btk[18,7],"m2W"=Btk[19,7],"m2Y"=Btk[20,7],"m2O"=1,
|
|
327 "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],
|
|
328 "m1S"=Btk[16,8],"m1T"=Btk[17,8],"m1V"=Btk[18,8],"m1W"=Btk[19,8],"m1Y"=Btk[20,8],"m1O"=1,
|
|
329 "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],
|
|
330 "d0S"=Btk[16,9],"d0T"=Btk[17,9],"d0V"=Btk[18,9],"d0W"=Btk[19,9],"d0Y"=Btk[20,9],"d0O"=1,
|
|
331 "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],
|
|
332 "p1S"=Btk[16,10],"p1T"=Btk[17,10],"p1V"=Btk[18,10],"p1W"=Btk[19,10],"p1Y"=Btk[20,10],"p1O"=1,
|
|
333 "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],
|
|
334 "p2S"=Btk[16,11],"p2T"=Btk[17,11],"p2V"=Btk[18,11],"p2W"=Btk[19,11],"p2Y"=Btk[20,11],"p2O"=1,
|
|
335 "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],
|
|
336 "p3S"=Btk[16,12],"p3T"=Btk[17,12],"p3V"=Btk[18,12],"p3W"=Btk[19,12],"p3Y"=Btk[20,12],"p3O"=1,
|
|
337 "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],
|
|
338 "p4S"=Btk[16,13],"p4T"=Btk[17,13],"p4V"=Btk[18,13],"p4W"=Btk[19,13],"p4Y"=Btk[20,13],"p4O"=1,
|
|
339 "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],
|
|
340 "p5S"=Btk[16,14],"p5T"=Btk[17,14],"p5V"=Btk[18,14],"p5W"=Btk[19,14],"p5Y"=Btk[20,14],"p5O"=1,
|
|
341 "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],
|
|
342 "p6S"=Btk[16,15],"p6T"=Btk[17,15],"p6V"=Btk[18,15],"p6W"=Btk[19,15],"p6Y"=Btk[20,15],"p6O"=1)
|
|
343
|
|
344 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],
|
|
345 "m6S"=Csk[16,3],"m6T"=Csk[17,3],"m6V"=Csk[18,3],"m6W"=Csk[19,3],"m6Y"=Csk[20,3],"m6O"=1,
|
|
346 "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],
|
|
347 "m5S"=Csk[16,4],"m5T"=Csk[17,4],"m5V"=Csk[18,4],"m5W"=Csk[19,4],"m5Y"=Csk[20,4],"m5O"=1,
|
|
348 "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],
|
|
349 "m4S"=Csk[16,5],"m4T"=Csk[17,5],"m4V"=Csk[18,5],"m4W"=Csk[19,5],"m4Y"=Csk[20,5],"m4O"=1,
|
|
350 "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],
|
|
351 "m3S"=Csk[16,6],"m3T"=Csk[17,6],"m3V"=Csk[18,6],"m3W"=Csk[19,6],"m3Y"=Csk[20,6],"m3O"=1,
|
|
352 "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],
|
|
353 "m2S"=Csk[16,7],"m2T"=Csk[17,7],"m2V"=Csk[18,7],"m2W"=Csk[19,7],"m2Y"=Csk[20,7],"m2O"=1,
|
|
354 "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],
|
|
355 "m1S"=Csk[16,8],"m1T"=Csk[17,8],"m1V"=Csk[18,8],"m1W"=Csk[19,8],"m1Y"=Csk[20,8],"m1O"=1,
|
|
356 "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],
|
|
357 "d0S"=Csk[16,9],"d0T"=Csk[17,9],"d0V"=Csk[18,9],"d0W"=Csk[19,9],"d0Y"=Csk[20,9],"d0O"=1,
|
|
358 "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],
|
|
359 "p1S"=Csk[16,10],"p1T"=Csk[17,10],"p1V"=Csk[18,10],"p1W"=Csk[19,10],"p1Y"=Csk[20,10],"p1O"=1,
|
|
360 "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],
|
|
361 "p2S"=Csk[16,11],"p2T"=Csk[17,11],"p2V"=Csk[18,11],"p2W"=Csk[19,11],"p2Y"=Csk[20,11],"p2O"=1,
|
|
362 "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],
|
|
363 "p3S"=Csk[16,12],"p3T"=Csk[17,12],"p3V"=Csk[18,12],"p3W"=Csk[19,12],"p3Y"=Csk[20,12],"p3O"=1,
|
|
364 "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],
|
|
365 "p4S"=Csk[16,13],"p4T"=Csk[17,13],"p4V"=Csk[18,13],"p4W"=Csk[19,13],"p4Y"=Csk[20,13],"p4O"=1,
|
|
366 "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],
|
|
367 "p5S"=Csk[16,14],"p5T"=Csk[17,14],"p5V"=Csk[18,14],"p5W"=Csk[19,14],"p5Y"=Csk[20,14],"p5O"=1,
|
|
368 "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],
|
|
369 "p6S"=Csk[16,15],"p6T"=Csk[17,15],"p6V"=Csk[18,15],"p6W"=Csk[19,15],"p6Y"=Csk[20,15],"p6O"=1)
|
|
370
|
|
371 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],
|
|
372 "m6S"=Fyn[16,3],"m6T"=Fyn[17,3],"m6V"=Fyn[18,3],"m6W"=Fyn[19,3],"m6Y"=Fyn[20,3],"m6O"=1,
|
|
373 "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],
|
|
374 "m5S"=Fyn[16,4],"m5T"=Fyn[17,4],"m5V"=Fyn[18,4],"m5W"=Fyn[19,4],"m5Y"=Fyn[20,4],"m5O"=1,
|
|
375 "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],
|
|
376 "m4S"=Fyn[16,5],"m4T"=Fyn[17,5],"m4V"=Fyn[18,5],"m4W"=Fyn[19,5],"m4Y"=Fyn[20,5],"m4O"=1,
|
|
377 "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],
|
|
378 "m3S"=Fyn[16,6],"m3T"=Fyn[17,6],"m3V"=Fyn[18,6],"m3W"=Fyn[19,6],"m3Y"=Fyn[20,6],"m3O"=1,
|
|
379 "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],
|
|
380 "m2S"=Fyn[16,7],"m2T"=Fyn[17,7],"m2V"=Fyn[18,7],"m2W"=Fyn[19,7],"m2Y"=Fyn[20,7],"m2O"=1,
|
|
381 "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],
|
|
382 "m1S"=Fyn[16,8],"m1T"=Fyn[17,8],"m1V"=Fyn[18,8],"m1W"=Fyn[19,8],"m1Y"=Fyn[20,8],"m1O"=1,
|
|
383 "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],
|
|
384 "d0S"=Fyn[16,9],"d0T"=Fyn[17,9],"d0V"=Fyn[18,9],"d0W"=Fyn[19,9],"d0Y"=Fyn[20,9],"d0O"=1,
|
|
385 "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],
|
|
386 "p1S"=Fyn[16,10],"p1T"=Fyn[17,10],"p1V"=Fyn[18,10],"p1W"=Fyn[19,10],"p1Y"=Fyn[20,10],"p1O"=1,
|
|
387 "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],
|
|
388 "p2S"=Fyn[16,11],"p2T"=Fyn[17,11],"p2V"=Fyn[18,11],"p2W"=Fyn[19,11],"p2Y"=Fyn[20,11],"p2O"=1,
|
|
389 "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],
|
|
390 "p3S"=Fyn[16,12],"p3T"=Fyn[17,12],"p3V"=Fyn[18,12],"p3W"=Fyn[19,12],"p3Y"=Fyn[20,12],"p3O"=1,
|
|
391 "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],
|
|
392 "p4S"=Fyn[16,13],"p4T"=Fyn[17,13],"p4V"=Fyn[18,13],"p4W"=Fyn[19,13],"p4Y"=Fyn[20,13],"p4O"=1,
|
|
393 "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],
|
|
394 "p5S"=Fyn[16,14],"p5T"=Fyn[17,14],"p5V"=Fyn[18,14],"p5W"=Fyn[19,14],"p5Y"=Fyn[20,14],"p5O"=1,
|
|
395 "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],
|
|
396 "p6S"=Fyn[16,15],"p6T"=Fyn[17,15],"p6V"=Fyn[18,15],"p6W"=Fyn[19,15],"p6Y"=Fyn[20,15],"p6O"=1)
|
|
397
|
|
398 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],
|
|
399 "m6S"=Hck[16,3],"m6T"=Hck[17,3],"m6V"=Hck[18,3],"m6W"=Hck[19,3],"m6Y"=Hck[20,3],"m6O"=1,
|
|
400 "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],
|
|
401 "m5S"=Hck[16,4],"m5T"=Hck[17,4],"m5V"=Hck[18,4],"m5W"=Hck[19,4],"m5Y"=Hck[20,4],"m5O"=1,
|
|
402 "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],
|
|
403 "m4S"=Hck[16,5],"m4T"=Hck[17,5],"m4V"=Hck[18,5],"m4W"=Hck[19,5],"m4Y"=Hck[20,5],"m4O"=1,
|
|
404 "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],
|
|
405 "m3S"=Hck[16,6],"m3T"=Hck[17,6],"m3V"=Hck[18,6],"m3W"=Hck[19,6],"m3Y"=Hck[20,6],"m3O"=1,
|
|
406 "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],
|
|
407 "m2S"=Hck[16,7],"m2T"=Hck[17,7],"m2V"=Hck[18,7],"m2W"=Hck[19,7],"m2Y"=Hck[20,7],"m2O"=1,
|
|
408 "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],
|
|
409 "m1S"=Hck[16,8],"m1T"=Hck[17,8],"m1V"=Hck[18,8],"m1W"=Hck[19,8],"m1Y"=Hck[20,8],"m1O"=1,
|
|
410 "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],
|
|
411 "d0S"=Hck[16,9],"d0T"=Hck[17,9],"d0V"=Hck[18,9],"d0W"=Hck[19,9],"d0Y"=Hck[20,9],"d0O"=1,
|
|
412 "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],
|
|
413 "p1S"=Hck[16,10],"p1T"=Hck[17,10],"p1V"=Hck[18,10],"p1W"=Hck[19,10],"p1Y"=Hck[20,10],"p1O"=1,
|
|
414 "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],
|
|
415 "p2S"=Hck[16,11],"p2T"=Hck[17,11],"p2V"=Hck[18,11],"p2W"=Hck[19,11],"p2Y"=Hck[20,11],"p2O"=1,
|
|
416 "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],
|
|
417 "p3S"=Hck[16,12],"p3T"=Hck[17,12],"p3V"=Hck[18,12],"p3W"=Hck[19,12],"p3Y"=Hck[20,12],"p3O"=1,
|
|
418 "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],
|
|
419 "p4S"=Hck[16,13],"p4T"=Hck[17,13],"p4V"=Hck[18,13],"p4W"=Hck[19,13],"p4Y"=Hck[20,13],"p4O"=1,
|
|
420 "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],
|
|
421 "p5S"=Hck[16,14],"p5T"=Hck[17,14],"p5V"=Hck[18,14],"p5W"=Hck[19,14],"p5Y"=Hck[20,14],"p5O"=1,
|
|
422 "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],
|
|
423 "p6S"=Hck[16,15],"p6T"=Hck[17,15],"p6V"=Hck[18,15],"p6W"=Hck[19,15],"p6Y"=Hck[20,15],"p6O"=1)
|
|
424
|
|
425 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],
|
|
426 "m6S"=JAK2[16,3],"m6T"=JAK2[17,3],"m6V"=JAK2[18,3],"m6W"=JAK2[19,3],"m6Y"=JAK2[20,3],"m6O"=1,
|
|
427 "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],
|
|
428 "m5S"=JAK2[16,4],"m5T"=JAK2[17,4],"m5V"=JAK2[18,4],"m5W"=JAK2[19,4],"m5Y"=JAK2[20,4],"m5O"=1,
|
|
429 "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],
|
|
430 "m4S"=JAK2[16,5],"m4T"=JAK2[17,5],"m4V"=JAK2[18,5],"m4W"=JAK2[19,5],"m4Y"=JAK2[20,5],"m4O"=1,
|
|
431 "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],
|
|
432 "m3S"=JAK2[16,6],"m3T"=JAK2[17,6],"m3V"=JAK2[18,6],"m3W"=JAK2[19,6],"m3Y"=JAK2[20,6],"m3O"=1,
|
|
433 "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],
|
|
434 "m2S"=JAK2[16,7],"m2T"=JAK2[17,7],"m2V"=JAK2[18,7],"m2W"=JAK2[19,7],"m2Y"=JAK2[20,7],"m2O"=1,
|
|
435 "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],
|
|
436 "m1S"=JAK2[16,8],"m1T"=JAK2[17,8],"m1V"=JAK2[18,8],"m1W"=JAK2[19,8],"m1Y"=JAK2[20,8],"m1O"=1,
|
|
437 "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],
|
|
438 "d0S"=JAK2[16,9],"d0T"=JAK2[17,9],"d0V"=JAK2[18,9],"d0W"=JAK2[19,9],"d0Y"=JAK2[20,9],"d0O"=1,
|
|
439 "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],
|
|
440 "p1S"=JAK2[16,10],"p1T"=JAK2[17,10],"p1V"=JAK2[18,10],"p1W"=JAK2[19,10],"p1Y"=JAK2[20,10],"p1O"=1,
|
|
441 "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],
|
|
442 "p2S"=JAK2[16,11],"p2T"=JAK2[17,11],"p2V"=JAK2[18,11],"p2W"=JAK2[19,11],"p2Y"=JAK2[20,11],"p2O"=1,
|
|
443 "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],
|
|
444 "p3S"=JAK2[16,12],"p3T"=JAK2[17,12],"p3V"=JAK2[18,12],"p3W"=JAK2[19,12],"p3Y"=JAK2[20,12],"p3O"=1,
|
|
445 "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],
|
|
446 "p4S"=JAK2[16,13],"p4T"=JAK2[17,13],"p4V"=JAK2[18,13],"p4W"=JAK2[19,13],"p4Y"=JAK2[20,13],"p4O"=1,
|
|
447 "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],
|
|
448 "p5S"=JAK2[16,14],"p5T"=JAK2[17,14],"p5V"=JAK2[18,14],"p5W"=JAK2[19,14],"p5Y"=JAK2[20,14],"p5O"=1,
|
|
449 "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],
|
|
450 "p6S"=JAK2[16,15],"p6T"=JAK2[17,15],"p6V"=JAK2[18,15],"p6W"=JAK2[19,15],"p6Y"=JAK2[20,15],"p6O"=1)
|
|
451
|
|
452 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],
|
|
453 "m6S"=Lck[16,3],"m6T"=Lck[17,3],"m6V"=Lck[18,3],"m6W"=Lck[19,3],"m6Y"=Lck[20,3],"m6O"=1,
|
|
454 "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],
|
|
455 "m5S"=Lck[16,4],"m5T"=Lck[17,4],"m5V"=Lck[18,4],"m5W"=Lck[19,4],"m5Y"=Lck[20,4],"m5O"=1,
|
|
456 "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],
|
|
457 "m4S"=Lck[16,5],"m4T"=Lck[17,5],"m4V"=Lck[18,5],"m4W"=Lck[19,5],"m4Y"=Lck[20,5],"m4O"=1,
|
|
458 "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],
|
|
459 "m3S"=Lck[16,6],"m3T"=Lck[17,6],"m3V"=Lck[18,6],"m3W"=Lck[19,6],"m3Y"=Lck[20,6],"m3O"=1,
|
|
460 "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],
|
|
461 "m2S"=Lck[16,7],"m2T"=Lck[17,7],"m2V"=Lck[18,7],"m2W"=Lck[19,7],"m2Y"=Lck[20,7],"m2O"=1,
|
|
462 "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],
|
|
463 "m1S"=Lck[16,8],"m1T"=Lck[17,8],"m1V"=Lck[18,8],"m1W"=Lck[19,8],"m1Y"=Lck[20,8],"m1O"=1,
|
|
464 "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],
|
|
465 "d0S"=Lck[16,9],"d0T"=Lck[17,9],"d0V"=Lck[18,9],"d0W"=Lck[19,9],"d0Y"=Lck[20,9],"d0O"=1,
|
|
466 "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],
|
|
467 "p1S"=Lck[16,10],"p1T"=Lck[17,10],"p1V"=Lck[18,10],"p1W"=Lck[19,10],"p1Y"=Lck[20,10],"p1O"=1,
|
|
468 "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],
|
|
469 "p2S"=Lck[16,11],"p2T"=Lck[17,11],"p2V"=Lck[18,11],"p2W"=Lck[19,11],"p2Y"=Lck[20,11],"p2O"=1,
|
|
470 "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],
|
|
471 "p3S"=Lck[16,12],"p3T"=Lck[17,12],"p3V"=Lck[18,12],"p3W"=Lck[19,12],"p3Y"=Lck[20,12],"p3O"=1,
|
|
472 "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],
|
|
473 "p4S"=Lck[16,13],"p4T"=Lck[17,13],"p4V"=Lck[18,13],"p4W"=Lck[19,13],"p4Y"=Lck[20,13],"p4O"=1,
|
|
474 "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],
|
|
475 "p5S"=Lck[16,14],"p5T"=Lck[17,14],"p5V"=Lck[18,14],"p5W"=Lck[19,14],"p5Y"=Lck[20,14],"p5O"=1,
|
|
476 "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],
|
|
477 "p6S"=Lck[16,15],"p6T"=Lck[17,15],"p6V"=Lck[18,15],"p6W"=Lck[19,15],"p6Y"=Lck[20,15],"p6O"=1)
|
|
478
|
|
479 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],
|
|
480 "m6S"=Lyn[16,3],"m6T"=Lyn[17,3],"m6V"=Lyn[18,3],"m6W"=Lyn[19,3],"m6Y"=Lyn[20,3],"m6O"=1,
|
|
481 "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],
|
|
482 "m5S"=Lyn[16,4],"m5T"=Lyn[17,4],"m5V"=Lyn[18,4],"m5W"=Lyn[19,4],"m5Y"=Lyn[20,4],"m5O"=1,
|
|
483 "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],
|
|
484 "m4S"=Lyn[16,5],"m4T"=Lyn[17,5],"m4V"=Lyn[18,5],"m4W"=Lyn[19,5],"m4Y"=Lyn[20,5],"m4O"=1,
|
|
485 "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],
|
|
486 "m3S"=Lyn[16,6],"m3T"=Lyn[17,6],"m3V"=Lyn[18,6],"m3W"=Lyn[19,6],"m3Y"=Lyn[20,6],"m3O"=1,
|
|
487 "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],
|
|
488 "m2S"=Lyn[16,7],"m2T"=Lyn[17,7],"m2V"=Lyn[18,7],"m2W"=Lyn[19,7],"m2Y"=Lyn[20,7],"m2O"=1,
|
|
489 "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],
|
|
490 "m1S"=Lyn[16,8],"m1T"=Lyn[17,8],"m1V"=Lyn[18,8],"m1W"=Lyn[19,8],"m1Y"=Lyn[20,8],"m1O"=1,
|
|
491 "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],
|
|
492 "d0S"=Lyn[16,9],"d0T"=Lyn[17,9],"d0V"=Lyn[18,9],"d0W"=Lyn[19,9],"d0Y"=Lyn[20,9],"d0O"=1,
|
|
493 "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],
|
|
494 "p1S"=Lyn[16,10],"p1T"=Lyn[17,10],"p1V"=Lyn[18,10],"p1W"=Lyn[19,10],"p1Y"=Lyn[20,10],"p1O"=1,
|
|
495 "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],
|
|
496 "p2S"=Lyn[16,11],"p2T"=Lyn[17,11],"p2V"=Lyn[18,11],"p2W"=Lyn[19,11],"p2Y"=Lyn[20,11],"p2O"=1,
|
|
497 "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],
|
|
498 "p3S"=Lyn[16,12],"p3T"=Lyn[17,12],"p3V"=Lyn[18,12],"p3W"=Lyn[19,12],"p3Y"=Lyn[20,12],"p3O"=1,
|
|
499 "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],
|
|
500 "p4S"=Lyn[16,13],"p4T"=Lyn[17,13],"p4V"=Lyn[18,13],"p4W"=Lyn[19,13],"p4Y"=Lyn[20,13],"p4O"=1,
|
|
501 "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],
|
|
502 "p5S"=Lyn[16,14],"p5T"=Lyn[17,14],"p5V"=Lyn[18,14],"p5W"=Lyn[19,14],"p5Y"=Lyn[20,14],"p5O"=1,
|
|
503 "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],
|
|
504 "p6S"=Lyn[16,15],"p6T"=Lyn[17,15],"p6V"=Lyn[18,15],"p6W"=Lyn[19,15],"p6Y"=Lyn[20,15],"p6O"=1)
|
|
505
|
|
506 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],
|
|
507 "m6S"=Pyk2[16,3],"m6T"=Pyk2[17,3],"m6V"=Pyk2[18,3],"m6W"=Pyk2[19,3],"m6Y"=Pyk2[20,3],"m6O"=1,
|
|
508 "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],
|
|
509 "m5S"=Pyk2[16,4],"m5T"=Pyk2[17,4],"m5V"=Pyk2[18,4],"m5W"=Pyk2[19,4],"m5Y"=Pyk2[20,4],"m5O"=1,
|
|
510 "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],
|
|
511 "m4S"=Pyk2[16,5],"m4T"=Pyk2[17,5],"m4V"=Pyk2[18,5],"m4W"=Pyk2[19,5],"m4Y"=Pyk2[20,5],"m4O"=1,
|
|
512 "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],
|
|
513 "m3S"=Pyk2[16,6],"m3T"=Pyk2[17,6],"m3V"=Pyk2[18,6],"m3W"=Pyk2[19,6],"m3Y"=Pyk2[20,6],"m3O"=1,
|
|
514 "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],
|
|
515 "m2S"=Pyk2[16,7],"m2T"=Pyk2[17,7],"m2V"=Pyk2[18,7],"m2W"=Pyk2[19,7],"m2Y"=Pyk2[20,7],"m2O"=1,
|
|
516 "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],
|
|
517 "m1S"=Pyk2[16,8],"m1T"=Pyk2[17,8],"m1V"=Pyk2[18,8],"m1W"=Pyk2[19,8],"m1Y"=Pyk2[20,8],"m1O"=1,
|
|
518 "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],
|
|
519 "d0S"=Pyk2[16,9],"d0T"=Pyk2[17,9],"d0V"=Pyk2[18,9],"d0W"=Pyk2[19,9],"d0Y"=Pyk2[20,9],"d0O"=1,
|
|
520 "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],
|
|
521 "p1S"=Pyk2[16,10],"p1T"=Pyk2[17,10],"p1V"=Pyk2[18,10],"p1W"=Pyk2[19,10],"p1Y"=Pyk2[20,10],"p1O"=1,
|
|
522 "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],
|
|
523 "p2S"=Pyk2[16,11],"p2T"=Pyk2[17,11],"p2V"=Pyk2[18,11],"p2W"=Pyk2[19,11],"p2Y"=Pyk2[20,11],"p2O"=1,
|
|
524 "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],
|
|
525 "p3S"=Pyk2[16,12],"p3T"=Pyk2[17,12],"p3V"=Pyk2[18,12],"p3W"=Pyk2[19,12],"p3Y"=Pyk2[20,12],"p3O"=1,
|
|
526 "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],
|
|
527 "p4S"=Pyk2[16,13],"p4T"=Pyk2[17,13],"p4V"=Pyk2[18,13],"p4W"=Pyk2[19,13],"p4Y"=Pyk2[20,13],"p4O"=1,
|
|
528 "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],
|
|
529 "p5S"=Pyk2[16,14],"p5T"=Pyk2[17,14],"p5V"=Pyk2[18,14],"p5W"=Pyk2[19,14],"p5Y"=Pyk2[20,14],"p5O"=1,
|
|
530 "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],
|
|
531 "p6S"=Pyk2[16,15],"p6T"=Pyk2[17,15],"p6V"=Pyk2[18,15],"p6W"=Pyk2[19,15],"p6Y"=Pyk2[20,15],"p6O"=1)
|
|
532
|
|
533 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],
|
|
534 "m6S"=Src[16,3],"m6T"=Src[17,3],"m6V"=Src[18,3],"m6W"=Src[19,3],"m6Y"=Src[20,3],"m6O"=1,
|
|
535 "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],
|
|
536 "m5S"=Src[16,4],"m5T"=Src[17,4],"m5V"=Src[18,4],"m5W"=Src[19,4],"m5Y"=Src[20,4],"m5O"=1,
|
|
537 "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],
|
|
538 "m4S"=Src[16,5],"m4T"=Src[17,5],"m4V"=Src[18,5],"m4W"=Src[19,5],"m4Y"=Src[20,5],"m4O"=1,
|
|
539 "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],
|
|
540 "m3S"=Src[16,6],"m3T"=Src[17,6],"m3V"=Src[18,6],"m3W"=Src[19,6],"m3Y"=Src[20,6],"m3O"=1,
|
|
541 "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],
|
|
542 "m2S"=Src[16,7],"m2T"=Src[17,7],"m2V"=Src[18,7],"m2W"=Src[19,7],"m2Y"=Src[20,7],"m2O"=1,
|
|
543 "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],
|
|
544 "m1S"=Src[16,8],"m1T"=Src[17,8],"m1V"=Src[18,8],"m1W"=Src[19,8],"m1Y"=Src[20,8],"m1O"=1,
|
|
545 "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],
|
|
546 "d0S"=Src[16,9],"d0T"=Src[17,9],"d0V"=Src[18,9],"d0W"=Src[19,9],"d0Y"=Src[20,9],"d0O"=1,
|
|
547 "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],
|
|
548 "p1S"=Src[16,10],"p1T"=Src[17,10],"p1V"=Src[18,10],"p1W"=Src[19,10],"p1Y"=Src[20,10],"p1O"=1,
|
|
549 "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],
|
|
550 "p2S"=Src[16,11],"p2T"=Src[17,11],"p2V"=Src[18,11],"p2W"=Src[19,11],"p2Y"=Src[20,11],"p2O"=1,
|
|
551 "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],
|
|
552 "p3S"=Src[16,12],"p3T"=Src[17,12],"p3V"=Src[18,12],"p3W"=Src[19,12],"p3Y"=Src[20,12],"p3O"=1,
|
|
553 "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],
|
|
554 "p4S"=Src[16,13],"p4T"=Src[17,13],"p4V"=Src[18,13],"p4W"=Src[19,13],"p4Y"=Src[20,13],"p4O"=1,
|
|
555 "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],
|
|
556 "p5S"=Src[16,14],"p5T"=Src[17,14],"p5V"=Src[18,14],"p5W"=Src[19,14],"p5Y"=Src[20,14],"p5O"=1,
|
|
557 "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],
|
|
558 "p6S"=Src[16,15],"p6T"=Src[17,15],"p6V"=Src[18,15],"p6W"=Src[19,15],"p6Y"=Src[20,15],"p6O"=1)
|
|
559
|
|
560 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],
|
|
561 "m6S"=Syk[16,3],"m6T"=Syk[17,3],"m6V"=Syk[18,3],"m6W"=Syk[19,3],"m6Y"=Syk[20,3],"m6O"=1,
|
|
562 "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],
|
|
563 "m5S"=Syk[16,4],"m5T"=Syk[17,4],"m5V"=Syk[18,4],"m5W"=Syk[19,4],"m5Y"=Syk[20,4],"m5O"=1,
|
|
564 "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],
|
|
565 "m4S"=Syk[16,5],"m4T"=Syk[17,5],"m4V"=Syk[18,5],"m4W"=Syk[19,5],"m4Y"=Syk[20,5],"m4O"=1,
|
|
566 "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],
|
|
567 "m3S"=Syk[16,6],"m3T"=Syk[17,6],"m3V"=Syk[18,6],"m3W"=Syk[19,6],"m3Y"=Syk[20,6],"m3O"=1,
|
|
568 "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],
|
|
569 "m2S"=Syk[16,7],"m2T"=Syk[17,7],"m2V"=Syk[18,7],"m2W"=Syk[19,7],"m2Y"=Syk[20,7],"m2O"=1,
|
|
570 "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],
|
|
571 "m1S"=Syk[16,8],"m1T"=Syk[17,8],"m1V"=Syk[18,8],"m1W"=Syk[19,8],"m1Y"=Syk[20,8],"m1O"=1,
|
|
572 "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],
|
|
573 "d0S"=Syk[16,9],"d0T"=Syk[17,9],"d0V"=Syk[18,9],"d0W"=Syk[19,9],"d0Y"=Syk[20,9],"d0O"=1,
|
|
574 "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],
|
|
575 "p1S"=Syk[16,10],"p1T"=Syk[17,10],"p1V"=Syk[18,10],"p1W"=Syk[19,10],"p1Y"=Syk[20,10],"p1O"=1,
|
|
576 "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],
|
|
577 "p2S"=Syk[16,11],"p2T"=Syk[17,11],"p2V"=Syk[18,11],"p2W"=Syk[19,11],"p2Y"=Syk[20,11],"p2O"=1,
|
|
578 "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],
|
|
579 "p3S"=Syk[16,12],"p3T"=Syk[17,12],"p3V"=Syk[18,12],"p3W"=Syk[19,12],"p3Y"=Syk[20,12],"p3O"=1,
|
|
580 "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],
|
|
581 "p4S"=Syk[16,13],"p4T"=Syk[17,13],"p4V"=Syk[18,13],"p4W"=Syk[19,13],"p4Y"=Syk[20,13],"p4O"=1,
|
|
582 "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],
|
|
583 "p5S"=Syk[16,14],"p5T"=Syk[17,14],"p5V"=Syk[18,14],"p5W"=Syk[19,14],"p5Y"=Syk[20,14],"p5O"=1,
|
|
584 "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],
|
|
585 "p6S"=Syk[16,15],"p6T"=Syk[17,15],"p6V"=Syk[18,15],"p6W"=Syk[19,15],"p6Y"=Syk[20,15],"p6O"=1)
|
|
586
|
|
587 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],
|
|
588 "m6S"=Yes[16,3],"m6T"=Yes[17,3],"m6V"=Yes[18,3],"m6W"=Yes[19,3],"m6Y"=Yes[20,3],"m6O"=1,
|
|
589 "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],
|
|
590 "m5S"=Yes[16,4],"m5T"=Yes[17,4],"m5V"=Yes[18,4],"m5W"=Yes[19,4],"m5Y"=Yes[20,4],"m5O"=1,
|
|
591 "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],
|
|
592 "m4S"=Yes[16,5],"m4T"=Yes[17,5],"m4V"=Yes[18,5],"m4W"=Yes[19,5],"m4Y"=Yes[20,5],"m4O"=1,
|
|
593 "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],
|
|
594 "m3S"=Yes[16,6],"m3T"=Yes[17,6],"m3V"=Yes[18,6],"m3W"=Yes[19,6],"m3Y"=Yes[20,6],"m3O"=1,
|
|
595 "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],
|
|
596 "m2S"=Yes[16,7],"m2T"=Yes[17,7],"m2V"=Yes[18,7],"m2W"=Yes[19,7],"m2Y"=Yes[20,7],"m2O"=1,
|
|
597 "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],
|
|
598 "m1S"=Yes[16,8],"m1T"=Yes[17,8],"m1V"=Yes[18,8],"m1W"=Yes[19,8],"m1Y"=Yes[20,8],"m1O"=1,
|
|
599 "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],
|
|
600 "d0S"=Yes[16,9],"d0T"=Yes[17,9],"d0V"=Yes[18,9],"d0W"=Yes[19,9],"d0Y"=Yes[20,9],"d0O"=1,
|
|
601 "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],
|
|
602 "p1S"=Yes[16,10],"p1T"=Yes[17,10],"p1V"=Yes[18,10],"p1W"=Yes[19,10],"p1Y"=Yes[20,10],"p1O"=1,
|
|
603 "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],
|
|
604 "p2S"=Yes[16,11],"p2T"=Yes[17,11],"p2V"=Yes[18,11],"p2W"=Yes[19,11],"p2Y"=Yes[20,11],"p2O"=1,
|
|
605 "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],
|
|
606 "p3S"=Yes[16,12],"p3T"=Yes[17,12],"p3V"=Yes[18,12],"p3W"=Yes[19,12],"p3Y"=Yes[20,12],"p3O"=1,
|
|
607 "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],
|
|
608 "p4S"=Yes[16,13],"p4T"=Yes[17,13],"p4V"=Yes[18,13],"p4W"=Yes[19,13],"p4Y"=Yes[20,13],"p4O"=1,
|
|
609 "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],
|
|
610 "p5S"=Yes[16,14],"p5T"=Yes[17,14],"p5V"=Yes[18,14],"p5W"=Yes[19,14],"p5Y"=Yes[20,14],"p5O"=1,
|
|
611 "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],
|
|
612 "p6S"=Yes[16,15],"p6T"=Yes[17,15],"p6V"=Yes[18,15],"p6W"=Yes[19,15],"p6Y"=Yes[20,15],"p6O"=1)
|
|
613
|
|
614 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],
|
|
615 "m6S"=FLT3[16,3],"m6T"=FLT3[17,3],"m6V"=FLT3[18,3],"m6W"=FLT3[19,3],"m6Y"=FLT3[20,3],"m6O"=1,
|
|
616 "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],
|
|
617 "m5S"=FLT3[16,4],"m5T"=FLT3[17,4],"m5V"=FLT3[18,4],"m5W"=FLT3[19,4],"m5Y"=FLT3[20,4],"m5O"=1,
|
|
618 "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],
|
|
619 "m4S"=FLT3[16,5],"m4T"=FLT3[17,5],"m4V"=FLT3[18,5],"m4W"=FLT3[19,5],"m4Y"=FLT3[20,5],"m4O"=1,
|
|
620 "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],
|
|
621 "m3S"=FLT3[16,6],"m3T"=FLT3[17,6],"m3V"=FLT3[18,6],"m3W"=FLT3[19,6],"m3Y"=FLT3[20,6],"m3O"=1,
|
|
622 "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],
|
|
623 "m2S"=FLT3[16,7],"m2T"=FLT3[17,7],"m2V"=FLT3[18,7],"m2W"=FLT3[19,7],"m2Y"=FLT3[20,7],"m2O"=1,
|
|
624 "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],
|
|
625 "m1S"=FLT3[16,8],"m1T"=FLT3[17,8],"m1V"=FLT3[18,8],"m1W"=FLT3[19,8],"m1Y"=FLT3[20,8],"m1O"=1,
|
|
626 "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],
|
|
627 "d0S"=FLT3[16,9],"d0T"=FLT3[17,9],"d0V"=FLT3[18,9],"d0W"=FLT3[19,9],"d0Y"=FLT3[20,9],"d0O"=1,
|
|
628 "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],
|
|
629 "p1S"=FLT3[16,10],"p1T"=FLT3[17,10],"p1V"=FLT3[18,10],"p1W"=FLT3[19,10],"p1Y"=FLT3[20,10],"p1O"=1,
|
|
630 "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],
|
|
631 "p2S"=FLT3[16,11],"p2T"=FLT3[17,11],"p2V"=FLT3[18,11],"p2W"=FLT3[19,11],"p2Y"=FLT3[20,11],"p2O"=1,
|
|
632 "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],
|
|
633 "p3S"=FLT3[16,12],"p3T"=FLT3[17,12],"p3V"=FLT3[18,12],"p3W"=FLT3[19,12],"p3Y"=FLT3[20,12],"p3O"=1,
|
|
634 "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],
|
|
635 "p4S"=FLT3[16,13],"p4T"=FLT3[17,13],"p4V"=FLT3[18,13],"p4W"=FLT3[19,13],"p4Y"=FLT3[20,13],"p4O"=1,
|
|
636 "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],
|
|
637 "p5S"=FLT3[16,14],"p5T"=FLT3[17,14],"p5V"=FLT3[18,14],"p5W"=FLT3[19,14],"p5Y"=FLT3[20,14],"p5O"=1,
|
|
638 "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],
|
|
639 "p6S"=FLT3[16,15],"p6T"=FLT3[17,15],"p6V"=FLT3[18,15],"p6W"=FLT3[19,15],"p6Y"=FLT3[20,15],"p6O"=1)
|
|
640
|
|
641 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],
|
|
642 "m6S"=ALK[16,3],"m6T"=ALK[17,3],"m6V"=ALK[18,3],"m6W"=ALK[19,3],"m6Y"=ALK[20,3],"m6O"=1,
|
|
643 "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],
|
|
644 "m5S"=ALK[16,4],"m5T"=ALK[17,4],"m5V"=ALK[18,4],"m5W"=ALK[19,4],"m5Y"=ALK[20,4],"m5O"=1,
|
|
645 "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],
|
|
646 "m4S"=ALK[16,5],"m4T"=ALK[17,5],"m4V"=ALK[18,5],"m4W"=ALK[19,5],"m4Y"=ALK[20,5],"m4O"=1,
|
|
647 "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],
|
|
648 "m3S"=ALK[16,6],"m3T"=ALK[17,6],"m3V"=ALK[18,6],"m3W"=ALK[19,6],"m3Y"=ALK[20,6],"m3O"=1,
|
|
649 "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],
|
|
650 "m2S"=ALK[16,7],"m2T"=ALK[17,7],"m2V"=ALK[18,7],"m2W"=ALK[19,7],"m2Y"=ALK[20,7],"m2O"=1,
|
|
651 "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],
|
|
652 "m1S"=ALK[16,8],"m1T"=ALK[17,8],"m1V"=ALK[18,8],"m1W"=ALK[19,8],"m1Y"=ALK[20,8],"m1O"=1,
|
|
653 "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],
|
|
654 "d0S"=ALK[16,9],"d0T"=ALK[17,9],"d0V"=ALK[18,9],"d0W"=ALK[19,9],"d0Y"=ALK[20,9],"d0O"=1,
|
|
655 "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],
|
|
656 "p1S"=ALK[16,10],"p1T"=ALK[17,10],"p1V"=ALK[18,10],"p1W"=ALK[19,10],"p1Y"=ALK[20,10],"p1O"=1,
|
|
657 "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],
|
|
658 "p2S"=ALK[16,11],"p2T"=ALK[17,11],"p2V"=ALK[18,11],"p2W"=ALK[19,11],"p2Y"=ALK[20,11],"p2O"=1,
|
|
659 "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],
|
|
660 "p3S"=ALK[16,12],"p3T"=ALK[17,12],"p3V"=ALK[18,12],"p3W"=ALK[19,12],"p3Y"=ALK[20,12],"p3O"=1,
|
|
661 "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],
|
|
662 "p4S"=ALK[16,13],"p4T"=ALK[17,13],"p4V"=ALK[18,13],"p4W"=ALK[19,13],"p4Y"=ALK[20,13],"p4O"=1,
|
|
663 "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],
|
|
664 "p5S"=ALK[16,14],"p5T"=ALK[17,14],"p5V"=ALK[18,14],"p5W"=ALK[19,14],"p5Y"=ALK[20,14],"p5O"=1,
|
|
665 "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],
|
|
666 "p6S"=ALK[16,15],"p6T"=ALK[17,15],"p6V"=ALK[18,15],"p6W"=ALK[19,15],"p6Y"=ALK[20,15],"p6O"=1)
|
|
667
|
|
668 }
|
|
669 #this creates the apply functions that are used below. Each amino acid at each position gets assigned a value from an EPM table
|
|
670
|
|
671
|
|
672 #this part turns amino acids at positions into values from the EPM tables using the above apply function
|
|
673 if(1==1){
|
|
674 ThisKinTablem6<-sapply(pPositionm6, function(x) ThisKinTable_props[x])
|
|
675 ThisKinTablem5<-sapply(pPositionm5, function(x) ThisKinTable_props[x])
|
|
676 ThisKinTablem4<-sapply(pPositionm4, function(x) ThisKinTable_props[x])
|
|
677 ThisKinTablem3<-sapply(pPositionm3, function(x) ThisKinTable_props[x])
|
|
678 ThisKinTablem2<-sapply(pPositionm2, function(x) ThisKinTable_props[x])
|
|
679 ThisKinTablem1<-sapply(pPositionm1, function(x) ThisKinTable_props[x])
|
|
680 ThisKinTabled0<-sapply(pPositiond0, function(x) ThisKinTable_props[x])
|
|
681 ThisKinTablep1<-sapply(pPositionp1, function(x) ThisKinTable_props[x])
|
|
682 ThisKinTablep2<-sapply(pPositionp2, function(x) ThisKinTable_props[x])
|
|
683 ThisKinTablep3<-sapply(pPositionp3, function(x) ThisKinTable_props[x])
|
|
684 ThisKinTablep4<-sapply(pPositionp4, function(x) ThisKinTable_props[x])
|
|
685 ThisKinTablep5<-sapply(pPositionp5, function(x) ThisKinTable_props[x])
|
|
686 ThisKinTablep6<-sapply(pPositionp6, function(x) ThisKinTable_props[x])
|
|
687
|
|
688 Ablm6<-sapply(pPositionm6, function(x) Abl_props[x])
|
|
689 Ablm5<-sapply(pPositionm5, function(x) Abl_props[x])
|
|
690 Ablm4<-sapply(pPositionm4, function(x) Abl_props[x])
|
|
691 Ablm3<-sapply(pPositionm3, function(x) Abl_props[x])
|
|
692 Ablm2<-sapply(pPositionm2, function(x) Abl_props[x])
|
|
693 Ablm1<-sapply(pPositionm1, function(x) Abl_props[x])
|
|
694 Abld0<-sapply(pPositiond0, function(x) Abl_props[x])
|
|
695 Ablp1<-sapply(pPositionp1, function(x) Abl_props[x])
|
|
696 Ablp2<-sapply(pPositionp2, function(x) Abl_props[x])
|
|
697 Ablp3<-sapply(pPositionp3, function(x) Abl_props[x])
|
|
698 Ablp4<-sapply(pPositionp4, function(x) Abl_props[x])
|
|
699 Ablp5<-sapply(pPositionp5, function(x) Abl_props[x])
|
|
700 Ablp6<-sapply(pPositionp6, function(x) Abl_props[x])
|
|
701
|
|
702 Argm6<-sapply(pPositionm6, function(x) Arg_props[x])
|
|
703 Argm5<-sapply(pPositionm5, function(x) Arg_props[x])
|
|
704 Argm4<-sapply(pPositionm4, function(x) Arg_props[x])
|
|
705 Argm3<-sapply(pPositionm3, function(x) Arg_props[x])
|
|
706 Argm2<-sapply(pPositionm2, function(x) Arg_props[x])
|
|
707 Argm1<-sapply(pPositionm1, function(x) Arg_props[x])
|
|
708 Argd0<-sapply(pPositiond0, function(x) Arg_props[x])
|
|
709 Argp1<-sapply(pPositionp1, function(x) Arg_props[x])
|
|
710 Argp2<-sapply(pPositionp2, function(x) Arg_props[x])
|
|
711 Argp3<-sapply(pPositionp3, function(x) Arg_props[x])
|
|
712 Argp4<-sapply(pPositionp4, function(x) Arg_props[x])
|
|
713 Argp5<-sapply(pPositionp5, function(x) Arg_props[x])
|
|
714 Argp6<-sapply(pPositionp6, function(x) Arg_props[x])
|
|
715
|
|
716 Btkm6<-sapply(pPositionm6, function(x) Btk_props[x])
|
|
717 Btkm5<-sapply(pPositionm5, function(x) Btk_props[x])
|
|
718 Btkm4<-sapply(pPositionm4, function(x) Btk_props[x])
|
|
719 Btkm3<-sapply(pPositionm3, function(x) Btk_props[x])
|
|
720 Btkm2<-sapply(pPositionm2, function(x) Btk_props[x])
|
|
721 Btkm1<-sapply(pPositionm1, function(x) Btk_props[x])
|
|
722 Btkd0<-sapply(pPositiond0, function(x) Btk_props[x])
|
|
723 Btkp1<-sapply(pPositionp1, function(x) Btk_props[x])
|
|
724 Btkp2<-sapply(pPositionp2, function(x) Btk_props[x])
|
|
725 Btkp3<-sapply(pPositionp3, function(x) Btk_props[x])
|
|
726 Btkp4<-sapply(pPositionp4, function(x) Btk_props[x])
|
|
727 Btkp5<-sapply(pPositionp5, function(x) Btk_props[x])
|
|
728 Btkp6<-sapply(pPositionp6, function(x) Btk_props[x])
|
|
729
|
|
730 Cskm6<-sapply(pPositionm6, function(x) Csk_props[x])
|
|
731 Cskm5<-sapply(pPositionm5, function(x) Csk_props[x])
|
|
732 Cskm4<-sapply(pPositionm4, function(x) Csk_props[x])
|
|
733 Cskm3<-sapply(pPositionm3, function(x) Csk_props[x])
|
|
734 Cskm2<-sapply(pPositionm2, function(x) Csk_props[x])
|
|
735 Cskm1<-sapply(pPositionm1, function(x) Csk_props[x])
|
|
736 Cskd0<-sapply(pPositiond0, function(x) Csk_props[x])
|
|
737 Cskp1<-sapply(pPositionp1, function(x) Csk_props[x])
|
|
738 Cskp2<-sapply(pPositionp2, function(x) Csk_props[x])
|
|
739 Cskp3<-sapply(pPositionp3, function(x) Csk_props[x])
|
|
740 Cskp4<-sapply(pPositionp4, function(x) Csk_props[x])
|
|
741 Cskp5<-sapply(pPositionp5, function(x) Csk_props[x])
|
|
742 Cskp6<-sapply(pPositionp6, function(x) Csk_props[x])
|
|
743
|
|
744 Fynm6<-sapply(pPositionm6, function(x) Fyn_props[x])
|
|
745 Fynm5<-sapply(pPositionm5, function(x) Fyn_props[x])
|
|
746 Fynm4<-sapply(pPositionm4, function(x) Fyn_props[x])
|
|
747 Fynm3<-sapply(pPositionm3, function(x) Fyn_props[x])
|
|
748 Fynm2<-sapply(pPositionm2, function(x) Fyn_props[x])
|
|
749 Fynm1<-sapply(pPositionm1, function(x) Fyn_props[x])
|
|
750 Fynd0<-sapply(pPositiond0, function(x) Fyn_props[x])
|
|
751 Fynp1<-sapply(pPositionp1, function(x) Fyn_props[x])
|
|
752 Fynp2<-sapply(pPositionp2, function(x) Fyn_props[x])
|
|
753 Fynp3<-sapply(pPositionp3, function(x) Fyn_props[x])
|
|
754 Fynp4<-sapply(pPositionp4, function(x) Fyn_props[x])
|
|
755 Fynp5<-sapply(pPositionp5, function(x) Fyn_props[x])
|
|
756 Fynp6<-sapply(pPositionp6, function(x) Fyn_props[x])
|
|
757
|
|
758 Hckm6<-sapply(pPositionm6, function(x) Hck_props[x])
|
|
759 Hckm5<-sapply(pPositionm5, function(x) Hck_props[x])
|
|
760 Hckm4<-sapply(pPositionm4, function(x) Hck_props[x])
|
|
761 Hckm3<-sapply(pPositionm3, function(x) Hck_props[x])
|
|
762 Hckm2<-sapply(pPositionm2, function(x) Hck_props[x])
|
|
763 Hckm1<-sapply(pPositionm1, function(x) Hck_props[x])
|
|
764 Hckd0<-sapply(pPositiond0, function(x) Hck_props[x])
|
|
765 Hckp1<-sapply(pPositionp1, function(x) Hck_props[x])
|
|
766 Hckp2<-sapply(pPositionp2, function(x) Hck_props[x])
|
|
767 Hckp3<-sapply(pPositionp3, function(x) Hck_props[x])
|
|
768 Hckp4<-sapply(pPositionp4, function(x) Hck_props[x])
|
|
769 Hckp5<-sapply(pPositionp5, function(x) Hck_props[x])
|
|
770 Hckp6<-sapply(pPositionp6, function(x) Hck_props[x])
|
|
771
|
|
772 JAK2m6<-sapply(pPositionm6, function(x) JAK2_props[x])
|
|
773 JAK2m5<-sapply(pPositionm5, function(x) JAK2_props[x])
|
|
774 JAK2m4<-sapply(pPositionm4, function(x) JAK2_props[x])
|
|
775 JAK2m3<-sapply(pPositionm3, function(x) JAK2_props[x])
|
|
776 JAK2m2<-sapply(pPositionm2, function(x) JAK2_props[x])
|
|
777 JAK2m1<-sapply(pPositionm1, function(x) JAK2_props[x])
|
|
778 JAK2d0<-sapply(pPositiond0, function(x) JAK2_props[x])
|
|
779 JAK2p1<-sapply(pPositionp1, function(x) JAK2_props[x])
|
|
780 JAK2p2<-sapply(pPositionp2, function(x) JAK2_props[x])
|
|
781 JAK2p3<-sapply(pPositionp3, function(x) JAK2_props[x])
|
|
782 JAK2p4<-sapply(pPositionp4, function(x) JAK2_props[x])
|
|
783 JAK2p5<-sapply(pPositionp5, function(x) JAK2_props[x])
|
|
784 JAK2p6<-sapply(pPositionp6, function(x) JAK2_props[x])
|
|
785
|
|
786 Lckm6<-sapply(pPositionm6, function(x) Lck_props[x])
|
|
787 Lckm5<-sapply(pPositionm5, function(x) Lck_props[x])
|
|
788 Lckm4<-sapply(pPositionm4, function(x) Lck_props[x])
|
|
789 Lckm3<-sapply(pPositionm3, function(x) Lck_props[x])
|
|
790 Lckm2<-sapply(pPositionm2, function(x) Lck_props[x])
|
|
791 Lckm1<-sapply(pPositionm1, function(x) Lck_props[x])
|
|
792 Lckd0<-sapply(pPositiond0, function(x) Lck_props[x])
|
|
793 Lckp1<-sapply(pPositionp1, function(x) Lck_props[x])
|
|
794 Lckp2<-sapply(pPositionp2, function(x) Lck_props[x])
|
|
795 Lckp3<-sapply(pPositionp3, function(x) Lck_props[x])
|
|
796 Lckp4<-sapply(pPositionp4, function(x) Lck_props[x])
|
|
797 Lckp5<-sapply(pPositionp5, function(x) Lck_props[x])
|
|
798 Lckp6<-sapply(pPositionp6, function(x) Lck_props[x])
|
|
799
|
|
800 Lynm6<-sapply(pPositionm6, function(x) Lyn_props[x])
|
|
801 Lynm5<-sapply(pPositionm5, function(x) Lyn_props[x])
|
|
802 Lynm4<-sapply(pPositionm4, function(x) Lyn_props[x])
|
|
803 Lynm3<-sapply(pPositionm3, function(x) Lyn_props[x])
|
|
804 Lynm2<-sapply(pPositionm2, function(x) Lyn_props[x])
|
|
805 Lynm1<-sapply(pPositionm1, function(x) Lyn_props[x])
|
|
806 Lynd0<-sapply(pPositiond0, function(x) Lyn_props[x])
|
|
807 Lynp1<-sapply(pPositionp1, function(x) Lyn_props[x])
|
|
808 Lynp2<-sapply(pPositionp2, function(x) Lyn_props[x])
|
|
809 Lynp3<-sapply(pPositionp3, function(x) Lyn_props[x])
|
|
810 Lynp4<-sapply(pPositionp4, function(x) Lyn_props[x])
|
|
811 Lynp5<-sapply(pPositionp5, function(x) Lyn_props[x])
|
|
812 Lynp6<-sapply(pPositionp6, function(x) Lyn_props[x])
|
|
813
|
|
814 Pyk2m6<-sapply(pPositionm6, function(x) Pyk2_props[x])
|
|
815 Pyk2m5<-sapply(pPositionm5, function(x) Pyk2_props[x])
|
|
816 Pyk2m4<-sapply(pPositionm4, function(x) Pyk2_props[x])
|
|
817 Pyk2m3<-sapply(pPositionm3, function(x) Pyk2_props[x])
|
|
818 Pyk2m2<-sapply(pPositionm2, function(x) Pyk2_props[x])
|
|
819 Pyk2m1<-sapply(pPositionm1, function(x) Pyk2_props[x])
|
|
820 Pyk2d0<-sapply(pPositiond0, function(x) Pyk2_props[x])
|
|
821 Pyk2p1<-sapply(pPositionp1, function(x) Pyk2_props[x])
|
|
822 Pyk2p2<-sapply(pPositionp2, function(x) Pyk2_props[x])
|
|
823 Pyk2p3<-sapply(pPositionp3, function(x) Pyk2_props[x])
|
|
824 Pyk2p4<-sapply(pPositionp4, function(x) Pyk2_props[x])
|
|
825 Pyk2p5<-sapply(pPositionp5, function(x) Pyk2_props[x])
|
|
826 Pyk2p6<-sapply(pPositionp6, function(x) Pyk2_props[x])
|
|
827
|
|
828 Srcm6<-sapply(pPositionm6, function(x) Src_props[x])
|
|
829 Srcm5<-sapply(pPositionm5, function(x) Src_props[x])
|
|
830 Srcm4<-sapply(pPositionm4, function(x) Src_props[x])
|
|
831 Srcm3<-sapply(pPositionm3, function(x) Src_props[x])
|
|
832 Srcm2<-sapply(pPositionm2, function(x) Src_props[x])
|
|
833 Srcm1<-sapply(pPositionm1, function(x) Src_props[x])
|
|
834 Srcd0<-sapply(pPositiond0, function(x) Src_props[x])
|
|
835 Srcp1<-sapply(pPositionp1, function(x) Src_props[x])
|
|
836 Srcp2<-sapply(pPositionp2, function(x) Src_props[x])
|
|
837 Srcp3<-sapply(pPositionp3, function(x) Src_props[x])
|
|
838 Srcp4<-sapply(pPositionp4, function(x) Src_props[x])
|
|
839 Srcp5<-sapply(pPositionp5, function(x) Src_props[x])
|
|
840 Srcp6<-sapply(pPositionp6, function(x) Src_props[x])
|
|
841
|
|
842 Sykm6<-sapply(pPositionm6, function(x) Syk_props[x])
|
|
843 Sykm5<-sapply(pPositionm5, function(x) Syk_props[x])
|
|
844 Sykm4<-sapply(pPositionm4, function(x) Syk_props[x])
|
|
845 Sykm3<-sapply(pPositionm3, function(x) Syk_props[x])
|
|
846 Sykm2<-sapply(pPositionm2, function(x) Syk_props[x])
|
|
847 Sykm1<-sapply(pPositionm1, function(x) Syk_props[x])
|
|
848 Sykd0<-sapply(pPositiond0, function(x) Syk_props[x])
|
|
849 Sykp1<-sapply(pPositionp1, function(x) Syk_props[x])
|
|
850 Sykp2<-sapply(pPositionp2, function(x) Syk_props[x])
|
|
851 Sykp3<-sapply(pPositionp3, function(x) Syk_props[x])
|
|
852 Sykp4<-sapply(pPositionp4, function(x) Syk_props[x])
|
|
853 Sykp5<-sapply(pPositionp5, function(x) Syk_props[x])
|
|
854 Sykp6<-sapply(pPositionp6, function(x) Syk_props[x])
|
|
855
|
|
856 Yesm6<-sapply(pPositionm6, function(x) Yes_props[x])
|
|
857 Yesm5<-sapply(pPositionm5, function(x) Yes_props[x])
|
|
858 Yesm4<-sapply(pPositionm4, function(x) Yes_props[x])
|
|
859 Yesm3<-sapply(pPositionm3, function(x) Yes_props[x])
|
|
860 Yesm2<-sapply(pPositionm2, function(x) Yes_props[x])
|
|
861 Yesm1<-sapply(pPositionm1, function(x) Yes_props[x])
|
|
862 Yesd0<-sapply(pPositiond0, function(x) Yes_props[x])
|
|
863 Yesp1<-sapply(pPositionp1, function(x) Yes_props[x])
|
|
864 Yesp2<-sapply(pPositionp2, function(x) Yes_props[x])
|
|
865 Yesp3<-sapply(pPositionp3, function(x) Yes_props[x])
|
|
866 Yesp4<-sapply(pPositionp4, function(x) Yes_props[x])
|
|
867 Yesp5<-sapply(pPositionp5, function(x) Yes_props[x])
|
|
868 Yesp6<-sapply(pPositionp6, function(x) Yes_props[x])
|
|
869
|
|
870 FLT3m6<-sapply(pPositionm6, function(x) FLT3_props[x])
|
|
871 FLT3m5<-sapply(pPositionm5, function(x) FLT3_props[x])
|
|
872 FLT3m4<-sapply(pPositionm4, function(x) FLT3_props[x])
|
|
873 FLT3m3<-sapply(pPositionm3, function(x) FLT3_props[x])
|
|
874 FLT3m2<-sapply(pPositionm2, function(x) FLT3_props[x])
|
|
875 FLT3m1<-sapply(pPositionm1, function(x) FLT3_props[x])
|
|
876 FLT3d0<-sapply(pPositiond0, function(x) FLT3_props[x])
|
|
877 FLT3p1<-sapply(pPositionp1, function(x) FLT3_props[x])
|
|
878 FLT3p2<-sapply(pPositionp2, function(x) FLT3_props[x])
|
|
879 FLT3p3<-sapply(pPositionp3, function(x) FLT3_props[x])
|
|
880 FLT3p4<-sapply(pPositionp4, function(x) FLT3_props[x])
|
|
881 FLT3p5<-sapply(pPositionp5, function(x) FLT3_props[x])
|
|
882 FLT3p6<-sapply(pPositionp6, function(x) FLT3_props[x])
|
|
883
|
|
884 ALKm6<-sapply(pPositionm6, function(x) ALK_props[x])
|
|
885 ALKm5<-sapply(pPositionm5, function(x) ALK_props[x])
|
|
886 ALKm4<-sapply(pPositionm4, function(x) ALK_props[x])
|
|
887 ALKm3<-sapply(pPositionm3, function(x) ALK_props[x])
|
|
888 ALKm2<-sapply(pPositionm2, function(x) ALK_props[x])
|
|
889 ALKm1<-sapply(pPositionm1, function(x) ALK_props[x])
|
|
890 ALKd0<-sapply(pPositiond0, function(x) ALK_props[x])
|
|
891 ALKp1<-sapply(pPositionp1, function(x) ALK_props[x])
|
|
892 ALKp2<-sapply(pPositionp2, function(x) ALK_props[x])
|
|
893 ALKp3<-sapply(pPositionp3, function(x) ALK_props[x])
|
|
894 ALKp4<-sapply(pPositionp4, function(x) ALK_props[x])
|
|
895 ALKp5<-sapply(pPositionp5, function(x) ALK_props[x])
|
|
896 ALKp6<-sapply(pPositionp6, function(x) ALK_props[x])
|
|
897 }
|
|
898 #############################
|
|
899 #here's my plan: EITHER
|
|
900 #make each peptide with different values for each kinase's EPM
|
|
901 #OR
|
|
902 #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
|
|
903 #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
|
|
904 #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
|
|
905 #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
|
|
906 #then I make the combinatorial library
|
|
907 #then I use an apply function to multiply all values together, which is how fisher does the scoring
|
|
908 ########################################
|
|
909
|
|
910
|
|
911 total=length(Positionp6)*length(Positionp5)*length(Positionp4)*length(Positionp3)*(length(Positionp2))*length(Positionp1)*
|
|
912 length(Positiond0)*length(Positionm1)*length(Positionm2)*length(Positionm3)*length(Positionm4)*length(Positionm5)*length(Positionm6)
|
|
913 #this is just a way to doublecheck that the length of the generated peptides vector is correct
|
|
914
|
|
915 GeneratedPeptides<-rep(NA, times=total*13)
|
|
916 GeneratedPeptides<-matrix(data = GeneratedPeptides,ncol = 13)
|
|
917
|
|
918
|
|
919 ThisKinPeptides<-GeneratedPeptides
|
|
920
|
|
921 AblPeptides<-GeneratedPeptides
|
|
922 ArgPeptides<-GeneratedPeptides
|
|
923 BtkPeptides<-GeneratedPeptides
|
|
924 CskPeptides<-GeneratedPeptides
|
|
925 FynPeptides<-GeneratedPeptides
|
|
926 HckPeptides<-GeneratedPeptides
|
|
927 JAK2Peptides<-GeneratedPeptides
|
|
928 LckPeptides<-GeneratedPeptides
|
|
929 LynPeptides<-GeneratedPeptides
|
|
930 Pyk2Peptides<-GeneratedPeptides
|
|
931 SrcPeptides<-GeneratedPeptides
|
|
932 SykPeptides<-GeneratedPeptides
|
|
933 YesPeptides<-GeneratedPeptides
|
|
934 FLT3Peptides<-GeneratedPeptides
|
|
935 ALKPeptides<-GeneratedPeptides
|
|
936
|
|
937 #create an empty vector of correct length by finding the number of each AAs per position and multiplying them
|
|
938 count<-0
|
|
939
|
|
940 # for (t in 1:length(Positionm7)) {
|
|
941 for (s in 1:length(Positionm6)) {
|
|
942 for (r in 1:length(Positionm5)) {
|
|
943 for (i in 1:length(Positionm4)) {
|
|
944 for (j in 1:length(Positionm3)) {
|
|
945 for (k in 1:length(Positionm2)) {
|
|
946 for (l in 1:length(Positionm1)) {
|
|
947 for (m in 1:length(Positiond0)) {
|
|
948 for (n in 1:length(Positionp1)) {
|
|
949 for (o in 1:length(Positionp2)) {
|
|
950 for (p in 1:length(Positionp3)) {
|
|
951 for (q in 1:length(Positionp4)) {
|
|
952 for (u in 1:length(Positionp5)) {
|
|
953 for (v in 1:length(Positionp6)) {
|
|
954 # for (w in 1:length(Positionp7)) {
|
|
955 # i=1
|
|
956 # j=1
|
|
957 # k=1
|
|
958 # l=1
|
|
959 # m=1
|
|
960 # n=1
|
|
961 # o=1
|
|
962 # p=1
|
|
963 # q=1
|
|
964 #
|
|
965 #for every single position, increment the count number, create a peptide using the AAs at that position
|
|
966 #then put them together into the generated peptides sequencex
|
|
967 count<-count+1
|
|
968 tabulation<-c(Positionm6[s],Positionm5[r],Positionm4[i],Positionm3[j],Positionm2[k],Positionm1[l],Positiond0[m],Positionp1[n],
|
|
969 Positionp2[o],Positionp3[p],Positionp4[q],Positionp5[u],Positionp6[v])
|
|
970 # 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])
|
|
971 #tabulation<-paste(tabulation, sep="", collapse="")
|
|
972 GeneratedPeptides[count,1:13]<-tabulation
|
|
973 # NumeratedPeptides[count,1:13]<-numeration
|
|
974
|
|
975 ThisKinTableNumeration<-c(ThisKinTablem6[s],ThisKinTablem5[r],ThisKinTablem4[i],ThisKinTablem3[j],ThisKinTablem2[k],ThisKinTablem1[l],ThisKinTabled0[m],ThisKinTablep1[n],
|
|
976 ThisKinTablep2[o],ThisKinTablep3[p],ThisKinTablep4[q],ThisKinTablep5[u],ThisKinTablep6[v])
|
|
977
|
|
978 AblNumeration<-c(Ablm6[s],Ablm5[r],Ablm4[i],Ablm3[j],Ablm2[k],Ablm1[l],Abld0[m],Ablp1[n],
|
|
979 Ablp2[o],Ablp3[p],Ablp4[q],Ablp5[u],Ablp6[v])
|
|
980 ArgNumeration<-c(Argm6[s],Argm5[r],Argm4[i],Argm3[j],Argm2[k],Argm1[l],Argd0[m],Argp1[n],
|
|
981 Argp2[o],Argp3[p],Argp4[q],Argp5[u],Argp6[v])
|
|
982 BtkNumeration<-c(Btkm6[s],Btkm5[r],Btkm4[i],Btkm3[j],Btkm2[k],Btkm1[l],Btkd0[m],Btkp1[n],
|
|
983 Btkp2[o],Btkp3[p],Btkp4[q],Btkp5[u],Btkp6[v])
|
|
984 CskNumeration<-c(Cskm6[s],Cskm5[r],Cskm4[i],Cskm3[j],Cskm2[k],Cskm1[l],Cskd0[m],Cskp1[n],
|
|
985 Cskp2[o],Cskp3[p],Cskp4[q],Cskp5[u],Cskp6[v])
|
|
986 FynNumeration<-c(Fynm6[s],Fynm5[r],Fynm4[i],Fynm3[j],Fynm2[k],Fynm1[l],Fynd0[m],Fynp1[n],
|
|
987 Fynp2[o],Fynp3[p],Fynp4[q],Fynp5[u],Fynp6[v])
|
|
988 HckNumeration<-c(Hckm6[s],Hckm5[r],Hckm4[i],Hckm3[j],Hckm2[k],Hckm1[l],Hckd0[m],Hckp1[n],
|
|
989 Hckp2[o],Hckp3[p],Hckp4[q],Hckp5[u],Hckp6[v])
|
|
990 JAK2Numeration<-c(JAK2m6[s],JAK2m5[r],JAK2m4[i],JAK2m3[j],JAK2m2[k],JAK2m1[l],JAK2d0[m],JAK2p1[n],
|
|
991 JAK2p2[o],JAK2p3[p],JAK2p4[q],JAK2p5[u],JAK2p6[v])
|
|
992 LckNumeration<-c(Lckm6[s],Lckm5[r],Lckm4[i],Lckm3[j],Lckm2[k],Lckm1[l],Lckd0[m],Lckp1[n],
|
|
993 Lckp2[o],Lckp3[p],Lckp4[q],Lckp5[u],Lckp6[v])
|
|
994 LynNumeration<-c(Lynm6[s],Lynm5[r],Lynm4[i],Lynm3[j],Lynm2[k],Lynm1[l],Lynd0[m],Lynp1[n],
|
|
995 Lynp2[o],Lynp3[p],Lynp4[q],Lynp5[u],Lynp6[v])
|
|
996 Pyk2Numeration<-c(Pyk2m6[s],Pyk2m5[r],Pyk2m4[i],Pyk2m3[j],Pyk2m2[k],Pyk2m1[l],Pyk2d0[m],Pyk2p1[n],
|
|
997 Pyk2p2[o],Pyk2p3[p],Pyk2p4[q],Pyk2p5[u],Pyk2p6[v])
|
|
998 SrcNumeration<-c(Srcm6[s],Srcm5[r],Srcm4[i],Srcm3[j],Srcm2[k],Srcm1[l],Srcd0[m],Srcp1[n],
|
|
999 Srcp2[o],Srcp3[p],Srcp4[q],Srcp5[u],Srcp6[v])
|
|
1000 SykNumeration<-c(Sykm6[s],Sykm5[r],Sykm4[i],Sykm3[j],Sykm2[k],Sykm1[l],Sykd0[m],Sykp1[n],
|
|
1001 Sykp2[o],Sykp3[p],Sykp4[q],Sykp5[u],Sykp6[v])
|
|
1002 YesNumeration<-c(Yesm6[s],Yesm5[r],Yesm4[i],Yesm3[j],Yesm2[k],Yesm1[l],Yesd0[m],Yesp1[n],
|
|
1003 Yesp2[o],Yesp3[p],Yesp4[q],Yesp5[u],Yesp6[v])
|
|
1004 FLT3Numeration<-c(FLT3m6[s],FLT3m5[r],FLT3m4[i],FLT3m3[j],FLT3m2[k],FLT3m1[l],FLT3d0[m],FLT3p1[n],
|
|
1005 FLT3p2[o],FLT3p3[p],FLT3p4[q],FLT3p5[u],FLT3p6[v])
|
|
1006 ALKNumeration<-c(ALKm6[s],ALKm5[r],ALKm4[i],ALKm3[j],ALKm2[k],ALKm1[l],ALKd0[m],ALKp1[n],
|
|
1007 ALKp2[o],ALKp3[p],ALKp4[q],ALKp5[u],ALKp6[v])
|
|
1008
|
|
1009 ThisKinPeptides[count,1:13]<-ThisKinTableNumeration
|
|
1010
|
|
1011 AblPeptides[count,1:13]<-AblNumeration
|
|
1012 ArgPeptides[count,1:13]<-ArgNumeration
|
|
1013 BtkPeptides[count,1:13]<-BtkNumeration
|
|
1014 CskPeptides[count,1:13]<-CskNumeration
|
|
1015 FynPeptides[count,1:13]<-FynNumeration
|
|
1016 HckPeptides[count,1:13]<-HckNumeration
|
|
1017 JAK2Peptides[count,1:13]<-JAK2Numeration
|
|
1018 LckPeptides[count,1:13]<-LckNumeration
|
|
1019 LynPeptides[count,1:13]<-LynNumeration
|
|
1020 Pyk2Peptides[count,1:13]<-Pyk2Numeration
|
|
1021 SrcPeptides[count,1:13]<-SrcNumeration
|
|
1022 SykPeptides[count,1:13]<-SykNumeration
|
|
1023 YesPeptides[count,1:13]<-YesNumeration
|
|
1024 FLT3Peptides[count,1:13]<-FLT3Numeration
|
|
1025 ALKPeptides[count,1:13]<-ALKNumeration
|
|
1026
|
|
1027
|
|
1028
|
|
1029 # }
|
|
1030 }
|
|
1031 }
|
|
1032 }
|
|
1033 }
|
|
1034 }
|
|
1035 }
|
|
1036 }
|
|
1037 }
|
|
1038 }
|
|
1039 }
|
|
1040 }
|
|
1041 }
|
|
1042 }
|
|
1043 # }
|
|
1044
|
|
1045 ThisKinPeptides<-as.data.frame(ThisKinPeptides)
|
|
1046 ThisKinPeptides$V1<-as.numeric(ThisKinPeptides$V1)
|
|
1047 ThisKinPeptides$V2<-as.numeric(ThisKinPeptides$V2)
|
|
1048 ThisKinPeptides$V3<-as.numeric(ThisKinPeptides$V3)
|
|
1049 ThisKinPeptides$V4<-as.numeric(ThisKinPeptides$V4)
|
|
1050 ThisKinPeptides$V5<-as.numeric(ThisKinPeptides$V5)
|
|
1051 ThisKinPeptides$V6<-as.numeric(ThisKinPeptides$V6)
|
|
1052 ThisKinPeptides$V7<-as.numeric(ThisKinPeptides$V7)
|
|
1053 ThisKinPeptides$V8<-as.numeric(ThisKinPeptides$V8)
|
|
1054 ThisKinPeptides$V9<-as.numeric(ThisKinPeptides$V9)
|
|
1055 ThisKinPeptides$V10<-as.numeric(ThisKinPeptides$V10)
|
|
1056 ThisKinPeptides$V11<-as.numeric(ThisKinPeptides$V11)
|
|
1057 ThisKinPeptides$V12<-as.numeric(ThisKinPeptides$V12)
|
|
1058 ThisKinPeptides$V13<-as.numeric(ThisKinPeptides$V13)
|
|
1059 ThisKinGeneratedScores<-with(ThisKinPeptides, V1*V2*V3*V4*V5*V6*V7*V8*V9*V10*V11*V12*V13)
|
|
1060
|
|
1061 AblPeptides<-as.data.frame(AblPeptides)
|
|
1062 AblGeneratedScores<-with(AblPeptides, V1*V2*V3*V4*V5*V6*V7*V8*V9*V10*V11*V12*V13)
|
|
1063 ArgPeptides<-as.data.frame(ArgPeptides)
|
|
1064 ArgGeneratedScores<-with(ArgPeptides, V1*V2*V3*V4*V5*V6*V7*V8*V9*V10*V11*V12*V13)
|
|
1065 BtkPeptides<-as.data.frame(BtkPeptides)
|
|
1066 BtkGeneratedScores<-with(BtkPeptides, V1*V2*V3*V4*V5*V6*V7*V8*V9*V10*V11*V12*V13)
|
|
1067 CskPeptides<-as.data.frame(CskPeptides)
|
|
1068 CskGeneratedScores<-with(CskPeptides, V1*V2*V3*V4*V5*V6*V7*V8*V9*V10*V11*V12*V13)
|
|
1069 FynPeptides<-as.data.frame(FynPeptides)
|
|
1070 FynGeneratedScores<-with(FynPeptides, V1*V2*V3*V4*V5*V6*V7*V8*V9*V10*V11*V12*V13)
|
|
1071 HckPeptides<-as.data.frame(HckPeptides)
|
|
1072 HckGeneratedScores<-with(HckPeptides, V1*V2*V3*V4*V5*V6*V7*V8*V9*V10*V11*V12*V13)
|
|
1073 JAK2Peptides<-as.data.frame(JAK2Peptides)
|
|
1074 JAK2GeneratedScores<-with(JAK2Peptides, V1*V2*V3*V4*V5*V6*V7*V8*V9*V10*V11*V12*V13)
|
|
1075 LckPeptides<-as.data.frame(LckPeptides)
|
|
1076 LckGeneratedScores<-with(LckPeptides, V1*V2*V3*V4*V5*V6*V7*V8*V9*V10*V11*V12*V13)
|
|
1077 LynPeptides<-as.data.frame(LynPeptides)
|
|
1078 LynGeneratedScores<-with(LynPeptides, V1*V2*V3*V4*V5*V6*V7*V8*V9*V10*V11*V12*V13)
|
|
1079 Pyk2Peptides<-as.data.frame(Pyk2Peptides)
|
|
1080 Pyk2GeneratedScores<-with(Pyk2Peptides, V1*V2*V3*V4*V5*V6*V7*V8*V9*V10*V11*V12*V13)
|
|
1081 SrcPeptides<-as.data.frame(SrcPeptides)
|
|
1082 SrcGeneratedScores<-with(SrcPeptides, V1*V2*V3*V4*V5*V6*V7*V8*V9*V10*V11*V12*V13)
|
|
1083 SykPeptides<-as.data.frame(SykPeptides)
|
|
1084 SykGeneratedScores<-with(SykPeptides, V1*V2*V3*V4*V5*V6*V7*V8*V9*V10*V11*V12*V13)
|
|
1085 YesPeptides<-as.data.frame(YesPeptides)
|
|
1086 YesGeneratedScores<-with(YesPeptides, V1*V2*V3*V4*V5*V6*V7*V8*V9*V10*V11*V12*V13)
|
|
1087 FLT3Peptides<-as.data.frame(FLT3Peptides)
|
|
1088 FLT3GeneratedScores<-with(FLT3Peptides, V1*V2*V3*V4*V5*V6*V7*V8*V9*V10*V11*V12*V13)
|
|
1089 ALKPeptides<-as.data.frame(ALKPeptides)
|
|
1090 ALKGeneratedScores<-with(ALKPeptides, V1*V2*V3*V4*V5*V6*V7*V8*V9*V10*V11*V12*V13)
|
|
1091
|
|
1092 AblGeneratedScores<-as.matrix(AblGeneratedScores, ncol=1)
|
|
1093 AblNorm<-AblNorm<-1/as.numeric(Abl[22,1])
|
|
1094 AblNormColumn<-as.matrix(rep(AblNorm,times=nrow(AblGeneratedScores)),ncol=1)
|
|
1095 AblGeneratedScores<-cbind(AblGeneratedScores,AblNormColumn)
|
|
1096 AblGeneratedScores<-as.data.frame(AblGeneratedScores)
|
|
1097 AblGeneratedScores<-with(AblGeneratedScores, V1/(V1+V2))
|
|
1098
|
|
1099 ArgGeneratedScores<-as.matrix(ArgGeneratedScores, ncol=1)
|
|
1100 ArgNorm<-ArgNorm<-1/as.numeric(Arg[22,1])
|
|
1101 ArgNormColumn<-as.matrix(rep(ArgNorm,times=nrow(ArgGeneratedScores)),ncol=1)
|
|
1102 ArgGeneratedScores<-cbind(ArgGeneratedScores,ArgNormColumn)
|
|
1103 ArgGeneratedScores<-as.data.frame(ArgGeneratedScores)
|
|
1104 ArgGeneratedScores<-with(ArgGeneratedScores, V1/(V1+V2))
|
|
1105
|
|
1106 BtkGeneratedScores<-as.matrix(BtkGeneratedScores, ncol=1)
|
|
1107 BtkNorm<-BtkNorm<-1/as.numeric(Btk[22,1])
|
|
1108 BtkNormColumn<-as.matrix(rep(BtkNorm,times=nrow(BtkGeneratedScores)),ncol=1)
|
|
1109 BtkGeneratedScores<-cbind(BtkGeneratedScores,BtkNormColumn)
|
|
1110 BtkGeneratedScores<-as.data.frame(BtkGeneratedScores)
|
|
1111 BtkGeneratedScores<-with(BtkGeneratedScores, V1/(V1+V2))
|
|
1112
|
|
1113 CskGeneratedScores<-as.matrix(CskGeneratedScores, ncol=1)
|
|
1114 CskNorm<-CskNorm<-1/as.numeric(Csk[22,1])
|
|
1115 CskNormColumn<-as.matrix(rep(CskNorm,times=nrow(CskGeneratedScores)),ncol=1)
|
|
1116 CskGeneratedScores<-cbind(CskGeneratedScores,CskNormColumn)
|
|
1117 CskGeneratedScores<-as.data.frame(CskGeneratedScores)
|
|
1118 CskGeneratedScores<-with(CskGeneratedScores, V1/(V1+V2))
|
|
1119
|
|
1120 FynGeneratedScores<-as.matrix(FynGeneratedScores, ncol=1)
|
|
1121 FynNorm<-FynNorm<-1/as.numeric(Fyn[22,1])
|
|
1122 FynNormColumn<-as.matrix(rep(FynNorm,times=nrow(FynGeneratedScores)),ncol=1)
|
|
1123 FynGeneratedScores<-cbind(FynGeneratedScores,FynNormColumn)
|
|
1124 FynGeneratedScores<-as.data.frame(FynGeneratedScores)
|
|
1125 FynGeneratedScores<-with(FynGeneratedScores, V1/(V1+V2))
|
|
1126
|
|
1127 HckGeneratedScores<-as.matrix(HckGeneratedScores, ncol=1)
|
|
1128 HckNorm<-HckNorm<-1/as.numeric(Hck[22,1])
|
|
1129 HckNormColumn<-as.matrix(rep(HckNorm,times=nrow(HckGeneratedScores)),ncol=1)
|
|
1130 HckGeneratedScores<-cbind(HckGeneratedScores,HckNormColumn)
|
|
1131 HckGeneratedScores<-as.data.frame(HckGeneratedScores)
|
|
1132 HckGeneratedScores<-with(HckGeneratedScores, 100*V1/(V1+V2))
|
|
1133
|
|
1134 JAK2GeneratedScores<-as.matrix(JAK2GeneratedScores, ncol=1)
|
|
1135 JAK2Norm<-JAK2Norm<-1/as.numeric(JAK2[22,1])
|
|
1136 JAK2NormColumn<-as.matrix(rep(JAK2Norm,times=nrow(JAK2GeneratedScores)),ncol=1)
|
|
1137 JAK2GeneratedScores<-cbind(JAK2GeneratedScores,JAK2NormColumn)
|
|
1138 JAK2GeneratedScores<-as.data.frame(JAK2GeneratedScores)
|
|
1139 JAK2GeneratedScores<-with(JAK2GeneratedScores, 100*V1/(V1+V2))
|
|
1140
|
|
1141 LckGeneratedScores<-as.matrix(LckGeneratedScores, ncol=1)
|
|
1142 LckNorm<-LckNorm<-1/as.numeric(Lck[22,1])
|
|
1143 LckNormColumn<-as.matrix(rep(LckNorm,times=nrow(LckGeneratedScores)),ncol=1)
|
|
1144 LckGeneratedScores<-cbind(LckGeneratedScores,LckNormColumn)
|
|
1145 LckGeneratedScores<-as.data.frame(LckGeneratedScores)
|
|
1146 LckGeneratedScores<-with(LckGeneratedScores, 100*V1/(V1+V2))
|
|
1147
|
|
1148 LynGeneratedScores<-as.matrix(LynGeneratedScores, ncol=1)
|
|
1149 LynNorm<-LynNorm<-1/as.numeric(Lyn[22,1])
|
|
1150 LynNormColumn<-as.matrix(rep(LynNorm,times=nrow(LynGeneratedScores)),ncol=1)
|
|
1151 LynGeneratedScores<-cbind(LynGeneratedScores,LynNormColumn)
|
|
1152 LynGeneratedScores<-as.data.frame(LynGeneratedScores)
|
|
1153 LynGeneratedScores<-with(LynGeneratedScores, 100*V1/(V1+V2))
|
|
1154
|
|
1155 Pyk2GeneratedScores<-as.matrix(Pyk2GeneratedScores, ncol=1)
|
|
1156 Pyk2Norm<-Pyk2Norm<-1/as.numeric(Pyk2[22,1])
|
|
1157 Pyk2NormColumn<-as.matrix(rep(Pyk2Norm,times=nrow(Pyk2GeneratedScores)),ncol=1)
|
|
1158 Pyk2GeneratedScores<-cbind(Pyk2GeneratedScores,Pyk2NormColumn)
|
|
1159 Pyk2GeneratedScores<-as.data.frame(Pyk2GeneratedScores)
|
|
1160 Pyk2GeneratedScores<-with(Pyk2GeneratedScores, 100*V1/(V1+V2))
|
|
1161
|
|
1162 SrcGeneratedScores<-as.matrix(SrcGeneratedScores, ncol=1)
|
|
1163 SrcNorm<-SrcNorm<-1/as.numeric(Src[22,1])
|
|
1164 SrcNormColumn<-as.matrix(rep(SrcNorm,times=nrow(SrcGeneratedScores)),ncol=1)
|
|
1165 SrcGeneratedScores<-cbind(SrcGeneratedScores,SrcNormColumn)
|
|
1166 SrcGeneratedScores<-as.data.frame(SrcGeneratedScores)
|
|
1167 SrcGeneratedScores<-with(SrcGeneratedScores, 100*V1/(V1+V2))
|
|
1168
|
|
1169 SykGeneratedScores<-as.matrix(SykGeneratedScores, ncol=1)
|
|
1170 SykNorm<-SykNorm<-1/as.numeric(Syk[22,1])
|
|
1171 SykNormColumn<-as.matrix(rep(SykNorm,times=nrow(SykGeneratedScores)),ncol=1)
|
|
1172 SykGeneratedScores<-cbind(SykGeneratedScores,SykNormColumn)
|
|
1173 SykGeneratedScores<-as.data.frame(SykGeneratedScores)
|
|
1174 SykGeneratedScores<-with(SykGeneratedScores, 100*V1/(V1+V2))
|
|
1175
|
|
1176 YesGeneratedScores<-as.matrix(YesGeneratedScores, ncol=1)
|
|
1177 YesNorm<-YesNorm<-1/as.numeric(Yes[22,1])
|
|
1178 YesNormColumn<-as.matrix(rep(YesNorm,times=nrow(YesGeneratedScores)),ncol=1)
|
|
1179 YesGeneratedScores<-cbind(YesGeneratedScores,YesNormColumn)
|
|
1180 YesGeneratedScores<-as.data.frame(YesGeneratedScores)
|
|
1181 YesGeneratedScores<-with(YesGeneratedScores, 100*V1/(V1+V2))
|
|
1182
|
|
1183 FLT3GeneratedScores<-as.matrix(FLT3GeneratedScores, ncol=1)
|
|
1184 FLT3Norm<-FLT3Norm<-1/as.numeric(FLT3[22,1])
|
|
1185 FLT3NormColumn<-as.matrix(rep(FLT3Norm,times=nrow(FLT3GeneratedScores)),ncol=1)
|
|
1186 FLT3GeneratedScores<-cbind(FLT3GeneratedScores,FLT3NormColumn)
|
|
1187 FLT3GeneratedScores<-as.data.frame(FLT3GeneratedScores)
|
|
1188 FLT3GeneratedScores<-with(FLT3GeneratedScores, 100*V1/(V1+V2))
|
|
1189
|
|
1190 ALKGeneratedScores<-as.matrix(ALKGeneratedScores, ncol=1)
|
|
1191 ALKNorm<-ALKNorm<-1/as.numeric(ALK[22,1])
|
|
1192 ALKNormColumn<-as.matrix(rep(ALKNorm,times=nrow(ALKGeneratedScores)),ncol=1)
|
|
1193 ALKGeneratedScores<-cbind(ALKGeneratedScores,ALKNormColumn)
|
|
1194 ALKGeneratedScores<-as.data.frame(ALKGeneratedScores)
|
|
1195 ALKGeneratedScores<-with(ALKGeneratedScores, 100*V1/(V1+V2))
|
|
1196
|
|
1197 ThisKinGeneratedScores1<-as.matrix(ThisKinGeneratedScores, ncol=1)
|
|
1198 ThisKinNorm<-ThisKinNorm<-1/as.numeric(NormalizationScore)
|
|
1199 ThisKinNormColumn<-as.matrix(rep(ThisKinNorm,times=nrow(ThisKinGeneratedScores1)),ncol=1)
|
|
1200 ThisKinGeneratedScores1<-cbind(ThisKinGeneratedScores1,ThisKinNormColumn)
|
|
1201 ThisKinGeneratedScores1<-as.data.frame(ThisKinGeneratedScores1)
|
|
1202 ThisKinGeneratedScores1<-with(ThisKinGeneratedScores1, V1/(V1+V2))
|
|
1203
|
|
1204 AblThresh<-as.numeric(Abl[24,1])
|
|
1205 AblActive<-unlist(AblGeneratedScores)>AblThresh/100
|
|
1206 ArgThresh<-as.numeric(Arg[24,1])
|
|
1207 ArgActive<-unlist(ArgGeneratedScores)>ArgThresh/100
|
|
1208 BtkThresh<-as.numeric(Btk[24,1])
|
|
1209 BtkActive<-unlist(BtkGeneratedScores)>BtkThresh/100
|
|
1210 CskThresh<-as.numeric(Csk[24,1])
|
|
1211 CskActive<-(CskGeneratedScores)>CskThresh/100
|
|
1212 FynThresh<-as.numeric(Fyn[24,1])
|
|
1213 FynActive<-unlist(FynGeneratedScores)>FynThresh/100
|
|
1214 HckThresh<-as.numeric(Hck[24,1])
|
|
1215 HckActive<-unlist(HckGeneratedScores)>HckThresh/100
|
|
1216 JAK2Thresh<-as.numeric(JAK2[24,1])
|
|
1217 JAk2Active<-unlist(JAK2GeneratedScores)>JAK2Thresh/100
|
|
1218 LckThresh<-as.numeric(Lck[24,1])
|
|
1219 LckActive<-unlist(LckGeneratedScores)>LckThresh/100
|
|
1220 LynThresh<-as.numeric(Lyn[24,1])
|
|
1221 LynActive<-unlist(LynGeneratedScores)>LynThresh/100
|
|
1222 Pyk2Thresh<-as.numeric(Pyk2[24,1])
|
|
1223 Pyk2Active<-unlist(Pyk2GeneratedScores)>Pyk2Thresh/100
|
|
1224 SrcThresh<-as.numeric(Src[24,1])
|
|
1225 SrcActive<-unlist(SrcGeneratedScores)>SrcThresh/100
|
|
1226 SykThresh<-as.numeric(Syk[24,1])
|
|
1227 SykActive<-unlist(SykGeneratedScores)>SykThresh/100
|
|
1228 YesThresh<-as.numeric(Yes[24,1])
|
|
1229 YesActive<-unlist(YesGeneratedScores)>YesThresh/100
|
|
1230 ALKThresh<-as.numeric(ALK[24,1])
|
|
1231 ALKActive<-unlist(ALKGeneratedScores)>ALKThresh/100
|
|
1232 FLT3Thresh<-as.numeric(FLT3[24,1])
|
|
1233 FLT3Active<-unlist(FLT3GeneratedScores)>FLT3Thresh/100
|
|
1234
|
|
1235
|
|
1236 AllActive<-AblActive+ArgActive+BtkActive+CskActive+FynActive+HckActive+JAk2Active+LckActive+LynActive+Pyk2Active+SrcActive+SykActive+YesActive+ALKActive+FLT3Active
|
|
1237 Scores<-ThisKinGeneratedScores
|
|
1238 ThresholdValues<-ThisKinGeneratedScores1
|
|
1239
|
|
1240 FullMotifs<-rep("Z",times=nrow(GeneratedPeptides))
|
|
1241 for (i in 1:nrow(GeneratedPeptides)) {
|
|
1242 motif<-GeneratedPeptides[i,1:13]
|
|
1243 motif<-paste(motif,sep = "", collapse = "")
|
|
1244 FullMotifs[i]<-motif
|
|
1245 }
|
|
1246
|
|
1247 Scores<-1/Scores
|
|
1248
|
|
1249 PeptidesWithRanks<-cbind.data.frame(FullMotifs,GeneratedPeptides,Scores,ThresholdValues)
|
|
1250
|
|
1251 PeptidesWithRanks<-cbind.data.frame(PeptidesWithRanks,AllActive,AblActive,ArgActive,BtkActive,CskActive,FynActive,HckActive,JAk2Active,LckActive,LynActive,Pyk2Active,SrcActive,SykActive,YesActive)
|
|
1252 RanksPeptides<-PeptidesWithRanks[order(PeptidesWithRanks$AllActive,
|
|
1253 PeptidesWithRanks$Scores,
|
|
1254 decreasing = FALSE),]
|
|
1255
|
|
1256 RanksPeptides$Scores<-1/RanksPeptides$Scores
|
|
1257
|
|
1258 #RanksPeptides$Scores<-1/RanksPeptides$Scores
|
|
1259
|
|
1260 #I want to see how much of the sequence space we still have available. I should do so
|
|
1261
|
|
1262 write.table(RanksPeptides,file = "output.csv",append = FALSE,row.names = FALSE,col.names = TRUE,sep = ",")
|
|
1263
|
6
|
1264 # })
|
|
1265
|
|
1266 options(warn=0) |