comparison kinatestid_r/Kinatest-R.R @ 19:9e520c365624 draft

Uploaded

18:e16ca3e9fc49

}


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:9)

for (i in 1:9){
  SumOfSigmasValue<-0
  for (j in 1:20){
    value<-0
    if (SDtable[j,i]>2){
      q=j+1

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:9)
for (i in 1:9){
  ExpectedValue<-0
  for (j in 1:20){
    value<-0
    if (SDtable[j,i]>2){
      value<-AllMeans[j]

  SumOfExpectedSigmaAAs[i]<-ExpectedValue*(length(substrates[,i])-sum(substrates[,i]%in% ""))/100
}

#SelectivityRow<-SumOfSigmaAAs/SumOfExpectedSigmaAAs

SelectivityRow<-rep(0,times=9)

for (q in 1:9) {
  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)
row.names(SDtable)<-NULL
SDtable<-data.frame(SetOfAAs,SDtable)

19:9e520c365624

}


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)

for (i in 1:15){
  SumOfSigmasValue<-0
  for (j in 1:20){
    value<-0
    if (SDtable[j,i]>2){
      q=j+1

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){
  ExpectedValue<-0
  for (j in 1:20){
    value<-0
    if (SDtable[j,i]>2){
      value<-AllMeans[j]

  SumOfExpectedSigmaAAs[i]<-ExpectedValue*(length(substrates[,i])-sum(substrates[,i]%in% ""))/100
}

#SelectivityRow<-SumOfSigmaAAs/SumOfExpectedSigmaAAs

SelectivityRow<-rep(0,times=15)

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)
row.names(SDtable)<-NULL
SDtable<-data.frame(SetOfAAs,SDtable)