stata - 从相关图存储 AC 结果
问题描述
我在存储自相关结果时遇到问题:
sysuse sp500.dta
tsset date
corrgram open
di `r(ac10)'
di `r(ac11)'
如您所见,该命令使用、等corrgram
打开一个表。我希望存储数据,但不知何故我无法访问所有信息。AC
PAC
Q
为什么我可以从中获取数据r(ac10)
而不能r(ac11)
?我的意思是,信息已经存在,有人可以启发我吗?
解决方案
该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 |
+-----------+