Hello, I am trying to calculate the autocovariance matrix for any general farima(p,d,q) with p,q > 1. Could anyone give an idea how to implement in R or if there is any package for this? thank you beforehand. Jose.