standard errors in tsayp591

[wardb]

I used the tsayp591.prg to reproduce the results for Markov switching GARCH example from Tsay's Analysis of Financial Time Series, Example 12.6, but cannot obtain the estimated standard errors of the estimated parameters for the two regimes and their differences. I'd be most grateful if anyone could assist me with code. Thanks in advance.

[TomDoan]

Just add a STATISTICS instruction inside that final loop. That will compute the mean and variance of the MCMC draws for each of the parameters.

Code: Select all

do i=1,10
   stats stats(i) 1 nkeep
   density(type=histogram,counts,maxgrid=12) stats(i) 1 nkeep xx fx
   scatter(style=bar,hlabel=vlabels(i))
   # xx fx
end do i

[wardb]

Thanks for the extra code, which was most helpful.