Following [1]}, an alternative approach is Quasi Maximum Likelihood (QML): this is obtained assuming \(\varepsilon _{t}\)  as Gamma distributed and then estimating \(\mathbf {\theta }\)  by Maximum Likelihood.
