changeset 11:a66a929246e4 draft default tip

Uploaded
author jfb
date Mon, 14 May 2018 17:32:18 -0400
parents 493658aefce8
children
files kinatestid_r/Kinatest-R.R
diffstat 1 files changed, 16 insertions(+), 17 deletions(-) [+]
line wrap: on
line diff
--- a/kinatestid_r/Kinatest-R.R	Wed Feb 28 13:57:02 2018 -0500
+++ b/kinatestid_r/Kinatest-R.R	Mon May 14 17:32:18 2018 -0400
@@ -535,29 +535,14 @@
   for (j in 1:20){
         value<-0
     if (SDtable[j,i]>2){
-      value<-sum(substrates[,i]==SetOfAAs[j])
+      k<-j+1
+      value<-sum(substrates[,i]==SetOfAAs[k])
     }
     SumOfSigmasValue<-SumOfSigmasValue+value
   }
   SumOfSigmaAAs[i]<-SumOfSigmasValue
 }
 
-# AAs1<-length(substrates[,1])-sum(substrates[,1]=="")
-# AAs2<-length(substrates[,2])-sum(substrates[,2]=="")
-# AAs3<-length(substrates[,3])-sum(substrates[,3]=="")
-# AAs4<-length(substrates[,4])-sum(substrates[,4]=="")
-# AAs5<-length(substrates[,5])-sum(substrates[,5]=="")
-# AAs6<-length(substrates[,6])-sum(substrates[,6]=="")
-# AAs7<-length(substrates[,7])-sum(substrates[,7]=="")
-# AAs8<-length(substrates[,8])-sum(substrates[,8]=="")
-# AAs9<-length(substrates[,9])-sum(substrates[,9]=="")
-# 
-# 
-# 
-# #AAsAtPositions<-c(AAs1,AAs2,AAs3,AAs4,AAs5,AAs6,AAs7,AAs8,AAs9)
-# AAsAtPositions<-c(length(substrates[,1]),length(substrates[,2]),length(substrates[,3]),length(substrates[,4]),
-#                   length(substrates[,5]),length(substrates[,6]),length(substrates[,7]),length(substrates[,8]),
-#                   length(substrates[,9]))
 
 SumOfExpectedSigmaAAs<-c(1:15)
 for (i in 1:15){
@@ -573,8 +558,22 @@
 }
 
 SelectivityRow<-SumOfSigmaAAs/SumOfExpectedSigmaAAs
+#SelectivitySheet<-rbind(SumOfSigmaAAs,SumOfExpectedSigmaAAs,SelectivityRow)
+
+SuperRow<-SelectivityRow
+
+
+for (x in 1:15) {
+  Numerator<-SumOfSigmaAAs[x]
+  Denomenator<-length(substrates[,x])-sum(substrates[,x]%in% "")
+  mult<-Numerator/Denomenator
+  value<-SuperRow[x]*mult
+  SelectivityRow[x]<-value
+}
+
 SelectivitySheet<-rbind(SumOfSigmaAAs,SumOfExpectedSigmaAAs,SelectivityRow)
 
+
 SetOfAAs<-matrix(data = SetOfAAs,ncol = 1)
 
 SDtableu<-SDtable