Mercurial > repos > jfb > kinatest_r_7_7testing

changeset 19:9e520c365624 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author jfb
date Tue, 17 Apr 2018 17:03:39 -0400
parents e16ca3e9fc49
children f7a6a3ec7710
files kinatestid_r/Kinatest-R.R
diffstat 1 files changed, 6 insertions(+), 8 deletions(-) [+]
line wrap: on
line diff
--- 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)