# NOT RUN {
if (require(incidence) && require(ape)) {
i <- incidence(rabies_car_2003$linelist$date, 28L)
plot(i)
tre <- nj(dist.dna(rabies_car_2003$dna))
plot(tre, main = "Neighbour-Joining tree")
}
# }
Run the code above in your browser using DataLab