Commit 7f40b9ac by Mustafa Tekpinar

Commented the lines containing pssm output.

parent 8e2a131c
...@@ -38,9 +38,9 @@ aliVeryCons = ali[pId>0.8,] ...@@ -38,9 +38,9 @@ aliVeryCons = ali[pId>0.8,]
N = c(dim(ali)[[1]],dim(aliCons)[[1]],dim(aliVeryCons)[[1]]) N = c(dim(ali)[[1]],dim(aliCons)[[1]],dim(aliVeryCons)[[1]])
#resAliCons = computePSSM(aliCons) #resAliCons = computePSSM(aliCons)
res = list(computePSSM(ali,N[1],npos),computePSSM(aliCons,N[2],npos),computePSSM(aliVeryCons,N[3],npos)) res = list(computePSSM(ali,N[1],npos),computePSSM(aliCons,N[2],npos),computePSSM(aliVeryCons,N[3],npos))
write.table(res[[1]][[3]],paste0(prot,"_pssm.txt")) # write.table(res[[1]][[3]],paste0(prot,"_pssm.txt"))
write.table(res[[2]][[3]],paste0(prot,"_pssm60.txt")) # write.table(res[[2]][[3]],paste0(prot,"_pssm60.txt"))
write.table(res[[3]][[3]],paste0(prot,"_pssm80.txt")) # write.table(res[[3]][[3]],paste0(prot,"_pssm80.txt"))
# read evolutionary traces computed by JET # read evolutionary traces computed by JET
jet=read.table(paste(prot,"_jet.res",sep=""),head=TRUE) jet=read.table(paste(prot,"_jet.res",sep=""),head=TRUE)
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment