# NOT RUN {
# }
# NOT RUN {
# read DNA from fasta file
fasta <- fasta.read("virus.fasta", "host.fasta")
fasta.v <- fasta[[1]]
fasta.h <- fasta[[2]]
RSCU.H <- RSCU.values(fasta.h)
RSCU.V <- RSCU.values(fasta.v)
# }
# NOT RUN {
# Calculate SiD
SiD <- SiD.value(RSCU.host, RSCU.virus)
# }
Run the code above in your browser using DataLab