Hello, Am trying to estimate a seasonal Arima by calling the R forecast package in Rpy2: fit = forecast.Arima(x = h02, order = order, seasonal = seasonal) Strangely I get the estimates of the non-seasonal part, only, even if the model is specified as (3,0,1)x(0,1,2). Can anyone tell me where the issue is? Regards, P. [[alternative HTML version deleted]]