example(NJ)
fit = pml(tree, Laurasiatherian)
anc.ml = ancestral.pml(fit, type = "ml")
anc.p = ancestral.pars(tree, Laurasiatherian)
require(seqLogo)
seqLogo(t(anc.ml[1:10,]), ic.scale=FALSE)
seqLogo(t(anc.p[1:10,]), ic.scale=FALSE)
Run the code above in your browser using DataLab