stata中对一列数据怎么依次取前一行、前两行……的最大值？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

如题。
比如说
set obs 10
gen x=uniform()
然后生成一列数，第一行的数据是x中前一行的最大值，第二行的数据是前两行的最大值，依次类推。请问这个怎么做?


再如果我要这样处理1*10^7个数，怎么做会快点？


谢谢！！

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏1 回帖

关键词：Stata tata 最大值 Uniform form 最大值

自己顶！

看不明白你到底想做什么

蓝色 发表于 2015-6-27 19:05
看不明白你到底想做什么

比如说我有一列10个数：
2
3
4
2
4
6
8
9
2
4
如果我想生成另外一列数，第一行是上述这一列数的前一个的最大值，第二行是上述这一列数的前两个中的最大值，第三行是上述这一列数的前三个中的最大值……
这个要怎么做？

junstata 发表于 2015-6-27 19:49
比如说我有一列10个数：
2
3

1. clear
   
2. 
   
3. . input  a
   
4. 
   
5.              a
   
6.   1.
   
7. . 2
   
8.   2. 3
   
9.   3. 4
   
10.   4. 2
    
11.   5. 4
    
12.   6. 6
    
13.   7. 8
    
14.   8. 9
    
15.   9. 2
    
16. 10. 4
    
17. 11.
    
18. . end
    
19. 
    
20. .
    
21. . list
    
22. 
    
23.      +---+
    
24.      | a |
    
25.      |---|
    
26.   1. | 2 |
    
27.   2. | 3 |
    
28.   3. | 4 |
    
29.   4. | 2 |
    
30.   5. | 4 |
    
31.      |---|
    
32.   6. | 6 |
    
33.   7. | 8 |
    
34.   8. | 9 |
    
35.   9. | 2 |
    
36. 10. | 4 |
    
37.      +---+
    
38. 
    
39. . count
    
40.    10
    
41. 
    
42. . local n = r(N)
    
43. 
    
44. . gen b = .
    
45. (10 missing values generated)
    
46. 
    
47. . list
    
48. 
    
49.      +-------+
    
50.      | a   b |
    
51.      |-------|
    
52.   1. | 2   . |
    
53.   2. | 3   . |
    
54.   3. | 4   . |
    
55.   4. | 2   . |
    
56.   5. | 4   . |
    
57.      |-------|
    
58.   6. | 6   . |
    
59.   7. | 8   . |
    
60.   8. | 9   . |
    
61.   9. | 2   . |
    
62. 10. | 4   . |
    
63.      +-------+
    
64. 
    
65. . forval i = 1/`n'{
    
66.   2.  egen c = max(a) if _n <= `i'
    
67.   3.  replace b = c in `i'
    
68.   4.  drop c
    
69.   5.  }
    
70. (9 missing values generated)
    
71. (1 real change made)
    
72. (8 missing values generated)
    
73. (1 real change made)
    
74. (7 missing values generated)
    
75. (1 real change made)
    
76. (6 missing values generated)
    
77. (1 real change made)
    
78. (5 missing values generated)
    
79. (1 real change made)
    
80. (4 missing values generated)
    
81. (1 real change made)
    
82. (3 missing values generated)
    
83. (1 real change made)
    
84. (2 missing values generated)
    
85. (1 real change made)
    
86. (1 missing value generated)
    
87. (1 real change made)
    
88. (1 real change made)
    
89. 
    
90. . list
    
91. 
    
92.      +-------+
    
93.      | a   b |
    
94.      |-------|
    
95.   1. | 2   2 |
    
96.   2. | 3   3 |
    
97.   3. | 4   4 |
    
98.   4. | 2   4 |
    
99.   5. | 4   4 |
    
100.      |-------|
     
101.   6. | 6   6 |
     
102.   7. | 8   8 |
     
103.   8. | 9   9 |
     
104.   9. | 2   9 |
     
105. 10. | 4   9 |
     
106.      +-------+

复制代码

andruw 发表于 2015-6-27 20:17

1. clear
   
2. input  a
   
3. 
   
4. 2
   
5. 3
   
6. 4
   
7. 2
   
8. 4
   
9. 6
   
10. 8
    
11. 9
    
12. 2
    
13. 4
    
14. 
    
15. end
    
16. 
    
17. list
    
18. count
    
19. local n = r(N)
    
20. gen b = .
    
21. list
    
22. forval i = 1/`n'{
    
23. egen c = max(a) if _n <= `i'
    
24. replace b = c in `i'
    
25. drop c
    
26. }
    
27. list

复制代码

andruw 发表于 2015-6-27 20:18

非常感谢！又知道了新的东西！

andruw 发表于 2015-6-27 20:18

你好！如果有10000000个数据，按照这个程序走实在太慢。。请问有没有更快的方法能处理？谢谢！

junstata 发表于 2015-7-9 10:05
你好！如果有10000000个数据，按照这个程序走实在太慢。。请问有没有更快的方法能处理？谢谢！

有人问类似的问题：
https://bbs.pinggu.org/thread-3790020-1-1.html

1. clear
   
2. set obs 10000000
   
3. gen x=uniform()
   
4. gen xm = x[1]
   
5. replace xm = cond(x > xm[_n-1], x, xm[_n-1]) if _n >= 2

复制代码

andruw 发表于 2015-7-9 16:52
有人问类似的问题：
https://bbs.pinggu.org/thread-3790020-1-1.html

好的谢谢！发现是我的基本功太差。。