# HG changeset patch
# User jfb
# Date 1523999019 14400
# Node ID 9e520c365624d4c9944d5bac24975a9e50342f29
# Parent  e16ca3e9fc4930c452e2e1048b74e4e6ccf59012
Uploaded

diff -r e16ca3e9fc49 -r 9e520c365624 kinatestid_r/Kinatest-R.R
--- a/kinatestid_r/Kinatest-R.R	Tue Apr 17 13:07:41 2018 -0400
+++ b/kinatestid_r/Kinatest-R.R	Tue Apr 17 17:03:39 2018 -0400
@@ -527,10 +527,9 @@
 
 SetOfAAs<-c("Letter","A","C","D","E","F","G","H","I","K","L","M","N","P","Q","R","S","T","V","W","Y")
 
+SumOfSigmaAAs<-c(1:15)
 
-SumOfSigmaAAs<-c(1:9)
-
-for (i in 1:9){
+for (i in 1:15){
   SumOfSigmasValue<-0
   for (j in 1:20){
     value<-0
@@ -557,8 +556,8 @@
 #                   length(substrates[,5]),length(substrates[,6]),length(substrates[,7]),length(substrates[,8]),
 #                   length(substrates[,9]))
 
-SumOfExpectedSigmaAAs<-c(1:9)
-for (i in 1:9){
+SumOfExpectedSigmaAAs<-c(1:15)
+for (i in 1:15){
   ExpectedValue<-0
   for (j in 1:20){
     value<-0
@@ -572,16 +571,15 @@
 
 #SelectivityRow<-SumOfSigmaAAs/SumOfExpectedSigmaAAs
 
-SelectivityRow<-rep(0,times=9)
+SelectivityRow<-rep(0,times=15)
 
-for (q in 1:9) {
+for (q in 1:15) {
   SelectivityRow[q]<-(SumOfSigmaAAs[q]/SumOfExpectedSigmaAAs[q])*(SumOfSigmaAAs[q]/AAsAtPositions[q])
 }
 
 SelectivitySheet<-rbind(SumOfSigmaAAs,SumOfExpectedSigmaAAs,SelectivityRow)
 
 SetOfAAs<-matrix(data = SetOfAAs,ncol = 1)
-
 SDtableu<-SDtable
 HeaderSD<-c(-7:7)
 SDtable<-rbind(HeaderSD,SDtableu)