一个用stata 计算股票成本的小程序

awer

# k9 ~: T2 [# _. a+ D/ S# [) A
2 @' N( c) I' L5 B2 l

   我是从西西河一直追虎大追到爱坛来的，虎大曾用stata 分析过股票的长期成本和一年的短期成本，但是直接用kdensity命令无法得出筹码峰的具体数值，因此我写了一个小程序来计算股票的筹码峰值、长期平均成本和短期平均成本。下面就给出具体实现方法。

计算之前需要有 ：

1） 招商证券系统（通达信其它客户端也行，功能都差不多）

2） Stata软件 （我用的版本是Stata/SE 11.0）

然后

1．新建一个文件夹储存股票数据，比如:

e:\stock\source

2. 股票数据导出

a 进入招商证券, 系统-盘后数据下载，更新数据

b快捷键“34”进行数据导出, 在弹出窗口选择
' t- E# p/ r* c* B: c( N“高级导出”,

在高级导出窗口选择 “日线导出”,

将导出目录选择为你新建的储存股票数据文件夹的目录（e:\stock\source）! 切记切记

文件名 为 XXXXXX.TXT

勾选  数据精确复权

添加品种中选择想要导出的品种

开始导出

3下载 stata, 将本贴附件中的hello. do文件copy至stata根目录下

用Stata中的 Do file editor 打开hello.do 文件， 确保hello.do 文件中 e:\stock\source 和你保存股票数据的文件夹目录一致

现在就可以计算了，以露天煤业（002128）为例

在Stata中 输入如下命令

do hello

hello 002128

- U5 B) |5 b* P# w( f/ g(002128为你想查询的任意一只股票的代码)

几秒钟之后应该会出现下面一张图

( P( t  b" A' Q9 d+ k$ h  ^

hello.do (5.25 KB, 下载次数: 12295)

2011-8-13 18:02 上传

点击文件名下载附件

其中蓝线为股票上市交易以来的平均成本值（注意!此值没有去掉一级市场发行和国家的0成本）

红线为一年内交易的平均成本值

黑实线和虚线为一年平均成本值附近两个筹码峰的值

黑方框代表的事20内每天股票的收盘价格

红方框代表股票最新收盘价

山寨完毕， 收工!!

awer

这点我也注意到了，所以我用 成交额/复权后的价格  作为成交量 用于权重计算

awer

沉默的前哨 发表于 2011-10-19 20:49
9 s. ?$ H) k* m能不能帮忙看看：

do hello

你算得哪个股票？

awer

沉默的前哨 发表于 2011-10-20 10:32
do hello 后变成
+ d7 k; B8 ~  \5 b......

awer

沉默的前哨 发表于 2011-10-19 20:45
/ y1 |/ {5 g+ k1 n为什么老是提示：type mismatch

注意 日期格式要选成  yyyy/mm/dd，要不然时间变量格式不对

awer

琴心三叠 发表于 2011-10-20 21:37
谢谢答复  还有一个问题  我导出的文件为何只有当天的数据  而不是历史全部数据呢 ...

# ~6 z: P4 v8 j* G有很多问题是数据不全造成的
$ g: j" F2 O/ i) @- n8 c1 \如果想要全部数据，
( W, W; L% N) P4 C; t系统--盘后数据下载--日线数据--日线和实时行情数据

  R! m, e) {! L: v* I把开始时间设成中国股市开盘的那天，下载全部数据

下载完全部数据后，在屏幕上使用快捷键 34

34--高级导出--

3 R0 t) s5 b8 }2 d然后选中   日线导出；选定导出目录；文件名设为：XXXXXX.TXT;勾选精确复权；不要生成导出头部；分隔格式为tab（）；日期格式为YYYMMDD，然后选择添加品种

" `& c, Y" B; m) [! y4 G  }我一般都是把所有数据下载到指定目录下，然后进行hello运算
! G) s: g/ V+ Q3 U# m

awer

改了一下，这次运算应该没啥问题了

