RT。
这是我在上的一个课要做的一个research project中的一步,大概是这样的:
已知所有parameter的值,用data step generate了一个AR(1)-GARCH(3,1)的model(即5个参数都知道,n=500的time series,burn-in period也是500),这样可以有一组Y,然后用procautoreg去算MLE的值,这些都在SAS output里可以找到,并且也会output出来MLE的standarderror。但是因为理论上MLE的asymptoticdistribution是normal的,也就是说如果我generate出来100000组Y,每一组都用autoreg算一个MLE,然后把这些MLE plot出来,会得到一个mean为我已知的那些parameter的normal distribution。
现在的问题是,如何让SAS算出来这个normal distribution的standard error是什么,这个是我现在卡住了的一个地方。我这个project的目的就是要看MLE的分布是不是和理论上的那个normal很一致。
我查了一些SAS的帮助,但是没有什么进展,所以到论坛上求助一下~~多谢各位的帮助!!!!!!
PS:要算的这些standard error(也是5个,每个parameter都会对应有一个)完全是通过已知的parameter的值和samplesize(这里是500)算出来的。SAS可以output出来sample covariancematrix,但是这个是和每一组Y有关系的,并不是我想要的那个结果。
再谢!!!!!
如果有知道如何解决这个问题的高手,可以在下面留言,我可以把我的code发给你看下。。