##simulate series with 500 observations:
set.seed(123)
y <- tegarch.sim(500, omega=0.01, phi1=0.9, kappa1=0.1, kappastar=0.05, df=10, skew=0.8)
##estimate a 1st. order Beta-skew-t-EGARCH model:
mymod <- tegarch.est(y)
##compute fitted values:
fitvals <- tegarch.fit(y, mymod)
##make plot of the conditional standard deviaiton:
plot(fitvals[,"stdev"])
##make histogram of the standardised residuals:
hist(fitvals[,"residstd"])
Run the code above in your browser using DataLab