Likelihood evaluations for stationary Gaussian time series are typically obtained via the DurbinLevinson algorithm, which scales as O(n^2) in the number of time series observations. This package provides a "superfast" O(n log^2 n) algorithm written in C++, crossing over with DurbinLevinson around n = 300. Efficient implementations of the score and Hessian functions are also provided, leading to superfast versions of inference algorithms such as NewtonRaphson and Hamiltonian Monte Carlo. The C++ code provides a Toeplitz matrix class packaged as a headeronly library, to simplify lowlevel usage in other packages and outside of R.
Package details 


Author  Yun Ling [aut], Martin Lysy [aut, cre] 
Maintainer  Martin Lysy <mlysy@uwaterloo.ca> 
License  GPL3 
Version  2.0.3 
Package repository  View on CRAN 
Installation 
Install the latest version of this package by entering the following in R:

Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.