Mercurial > repos > jfb > kinatest_r_7_7testing
changeset 19:9e520c365624 draft
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)