Mercurial > repos > jfb > kinatest_7_7
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