Restricted Maximum Likelihood (REML) in ARMA model

Greetings,

I'm trying to fit a generalized linear model (glm) with a ARMA errors. R using the nlme package with the gls() function can fit REML approach. I think matlab usi ML estimate and I was wondering if it's possible to fit a glm with REML for a ARIMA error structure