From Ape's log, `chronogram`

has indeed been deprecated.

CHANGES IN APE VERSION 2.5-2

DEPRECATED & DEFUNCT

o evolve.phylo() and plot.ancestral() have been removed.

o chronogram(), ratogram(), and NPRS.criterion() have been removed.

But as mentionned in this page (thanks @cpad0112), using the parameter `lambda=0`

in `chronos`

is similar to applying the NPRS.

**[EDIT]** After looking at the function's code and confirming by reading Sanderson 2002, lambda is the level of penalization of the loglikelihood. So the higher the lambda, the closer you are from NPRS. I am pretty sure the page I quoted here is wrong. Here is the function of interest (penalized likelihood with correlated rate), in chronos:

```
correlated = function(rate,
node.time) {
loglik <- log.lik.poisson(rate, node.time)
if (!is.finite(loglik)) return(-1e+100)
loglik - lambda * (sum((rate[ind1] - rate[ind2])^2) +
var(rate[basal]))
}
```

This page may help: http://grokbase.com/t/r/r-sig-phylo/112as2y3m8/replacement-for-chronogram

1.7k