Commit 5647766c by Mustafa Tekpinar

The program won't write *_pred_evolEpi.txt and *_pred_evolInd.txt files anymore.

parent f5de1a26
...@@ -129,7 +129,7 @@ predInd = computePredNbSeqs(wt,nbSeqs) ...@@ -129,7 +129,7 @@ predInd = computePredNbSeqs(wt,nbSeqs)
rownames(predInd)=aa rownames(predInd)=aa
# output the sequence counts log-odd ratios # output the sequence counts log-odd ratios
write.table(predInd,paste0(prot,"_pred_evolInd.txt")) # write.table(predInd,paste0(prot,"_pred_evolInd.txt"))
print("done") print("done")
print("running global epistatic model...") print("running global epistatic model...")
...@@ -139,7 +139,7 @@ rownames(pred)=aa ...@@ -139,7 +139,7 @@ rownames(pred)=aa
# output the evolutionary distances between the query and the closest variants # output the evolutionary distances between the query and the closest variants
evolDist = pred/sum(trace^2) evolDist = pred/sum(trace^2)
evolDist[is.na(evolDist)] = 1 evolDist[is.na(evolDist)] = 1
write.table(evolDist,paste0(prot,"_pred_evolEpi.txt")) # write.table(evolDist,paste0(prot,"_pred_evolEpi.txt"))
print("done") print("done")
print("running normalization...") print("running normalization...")
......
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