( [' \, Q% N! o+ w$ @# B
. k# b" T; V, I+ Dcapture program drop hello
" z* H4 a- Q' G# [
1 j. y: R' Q7 Lprogram hello
7 Y2 u. f# L( I, N, @" A1 t5 i
quietly {
& D- w+ H+ O3 x3 U: Q& E
* M8 a  j8 n& w+ B& q1 aclear
$ J+ L- _; l3 Q5 l. b; r# ^
% u1 P4 Q2 E& [; U% Z2 s* Ocd e:\finance\stock\20110101\source

insheet using `1'.txt
$ i" }5 x+ b2 @! Y+ Q
   keep if v5>0
   
# L+ ^- p! V5 W   gen vt=_n
( a9 b8 P# q7 q& r+ R+ ^1 ~   
  _5 D, A& T+ U2 H$ ?   tsset vt
9 `3 k+ ~. L$ @; Y( d# r2 @
. E1 F3 N6 N5 g+ d3 ~- i   gen zf=D1.v5

" |& o! ~0 R- \7 q1 P, W5 J, R; j   gsort -v1
) S  k: o$ \$ }- r! R( R# S0 y9 K
   gen avol=v7/v5
; g0 Z: u1 r" ]! F: e3 K" k+ {   
   keep if avol>0
   
7 D8 C, ^* J, t; N( k$ c+ r& d   drop vt
  
   gen vt=_n
   
" v2 h: z' ^, v+ O# i   tsset vt
8 l. T" }# A5 q" d   
    sum v5 [aw =avol]

7 H3 U/ ?1 u" E3 w  L   gen amean=r(mean) in 1
4 i( |. v- o* K7 F1 y: {   
   * 总加权平均值 amean
* a# u2 }; S) g0 k: W; u0 ^   
, K  D- k) D3 h* a" H' U   
$ n: [+ `8 x/ v7 U' I, s& P3 |) o   keep in 1/240

   sum v5 [aw =avol]

, `+ h9 |+ p8 A0 j1 n4 k2 A6 x   gen mean240=r(mean) in 1
' B& I5 g/ I% q* t  A; w   
   *240加权平均值 mean240
   
   set obs 500
   
   kdensity v5 [aw =avol],gaussian n(500) generate(dx240 dy240) nograph
/ ^7 f- l. T. U   
   *240 kdensity(dx240 dy240)
   
   sum dy240
   
4 G9 {8 g, _6 J* R4 S/ |   gen yxis=r(max) in 1
   
5 r, t; F8 ?. A3 g# `) q, j   gen yxis2=yxis/2
' _6 q- H$ X! `, j8 g# f3 K, _& H. V   
0 k+ B4 r& S# n8 s4 E. ?8 C   *yxis
   
% N' M$ t( m& @, k4 {% _4 i   gen dy240max=r(max)
' W6 F1 S% \4 n5 h: J' F& l/ w7 ^   
/ w2 @( [; p( F# y; v7 J5 U& d# o   gen dy240minp=r(min)

; R3 }3 E: I/ ?    gen vtt=_n
1 h8 Z' y* z8 s* a
    tsset vtt

   gen ddy240=D1.dy240/D1.dx240
   
   sum ddy240
   
   gen ddy240max=r(max)
   
! U& H( Z6 v/ `7 u- B   gen ddy240min=r(min)
1 g, w; ^7 M* l* f3 K     
   gen ddy240g=ddy240*(dy240max-dy240min)/(ddy240max-ddy240min)
   
   sum ddy240g
- B6 H5 M9 m6 d/ V% }! Z   
   gen yxis3=r(min) in 1
8 C' I1 v/ F6 l7 q; W- p   
   gen yxis4=0 in 1
% v! \+ k1 A% y  
: n( N' t, u) u  *ddy240g
   
   
7 Z: ^+ ^8 R0 i" a, g  preserve
  
. J* ^0 h- @! F; X/ O' s  sum v5 [aw =avol]
7 c$ S! ~, P# M" T  
  keep if dx240<r(mean)
  
  sum ddy240
) R$ |, u: s  d' X5 j: ?" e  
  if r(min)<0 {
3 r6 V) K' \, H! t" u/ T6 V        
" ^5 o& T/ _4 F4 S8 J                keep if ddy240<=0
               
                count
               
0 [6 P  w+ t+ s- D                if r(N)>2 {
& W7 B% Q, s/ w4 d" j                 
                 pkexamine dx240 dy240
               
                local tempminp=r(tomc)
( }; `" B1 H* x                restore
                 }
( ^) q- t7 C( |# A4 M7 j                 
' U* V" k  L% t1 f5 l        else {
5 D2 t5 i# {; R/ w" s3 h. ?, V               
                 di "cannot find minp"
                 
                restore
3 O% h" E; k3 o4 k! Z- L                 sum dx240
                 local tempminp=r(mean)
2 E+ L: X/ g9 S' z8 S7 A                 
                 }
                }
   else {
3 V1 [0 I: N/ Q3 N+ a- O* n3 \0 O                 
                 keep if ddy240<=r(min)
                 
                 sum dx240
                 
                 local tempminp=r(min)
                 restore
# a! d/ G) k, ~0 C' ^5 A2 @& T9 }                 }



   gen minp=`tempminp' in 1
   *240 minp
0 b% K2 X  X6 g- U/ Z

& R5 F' b5 R8 N4 D1 \   
   preserve
4 m* X- ?; S) v- ^6 R   
   sum v5 [aw =avol]
; @0 j* E7 o0 ]( W3 y$ m, [& v: K   
0 z2 z. S2 I& g( N* l$ B0 a   keep if dx240>r(mean)
4 w4 Q9 U$ N2 [, F4 _) Q7 Y   
   sum ddy240
   
   if r(max)>0 {
! B0 v: I. ?2 u9 h1 p               
3 r# U% d. s6 H% i9 f                keep if ddy240>=0
6 t; Q. T5 r7 |4 g8 g( r               
/ p9 Z- o& i) e  P0 t- v) b                count
               
8 h) T3 w1 ~2 ]8 f+ k                if r(N)>2 {
                 
% p% j, l  J" w: V                 pkexamine dx240 dy240
: N; M4 w8 |4 b3 [- [* H               
                local tempmaxp=r(tomc)
                restore
                 }
                 
                 
" A/ M9 ^. P! A  Z5 V  [        else {
                restore
; X+ B/ b4 y& r                 di "cannot find maxp"
               
1 A' ^, |  o, U  a0 E3 {2 @+ P7 a* K               
* c# t8 ?7 H' V  H8 D, C                 sum dx240
1 r1 w* G+ d+ P/ P/ [7 t                 local tempmaxp=r(mean)
                 
                 }
                 
                 }
                 
    else {
       
           keep if ddy240>=r(max)
          
           sum dx240
          
           local tempmaxp=r(max)
$ O8 [( w  }; \8 j/ N+ a" C           restore
) v: u! O: y/ O! h* U           }
       
' O7 @1 b' M8 B9 I) F
  

2 N7 i9 }4 O5 _7 z   gen maxp=`tempmaxp' in 1
1 ~' Q' S# b2 d& }! q+ g; e& @   
- U) ]+ ^3 @- ~. T' Z: U' |   *240maxp
   
   
" E* _+ A# C, T" ]. J6 R6 a$ y6 }   *-------------ddy240pminp左侧高峰值-----
$ N- ~0 g: z, E   
/ L! i! E+ e5 {  U4 g    preserve
   
; f; o# @0 A$ f   keep if dx240<`tempminp'
# |' e) \! O- a: z6 M' R% ~: S   
" p) w, W5 J7 p( l+ R% l   keep if ddy240>0
4 ~+ ]9 b: G7 U/ `0 Y/ Z- S' f   
   count
# v* u9 j' D( \# K               
$ u* j, D$ R; j% E4 x8 |# T                if r(N)>2 {
9 b' l& S# B0 h3 {! X( J9 I                 
+ X$ O  {/ R; P2 M  w                pkexamine dx240 ddy240
8 r6 G. i* c% p5 L" q" k/ |               
0 \- u. q' S! B+ b; ]                local tempminphddy240=r(tomc)
                restore
                 }
/ l; ^+ _  ?. R; r) e/ m% ]                 
; l; {1 L3 ~0 E5 x, @9 I7 l9 P        else {
/ b4 u4 r/ G0 F                restore
                 di "cannot find pminphddy240"
! I6 j+ J! R0 ?# x; Q5 o- Z- ?: O                 
               
                 sum dx240
2 i& d4 D; _$ D/ `# t5 m8 |                 local tempminphddy240=r(mean)
               
                 
- D8 j; d; f7 ^( a$ ^2 ]1 A" c, h                 }
4 T) ?9 {5 I; o                 
* ^7 U' Q  [9 }& t  g1 V   
( _9 u4 K, [2 f4 c) y! @4 J+ a& W   
6 Y& N; u$ i; b   
/ x+ y& w8 f: _1 p7 `; g3 G! {; \6 I, L   
   
% A" Z/ d- w$ {  x4 ]. J1 r        gen pminphddy240=`tempminphddy240' in 1
       
4 R8 d! ^( U3 z1 ?; G1 F       
        *-------------ddy240pminp右侧低峰值 -----------------
    preserve

   keep if dx240>`tempminp'
  v: K% W! g) q3 T6 a& U   
* r  I3 S& l, F8 u   keep if dx240<`tempmaxp'
4 ?& q( M; q6 H( }1 S4 c   
   keep if ddy240<0
9 E7 e/ d' R2 j   
   gen temddy240=-ddy240
   
; x! U. i3 @, a0 B4 I1 J4 a    count
               
                if r(N)>2 {
: l. l8 y( I1 w: V) \/ l                 
                 pkexamine dx240 temddy240
8 e0 L& o% T0 k. i2 \6 W   
' i/ C" ^; l2 M2 G( S9 \  local tempminplddy240=r(tomc)
  L; i1 U1 k. I0 H9 q& }: Q                restore
1 a( @0 T# W" Q% \8 `                 }
+ e4 w5 S+ a" S, {, x! C                 
        else {
               
8 L' C2 M- I9 X% M% D; w                 di "cannot find pminplddy240"
                restore
                 
+ Q0 k4 a; j' n                  sum dx240
: n' a; u' ?  h4 s                 local tempminplddy240=r(mean)
( Z' }% e' [5 C$ @& _1 k) ^                 
) y: G# w" x: |                 
                 }
   
* Q0 k2 Z- s  e1 c  [7 s   
   
" n" K# Y: ?4 x. v! n   
/ Y7 v) Y3 y) @% E: @   
        gen pminplddy240=`tempminplddy240' in 1
       
        *-------------ddy240pmaxp左侧高峰值 -------------------
9 m$ q' A: r2 L% m- ?+ g/ y       
& {2 ?( `) z6 b" z% {0 z9 M         preserve
   
* W% s0 B3 v* g7 J# s   keep if dx240>`tempminp'
   
   keep if dx240<`tempmaxp'
. V" w0 \/ |8 u! \0 U' M8 a' u9 I   
1 s9 L% e. f; X! g2 q$ k) P   keep if ddy240>0
, e5 H1 Q! E! Q. ~# @   count
               
; `" V3 r) Z$ R! R4 C                if r(N)>2 {
1 M2 L( Z6 P2 a5 Y                 
                 pkexamine dx240 ddy240
   
   local tempmaxphddy240=r(tomc)
   
   restore
2 l& @- N' U6 h1 ?                 }
                 
        else {
9 _* s. i* H+ v: x- w, K$ [               
: K3 A: @4 p2 j2 q                 di "cannot find pmaxphddy240"
                 restore
               
                 sum dx240
                 local tempmaxphddy240=r(mean)
: e* y( K1 y& C4 s% _$ V                 
                 
# E# _, R8 E% A2 D& H7 z0 @                 }
   
   
( P) m! D7 N; m* L' D+ f% T- U   
0 t: L' h8 U, ~- F% I8 n& H" S$ m   
   
# i) z$ B, N$ r  k' a# w   
        gen pmaxphddy240=`tempmaxphddy240' in 1
6 {* q) y! P6 w+ X" }& }8 @4 M* a5 a       
8 ]# ~9 P3 K( j" M       
        *-------------ddy240pmaxp右侧低峰值 -------------------
& _  X9 O# C: p3 _" H: |       
         preserve
  
   keep if dx240>`tempmaxp'
   
1 Z6 N5 e- h( U1 L" X5 {   gen temddy240=-ddy240
9 b* \! n% ?2 A9 i0 \   
   keep if temddy240>0
4 t5 }, n) w* U9 M1 C) C   
% c! E( h$ g& F   count
               
/ f: J  q) }6 M+ O! b                if r(N)>2 {
                 
                pkexamine dx240 temddy240
1 S  e0 [4 ^- @% G0 `   
/ k: ^: s( n& p% k) T$ I( V3 R* z* s   local tempmaxplddy240=r(tomc)
   restore
7 U) C: i, B3 P( @+ w" O% Q                 }
; f7 z0 {, y7 e4 s  i( \$ J4 N                 
# S" \9 `2 p) R$ B9 n+ e9 Y        else {
. }7 p2 J2 A8 d0 r                restore
                 di "cannot find pmaxlhddy240"
7 s9 H0 K, R* N                 
) ^' M; K. j- _4 q/ r- f; U7 d& `               
) p) b: a- E9 A! P3 g7 N                           sum dx240
                 local tempmaxplddy240=r(mean) in 1
6 H+ g# O7 m9 `               
; |9 O( N' W6 t                 
4 \# n. b. K& x% q7 G9 z* t0 i# y  j                 }
   
' J' i. A& W" c" n   
/ F+ Q) B# z6 ~6 O' I, r1 k$ m   
- W( U' H, J" j- l+ V& P: j   
   
   
        gen pmaxplddy240=`tempmaxphddy240' in 1
7 u) i+ t2 P/ T7 {6 h6 C       
        *-------------------------------------
   
   gen price5d=v5 in 1/20
      
   gen price5yy=_n in 1/20
   
   gen price5y=price5yy*(dy240max-dy240min)/20
   
/ i; u2 w: e# s$ o* ?9 N   *price5y &price5d
   
' X7 B. d4 Y# t3 z! L; ?   gen priceny=price5y in 1
   
   gen pricen=v5 in 1
6 k) y7 }3 H. }9 o
     
   *priceny &pricen
   
) _+ T" P  X* Q8 R   
3 o3 n  n1 c' z& m$ z" {   *-------------------------graph------------------
  
" j- p0 X( b; n" m- o0 K" S4 j  tw(line dy240 ddy240g dx240)(spike yxis amean,lc(blue))(spike yxis mean240, lc(red))(spike yxis maxp, lp(-) lc(black))(spike yxis minp, lc(black))(scatter yxis2 amean, lc(blue) msymbol(+) mlabel(amean)) (scatter yxis2 mean240,lc(red) msymbol(+) mlabel(mean240) lc(red))(scatter yxis maxp,msymbol(+) mlabel(maxp) lc(black))(scatter yxis minp, msymbol(+) mlabel(minp) lc(black))(scatter yxis4 pminphddy240, msymbol(O) c(l[-]) mlabel(pminphddy240))(scatter yxis3 pminplddy240, msymbol(O) c(l[-]) mlabel(pminplddy240))(scatter yxis4 pmaxphddy240, msymbol(O) c(l[-]) mlabel(pmaxphddy240))(scatter yxis3 pmaxplddy240, msymbol(O) c(l[-]) mlabel(pmaxplddy240))(scatter price5y price5d, msymbol(O) mfc(black) mlc(black) c(l[-])) (scatter priceny pricen, msymbol(O) mfc(red) mlc(red)), title("`f'") xtitle("`1'") ytitle("density") legend(off)
6 ?0 d. `; j& w+ B# S+ O
* ---------可选项目----------------------
1 V  p, K" Z# G+ h
$ m4 x$ R' W! y9 ]! Q  {9 Q}
end

awer

旅途愉快 发表于 2012-3-28 23:23
$ D5 R. k8 H. R- F9 w4 j这段代码是重新复制粘贴到hello.do里面的？

对的，就是hello里面的那个

awer

maodeng 发表于 2013-1-28 13:13
3 M8 q4 f3 Z+ M% p0 q9 ~! M. hello 300236
7 }- q  c: T& _) [1 cstatement out of context
r(119);

抱歉，最近太忙了
$ f9 M) Q- T! L5 l# b: h
300236 我用下面的程序试过，没问题

你把 D:\stock\test\source 改成你导出的数据存放的目录，再试一试

2 e4 `2 {; Z& d* `% ccapture program drop hello

# G2 x6 ~" T+ a$ e: Z/ dprogram hello

5 b1 R$ L8 ^: ?) o quietly {
5 D" V! i* d, b  }, b* M
clear

2 X3 i# w4 s0 t3 w0 L: z% Wcd D:\stock\test\source

; e5 f- g, w. \) [* V4 X8 ~  k1 Oinsheet using `1'.txt

   keep if v5>0
   
   gen vt=_n
   
   tsset vt
) P: y7 W- e) M6 |
9 m$ f, G9 A* L& T+ x4 w( l   gen zf=D1.v5
8 D/ q& h/ A$ P9 b) Z( u
* a4 j& u2 g7 A9 Z' c2 r- M* h7 n1 A   gsort -v1
+ x( d' E: V, H4 e
& Z7 D/ C+ H3 s! @, p/ i) n: Z   gen avol=v7/v5
' t; V& Z9 D# ^/ m   
- U9 h% q! u, r( V0 N( U/ W. D   keep if avol>0
8 D2 F$ A! ?8 V1 u   
& R6 |, N6 ^" L8 P/ B4 p1 M( ~   drop vt
  
0 A" V8 t, y- d1 [4 I   gen vt=_n
; h' V& {" m7 X' p5 i) ~   
% s& N% A  r3 b5 g) K4 f   tsset vt
+ O9 A' B- i1 E5 W1 B: C7 i   
5 b+ e3 s( a' H7 b! a' W  {) T9 S    sum v5 [aw =avol]

   gen amean=r(mean) in 1
   
5 q7 h) G9 L* h8 q" m6 a/ N   * 总加权平均值 amean
/ }' u  {2 u. q   
3 [5 [( Q- r- h7 V& z" B   
: ~& I! s. c' d7 I1 {$ E2 e   keep in 1/240
. G4 m: {* @, L4 X9 q2 R6 F/ D4 A; q
   sum v5 [aw =avol]
# h3 i: \8 s' ^8 E) n! U3 w2 O) R5 L  u
   gen mean240=r(mean) in 1
4 y( h7 Y) U' ?; v: B0 I; a   
   *240加权平均值 mean240
- g) f5 ~9 g( J/ D4 m   
! k0 L7 \- M6 q   set obs 500
4 Y9 g- @! _1 \   
  `; R/ z5 b9 G% h/ J   kdensity v5 [aw =avol],gaussian n(500) generate(dx240 dy240) nograph
   
   *240 kdensity(dx240 dy240)
   
0 e: d1 x6 O' `  g  i  P6 @   sum dy240
   
   gen yxis=r(max) in 1
+ o4 a4 ^9 h1 [, B8 w   
; b* j7 X/ @  Y% _, F: [   gen yxis2=yxis/2
0 b9 k2 l9 G8 w* E' c   
5 ~9 ~1 E, K/ X7 A0 A: E; |; \   *yxis
6 I' Z7 q5 t. y! ?' i' v   
   gen dy240max=r(max)
   
   gen dy240minp=r(min)

    gen vtt=_n

! j/ ?" }8 ~( X4 ?/ E+ m" h    tsset vtt
) A' t" @1 f' a7 B( `" o$ P
8 r2 C+ A" k: p( S" g" W   gen ddy240=D1.dy240/D1.dx240
   
  U' I* z& s+ g   sum ddy240
8 G8 V6 L+ I6 _" j5 }% }+ ^. t$ m   
; }) O. F, g* r& J   gen ddy240max=r(max)
3 F3 S. P. E2 ]- F: C   
: }/ Z' A' S$ ?) x& a; ]# h% ]) w- \   gen ddy240min=r(min)
4 i+ F& t! P3 b1 x# f. u     
   gen ddy240g=ddy240*(dy240max-dy240min)/(ddy240max-ddy240min)
   
/ O7 h( w6 I6 B- b8 S   sum ddy240g
7 i0 x: A4 ]$ L' J( t. W8 l; h   
4 q0 b! N" \- Y: ^   gen yxis3=r(min) in 1
   
   gen yxis4=0 in 1
; ~: R$ |! c. i/ B: V' b  
  *ddy240g
   
   
  preserve
  
! ?: h7 W% o' K4 a/ f  sum v5 [aw =avol]
; e: Y! T. l0 ]" Z( C( \6 ^  
7 x" q0 y" @; a! C  keep if dx240<r(mean)
! t) ~0 b9 h; F  
2 ?$ J2 w8 i# j) n9 g( R% F  sum ddy240
  
. |& g1 W, @. E  if r(min)<0 {
7 E: s, [% a; }' X7 S  ]        
                keep if ddy240<=0
               
                count
               
% `7 \, ^7 ?! M: \7 j                if r(N)>2 {
                 
                 pkexamine dx240 dy240
( Y& d8 c2 _, t7 z- ?; f               
! q# S( b7 _2 |) [7 h5 P                local tempminp=r(tomc)
; \0 k/ c# G9 v* q                restore
# k# x, Y- F. `$ F9 j                 }
                 
2 Z+ C% q& p4 C" o        else {
6 @* j9 r1 Z' |; w: P) E- U1 b               
4 Y: e  h* n( E1 E                 di "cannot find minp"
                 
                restore
8 @: W& o' o6 s0 L: k7 `) w: W                 sum dx240
                 local tempminp=r(mean)
                 
                 }
                }
   else {
2 S, H; H! K3 c, P                 
                 keep if ddy240<=r(min)
                 
                 sum dx240
                 
  j! A0 Y# @) A6 t$ a) M$ T                 local tempminp=r(min)
+ {$ Y, k" M* }5 H) Y. t                 restore
1 K2 X* C" l$ E! y: ]                 }
! z6 S3 I4 k( p6 z9 f3 S
# e0 ^) B0 W& u" S
5 j. y1 o1 E/ b$ l
   gen minp=`tempminp' in 1
   *240 minp

& |; H* u2 A' z5 C
, a' W' o9 ?& k- p# R6 S   
   preserve
3 m# o$ T/ {3 V   
   sum v5 [aw =avol]
   
   keep if dx240>r(mean)
' ^+ c4 t. i2 S   
   sum ddy240
/ ]+ \- D6 w7 M- m/ X% A   
   if r(max)>0 {
9 v+ S7 S6 Z& S/ G               
" L9 N+ s6 g' J9 a0 M: d& g                keep if ddy240>=0
               
+ l8 R. ~5 e, x- H3 Y, L4 n                count
  }  `6 U5 D2 \2 j+ ?+ i               
                if r(N)>2 {
                 
                 pkexamine dx240 dy240
               
3 ?  B+ c! S1 t( g4 v% a                local tempmaxp=r(tomc)
                restore
( V9 M5 f. p% m, Y0 S3 C/ c3 W- M) j                 }
                 
                 
        else {
                restore
2 T) F9 G9 G! A9 B$ R" Q6 @                 di "cannot find maxp"
               
               
                 sum dx240
6 X7 v, K/ [, x: B                 local tempmaxp=r(mean)
/ Z) p' o# P! h0 W: @! H$ h                 
                 }
                 
$ f! a! _& z) u2 c1 W                 }
                 
" y3 p, l) |) a$ l    else {
7 ^6 r' }; }. J& B       
           keep if ddy240>=r(max)
, c% d, T1 {' ?, N! T          
, s7 |" f8 ?% l# D/ Y. S           sum dx240
          
           local tempmaxp=r(max)
/ q: ?6 J# n* q4 [           restore
           }
       
5 T1 m* \1 s* d% }) I' R9 C
  
. a4 e- a4 l1 Q* T, A2 K+ U
; b2 f* K! \7 R! g/ J   gen maxp=`tempmaxp' in 1
  W, z. G* R# ]   
& \& d1 m; i2 t) F' k2 p1 i   *240maxp
   
' U# w. t. i- J1 R- o7 Y   
   *-------------ddy240pminp左侧高峰值-----
0 M1 V, k* U5 p8 L: u- s1 J4 Z   
    preserve
/ w3 @" o" T' k! ]) |$ f8 ~- d   
   keep if dx240<`tempminp'
) |4 n8 w8 [% u) M! |3 F0 y. K   
6 @, B! N8 E2 J2 k& S   keep if ddy240>0
   
, l+ }( ?3 w1 A; d   count
& F" K' E. c* p, C0 G               
" z2 `' P% }8 `* a7 w( t0 U* O                if r(N)>2 {
4 E+ |- F, X( t4 Q: U! Z4 D& c                 
                pkexamine dx240 ddy240
- c2 `: X* H" D9 B* a1 j2 Z               
                local tempminphddy240=r(tomc)
                restore
                 }
; A, B4 k+ `( {7 f                 
; M5 t: ~- C# q: x2 \% j$ q9 M7 h        else {
9 ?  G. w: N  h                restore
; u$ B* ?: S1 ?) }/ z; p: V                 di "cannot find pminphddy240"
( a4 ~+ D% d0 X& C5 b                 
! D+ D, D3 e+ v$ p/ A               
                 sum dx240
                 local tempminphddy240=r(mean)
$ I9 g( F2 Y! L               
' _5 h! j9 Z# }5 c- Z                 
                 }
. E: c2 O+ M4 N; \                 
- f5 O- e% q. m8 P   
* n$ [: G! U& i" x4 _   
6 ~: J# X5 `9 ]8 v) }* S# X7 t   
   
7 v: a7 m# u# f9 }, D   
        gen pminphddy240=`tempminphddy240' in 1
       
       
        *-------------ddy240pminp右侧低峰值 -----------------
    preserve
# R# F3 p5 N9 a: D/ R$ h  W* u
   keep if dx240>`tempminp'
  ~; t' W( t" S" H0 ^: l   
   keep if dx240<`tempmaxp'
   
4 `, g! e3 ]" K7 A% @1 k   keep if ddy240<0
2 U( D$ @+ _- l4 ~) y   
& @! f# C" l% G( ~" A; b. n   gen temddy240=-ddy240
   
    count
               
. h2 c5 Q! |  q1 \4 n" d                if r(N)>2 {
                 
                 pkexamine dx240 temddy240
   
  local tempminplddy240=r(tomc)
                restore
                 }
                 
        else {
- _& z' h3 b- r8 ]! B8 Y: j               
2 @: r8 T7 G! }% Q2 k4 c* D( v                 di "cannot find pminplddy240"
                restore
                 
& o2 G/ [* ]; Q% r- B                  sum dx240
                 local tempminplddy240=r(mean)
) Y3 s+ d. u; }9 ~* F                 
: X8 T$ }9 N  \                 
1 c+ ?/ s  y) m0 @                 }
   
   
   
   
   
& k/ W' |6 t' \        gen pminplddy240=`tempminplddy240' in 1
       
+ Q/ V+ y7 t' }# l1 z6 [        *-------------ddy240pmaxp左侧高峰值 -------------------
       
6 N& b7 |0 h6 o4 [7 ~         preserve
   
$ k7 w6 |5 O" M- R% p- v   keep if dx240>`tempminp'
   
   keep if dx240<`tempmaxp'
5 F- z! y4 T, y9 S: E+ I5 x0 Q   
6 h% R- S. ?- c( N5 C   keep if ddy240>0
   count
               
" p5 M/ b2 e* m                if r(N)>2 {
5 r9 @1 E5 X% `7 I4 O! y1 \- s                 
                 pkexamine dx240 ddy240
   
   local tempmaxphddy240=r(tomc)
   
   restore
+ d" f% c8 E8 v2 }9 B! N( m6 ?                 }
                 
        else {
$ b; o* k& i* A* {/ e8 g% {               
                 di "cannot find pmaxphddy240"
                 restore
               
+ O  _; U$ a6 j% k- U# Q1 h$ H                 sum dx240
                 local tempmaxphddy240=r(mean)
" h& A5 r  \3 T, d0 w% _                 
  i1 X+ I# R0 Q$ F. q                 
                 }
" w8 l6 u. U3 V' u( c   
( n; }3 x, Y8 v2 V* R+ |( Q   
5 r' |- s) ]) h/ c7 j/ t+ E   
   
( S2 e, B3 x- q2 O* l: p   
% S+ L" C: ^- }- v! j: ?   
        gen pmaxphddy240=`tempmaxphddy240' in 1
, D) ?9 b& o. K$ S2 z       
       
        *-------------ddy240pmaxp右侧低峰值 -------------------
       
         preserve
5 W2 G! h; g0 s% Z* |3 u5 N' ^( S  
   keep if dx240>`tempmaxp'
2 c4 l4 ?, ~& C1 J& r4 s( N/ u6 n   
  p) F0 v8 @7 e2 ^6 x   gen temddy240=-ddy240
   
1 }; ~1 b3 Q" S) M   keep if temddy240>0
0 Z: z& g, u% [   
' [" _6 [% i) s6 h( P+ q" h' p   count
' C7 r* C' [: ?# r. }8 X               
8 {  c6 Q! q  G) v% J# j% J: g- k, B                if r(N)>2 {
                 
                pkexamine dx240 temddy240
5 q- Z5 j9 u/ i* O   
3 R1 t9 Q- o8 i* e- j   local tempmaxplddy240=r(tomc)
   restore
                 }
6 ]$ |! ?. P! b# F                 
        else {
                restore
                 di "cannot find pmaxlhddy240"
2 B1 `4 ~  j% P- f; {                 
               
+ z8 Z3 _, y) {                           sum dx240
' n1 f& z6 H9 \                 local tempmaxplddy240=r(mean) in 1
               
; c& f7 p% {; x% @                 
                 }
   
   
   
   
" Z% e2 S7 z$ c9 \5 C. K* c6 w   
   
        gen pmaxplddy240=`tempmaxphddy240' in 1
       
0 I  s$ k- e2 ^  ]! S9 n0 G        *-------------------------------------
   
0 \: P# Q& h) _   gen price5d=v5 in 1/20
$ [; I% {% H7 R      
2 j0 f" x0 n$ E; r( S6 }6 a   gen price5yy=_n in 1/20
& E! y: n( H$ q0 t  g# a; Z   
. s+ x9 l6 S+ P; A  r3 L   gen price5y=price5yy*(dy240max-dy240min)/20
   
' r- g- p1 ^& u$ X  q/ {5 d$ I$ v   *price5y &price5d
9 u6 B# A! [( @; I* ]- C- ?   
   gen priceny=price5y in 1
   
   gen pricen=v5 in 1
- Y: \6 d7 [6 f, S2 y; f9 l
     
: M5 @- h0 z" b6 j% {9 ~* ~   *priceny &pricen
. E0 Q  ~' e& T& i5 |) P5 u- ?4 ]' p   
   
* T: ^5 K% @, W+ G% r8 ^/ L   *-------------------------graph------------------
/ [1 p$ P4 W0 Q& X  
  tw(line dy240 ddy240g dx240)(spike yxis mean240, lc(red))(spike yxis maxp, lp(-) lc(black))(spike yxis minp, lc(black)) (scatter yxis2 mean240,lc(red) msymbol(+) mlabel(mean240) lc(red))(scatter yxis maxp,msymbol(+) mlabel(maxp) lc(black))(scatter yxis minp, msymbol(+) mlabel(minp) lc(black))(scatter yxis4 pminphddy240, msymbol(O) c(l[-]) mlabel(pminphddy240))(scatter yxis3 pminplddy240, msymbol(O) c(l[-]) mlabel(pminplddy240))(scatter yxis4 pmaxphddy240, msymbol(O) c(l[-]) mlabel(pmaxphddy240))(scatter yxis3 pmaxplddy240, msymbol(O) c(l[-]) mlabel(pmaxplddy240))(scatter price5y price5d, msymbol(O) mfc(black) mlc(black) c(l[-])) (scatter priceny pricen, msymbol(O) mfc(red) mlc(red)), title("`f'") xtitle("`1'") ytitle("density") legend(off)

  d/ I) I% r: t* ---------可选项目----------------------
# |, [4 \' g! }! p*(spike yxis amean,lc(blue))(scatter yxis2 amean, lc(blue) msymbol(+) mlabel(amean))
9 v* a3 R: M+ c}
end
% Y- _, t1 c9 r