stata 编程问题求助

各位大侠，您好，我在stata应用过程中遇到了两个问题，现特向各位stata高手请教求助. dataex

----------------------- copy starting from the next line -----------------------

1. * Example generated by -dataex-. For more info, type help dataex
   
2. clear
   
3. input str1 keys float group
   
4. "a" 1
   
5. "b" 1
   
6. "b" 1
   
7. "d" 1
   
8. "g" 1
   
9. "g" 1
   
10. "h" 1
    
11. "m" 1
    
12. "a" 2
    
13. "b" 2
    
14. "b" 2
    
15. "b" 2
    
16. "c" 2
    
17. "f" 2
    
18. "g" 2
    
19. "h" 2
    
20. "m" 2
    
21. "a" 3
    
22. "b" 3
    
23. "b" 3
    
24. "b" 3
    
25. "c" 3
    
26. "c" 3
    
27. "c" 3
    
28. "d" 3
    
29. "d" 3
    
30. "m" 3
    
31. "a" 4
    
32. "b" 4
    
33. "b" 4
    
34. "b" 4
    
35. "b" 4
    
36. "b" 4
    
37. "c" 4
    
38. "d" 4
    
39. "d" 4
    
40. "m" 4
    
41. end

复制代码

------------------ copy up to and including the previous line ------------------



问题1： 提取变量keys中的某特定值的序号并存储于暂元中，如m,在以上数据的序号分别是8，17，27，37，当然实际他们的分布没有特定的规律，每组有一个，但位置不固定，想将序号提取出来，然后应用于循环中。
set obs _n
levelsof obs if keys == "m", gen(num)
ret list
这种方法我用过了，想看看还有没有别的好方法

问题2：变量keys中的某特定值如b在各组中连续出现（注：一定要连续出现，如不连续出现则不计）的次数最大值，如b在第1、2、3、4组中分别出现了2，3，3，5次，最大值为5，需要提取这个最大值和存储于暂元中。
我用tabstat能求出连续值的个数列表，没能继续求结果，请各位大侠出手 帮助

以上两问题的各种解答方法（第一问我已经使用并呈现的方法除外）按方法计都奖励30积分，不同方法累积积分，如果相同方法则第一个提供者得到奖励。（注：用不同的命令属于不同的方法，用矩阵或者mata也是不同的方法，用不同的程序思路解决问题也属于不同的方法……但是如第一问中自己提供的方法，如有另外提供如下程序:set a =_n,levelsof a if keys == "m", gen(b),则属于同一种方法。请各位大侠积极出手相助。
谢谢大家，再次感谢！！！

第一个问题你已经解决了我就不回答了。回答第二个问题，为了满足你的要求，我在group=4时加了一个不连续的b,让你看得更清楚些。

1. * Example generated by -dataex-. For more info, type help dataex
   
2. clear
   
3. input str1 keys float group
   
4. "a" 1
   
5. "b" 1
   
6. "b" 1
   
7. "d" 1
   
8. "g" 1
   
9. "g" 1
   
10. "h" 1
    
11. "m" 1
    
12. "a" 2
    
13. "b" 2
    
14. "b" 2
    
15. "b" 2
    
16. "c" 2
    
17. "f" 2
    
18. "g" 2
    
19. "h" 2
    
20. "m" 2
    
21. "a" 3
    
22. "b" 3
    
23. "b" 3
    
24. "b" 3
    
25. "c" 3
    
26. "c" 3
    
27. "c" 3
    
28. "d" 3
    
29. "d" 3
    
30. "m" 3
    
31. "a" 4
    
32. "b" 4
    
33. "b" 4
    
34. "b" 4
    
35. "b" 4
    
36. "b" 4
    
37. "c" 4
    
38. "b" 4
    
39. "d" 4
    
40. "m" 4
    
41. end
    
42. 
    
43. 
    
44. spell keys,by(group)
    
45. gen d=_seq*_end if keys=="b"
    
46. egen wanted=max(d)
    
47. levelsof wanted
    
48. 
    
49. 5

复制代码

第二个也可以
bys group:egen count=count(keys) if keys=="b"
egen max=max(count)
local max=max

第二个只能针对你现有数据有效 还是楼上的正解
给个第一个的其他解法
gen id=_n
gen id2=id if keys=="m"
levelsof id2
dis r(levels)

pengxhan 发表于 2021-8-13 16:25
第二个也可以
bys group:egen count=count(keys) if keys=="b"
egen max=max(count)

这个代码不行哦，他要的是连续的b的个数最大值，举个例子：

1. clear
   
2. input str1 keys float group
   
3. "a" 1
   
4. "b" 1
   
5. "b" 1
   
6. "d" 1
   
7. "g" 1
   
8. "g" 1
   
9. "h" 1
   
