annotate kinatestid_r/.Rhistory @ 18:e16ca3e9fc49 draft

Uploaded
author jfb
date Tue, 17 Apr 2018 13:07:41 -0400
parents
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
18
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
1 SubstrateBackgroundFrequency<- read.csv("7to7PositiveSubBAckFreq_11-7-17.csv", stringsAsFactors=FALSE, header = FALSE)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
2 ScreenerFilename<-"screener7-7.csv"
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
3 FILENAME<-"output1.csv"
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
4 FILENAME2<-"output2.csv"
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
5 FILENAME3<-"output3.csv"
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
6 SubstrateBackgroundFrequency<-t(SubstrateBackgroundFrequency)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
7 # number<-nrow(SubstrateBackgroundFrequency)-1
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
8 SubstrateBackgroundFrequency<-SubstrateBackgroundFrequency[2:nrow(SubstrateBackgroundFrequency),]
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
9 Sub<-na.omit(SubstrateBackgroundFrequency)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
10 SubstrateBackgroundFrequency<-Sub
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
11 #
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
12 # args = commandArgs(trailingOnly=TRUE)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
13 # TodaysKinase<-as.character(args[1])
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
14 OutputMatrix<-"KinaseMatrix.csv"
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
15 CharacterizationTable<-"CharacterizationTableForThisKinase.csv"
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
16 SDtable<-"SDtableforthisKinase"
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
17 SiteSelectivityTable<-"SiteSelectivityForThisKinase"
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
18 substrates<-matrix(rep("A",times=((nrow(ImportedSubstrateList)-1)*15)),ncol = 15)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
19 #SeqsToBeScored<-"asdasd"
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
20 for (i in 2:nrow(ImportedSubstrateList))
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
21 {
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
22 substratemotif<-ImportedSubstrateList[i,4:18]
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
23 substratemotif[8]<-"Y"
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
24 #substratemotif<-paste(substratemotif,sep = "",collapse = "")
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
25 j=i-1
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
26 substratemotif<-unlist(substratemotif)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
27 substrates[j,1:15]<-substratemotif
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
28 }
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
29 # SpacesToOs<-c(""="O",)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
30 # substrates<-SpacesToOs[substrates]
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
31 SubstrateBackgroundFrequency[nrow(SubstrateBackgroundFrequency),2]
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
32 if(2==2){
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
33 Amean<-mean(as.numeric(SubstrateBackgroundFrequency[1:(nrow(SubstrateBackgroundFrequency)),2]), na.rm=TRUE)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
34 Cmean<-mean(as.numeric(SubstrateBackgroundFrequency[1:(nrow(SubstrateBackgroundFrequency)),3]), na.rm=TRUE)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
35 Dmean<-mean(as.numeric(SubstrateBackgroundFrequency[1:(nrow(SubstrateBackgroundFrequency)),4]), na.rm=TRUE)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
36 Emean<-mean(as.numeric(SubstrateBackgroundFrequency[1:(nrow(SubstrateBackgroundFrequency)),5]), na.rm=TRUE)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
37 Fmean<-mean(as.numeric(SubstrateBackgroundFrequency[1:(nrow(SubstrateBackgroundFrequency)),6]), na.rm=TRUE)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
38 Gmean<-mean(as.numeric(SubstrateBackgroundFrequency[1:(nrow(SubstrateBackgroundFrequency)),7]), na.rm=TRUE)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
39 Hmean<-mean(as.numeric(SubstrateBackgroundFrequency[1:(nrow(SubstrateBackgroundFrequency)),8]), na.rm=TRUE)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
40 Imean<-mean(as.numeric(SubstrateBackgroundFrequency[1:(nrow(SubstrateBackgroundFrequency)),9]), na.rm=TRUE)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
41 Kmean<-mean(as.numeric(SubstrateBackgroundFrequency[1:(nrow(SubstrateBackgroundFrequency)),10]), na.rm=TRUE)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
42 Lmean<-mean(as.numeric(SubstrateBackgroundFrequency[1:(nrow(SubstrateBackgroundFrequency)),11]), na.rm=TRUE)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
43 Mmean<-mean(as.numeric(SubstrateBackgroundFrequency[1:(nrow(SubstrateBackgroundFrequency)),12]), na.rm=TRUE)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
44 Nmean<-mean(as.numeric(SubstrateBackgroundFrequency[1:(nrow(SubstrateBackgroundFrequency)),13]), na.rm=TRUE)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
45 Pmean<-mean(as.numeric(SubstrateBackgroundFrequency[1:(nrow(SubstrateBackgroundFrequency)),14]), na.rm=TRUE)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
46 Qmean<-mean(as.numeric(SubstrateBackgroundFrequency[1:(nrow(SubstrateBackgroundFrequency)),15]), na.rm=TRUE)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
47 Rmean<-mean(as.numeric(SubstrateBackgroundFrequency[1:(nrow(SubstrateBackgroundFrequency)),16]), na.rm=TRUE)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
48 Smean<-mean(as.numeric(SubstrateBackgroundFrequency[1:(nrow(SubstrateBackgroundFrequency)),17]), na.rm=TRUE)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
49 Tmean<-mean(as.numeric(SubstrateBackgroundFrequency[1:(nrow(SubstrateBackgroundFrequency)),18]), na.rm=TRUE)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
50 Vmean<-mean(as.numeric(SubstrateBackgroundFrequency[1:(nrow(SubstrateBackgroundFrequency)),19]), na.rm=TRUE)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
51 Wmean<-mean(as.numeric(SubstrateBackgroundFrequency[1:(nrow(SubstrateBackgroundFrequency)),20]), na.rm=TRUE)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
52 Ymean<-mean(as.numeric(SubstrateBackgroundFrequency[1:(nrow(SubstrateBackgroundFrequency)),21]), na.rm=TRUE)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
53 AllMeans<-c(Amean,Cmean,Dmean,Emean,Fmean,Gmean,Hmean,Imean,Kmean,Lmean,Mmean,Nmean,Pmean,Qmean,Rmean,Smean,Tmean,Vmean,Wmean,Ymean)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
54 Asd<-sd(as.numeric(SubstrateBackgroundFrequency[1:(nrow(SubstrateBackgroundFrequency)),2]), na.rm=TRUE)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
55 Csd<-sd(as.numeric(SubstrateBackgroundFrequency[1:(nrow(SubstrateBackgroundFrequency)),3]), na.rm=TRUE)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
56 Dsd<-sd(as.numeric(SubstrateBackgroundFrequency[1:(nrow(SubstrateBackgroundFrequency)),4]), na.rm=TRUE)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
57 Esd<-sd(as.numeric(SubstrateBackgroundFrequency[1:(nrow(SubstrateBackgroundFrequency)),5]), na.rm=TRUE)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
58 Fsd<-sd(as.numeric(SubstrateBackgroundFrequency[1:(nrow(SubstrateBackgroundFrequency)),6]), na.rm=TRUE)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
59 Gsd<-sd(as.numeric(SubstrateBackgroundFrequency[1:(nrow(SubstrateBackgroundFrequency)),7]), na.rm=TRUE)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
60 Hsd<-sd(as.numeric(SubstrateBackgroundFrequency[1:(nrow(SubstrateBackgroundFrequency)),8]), na.rm=TRUE)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
61 Isd<-sd(as.numeric(SubstrateBackgroundFrequency[1:(nrow(SubstrateBackgroundFrequency)),9]), na.rm=TRUE)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
62 Ksd<-sd(as.numeric(SubstrateBackgroundFrequency[1:(nrow(SubstrateBackgroundFrequency)),10]), na.rm=TRUE)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
63 Lsd<-sd(as.numeric(SubstrateBackgroundFrequency[1:(nrow(SubstrateBackgroundFrequency)),11]), na.rm=TRUE)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
64 Msd<-sd(as.numeric(SubstrateBackgroundFrequency[1:(nrow(SubstrateBackgroundFrequency)),12]), na.rm=TRUE)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
65 Nsd<-sd(as.numeric(SubstrateBackgroundFrequency[1:(nrow(SubstrateBackgroundFrequency)),13]), na.rm=TRUE)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
66 Psd<-sd(as.numeric(SubstrateBackgroundFrequency[1:(nrow(SubstrateBackgroundFrequency)),14]), na.rm=TRUE)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
67 Qsd<-sd(as.numeric(SubstrateBackgroundFrequency[1:(nrow(SubstrateBackgroundFrequency)),15]), na.rm=TRUE)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
68 Rsd<-sd(as.numeric(SubstrateBackgroundFrequency[1:(nrow(SubstrateBackgroundFrequency)),16]), na.rm=TRUE)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
69 Ssd<-sd(as.numeric(SubstrateBackgroundFrequency[1:(nrow(SubstrateBackgroundFrequency)),17]), na.rm=TRUE)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
70 Tsd<-sd(as.numeric(SubstrateBackgroundFrequency[1:(nrow(SubstrateBackgroundFrequency)),18]), na.rm=TRUE)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
71 Vsd<-sd(as.numeric(SubstrateBackgroundFrequency[1:(nrow(SubstrateBackgroundFrequency)),19]), na.rm=TRUE)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
72 Wsd<-sd(as.numeric(SubstrateBackgroundFrequency[1:(nrow(SubstrateBackgroundFrequency)),20]), na.rm=TRUE)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
73 Ysd<-sd(as.numeric(SubstrateBackgroundFrequency[1:(nrow(SubstrateBackgroundFrequency)),21]), na.rm=TRUE)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
74 }
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
75 AllSDs<-c(Asd,Csd,Dsd,Esd,Fsd,Gsd,Hsd,Isd,Ksd,Lsd,Msd,Nsd,Psd,Qsd,Rsd,Ssd,Tsd,Vsd,Wsd,Ysd)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
76 #this is subbackfreq SDs
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
77 SBF_statisticalvalues<-cbind(AllMeans,AllSDs)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
78 #create the percent table
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
79 if (1==1){
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
80 Column1<-substrates[,1]
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
81 Column2<-substrates[,2]
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
82 Column3<-substrates[,3]
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
83 Column4<-substrates[,4]
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
84 Column5<-substrates[,5]
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
85 Column6<-substrates[,6]
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
86 Column7<-substrates[,7]
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
87 Column8<-substrates[,8]
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
88 Column9<-substrates[,9]
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
89 Column10<-substrates[,10]
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
90 Column11<-substrates[,11]
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
91 Column12<-substrates[,12]
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
92 Column13<-substrates[,13]
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
93 Column14<-substrates[,14]
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
94 Column15<-substrates[,15]
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
95 spaces1<-sum((Column1%in% ""))
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
96 spaces2<-sum(Column2%in% "")
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
97 spaces3<-sum(Column3%in% "")
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
98 spaces4<-sum(Column4%in% "")
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
99 spaces5<-sum(Column5%in% "")
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
100 spaces6<-sum(Column6%in% "")
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
101 spaces7<-sum(Column7%in% "")
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
102 spaces8<-sum(Column8%in% "")
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
103 spaces9<-sum(Column9%in% "")
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
104 spaces10<-sum(Column10%in% "")
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
105 spaces11<-sum(Column11%in% "")
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
106 spaces12<-sum(Column12%in% "")
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
107 spaces13<-sum(Column13%in% "")
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
108 spaces14<-sum(Column14%in% "")
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
109 spaces15<-sum(Column15%in% "")
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
110 A1<-sum(Column1 %in% "A")/(length(Column1)-spaces1)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
111 A2<-sum(Column2 %in% "A")/(length(Column2)-spaces2)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
112 A3<-sum(Column3 %in% "A")/(length(Column3)-spaces3)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
113 A4<-sum(Column4 %in% "A")/(length(Column4)-spaces4)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
114 A5<-sum(Column5 %in% "A")/(length(Column5)-spaces5)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
115 A6<-sum(Column6 %in% "A")/(length(Column6)-spaces6)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
116 A7<-sum(Column7 %in% "A")/(length(Column7)-spaces7)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
117 A8<-sum(Column8 %in% "A")/(length(Column8)-spaces8)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
118 A9<-sum(Column9 %in% "A")/(length(Column9)-spaces9)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
119 A10<-sum(Column10 %in% "A")/(length(Column10)-spaces10)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
120 A11<-sum(Column11 %in% "A")/(length(Column11)-spaces11)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
121 A12<-sum(Column12 %in% "A")/(length(Column12)-spaces12)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
122 A13<-sum(Column13 %in% "A")/(length(Column13)-spaces13)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
123 A14<-sum(Column14 %in% "A")/(length(Column14)-spaces14)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
124 A15<-sum(Column15 %in% "A")/(length(Column15)-spaces15)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
125 AllAs<-cbind(A1,A2,A3,A4,A5,A6,A7,A8,A9,A10,A11,A12,A13,A14,A15)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
126 C1<-sum(Column1 %in% "C")/(length(Column1)-spaces1)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
127 C2<-sum(Column2 %in% "C")/(length(Column2)-spaces2)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
128 C3<-sum(Column3 %in% "C")/(length(Column3)-spaces3)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
129 C4<-sum(Column4 %in% "C")/(length(Column4)-spaces4)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
130 C5<-sum(Column5 %in% "C")/(length(Column5)-spaces5)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
131 C6<-sum(Column6 %in% "C")/(length(Column6)-spaces6)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
132 C7<-sum(Column7 %in% "C")/(length(Column7)-spaces7)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
133 C8<-sum(Column8 %in% "C")/(length(Column8)-spaces8)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
134 C9<-sum(Column9 %in% "C")/(length(Column9)-spaces9)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
135 C10<-sum(Column10 %in% "C")/(length(Column10)-spaces10)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
136 C11<-sum(Column11 %in% "C")/(length(Column11)-spaces11)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
137 C12<-sum(Column12 %in% "C")/(length(Column12)-spaces12)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
138 C13<-sum(Column13 %in% "C")/(length(Column13)-spaces13)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
139 C14<-sum(Column14 %in% "C")/(length(Column14)-spaces14)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
140 C15<-sum(Column15 %in% "C")/(length(Column15)-spaces15)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
141 CllCs<-cbind(C1,C2,C3,C4,C5,C6,C7,C8,C9,C10,C11,C12,C13,C14,C15)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
142 D1<-sum(Column1 %in% "D")/(length(Column1)-spaces1)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
143 D2<-sum(Column2 %in% "D")/(length(Column2)-spaces2)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
144 D3<-sum(Column3 %in% "D")/(length(Column3)-spaces3)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
145 D4<-sum(Column4 %in% "D")/(length(Column4)-spaces4)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
146 D5<-sum(Column5 %in% "D")/(length(Column5)-spaces5)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
147 D6<-sum(Column6 %in% "D")/(length(Column6)-spaces6)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
148 D7<-sum(Column7 %in% "D")/(length(Column7)-spaces7)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
149 D8<-sum(Column8 %in% "D")/(length(Column8)-spaces8)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
150 D9<-sum(Column9 %in% "D")/(length(Column9)-spaces9)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
151 D10<-sum(Column10 %in% "D")/(length(Column10)-spaces10)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
152 D11<-sum(Column11 %in% "D")/(length(Column11)-spaces11)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
153 D12<-sum(Column12 %in% "D")/(length(Column12)-spaces12)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
154 D13<-sum(Column13 %in% "D")/(length(Column13)-spaces13)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
155 D14<-sum(Column14 %in% "D")/(length(Column14)-spaces14)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
156 D15<-sum(Column15 %in% "D")/(length(Column15)-spaces15)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
157 DllDs<-cbind(D1,D2,D3,D4,D5,D6,D7,D8,D9,D10,D11,D12,D13,D14,D15)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
158 E1<-sum(Column1 %in% "E")/(length(Column1)-spaces1)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
159 E2<-sum(Column2 %in% "E")/(length(Column2)-spaces2)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
160 E3<-sum(Column3 %in% "E")/(length(Column3)-spaces3)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
161 E4<-sum(Column4 %in% "E")/(length(Column4)-spaces4)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
162 E5<-sum(Column5 %in% "E")/(length(Column5)-spaces5)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
163 E6<-sum(Column6 %in% "E")/(length(Column6)-spaces6)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
164 E7<-sum(Column7 %in% "E")/(length(Column7)-spaces7)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
165 E8<-sum(Column8 %in% "E")/(length(Column8)-spaces8)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
166 E9<-sum(Column9 %in% "E")/(length(Column9)-spaces9)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
167 E10<-sum(Column10 %in% "E")/(length(Column10)-spaces10)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
168 E11<-sum(Column11 %in% "E")/(length(Column11)-spaces11)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
169 E12<-sum(Column12 %in% "E")/(length(Column12)-spaces12)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
170 E13<-sum(Column13 %in% "E")/(length(Column13)-spaces13)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
171 E14<-sum(Column14 %in% "E")/(length(Column14)-spaces14)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
172 E15<-sum(Column15 %in% "E")/(length(Column15)-spaces15)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
173 EllEs<-cbind(E1,E2,E3,E4,E5,E6,E7,E8,E9,E10,E11,E12,E13,E14,E15)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
174 F1<-sum(Column1 %in% "F")/(length(Column1)-spaces1)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
175 F2<-sum(Column2 %in% "F")/(length(Column2)-spaces2)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
176 F3<-sum(Column3 %in% "F")/(length(Column3)-spaces3)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
177 F4<-sum(Column4 %in% "F")/(length(Column4)-spaces4)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
178 F5<-sum(Column5 %in% "F")/(length(Column5)-spaces5)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
179 F6<-sum(Column6 %in% "F")/(length(Column6)-spaces6)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
180 F7<-sum(Column7 %in% "F")/(length(Column7)-spaces7)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
181 F8<-sum(Column8 %in% "F")/(length(Column8)-spaces8)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
182 F9<-sum(Column9 %in% "F")/(length(Column9)-spaces9)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
183 F10<-sum(Column10 %in% "F")/(length(Column10)-spaces10)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
184 F11<-sum(Column11 %in% "F")/(length(Column11)-spaces11)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
185 F12<-sum(Column12 %in% "F")/(length(Column12)-spaces12)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
186 F13<-sum(Column13 %in% "F")/(length(Column13)-spaces13)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
187 F14<-sum(Column14 %in% "F")/(length(Column14)-spaces14)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
188 F15<-sum(Column15 %in% "F")/(length(Column15)-spaces15)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
189 FllFs<-cbind(F1,F2,F3,F4,F5,F6,F7,F8,F9,F10,F11,F12,F13,F14,F15)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
190 G1<-sum(Column1 %in% "G")/(length(Column1)-spaces1)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
191 G2<-sum(Column2 %in% "G")/(length(Column2)-spaces2)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
192 G3<-sum(Column3 %in% "G")/(length(Column3)-spaces3)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
193 G4<-sum(Column4 %in% "G")/(length(Column4)-spaces4)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
194 G5<-sum(Column5 %in% "G")/(length(Column5)-spaces5)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
195 G6<-sum(Column6 %in% "G")/(length(Column6)-spaces6)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
196 G7<-sum(Column7 %in% "G")/(length(Column7)-spaces7)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
197 G8<-sum(Column8 %in% "G")/(length(Column8)-spaces8)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
198 G9<-sum(Column9 %in% "G")/(length(Column9)-spaces9)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
199 G10<-sum(Column10 %in% "G")/(length(Column10)-spaces10)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
200 G11<-sum(Column11 %in% "G")/(length(Column11)-spaces11)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
201 G12<-sum(Column12 %in% "G")/(length(Column12)-spaces12)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
202 G13<-sum(Column13 %in% "G")/(length(Column13)-spaces13)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
203 G14<-sum(Column14 %in% "G")/(length(Column14)-spaces14)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
204 G15<-sum(Column15 %in% "G")/(length(Column15)-spaces15)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
205 GllGs<-cbind(G1,G2,G3,G4,G5,G6,G7,G8,G9,G10,G11,G12,G13,G14,G15)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
206 H1<-sum(Column1 %in% "H")/(length(Column1)-spaces1)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
207 H2<-sum(Column2 %in% "H")/(length(Column2)-spaces2)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
208 H3<-sum(Column3 %in% "H")/(length(Column3)-spaces3)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
209 H4<-sum(Column4 %in% "H")/(length(Column4)-spaces4)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
210 H5<-sum(Column5 %in% "H")/(length(Column5)-spaces5)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
211 H6<-sum(Column6 %in% "H")/(length(Column6)-spaces6)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
212 H7<-sum(Column7 %in% "H")/(length(Column7)-spaces7)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
213 H8<-sum(Column8 %in% "H")/(length(Column8)-spaces8)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
214 H9<-sum(Column9 %in% "H")/(length(Column9)-spaces9)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
215 H10<-sum(Column10 %in% "H")/(length(Column10)-spaces10)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
216 H11<-sum(Column11 %in% "H")/(length(Column11)-spaces11)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
217 H12<-sum(Column12 %in% "H")/(length(Column12)-spaces12)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
218 H13<-sum(Column13 %in% "H")/(length(Column13)-spaces13)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
219 H14<-sum(Column14 %in% "H")/(length(Column14)-spaces14)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
220 H15<-sum(Column15 %in% "H")/(length(Column15)-spaces15)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
221 HllHs<-cbind(H1,H2,H3,H4,H5,H6,H7,H8,H9,H10,H11,H12,H13,H14,H15)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
222 I1<-sum(Column1 %in% "I")/(length(Column1)-spaces1)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
223 I2<-sum(Column2 %in% "I")/(length(Column2)-spaces2)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
224 I3<-sum(Column3 %in% "I")/(length(Column3)-spaces3)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
225 I4<-sum(Column4 %in% "I")/(length(Column4)-spaces4)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
226 I5<-sum(Column5 %in% "I")/(length(Column5)-spaces5)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
227 I6<-sum(Column6 %in% "I")/(length(Column6)-spaces6)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
228 I7<-sum(Column7 %in% "I")/(length(Column7)-spaces7)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
229 I8<-sum(Column8 %in% "I")/(length(Column8)-spaces8)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
230 I9<-sum(Column9 %in% "I")/(length(Column9)-spaces9)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
231 I10<-sum(Column10 %in% "I")/(length(Column10)-spaces10)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
232 I11<-sum(Column11 %in% "I")/(length(Column11)-spaces11)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
233 I12<-sum(Column12 %in% "I")/(length(Column12)-spaces12)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
234 I13<-sum(Column13 %in% "I")/(length(Column13)-spaces13)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
235 I14<-sum(Column14 %in% "I")/(length(Column14)-spaces14)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
236 I15<-sum(Column15 %in% "I")/(length(Column15)-spaces15)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
237 IllIs<-cbind(I1,I2,I3,I4,I5,I6,I7,I8,I9,I10,I11,I12,I13,I14,I15)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
238 K1<-sum(Column1 %in% "K")/(length(Column1)-spaces1)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
239 K2<-sum(Column2 %in% "K")/(length(Column2)-spaces2)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
240 K3<-sum(Column3 %in% "K")/(length(Column3)-spaces3)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
241 K4<-sum(Column4 %in% "K")/(length(Column4)-spaces4)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
242 K5<-sum(Column5 %in% "K")/(length(Column5)-spaces5)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
243 K6<-sum(Column6 %in% "K")/(length(Column6)-spaces6)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
244 K7<-sum(Column7 %in% "K")/(length(Column7)-spaces7)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
245 K8<-sum(Column8 %in% "K")/(length(Column8)-spaces8)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
246 K9<-sum(Column9 %in% "K")/(length(Column9)-spaces9)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
247 K10<-sum(Column10 %in% "K")/(length(Column10)-spaces10)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
248 K11<-sum(Column11 %in% "K")/(length(Column11)-spaces11)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
249 K12<-sum(Column12 %in% "K")/(length(Column12)-spaces12)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
250 K13<-sum(Column13 %in% "K")/(length(Column13)-spaces13)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
251 K14<-sum(Column14 %in% "K")/(length(Column14)-spaces14)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
252 K15<-sum(Column15 %in% "K")/(length(Column15)-spaces15)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
253 KllKs<-cbind(K1,K2,K3,K4,K5,K6,K7,K8,K9,K10,K11,K12,K13,K14,K15)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
254 L1<-sum(Column1 %in% "L")/(length(Column1)-spaces1)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
255 L2<-sum(Column2 %in% "L")/(length(Column2)-spaces2)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
256 L3<-sum(Column3 %in% "L")/(length(Column3)-spaces3)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
257 L4<-sum(Column4 %in% "L")/(length(Column4)-spaces4)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
258 L5<-sum(Column5 %in% "L")/(length(Column5)-spaces5)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
259 L6<-sum(Column6 %in% "L")/(length(Column6)-spaces6)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
260 L7<-sum(Column7 %in% "L")/(length(Column7)-spaces7)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
261 L8<-sum(Column8 %in% "L")/(length(Column8)-spaces8)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
262 L9<-sum(Column9 %in% "L")/(length(Column9)-spaces9)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
263 L10<-sum(Column10 %in% "L")/(length(Column10)-spaces10)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
264 L11<-sum(Column11 %in% "L")/(length(Column11)-spaces11)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
265 L12<-sum(Column12 %in% "L")/(length(Column12)-spaces12)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
266 L13<-sum(Column13 %in% "L")/(length(Column13)-spaces13)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
267 L14<-sum(Column14 %in% "L")/(length(Column14)-spaces14)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
268 L15<-sum(Column15 %in% "L")/(length(Column15)-spaces15)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
269 LllLs<-cbind(L1,L2,L3,L4,L5,L6,L7,L8,L9,L10,L11,L12,L13,L14,L15)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
270 M1<-sum(Column1 %in% "M")/(length(Column1)-spaces1)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
271 M2<-sum(Column2 %in% "M")/(length(Column2)-spaces2)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
272 M3<-sum(Column3 %in% "M")/(length(Column3)-spaces3)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
273 M4<-sum(Column4 %in% "M")/(length(Column4)-spaces4)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
274 M5<-sum(Column5 %in% "M")/(length(Column5)-spaces5)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
275 M6<-sum(Column6 %in% "M")/(length(Column6)-spaces6)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
276 M7<-sum(Column7 %in% "M")/(length(Column7)-spaces7)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
277 M8<-sum(Column8 %in% "M")/(length(Column8)-spaces8)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
278 M9<-sum(Column9 %in% "M")/(length(Column9)-spaces9)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
279 M10<-sum(Column10 %in% "M")/(length(Column10)-spaces10)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
280 M11<-sum(Column11 %in% "M")/(length(Column11)-spaces11)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
281 M12<-sum(Column12 %in% "M")/(length(Column12)-spaces12)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
282 M13<-sum(Column13 %in% "M")/(length(Column13)-spaces13)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
283 M14<-sum(Column14 %in% "M")/(length(Column14)-spaces14)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
284 M15<-sum(Column15 %in% "M")/(length(Column15)-spaces15)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
285 MllMs<-cbind(M1,M2,M3,M4,M5,M6,M7,M8,M9,M10,M11,M12,M13,M14,M15)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
286 N1<-sum(Column1 %in% "N")/(length(Column1)-spaces1)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
287 N2<-sum(Column2 %in% "N")/(length(Column2)-spaces2)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
288 N3<-sum(Column3 %in% "N")/(length(Column3)-spaces3)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
289 N4<-sum(Column4 %in% "N")/(length(Column4)-spaces4)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
290 N5<-sum(Column5 %in% "N")/(length(Column5)-spaces5)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
291 N6<-sum(Column6 %in% "N")/(length(Column6)-spaces6)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
292 N7<-sum(Column7 %in% "N")/(length(Column7)-spaces7)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
293 N8<-sum(Column8 %in% "N")/(length(Column8)-spaces8)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
294 N9<-sum(Column9 %in% "N")/(length(Column9)-spaces9)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
295 N10<-sum(Column10 %in% "N")/(length(Column10)-spaces10)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
296 N11<-sum(Column11 %in% "N")/(length(Column11)-spaces11)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
297 N12<-sum(Column12 %in% "N")/(length(Column12)-spaces12)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
298 N13<-sum(Column13 %in% "N")/(length(Column13)-spaces13)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
299 N14<-sum(Column14 %in% "N")/(length(Column14)-spaces14)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
300 N15<-sum(Column15 %in% "N")/(length(Column15)-spaces15)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
301 NllNs<-cbind(N1,N2,N3,N4,N5,N6,N7,N8,N9,N10,N11,N12,N13,N14,N15)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
302 P1<-sum(Column1 %in% "P")/(length(Column1)-spaces1)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
303 P2<-sum(Column2 %in% "P")/(length(Column2)-spaces2)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
304 P3<-sum(Column3 %in% "P")/(length(Column3)-spaces3)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
305 P4<-sum(Column4 %in% "P")/(length(Column4)-spaces4)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
306 P5<-sum(Column5 %in% "P")/(length(Column5)-spaces5)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
307 P6<-sum(Column6 %in% "P")/(length(Column6)-spaces6)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
308 P7<-sum(Column7 %in% "P")/(length(Column7)-spaces7)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
309 P8<-sum(Column8 %in% "P")/(length(Column8)-spaces8)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
310 P9<-sum(Column9 %in% "P")/(length(Column9)-spaces9)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
311 P10<-sum(Column10 %in% "P")/(length(Column10)-spaces10)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
312 P11<-sum(Column11 %in% "P")/(length(Column11)-spaces11)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
313 P12<-sum(Column12 %in% "P")/(length(Column12)-spaces12)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
314 P13<-sum(Column13 %in% "P")/(length(Column13)-spaces13)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
315 P14<-sum(Column14 %in% "P")/(length(Column14)-spaces14)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
316 P15<-sum(Column15 %in% "P")/(length(Column15)-spaces15)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
317 PllPs<-cbind(P1,P2,P3,P4,P5,P6,P7,P8,P9,P10,P11,P12,P13,P14,P15)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
318 Q1<-sum(Column1 %in% "Q")/(length(Column1)-spaces1)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
319 Q2<-sum(Column2 %in% "Q")/(length(Column2)-spaces2)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
320 Q3<-sum(Column3 %in% "Q")/(length(Column3)-spaces3)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
321 Q4<-sum(Column4 %in% "Q")/(length(Column4)-spaces4)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
322 Q5<-sum(Column5 %in% "Q")/(length(Column5)-spaces5)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
323 Q6<-sum(Column6 %in% "Q")/(length(Column6)-spaces6)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
324 Q7<-sum(Column7 %in% "Q")/(length(Column7)-spaces7)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
325 Q8<-sum(Column8 %in% "Q")/(length(Column8)-spaces8)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
326 Q9<-sum(Column9 %in% "Q")/(length(Column9)-spaces9)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
327 Q10<-sum(Column10 %in% "Q")/(length(Column10)-spaces10)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
328 Q11<-sum(Column11 %in% "Q")/(length(Column11)-spaces11)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
329 Q12<-sum(Column12 %in% "Q")/(length(Column12)-spaces12)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
330 Q13<-sum(Column13 %in% "Q")/(length(Column13)-spaces13)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
331 Q14<-sum(Column14 %in% "Q")/(length(Column14)-spaces14)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
332 Q15<-sum(Column15 %in% "Q")/(length(Column15)-spaces15)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
333 QllQs<-cbind(Q1,Q2,Q3,Q4,Q5,Q6,Q7,Q8,Q9,Q10,Q11,Q12,Q13,Q14,Q15)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
334 R1<-sum(Column1 %in% "R")/(length(Column1)-spaces1)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
335 R2<-sum(Column2 %in% "R")/(length(Column2)-spaces2)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
336 R3<-sum(Column3 %in% "R")/(length(Column3)-spaces3)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
337 R4<-sum(Column4 %in% "R")/(length(Column4)-spaces4)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
338 R5<-sum(Column5 %in% "R")/(length(Column5)-spaces5)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
339 R6<-sum(Column6 %in% "R")/(length(Column6)-spaces6)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
340 R7<-sum(Column7 %in% "R")/(length(Column7)-spaces7)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
341 R8<-sum(Column8 %in% "R")/(length(Column8)-spaces8)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
342 R9<-sum(Column9 %in% "R")/(length(Column9)-spaces9)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
343 R10<-sum(Column10 %in% "R")/(length(Column10)-spaces10)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
344 R11<-sum(Column11 %in% "R")/(length(Column11)-spaces11)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
345 R12<-sum(Column12 %in% "R")/(length(Column12)-spaces12)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
346 R13<-sum(Column13 %in% "R")/(length(Column13)-spaces13)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
347 R14<-sum(Column14 %in% "R")/(length(Column14)-spaces14)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
348 R15<-sum(Column15 %in% "R")/(length(Column15)-spaces15)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
349 RllRs<-cbind(R1,R2,R3,R4,R5,R6,R7,R8,R9,R10,R11,R12,R13,R14,R15)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
350 S1<-sum(Column1 %in% "S")/(length(Column1)-spaces1)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
351 S2<-sum(Column2 %in% "S")/(length(Column2)-spaces2)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
352 S3<-sum(Column3 %in% "S")/(length(Column3)-spaces3)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
353 S4<-sum(Column4 %in% "S")/(length(Column4)-spaces4)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
354 S5<-sum(Column5 %in% "S")/(length(Column5)-spaces5)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
355 S6<-sum(Column6 %in% "S")/(length(Column6)-spaces6)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
356 S7<-sum(Column7 %in% "S")/(length(Column7)-spaces7)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
357 S8<-sum(Column8 %in% "S")/(length(Column8)-spaces8)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
358 S9<-sum(Column9 %in% "S")/(length(Column9)-spaces9)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
359 S10<-sum(Column10 %in% "S")/(length(Column10)-spaces10)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
360 S11<-sum(Column11 %in% "S")/(length(Column11)-spaces11)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
361 S12<-sum(Column12 %in% "S")/(length(Column12)-spaces12)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
362 S13<-sum(Column13 %in% "S")/(length(Column13)-spaces13)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
363 S14<-sum(Column14 %in% "S")/(length(Column14)-spaces14)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
364 S15<-sum(Column15 %in% "S")/(length(Column15)-spaces15)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
365 SllSs<-cbind(S1,S2,S3,S4,S5,S6,S7,S8,S9,S10,S11,S12,S13,S14,S15)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
366 T1<-sum(Column1 %in% "T")/(length(Column1)-spaces1)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
367 T2<-sum(Column2 %in% "T")/(length(Column2)-spaces2)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
368 T3<-sum(Column3 %in% "T")/(length(Column3)-spaces3)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
369 T4<-sum(Column4 %in% "T")/(length(Column4)-spaces4)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
370 T5<-sum(Column5 %in% "T")/(length(Column5)-spaces5)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
371 T6<-sum(Column6 %in% "T")/(length(Column6)-spaces6)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
372 T7<-sum(Column7 %in% "T")/(length(Column7)-spaces7)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
373 T8<-sum(Column8 %in% "T")/(length(Column8)-spaces8)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
374 T9<-sum(Column9 %in% "T")/(length(Column9)-spaces9)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
375 T10<-sum(Column10 %in% "T")/(length(Column10)-spaces10)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
376 T11<-sum(Column11 %in% "T")/(length(Column11)-spaces11)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
377 T12<-sum(Column12 %in% "T")/(length(Column12)-spaces12)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
378 T13<-sum(Column13 %in% "T")/(length(Column13)-spaces13)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
379 T14<-sum(Column14 %in% "T")/(length(Column14)-spaces14)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
380 T15<-sum(Column15 %in% "T")/(length(Column15)-spaces15)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
381 TllTs<-cbind(T1,T2,T3,T4,T5,T6,T7,T8,T9,T10,T11,T12,T13,T14,T15)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
382 V1<-sum(Column1 %in% "V")/(length(Column1)-spaces1)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
383 V2<-sum(Column2 %in% "V")/(length(Column2)-spaces2)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
384 V3<-sum(Column3 %in% "V")/(length(Column3)-spaces3)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
385 V4<-sum(Column4 %in% "V")/(length(Column4)-spaces4)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
386 V5<-sum(Column5 %in% "V")/(length(Column5)-spaces5)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
387 V6<-sum(Column6 %in% "V")/(length(Column6)-spaces6)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
388 V7<-sum(Column7 %in% "V")/(length(Column7)-spaces7)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
389 V8<-sum(Column8 %in% "V")/(length(Column8)-spaces8)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
390 V9<-sum(Column9 %in% "V")/(length(Column9)-spaces9)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
391 V10<-sum(Column10 %in% "V")/(length(Column10)-spaces10)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
392 V11<-sum(Column11 %in% "V")/(length(Column11)-spaces11)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
393 V12<-sum(Column12 %in% "V")/(length(Column12)-spaces12)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
394 V13<-sum(Column13 %in% "V")/(length(Column13)-spaces13)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
395 V14<-sum(Column14 %in% "V")/(length(Column14)-spaces14)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
396 V15<-sum(Column15 %in% "V")/(length(Column15)-spaces15)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
397 VllVs<-cbind(V1,V2,V3,V4,V5,V6,V7,V8,V9,V10,V11,V12,V13,V14,V15)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
398 W1<-sum(Column1 %in% "W")/(length(Column1)-spaces1)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
399 W2<-sum(Column2 %in% "W")/(length(Column2)-spaces2)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
400 W3<-sum(Column3 %in% "W")/(length(Column3)-spaces3)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
401 W4<-sum(Column4 %in% "W")/(length(Column4)-spaces4)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
402 W5<-sum(Column5 %in% "W")/(length(Column5)-spaces5)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
403 W6<-sum(Column6 %in% "W")/(length(Column6)-spaces6)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
404 W7<-sum(Column7 %in% "W")/(length(Column7)-spaces7)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
405 W8<-sum(Column8 %in% "W")/(length(Column8)-spaces8)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
406 W9<-sum(Column9 %in% "W")/(length(Column9)-spaces9)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
407 W10<-sum(Column10 %in% "W")/(length(Column10)-spaces10)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
408 W11<-sum(Column11 %in% "W")/(length(Column11)-spaces11)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
409 W12<-sum(Column12 %in% "W")/(length(Column12)-spaces12)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
410 W13<-sum(Column13 %in% "W")/(length(Column13)-spaces13)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
411 W14<-sum(Column14 %in% "W")/(length(Column14)-spaces14)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
412 W15<-sum(Column15 %in% "W")/(length(Column15)-spaces15)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
413 WllWs<-cbind(W1,W2,W3,W4,W5,W6,W7,W8,W9,W10,W11,W12,W13,W14,W15)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
414 Y1<-sum(Column1 %in% "Y")/(length(Column1)-spaces1)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
415 Y2<-sum(Column2 %in% "Y")/(length(Column2)-spaces2)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
416 Y3<-sum(Column3 %in% "Y")/(length(Column3)-spaces3)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
417 Y4<-sum(Column4 %in% "Y")/(length(Column4)-spaces4)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
418 Y5<-sum(Column5 %in% "Y")/(length(Column5)-spaces5)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
419 Y6<-sum(Column6 %in% "Y")/(length(Column6)-spaces6)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
420 Y7<-sum(Column7 %in% "Y")/(length(Column7)-spaces7)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
421 Y8<-sum(Column8 %in% "Y")/(length(Column8)-spaces8)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
422 Y9<-sum(Column9 %in% "Y")/(length(Column9)-spaces9)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
423 Y10<-sum(Column10 %in% "Y")/(length(Column10)-spaces10)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
424 Y11<-sum(Column11 %in% "Y")/(length(Column11)-spaces11)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
425 Y12<-sum(Column12 %in% "Y")/(length(Column12)-spaces12)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
426 Y13<-sum(Column13 %in% "Y")/(length(Column13)-spaces13)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
427 Y14<-sum(Column14 %in% "Y")/(length(Column14)-spaces14)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
428 Y15<-sum(Column15 %in% "Y")/(length(Column15)-spaces15)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
429 YllYs<-cbind(Y1,Y2,Y3,Y4,Y5,Y6,Y7,Y8,Y9,Y10,Y11,Y12,Y13,Y14,Y15)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
430 }
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
431 #this is substrate percents
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
432 #A C D E F G H I K L N P Q R S T V W Y
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
433 PercentTable<-rbind(AllAs,CllCs,DllDs,EllEs,FllFs,GllGs,HllHs,IllIs,KllKs,LllLs,MllMs,NllNs,PllPs,QllQs,RllRs,SllSs,TllTs,VllVs,WllWs,YllYs)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
434 PercentTable<-PercentTable*100
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
435 #create the SD table
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
436 SDtable<-matrix(data = rep(1,times=(nrow(PercentTable)*ncol(PercentTable))),nrow = nrow(PercentTable),ncol = ncol(PercentTable))
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
437 #for every row, a percertage minus the same mean over the same SD
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
438 if(1==1){
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
439 SDtable[1,]<-(PercentTable[1,]-Amean)/Asd
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
440 SDtable[2,]<-(PercentTable[2,]-Cmean)/Csd
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
441 SDtable[3,]<-(PercentTable[3,]-Dmean)/Dsd
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
442 SDtable[4,]<-(PercentTable[4,]-Emean)/Esd
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
443 SDtable[5,]<-(PercentTable[5,]-Fmean)/Fsd
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
444 SDtable[6,]<-(PercentTable[6,]-Gmean)/Gsd
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
445 SDtable[7,]<-(PercentTable[7,]-Hmean)/Hsd
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
446 SDtable[8,]<-(PercentTable[8,]-Imean)/Isd
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
447 SDtable[9,]<-(PercentTable[9,]-Kmean)/Ksd
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
448 SDtable[10,]<-(PercentTable[10,]-Lmean)/Lsd
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
449 SDtable[11,]<-(PercentTable[11,]-Mmean)/Msd
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
450 SDtable[12,]<-(PercentTable[12,]-Nmean)/Nsd
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
451 SDtable[13,]<-(PercentTable[13,]-Pmean)/Psd
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
452 SDtable[14,]<-(PercentTable[14,]-Qmean)/Qsd
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
453 SDtable[15,]<-(PercentTable[15,]-Rmean)/Rsd
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
454 SDtable[16,]<-(PercentTable[16,]-Smean)/Ssd
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
455 SDtable[17,]<-(PercentTable[17,]-Tmean)/Tsd
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
456 SDtable[18,]<-(PercentTable[18,]-Vmean)/Vsd
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
457 SDtable[19,]<-(PercentTable[19,]-Wmean)/Wsd
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
458 SDtable[20,]<-(PercentTable[20,]-Ymean)/Ysd
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
459 }
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
460 SetOfAAs<-c("Letter","A","C","D","E","F","G","H","I","K","L","M","N","P","Q","R","S","T","V","W","Y")
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
461 SumOfSigmaAAs<-c(1:15)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
462 for (i in 1:15){
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
463 SumOfSigmasValue<-0
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
464 for (j in 1:20){
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
465 value<-0
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
466 if (SDtable[j,i]>2){
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
467 value<-sum(substrates[,i]==SetOfAAs[j])
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
468 }
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
469 SumOfSigmasValue<-SumOfSigmasValue+value
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
470 }
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
471 SumOfSigmaAAs[i]<-SumOfSigmasValue
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
472 }
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
473 # AAs1<-length(substrates[,1])-sum(substrates[,1]=="")
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
474 # AAs2<-length(substrates[,2])-sum(substrates[,2]=="")
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
475 # AAs3<-length(substrates[,3])-sum(substrates[,3]=="")
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
476 # AAs4<-length(substrates[,4])-sum(substrates[,4]=="")
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
477 # AAs5<-length(substrates[,5])-sum(substrates[,5]=="")
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
478 # AAs6<-length(substrates[,6])-sum(substrates[,6]=="")
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
479 # AAs7<-length(substrates[,7])-sum(substrates[,7]=="")
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
480 # AAs8<-length(substrates[,8])-sum(substrates[,8]=="")
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
481 # AAs9<-length(substrates[,9])-sum(substrates[,9]=="")
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
482 #
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
483 #
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
484 #
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
485 # #AAsAtPositions<-c(AAs1,AAs2,AAs3,AAs4,AAs5,AAs6,AAs7,AAs8,AAs9)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
486 # AAsAtPositions<-c(length(substrates[,1]),length(substrates[,2]),length(substrates[,3]),length(substrates[,4]),
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
487 # length(substrates[,5]),length(substrates[,6]),length(substrates[,7]),length(substrates[,8]),
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
488 # length(substrates[,9]))
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
489 SumOfExpectedSigmaAAs<-c(1:15)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
490 for (i in 1:15){
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
491 ExpectedValue<-0
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
492 for (j in 1:20){
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
493 value<-0
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
494 if (SDtable[j,i]>2){
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
495 value<-AllMeans[j]
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
496 }
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
497 ExpectedValue<-ExpectedValue+value
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
498 }
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
499 SumOfExpectedSigmaAAs[i]<-ExpectedValue*(length(substrates[,i])-sum(substrates[,i]%in% ""))/100
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
500 }
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
501 SelectivityRow<-SumOfSigmaAAs/SumOfExpectedSigmaAAs
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
502 SelectivitySheet<-rbind(SumOfSigmaAAs,SumOfExpectedSigmaAAs,SelectivityRow)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
503 SetOfAAs<-matrix(data = SetOfAAs,ncol = 1)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
504 SDtable
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
505 Amean
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
506 Cmean
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
507 Wsd
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
508 Amean
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
509 head(SubstrateBackgroundFrequency)
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
510 source('C:/Users/User Name/Desktop/UMN 2017 Spring/Parker Lab/KinatestSwitcheroo.R')
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
511 source('C:/Users/User Name/Desktop/UMN 2017 Spring/Parker Lab/KinatestSwitcheroo.R')
e16ca3e9fc49 Uploaded
jfb
parents:
diff changeset
512 source('C:/Users/User Name/Desktop/UMN 2017 Spring/Parker Lab/KinatestSwitcheroo.R')