<div dir="ltr"><div>Hi Alexios</div><div><br></div><div>I am modelling the same data and getting vastly different estimates using "rugarch" package and "fGarch" package (all installs and packages have been recently downloaded and should be up to date). The data is 242 monthly log returns in raw log return format (i.e. not multiplied by 100 to get percent format). I have attached the file for reproducible results.</div><div><br></div><div>In estimating an ARMA(2,2)-GARCH(1,1) with "normal" errors, I get the following output using:</div><div><br></div><div>(1) "rugarch" package:</div><div>spec <- ugarchspec(variance.model=list(model="sGARCH",garchOrder=c(1,1),</div><div>          submodel=NULL,external.regressors=NULL,variance.targeting=F), </div><div>          mean.model=list(armaOrder=c(2,2),include.mean=T,external.regressors=NULL),</div><div>          distribution.model="norm")</div><div>tempgarch <- ugarchfit(spec=spec,data=ALSI.reg.log.ret,solver="hybrid")</div><div>show(tempgarch)</div><div><br></div><div>Output:</div><div>Optimal Parameters<br></div><div><div>------------------------------------</div><div>              Estimate     Std. Error     t value       Pr(>|t|)</div><div>mu         0.015107    0.003310     4.56480      0.000005</div><div>ar1         0.561145    0.577557     0.97158      0.331258</div><div>ar2        -0.303709    0.526705    -0.57662      0.564195</div><div>ma1      -0.629355    0.555266    -1.13343      0.257033</div><div>ma2       0.430343    0.511488     0.84135      0.400149</div><div>omega   0.000257    0.000343     0.75094      0.452687</div><div>alpha1   0.197622    0.148026     1.33505      0.181862</div><div>beta1     0.730256    0.228730     3.19265      0.001410</div></div><div><br></div><div>(2) "fGarch" package</div><div><div>garch.fit=garchFit(formula=~arma(2,2)+garch(1,1),data=ALSI.reg.log.ret,cond.dist="norm",trace=F,include.mean=T)</div><div>summary(garch.fit)</div></div><div><br></div><div>Output:</div><div>Error Analysis:<br></div><div><div>              Estimate      Std. Error      t value    Pr(>|t|)    </div><div>mu          0.0368275   0.0070622    5.215      0.000000184 ***</div><div>ar1        -0.4909910   0.0955405   -5.139      0.000000276 ***</div><div>ar2        -0.8424945   0.0538062  -15.658     < 2e-16 ***</div><div>ma1       0.4733331   0.0904955    5.230      0.000000169 ***</div><div>ma2       0.8844500   0.0550194   16.075     < 2e-16 ***</div><div>omega   0.0003106   0.0001838    1.690      0.0910 .  </div><div>alpha1   0.3696620   0.1496237    2.471      0.0135 *  </div><div>beta1     0.5808720   0.1434332    4.050      0.000051267 ***</div></div><div><br></div><div>I am confused as to why the output differs to such a great extent. Any idea why this is happening?</div><div><br></div><div>Thank you very much</div><div>Gareth</div><div><br></div></div>