Skip to content
Prev 2215 / 15274 Next

fitted fGarch model

I need to modify my earlier response.

library(FinTS)
data(sp500)
spFit00.11a <- garchFit(sp500~garch(1,1), data=as.numeric(sp500))
at <- spFit00.11a at residuals
quantile(at)
           0%          25%          50%          75%         100%
-0.306849724 -0.028099724  0.001250276  0.029775276  0.414750276

	  From this, it is clear that garchFit()@residuals are NOT standardized.

	
st <- spFit00.11a at sigma.t
qqnorm(at, datax=TRUE)
qqnorm(at/st, datax=TRUE)

	  From comparing these two normal probability plots, it is apparent 
that 'at/st' is more normal than 'at'.

	  However, sd(at/st) = 0.9979, so the standardization is not perfect.

all.equal(as.numeric(sp500), spFit00.11a at fitted+at)	
# TRUE

	  Your pseudocode suggests that 'vt' are the standardized residuals, 
but that's not the case in this example.  If you have a real example, 
that contradicts this case, could you please report it?  Otherwise, I 
will assume that the test you reported was different from what you 
actually did.  It's always best to test an answer before giving it -- 
something I failed to do in my initial response to this question.

	  Best Wishes,
	  Spencer
babel at centrum.sk wrote: