首页 > 解决方案 > 从相关图存储 AC 结果

问题描述

我在存储自相关结果时遇到问题:

sysuse sp500.dta
tsset date
corrgram open

di `r(ac10)'
di `r(ac11)'

如您所见,该命令使用、等corrgram打开一个表。我希望存储数据,但不知何故我无法访问所有信息。ACPACQ

为什么我可以从中获取数据r(ac10)而不能r(ac11)?我的意思是,信息已经存在,有人可以启发我吗?

标签: statastata-macrosautocorrelation

解决方案


corrgram命令返回某些结果,例如r(ac1)-r(ac10)作为标量:

. return list

scalars:
               r(lags) =  40
                r(q10) =  997.7976424554661
              r(pac10) =  .
               r(ac10) =  .4660537633851876
                 r(q9) =  941.2145380745385
               r(pac9) =  .
                r(ac9) =  .4926995142415946
                 r(q8) =  878.2410054272498
               r(pac8) =  .
                r(ac8) =  .6798652708950988
                 r(q7) =  758.8350022301894
               r(pac7) =  .
                r(ac7) =  .8497950586730207
                 r(q6) =  573.0532520717032
               r(pac6) =  .
                r(ac6) =  .7058394210650882
                 r(q5) =  445.4128558881126
               r(pac5) =  .
                r(ac5) =  .537862961940051
                 r(q4) =  371.600563372601
               r(pac4) =  -.123160137402293
                r(ac4) =  .5238127059274379
                 r(q3) =  301.8811184740619
               r(pac3) =  .0903652372655368
                r(ac3) =  .5378744937938151
                 r(q2) =  228.6682359982045
               r(pac2) =  -.0172170443544806
                r(ac2) =  .5633576757209979
                 r(q1) =  148.6802035217271
               r(pac1) =  .9912569913768637
                r(ac1) =  .769625068645877

但是,它也将所有内容作为矩阵返回:

matrices:
                  r(Q) :  1 x 40
                r(PAC) :  1 x 4
                 r(AC) :  1 x 40

您可以从相应返回的矩阵中访问您需要的内容,例如r(Q)r(PAC)r(AC)

具体来说ac11

. display el(r(AC), 1, 11)
.46352669

或者

. matrix A = r(AC)
. display A[1,11]
.46352669

要将返回的结果直接保存为具有名称的 Stata 变量ac

mata: st_store((1::40), st_addvar("double","ac"), colshape(st_matrix("r(AC)"), 1))

. list ac in 1 / 40

     +-----------+
     |        ac |
     |-----------|
  1. | .76962507 |
  2. | .56335768 |
  3. | .53787449 |
  4. | .52381271 |
  5. | .53786296 |
     |-----------|
  6. | .70583942 |
  7. | .84979506 |
  8. | .67986527 |
  9. | .49269951 |
 10. | .46605376 |
     |-----------|
 11. | .46352669 |
 12. | .47980565 |
 13. | .62691605 |
 14. | .76459033 |
 15. | .59986479 |
     |-----------|
 16. | .43064803 |
 17. | .41606848 |
 18. |  .4143769 |
 19. | .42979786 |
 20. | .55313291 |
     |-----------|
 21. |  .6518829 |
 22. | .50817027 |
 23. | .35552032 |
 24. | .34330162 |
 25. | .34298527 |
     |-----------|
 26. | .34803055 |
 27. | .44354181 |
 28. | .52455024 |
 29. | .40499719 |
 30. | .28776225 |
     |-----------|
 31. |  .2736146 |
 32. | .27286332 |
 33. | .27667982 |
 34. | .34373118 |
 35. | .40664714 |
     |-----------|
 36. | .30633385 |
 37. | .21286866 |
 38. | .19836486 |
 39. | .20569954 |
 40. | .20306142 |
     +-----------+

推荐阅读