Bonjour!

J'aimerais savoir s'il est possible sous R d'ajuster des données à un modèle linéaire généralisé dans lequel les résidus suivent un processus SARIMA(1, 0, 0)(1, 0, 0)_12.

En fait, mes données sont découpés en plusieurs sous groupes et pour chaque sous groupe, les résidus du modèle suivent le processus SARIMA que je viens de donner. Avec la libraire nlme, il y a moyen de faire ce genre de truc quand les résidus suivent un processus ARMA mais pas SARIMA
(avec la syntaxe : gls(Y~facteur1 + facteur2 + facteur3, correlation=corARMA(form=~temps|facteur1)) par exemple)

Quelqu'un sait comment je pourrais procéder?