# NOT RUN {
# read DNA from fasta file
fasta <- fasta.read("virus.fasta", "host.fasta")
fasta.v <- fasta[[1]]
fasta.h <- fasta[[2]]
# }
# NOT RUN {
# Calculate zscore using (codon model)
codon <- dinuc.codon(fasta.v, permutations = 10)
# }
Run the code above in your browser using DataLab