comparison kinatestid_r/Kinatest-R.R @ 11:a66a929246e4 draft default tip

Uploaded

10:493658aefce8

for (i in 1:15){
  SumOfSigmasValue<-0
  for (j in 1:20){
        value<-0
    if (SDtable[j,i]>2){
      value<-sum(substrates[,i]==SetOfAAs[j])
    }
    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){
  ExpectedValue<-0
  for (j in 1:20){

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

SelectivityRow<-SumOfSigmaAAs/SumOfExpectedSigmaAAs
SelectivitySheet<-rbind(SumOfSigmaAAs,SumOfExpectedSigmaAAs,SelectivityRow)

SetOfAAs<-matrix(data = SetOfAAs,ncol = 1)

SDtableu<-SDtable
HeaderSD<-c(-7:7)

11:a66a929246e4

for (i in 1:15){
  SumOfSigmasValue<-0
  for (j in 1:20){
        value<-0
    if (SDtable[j,i]>2){
      k<-j+1
      value<-sum(substrates[,i]==SetOfAAs[k])
    }
    SumOfSigmasValue<-SumOfSigmasValue+value
  }
  SumOfSigmaAAs[i]<-SumOfSigmasValue
}


SumOfExpectedSigmaAAs<-c(1:15)
for (i in 1:15){
  ExpectedValue<-0
  for (j in 1:20){

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

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
HeaderSD<-c(-7:7)