10. "m" 1
    
11. "a" 2
    
12. "b" 2
    
13. "b" 2
    
14. "b" 2
    
15. "c" 2
    
16. "f" 2
    
17. "g" 2
    
18. "h" 2
    
19. "m" 2
    
20. "a" 3
    
21. "b" 3
    
22. "b" 3
    
23. "b" 3
    
24. "c" 3
    
25. "c" 3
    
26. "c" 3
    
27. "d" 3
    
28. "d" 3
    
29. "m" 3
    
30. "a" 4
    
31. "b" 4
    
32. "b" 4
    
33. "f" 4
    
34. "b" 4
    
35. "b" 4
    
36. "c" 4
    
37. "b" 4
    
38. "d" 4
    
39. "m" 4
    
40. end
    
41. 
    
42. bys group:egen count=count(keys) if keys=="b"
    
43. egen max=max(count)
    
44. local max=max
    
45. 
    
46. spell keys,by(group)
    
47. gen d=_seq*_end if keys=="b"
    
48. egen wanted=max(d)
    
49. levelsof wanted
    
50. 
    
51. drop _* d
    
52. list
    
53. 
    
54.      +-------------------------------------+
    
55.      | keys   group   count   max   wanted |
    
56.      |-------------------------------------|
    
57.   1. |    a       1       .     5        3 |
    
58.   2. |    b       1       2     5        3 |
    
59.   3. |    b       1       2     5        3 |
    
60.   4. |    d       1       .     5        3 |
    
61.   5. |    g       1       .     5        3 |
    
62.      |-------------------------------------|
    
63.   6. |    g       1       .     5        3 |
    
64.   7. |    h       1       .     5        3 |
    
65.   8. |    m       1       .     5        3 |
    
66.   9. |    a       2       .     5        3 |
    
67. 10. |    b       2       3     5        3 |
    
68.      |-------------------------------------|
    
69. 11. |    b       2       3     5        3 |
    
70. 12. |    b       2       3     5        3 |
    
71. 13. |    c       2       .     5        3 |
    
72. 14. |    f       2       .     5        3 |
    
73. 15. |    g       2       .     5        3 |
    
74.      |-------------------------------------|
    
75. 16. |    h       2       .     5        3 |
    
76. 17. |    m       2       .     5        3 |
    
77. 18. |    a       3       .     5        3 |
    
78. 19. |    b       3       3     5        3 |
    
79. 20. |    b       3       3     5        3 |
    
80.      |-------------------------------------|
    
81. 21. |    b       3       3     5        3 |
    
82. 22. |    c       3       .     5        3 |
    
83. 23. |    c       3       .     5        3 |
    
84. 24. |    c       3       .     5        3 |
    
85. 25. |    d       3       .     5        3 |
    
86.      |-------------------------------------|
    
87. 26. |    d       3       .     5        3 |
    
88. 27. |    m       3       .     5        3 |
    
89. 28. |    a       4       .     5        3 |
    
90. 29. |    b       4       5     5        3 |
    
91. 30. |    b       4       5     5        3 |
    
92.      |-------------------------------------|
    
93. 31. |    f       4       .     5        3 |
    
94. 32. |    b       4       5     5        3 |
    
95. 33. |    b       4       5     5        3 |
    
96. 34. |    c       4       .     5        3 |
    
97. 35. |    b       4       5     5        3 |
    
98.      |-------------------------------------|
    
99. 36. |    d       4       .     5        3 |
    
100. 37. |    m       4       .     5        3 |
     
101.      +-------------------------------------+

复制代码

按照我的代码得到的是3，按你的代码得到的是5.他要求b一定要连续出现哦，如不连续则不计。

谢谢楼主的回复：很好的思路，您提供的两种方法，
第一种方法：我改进如下
spell keys,by(group)
gen d=_seq*_end if keys=="b"
sum d
local max `r(max) '

第二种方法:有些小瑕疵，因为它没有识别是否为连续出现的b,连续的最大计数为5，而不是6，该方法把所组内所有的结果都计篱

楼上正解，第一种方法可行，第二种方法不是我想要的方法，但是我错给第二种方法点赞了

fugangxx 发表于 2021-8-13 11:59
各位大侠，您好，我在stata应用过程中遇到了两个问题，现特向各位stata高手请教求助. dataex

---------- ...

我说了呀，只针对他目前已经连续数据有效。

pengxhan 发表于 2021-8-13 17:28
我说了呀，只针对他目前已经连续数据有效。

嗯，我在你发那个帖子之前发出去的，在审核中，等你解释的时候，已经撤不回来了

纯粹为了练习 我编写了下面命令可识别不连续 不用spell
bys group:egen seq=seq() if keys=="b"
egen  max=max(seq) if  seq[_n-1]!=.
levelsof  max