急求！如何使用循环语句计算不同观察值同一变量任意相减？

求教大神！我有2008和2009年两年患者的门诊就诊数据，想要生成一个变量timegap，等于2008年-2009年期间的任意2次门诊就诊时间之差，如果30≤timegap≤365，即认为患者为我的目标人群。
现在的问题是，不同患者2008年-2009年就诊次数不一样，而且有的次数多达上百次，一条一条手写不现实，请问如何用循环的语句实现我的目标呢？

举例：
序号  患者编号         就诊时间          就诊年份
1             1              2008.10.01      2008
2             1              2008.12.15      2008
3             1              2009.01.10      2009
4             1              2009.05.01      2009
5             2              2008.08.15      2008
6             2              2008.10.01      2008
7             2              2008.12.15      2008
8             2              2008.12.31      2008
9             2              2009.01.31      2009
10           2              2009.05.31      2009
11           2              2009.08.31      2009

对患者编号为1的人，我想实现的是序号1-4这4条记录中，任意两条记录的就诊时间相减，差值满足30≤timegap≤365即可
对患者编号为2的人，我想实现的是序号5-11这7条记录中，任意两条记录时间的差值满足30≤timegap≤365即可

急求循环语句~！谢谢各位大神

See below.

顶一个~

尔后建议请用 dataex (先 ssc install dataex 并见说明) 将原始 Stata 资料中具有”代表性”的一部分资料列出，以供有意回答者实验之用，并能提供具体操作指令。并请参考 http://www.jianshu.com/p/9870080fe769,  https://bbs.pinggu.org/thread-5048204-1-1.html, 与 https://bbs.pinggu.org/thread-5917273-1-1.html。

逍遥梦蝶 发表于 2018-2-27 07:01

大神，从34行开始能不能简单说一下思路~我看不太懂命令

tracychen1900 发表于 2018-2-27 13:07
大神，从34行开始能不能简单说一下思路~我看不太懂命令

稍微有一点错误。应该在while循环下再加一层循环，才能计算出所有可能组合的差值。这一步不难，你可以自己尝试一下。

1. *- Input data
   
2. cls
   
3. clear all
   
4. input patient str15 date year
   
5. 1 "2008.10.01" 2008
   
6. 1 "2008.12.15" 2008
   
7. 1 "2009.01.10" 2009
   
8. 1 "2009.05.01" 2009
   
9. 2 "2008.08.15" 2008
   
10. 2 "2008.10.01" 2008
    
11. 2 "2008.12.15" 2008
    
12. 2 "2008.12.31" 2008
    
13. 2 "2009.01.31" 2009
    
14. 2 "2009.05.31" 2009
    
15. 2 "2009.08.31" 2009
    
16. end
    
17. 
    
18. *--- 思路 ---*
    
19. *
    
20. * - 对某一病人任意两次就诊日期做差，依据本人经验，将日期置于一列之内
    
21. *   不方便计算。并且，由于有 anymatch() 拓展函数的需要，应该可以将这
    
22. *   一长数据转换为宽数据，进行变量间的做差运算。
    
23. *
    
24. * - 由于不同病人就诊的次数不相同，预期：转化为宽数据后，变宽的变量
    
25. *  （本例中是 date 变量）的个数由就诊次数最多的病人决定。就诊次数相对
    
26. *   少的病人的就诊日期会有被标记为缺失值。
    
27. *
    
28. * - 由于日期间隔均是一个较后的日期减去较前的日期，取所有可能的差值组合
    
29. *   的方法是：固定住一个较前的日期，用所有较后的日期去减它。比如：固定
    
30. *   住 date1，而后有
    
31. *      date2 - date1
    
32. *      date3 - date1
    
33. *      ...
    
34. *      date7 - date1
    
35. *   类似地，固定住 date2，而后有
    
36. *      date3 - date2
    
37. *      date4 - date2
    
38. *      ...
    
39. *      date7 - date2
    
40. *   这样的组合理论上应该 = 单一病人就诊次数最大值 * (单一病人就诊次数最大值 - 1)/2
    
41. *   需要采取合适的循环结构来解决。
    
42. 
    
43. *- Date
    
44. gen dt = date(date, "YMD")
    
45. format dt %tdCCYY-Mon-DD
    
46. drop date year
    
47. rename dt date
    
48. order date, after(patient)
    
49. 
    
50. *- Reshape data
    
51. sort patient date
    
52. by patient: gen rank = _n
    
53. order rank, after(patient)
    
54. quiet sum rank
    
55. local m = `r(max)'  // 提取出单一病人就诊次数的最大值
    
56. 
    
57. reshape wide date, i(patient) j(rank)
    
58. 
    
59. *- Calculate pairwise differences
    
60. local i = 1
    
61. local j = 1
    
62. local k = 1
    
63. while `i' < `m' {                         // 对减数循环，只要 i 小于单一病人就诊次数最大值就继续循环
    
64.    forval j = `=(`i'+1)'/`m' {            // 对被减数循环，被减总是更靠后的日期
    
65.        gen diff_`k++' = date`j' - date`i' // 生成日期间隔天数
    
66.    }
    
67.    local i = `i' + 1
    
68. }
    
69. drop date*
    
70. 
    
71. *- Create indicator
    
72. foreach v of varlist diff_* {
    
73.    replace `v' = inrange(`v', 30, 365) if !missing(`v')
    
74. }
    
75. *
    
76. egen target = anymatch(diff_*), value(1)
    
77. label var target "Enter the sample or not"

复制代码

逍遥梦蝶 发表于 2018-2-27 17:43
稍微有一点错误。应该在while循环下再加一层循环，才能计算出所有可能组合的差值。这一步不难，你可以自己 ...

还请大神指导一下。。。

请先安装 (ssc install) rangestat，然后：

1. * Example generated by -dataex-. To install: ssc install dataex
   
2. clear
   
3. input float id str12 ymd float n
   
4. 1 "2008.10.01" 1
   
5. 1 "2008.10.15" 2
   
6. 1 "2008.10.31" 3
   
7. 1 "2008.11.15" 4
   
8. 2 "2008.08.15" 1
   
9. 2 "2008.10.01" 2
   
10. 2 "2008.12.15" 3
    
11. 2 "2008.12.31" 4
    
12. 2 "2009.01.31" 5
    
13. 2 "2009.05.31" 6
    
14. 2 "2009.08.31" 7
    
15. 3 "2009.05.30" 1
    
16. 3 "2009.08.31" 2
    
17. 3 "XXXX.XX.XX" 3
    
18. 3 "2011.09.30" 4
    
19. end
    
20. gen date = daily(ymd, "YMD")
    
21. format date %td
    
22. 
    
23. // Robert Picard
    
24. rangestat (count)  over=id, by(id) interval(date 30 365)
    
25. rangestat (count) under=id, by(id) interval(date -365 -30)
    
26. 
    
27. gen inwindow = over < .  | under < .

复制代码

逍遥梦蝶 发表于 2018-2-27 22:18

谢谢大神！！！！！