# NOT RUN {
# create a sequence with some processes attached
s<-Sequence(
string="ATGC",
alphabets=list(NucleotideAlphabet()),
processes=list(list(JC69()),list(JC69(),GTR()))
)
# get total rates for positions 1 and 3
getTotalRatesFromRange(s,c(1,3))
# get all total rates via virtual field
s$totalRates # via the "getTotalRates.Sequence" method
# }
Run the code above in your browser using DataLab