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

awer

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

计算之前需要有 ：

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

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

然后

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

e:\stock\source

2. 股票数据导出

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

b快捷键“34”进行数据导出, 在弹出窗口选择
“高级导出”,

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

将导出目录选择为你新建的储存股票数据文件夹的目录（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

(002128为你想查询的任意一只股票的代码)

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

2 _( G* T! H+ r3 |  Z

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

2011-8-13 18:02 上传

点击文件名下载附件

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

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

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

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

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

山寨完毕， 收工!!

awer

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

awer

沉默的前哨 发表于 2011-10-19 20:49
0 M2 m, S; h* U能不能帮忙看看：

do hello

  Z  h! q" y. ^- }& S你算得哪个股票？

awer

沉默的前哨 发表于 2011-10-20 10:32
$ n  |$ c4 t' R) n7 m/ ^( Y3 Kdo hello 后变成
/ F+ B' q. c/ B# u& W......

awer

沉默的前哨 发表于 2011-10-19 20:45
$ @: j+ O" F6 |5 p# H+ ~9 v为什么老是提示：type mismatch

  U! O+ Q# I0 i: ]注意 日期格式要选成  yyyy/mm/dd，要不然时间变量格式不对

awer

琴心三叠 发表于 2011-10-20 21:37
0 m) s, `7 ?" R1 \+ M谢谢答复  还有一个问题  我导出的文件为何只有当天的数据  而不是历史全部数据呢 ...

, b/ i% j$ ]' G$ \/ Y8 r有很多问题是数据不全造成的
5 w% a5 g: ^% \8 R如果想要全部数据，
系统--盘后数据下载--日线数据--日线和实时行情数据

把开始时间设成中国股市开盘的那天，下载全部数据
3 t. d* k3 \3 z# j
$ q( C7 H$ m9 C  }+ W3 o% o5 a下载完全部数据后，在屏幕上使用快捷键 34
. o: H! f4 R* G/ `5 j  h3 ]
34--高级导出--
3 F+ }! H5 ]- d, m9 j; b0 A, x
3 k. X! Z" F' s8 o然后选中   日线导出；选定导出目录；文件名设为：XXXXXX.TXT;勾选精确复权；不要生成导出头部；分隔格式为tab（）；日期格式为YYYMMDD，然后选择添加品种
% P, `8 F$ |" K7 ]$ \8 }: m( i
我一般都是把所有数据下载到指定目录下，然后进行hello运算
5 X* e" ?. O0 a* K

awer

改了一下，这次运算应该没啥问题了
, d% I1 R0 b6 e0 ~
. A6 q% X3 ~, `, g
7 N! [' u; ?' Icapture program drop hello
: W2 G2 }5 p. X" s
program hello
' ]8 _0 g- B* c% K/ C; z- q
, b: y. G) o5 P quietly {
4 q; a4 f8 l( d9 D" j
clear

0 d: p1 \! b3 e! a+ k, X, K" N- ycd e:\finance\stock\20110101\source

insheet using `1'.txt

  K/ w6 J; [6 d$ t   keep if v5>0
% W' x4 h! _& I. n0 e- w0 d1 M   
. b9 t3 `3 L% Z0 s  I0 `) w   gen vt=_n
6 G9 m- U( Y4 C4 L   
. |' X+ f6 \6 p- X4 P; n& v' _   tsset vt
+ R$ E/ f4 o& H1 m/ U
; X2 F9 F* A' i0 d   gen zf=D1.v5

5 m3 U9 g" d# R9 \' S7 x   gsort -v1

& g% {3 |# }& A# j4 e7 O" _6 k, k   gen avol=v7/v5
8 [/ p2 N  |* H/ u9 c1 m+ X   
   keep if avol>0
3 Z& }& B/ R  Z4 A   
   drop vt
3 }# ?- L' B' ~/ \  
   gen vt=_n
8 E0 D5 g6 _% m1 p1 _   
: l) z: d0 ^9 D1 Q/ Z0 O   tsset vt
   
: p* B4 G, J! K( ]% j) }' t; I    sum v5 [aw =avol]
; u- C/ q: {+ E4 b5 p9 [
) ?1 x. i1 E9 d   gen amean=r(mean) in 1
   
   * 总加权平均值 amean
   
   
   keep in 1/240

   sum v5 [aw =avol]

" E) F! C+ C2 I   gen mean240=r(mean) in 1
   
$ ^  _& \; `* k, D; w# }   *240加权平均值 mean240
& w0 Z  ^5 w9 j. I% u! n* Q5 L! {   
   set obs 500
   
! x' q4 e, q) ~& J" x2 X) K   kdensity v5 [aw =avol],gaussian n(500) generate(dx240 dy240) nograph
3 V3 H2 n2 ^7 z1 N  J   
  r( H& F' N1 _$ c   *240 kdensity(dx240 dy240)
% R' D+ _2 r. x! r   
5 `* P0 e" r' n) K2 o1 Q, s   sum dy240
) p7 X0 ]  E- A. D5 {" }   
& F+ n8 v/ K5 g- L2 {   gen yxis=r(max) in 1
1 ?0 `6 {& f2 N. R, H   
/ X) \* c/ Q4 \4 z* u" B: A   gen yxis2=yxis/2
   
   *yxis
# f" L% H0 C, w5 X# p6 r   
   gen dy240max=r(max)
   
6 r6 |: J: r. x4 t5 j, p( X; q   gen dy240minp=r(min)
8 q6 R. h( q3 @) k
    gen vtt=_n
: ^+ }$ [5 R% I1 B
    tsset vtt
. b$ J) y% s3 C7 p& T/ N. j
* c: W/ v( D  y" i   gen ddy240=D1.dy240/D1.dx240
   
   sum ddy240
/ w2 o# c6 y8 }& t2 r   
   gen ddy240max=r(max)
' A( S! T5 n% h" V+ n   
   gen ddy240min=r(min)
     
   gen ddy240g=ddy240*(dy240max-dy240min)/(ddy240max-ddy240min)
. p3 P1 g0 K& g6 e   
   sum ddy240g
   
   gen yxis3=r(min) in 1
   
   gen yxis4=0 in 1
- b. U( C' i2 r; D3 G  
  *ddy240g
6 k5 O- I5 A8 ]) ]) P5 I$ V( j; M: V! w5 ?   
/ r4 v# O7 z0 Z9 y6 F   
! n2 b# o. _+ z1 P0 c  A" g$ M5 I4 Q  preserve
6 S( H6 F; r5 X' s- b! w  
9 H) \, \% T1 X9 S: y! ?6 y; _  sum v5 [aw =avol]
: T- b/ X# k: g2 S  
! ~5 V4 ~; U$ \% z3 K5 [  keep if dx240<r(mean)
  S2 Y$ u+ X" l2 {0 Z+ i  
  sum ddy240
  
1 q  S) I- M7 @' x/ |( k/ \/ F  if r(min)<0 {
        
' ^! I& d! }( F8 n: \% R                keep if ddy240<=0
% k6 F  O+ f, C+ m               
; v1 M; c+ `/ u0 g                count
  z" ?) {( b* ?1 z$ O               
& Y2 J- V. }9 A                if r(N)>2 {
2 n0 M; Y/ T5 B8 E5 P                 
                 pkexamine dx240 dy240
               
                local tempminp=r(tomc)
% Y7 S0 S7 L- _2 g1 `- D# ?                restore
                 }
& P- H. G& E* ^1 X7 y                 
" ^$ B! X# W8 n6 \        else {
* \0 N# g! v6 h$ g3 E               
                 di "cannot find minp"
2 Y' p4 }+ \1 U# i* Q) r+ l9 v                 
+ B( L9 \. J  ]# ?                restore
* k, [3 ?7 S; d( }" M                 sum dx240
                 local tempminp=r(mean)
                 
                 }
) Z5 I# f) s7 v; x) q! S0 ~1 P" [                }
, ?9 g! W4 h/ |- z+ r* J0 T   else {
' u; v. A' D( ?  z( S# ?* }5 X                 
* E; e6 H  h! H0 }0 L" W6 Y                 keep if ddy240<=r(min)
                 
1 a) G( B5 H( D) \2 Q2 `0 L+ `                 sum dx240
                 
                 local tempminp=r(min)
                 restore
                 }
' c) `9 K; ]3 N3 |( t2 O/ l/ \
* u6 S7 @: t4 T. ]' f$ H
( O2 y0 i, q+ Q* f" H
  n! N4 k9 y3 u4 c9 b   gen minp=`tempminp' in 1
: x5 s0 L1 A  E& m   *240 minp
5 |8 c# C* u. F9 o0 y
" L9 q4 r- W7 M8 l) Q/ q4 y3 R; C
! X  a; u3 A5 ^* ]* D; i   
- }" \* Z3 s& C- |   preserve
2 ]( G6 k$ \8 X0 ?' p! i   
   sum v5 [aw =avol]
- e$ r  U% B$ }4 m/ A7 U2 E   
" R, |2 x0 m# p" ^   keep if dx240>r(mean)
! O0 f* I# z/ B' I* C" v   
, e, z  p6 D6 g- r# K5 e2 d   sum ddy240
* ]7 `- F# y6 l" G   
+ O% Y2 M7 @! S7 S5 v9 t   if r(max)>0 {
               
                keep if ddy240>=0
- g) \! q( L$ ?/ g- A               
                count
               
                if r(N)>2 {
                 
                 pkexamine dx240 dy240
2 ]% ?% Y& b% x: G% X               
# e- t# b" o3 y& F' K# {, y                local tempmaxp=r(tomc)
: ~- O! }# L' A+ k1 q  t" d  g' C                restore
; `8 k' y5 V5 b/ e1 D+ o4 H* G                 }
5 R' D. N- M  B& l                 
                 
8 Y/ Z3 X7 v8 h9 m" e+ t# A        else {
                restore
; l) Z# @# u' u) Q4 ^                 di "cannot find maxp"
; I3 ]% y. F/ I6 y               
4 Z8 a  D- i/ E5 x  M, [* N  H. c               
                 sum dx240
  V+ R% W) L1 k" s  Z! u2 m% K                 local tempmaxp=r(mean)
                 
                 }
7 U- x6 p3 e+ [5 a                 
                 }
4 D; B2 ^6 m8 V4 w0 `9 j: F, K                 
  {) Z1 B/ [1 N- Z$ O6 l) D    else {
       
           keep if ddy240>=r(max)
          
           sum dx240
6 S  y, P, u1 P2 g: v3 O          
           local tempmaxp=r(max)
6 O! D" I% j- a( s$ g5 F           restore
           }
- T; d/ z$ r* L' r: Z- c       
- v# f. f% ?- ?: I# d5 x  a7 f
  

   gen maxp=`tempmaxp' in 1
. z$ D' D( |0 u' \- \2 H" A   
   *240maxp
   
   
6 F6 H( K+ [% w2 B0 Q/ b: j1 c   *-------------ddy240pminp左侧高峰值-----
1 r+ Y$ i6 j2 z" w9 x: N% p3 @   
    preserve
8 f2 y: v% d. |. H, L: b( ^   
) }/ T* r5 j0 A; S( o   keep if dx240<`tempminp'
   
6 G. J- z. c# h# [   keep if ddy240>0
   
   count
* L3 z$ i1 E$ S               
- _- a1 ?6 D* z4 `. c                if r(N)>2 {
                 
& {0 Q0 W1 x; j) T8 {                pkexamine dx240 ddy240
               
                local tempminphddy240=r(tomc)
                restore
+ L% R1 Y/ x8 z5 B                 }
                 
        else {
. ?" q) c' [4 n8 w) p. x                restore
8 S) c8 ?2 C0 F# i' w; j/ V                 di "cannot find pminphddy240"
7 B4 ~3 O  I, G" N; u                 
               
) {$ v% K$ d- S: l" `' B                 sum dx240
$ S, F7 u2 q7 S' p4 ~- i; T                 local tempminphddy240=r(mean)
               
                 
                 }
                 
   
   
( e0 c7 Z/ G9 @' s  A' c& H   
! _, P/ `6 D1 W" r# w3 M3 L; _   
2 F8 ]3 p2 d( \* C& P. {- B# s/ d   
        gen pminphddy240=`tempminphddy240' in 1
/ q" l/ h* Y0 k4 D: U       
       
7 ~3 I0 k/ ~# |9 @0 S        *-------------ddy240pminp右侧低峰值 -----------------
    preserve
/ z: S' |5 j1 N/ s# p
   keep if dx240>`tempminp'
0 |1 j/ p* c" ~8 c. N( F2 `   
& h/ F5 ?0 b1 w" ^* J  i. I   keep if dx240<`tempmaxp'
   
   keep if ddy240<0
   
1 B7 G6 I' H3 w, j, p   gen temddy240=-ddy240
   
    count
9 N4 C' V; D7 t' ?, y2 j1 z               
" o5 N. }! H1 \- M/ g                if r(N)>2 {
; B2 B7 b5 ?1 x% |) b  r4 r- D                 
; `5 ?; I. b; |" s6 U; h8 k                 pkexamine dx240 temddy240
2 B9 ^, M& N7 l" M   
1 H$ R7 J; v" r  local tempminplddy240=r(tomc)
                restore
                 }
                 
% x: }7 f& {7 {+ y8 m7 C& y        else {
& ?* B7 X( t4 ^2 z$ M               
" j; J. \  P2 T- D  ^                 di "cannot find pminplddy240"
2 W8 r+ U6 P7 H- ^1 W! I                restore
& A2 e2 J. d7 z* @! k( E                 
                  sum dx240
                 local tempminplddy240=r(mean)
, M& ?( i, _% J                 
% n' S6 o+ E" X                 
9 R0 d$ f4 `% _+ P; [                 }
& [3 H; a* x8 y2 Y7 v! D& Q# L   
   
* b6 M4 @: z, u0 g; ]4 n   
   
* T. ?/ F7 Q) }8 e   
        gen pminplddy240=`tempminplddy240' in 1
       
        *-------------ddy240pmaxp左侧高峰值 -------------------
9 R. l! i1 J9 M* P# m. p# W% A, q" i       
' j1 Q7 K5 F% [  f& e  b# s         preserve
1 G4 i# {4 O0 ^; h! p5 \   
   keep if dx240>`tempminp'
& i- q# _% x' j4 k+ P  I   
   keep if dx240<`tempmaxp'
: ]* y" {4 q; q) V+ O/ ~) V   
   keep if ddy240>0
   count
. r" `" h' v7 o) R, C3 E( }               
                if r(N)>2 {
                 
/ ]+ ^: M& U$ e. \% B6 d0 H8 L                 pkexamine dx240 ddy240
   
8 a' F0 n. h9 L5 Y; q   local tempmaxphddy240=r(tomc)
   
% L; Q$ j: ]5 U( @/ y9 Y   restore
: |2 e" X! Y# Z6 [                 }
: _% g3 E% f. A                 
# K0 D3 h3 Z6 R. O        else {
; p4 M* n8 a( e- l, D: {( a) s' {               
                 di "cannot find pmaxphddy240"
                 restore
               
                 sum dx240
                 local tempmaxphddy240=r(mean)
, K1 G! G. r& o" X7 `0 Y  i                 
                 
                 }
) B8 M5 f# C3 X5 |   
   
   
# D$ @. y. ~' Z7 {4 m   
5 b7 ^3 l1 W* M   
# i4 R+ n- A9 M& H- \8 |   
; u% L3 W% I" ]9 \' Z( G- F4 d" @        gen pmaxphddy240=`tempmaxphddy240' in 1
3 |+ ?6 @4 F) I2 J0 O- {       
       
        *-------------ddy240pmaxp右侧低峰值 -------------------
) C: T4 p% ?* r& K3 ~       
/ k8 B* U9 _6 D8 l# L  [2 |0 c         preserve
: h. N, {* J0 f; S  
- R/ `7 K+ r7 U# A! \   keep if dx240>`tempmaxp'
! P) E  C2 R6 y$ R, r) X9 m2 {, l' u   
' g" k9 R- i' W2 y' W+ b; C( K, M6 x   gen temddy240=-ddy240
   
& q( P$ U' j# O7 z1 @* Y   keep if temddy240>0
   
5 c; l; R" H! ?$ S% v9 E( q   count
               
                if r(N)>2 {
                 
                pkexamine dx240 temddy240
   
   local tempmaxplddy240=r(tomc)
   restore
                 }
                 
, A$ `. ~1 D) O- I' q+ b( g8 i# S        else {
                restore
  f$ T: k  y2 N8 P9 T2 \5 K                 di "cannot find pmaxlhddy240"
" m7 _: q( Z% R6 u9 x' i4 G! p                 
               
                           sum dx240
# a$ e! W  Z6 G                 local tempmaxplddy240=r(mean) in 1
0 A; i3 t. O) [+ r* x* v               
3 R4 ]/ ]' w# g2 E5 U/ u6 B                 
( h/ n1 \  f8 j: J! u                 }
   
   
# B3 r- s7 A: l1 ?' g' O   
   
0 k3 l: g5 v/ E0 `# ^7 A. ~" s2 _   
   
        gen pmaxplddy240=`tempmaxphddy240' in 1
       
" Z  Z- E; O7 ~. ]        *-------------------------------------
) |3 O7 _7 P9 N" g9 e   
   gen price5d=v5 in 1/20
/ Y& O8 H; D, Y9 L* F' l      
   gen price5yy=_n in 1/20
   
   gen price5y=price5yy*(dy240max-dy240min)/20
   
, }; Q: k! E+ r% l2 A   *price5y &price5d
0 G; Q' _+ o' ~2 g' ]   
   gen priceny=price5y in 1
, W: N; v7 c) A+ g; V   
2 Y8 D% R2 W$ W' @+ j9 o$ i! G5 d   gen pricen=v5 in 1
6 ]. `  m. m6 P2 E% \! V+ T
4 L5 E* e' [* Y, a! I     
' [: y+ k  t1 f8 a   *priceny &pricen
   
1 w# E% K/ v1 M7 f4 M   
+ }( K4 j6 [" e6 v   *-------------------------graph------------------
9 E( U. z+ e# D6 m: u4 g  \, t  
7 ^$ K5 m4 J# t3 O7 [6 ^  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)

* ---------可选项目----------------------
: |* L  h4 k3 l, l6 g, s
5 U' W  v5 r* J7 h! U6 _( S}
end

awer

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

; G4 |5 a- S6 L4 O) v# V; _6 s% p对的，就是hello里面的那个

awer

maodeng 发表于 2013-1-28 13:13
. hello 300236
statement out of context
. c8 H0 Y; A6 w5 g9 S1 q- s  v0 ir(119);

3 {; v9 E; R/ x/ Q5 H( Y抱歉，最近太忙了
  g  F  q+ I5 f+ B( @
300236 我用下面的程序试过，没问题

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

8 x, H9 D6 S/ ycapture program drop hello

5 {0 R( L/ T6 ^program hello

7 Q# ^" V$ p1 v$ w+ L* i quietly {
) N% J8 }2 V. Y, ]2 o1 g4 S
6 F) k2 `) j+ G2 n* f; Cclear
9 e; W& m% W, I+ H9 `$ D
cd D:\stock\test\source

insheet using `1'.txt
5 C( c1 ]1 Y& c# V) j4 j1 r2 S8 d
! i- w' B+ P+ D4 M8 t7 e$ Q0 \   keep if v5>0
8 `2 c* a# a, o2 M8 z5 a. L8 ~   
/ p/ h) {+ C  O/ Y; Q, K0 \2 [( @   gen vt=_n
" ]/ Y& C1 X1 w9 X) I5 S   
   tsset vt
% R+ e0 C7 q" ?# h# z, e
   gen zf=D1.v5

% w" O4 ^0 t/ C: y   gsort -v1

   gen avol=v7/v5
   
, ]( p% |, Y0 h2 @4 `  Q   keep if avol>0
   
   drop vt
  
! k" r- a* s* i) {" O   gen vt=_n
   
   tsset vt
# T, q1 [9 r# F' I" X   
    sum v5 [aw =avol]
& @* d( }! _. }( P
, Q) h, F) [8 K; y2 R   gen amean=r(mean) in 1
7 O, ?* i, w  W- A5 [" g   
: p9 }- L: K: R7 d6 r3 [   * 总加权平均值 amean
   
   
8 ]& l" n" @# f, e   keep in 1/240
9 {* X( G0 t# N# t
   sum v5 [aw =avol]
7 A( _) T5 v& U5 v0 t5 i7 _9 i$ [4 E
8 |4 E: r6 b4 t4 k- I/ A   gen mean240=r(mean) in 1
   
   *240加权平均值 mean240
& U. M  |/ T+ }! @8 T, t   
8 k1 p- j! F2 \4 F6 B   set obs 500
3 M0 s0 L2 ^! h0 ]$ O  W* }/ U   
   kdensity v5 [aw =avol],gaussian n(500) generate(dx240 dy240) nograph
% L& @% o  o+ I' d, W   
   *240 kdensity(dx240 dy240)
   
   sum dy240
" c/ f2 D0 I2 x0 _   
/ O0 n+ m0 a9 d4 b8 o   gen yxis=r(max) in 1
0 e0 V! d9 s% {/ X" U   
+ b" ^' }1 X' u" Z( ?   gen yxis2=yxis/2
' k, t' R' P; a- ?3 E   
/ ^6 Q5 O: K6 m' H* F! R7 }  J   *yxis
/ Z# e2 O; _' ^3 c   
- y- X* w+ |. ?   gen dy240max=r(max)
   
   gen dy240minp=r(min)
$ c8 s0 k0 r9 V+ h# ?
& m" `1 t9 g3 Y% D2 E' j    gen vtt=_n
! _& \  i) n0 d. n; n$ n& W$ U
    tsset vtt
+ N% \+ U  ^2 @& Y/ l
9 O$ f+ @4 k! v# K   gen ddy240=D1.dy240/D1.dx240
; d6 a+ {) r) f: R   
" O6 N% j: y2 T+ g" `. V   sum ddy240
   
   gen ddy240max=r(max)
   
   gen ddy240min=r(min)
     
# h9 M9 D1 A  a   gen ddy240g=ddy240*(dy240max-dy240min)/(ddy240max-ddy240min)
. R1 z$ r( F. P& f6 L   
   sum ddy240g
   
. z9 B- w1 n) U   gen yxis3=r(min) in 1
( S  a2 X- B: m- y3 H   
   gen yxis4=0 in 1
* f6 R2 U* E# {2 m1 I- q  
  *ddy240g
! s8 i- I8 D* F* I+ z8 Q" f( {   
1 Z4 ~/ h0 u4 O) y# B5 y5 f   
% ]5 ?* V2 g. w" j$ u( \9 [  preserve
  
$ _+ @3 c: f2 J& G. y  sum v5 [aw =avol]
  
  keep if dx240<r(mean)
  e7 Q& J9 o" D  
  sum ddy240
  
  if r(min)<0 {
! s: K4 W' e" H        
& K: q2 |; X' M                keep if ddy240<=0
               
                count
7 T- K/ ~- X) c; C1 o9 Z, v3 H7 \' I               
                if r(N)>2 {
. ?, g0 ]& J% Q                 
1 ?! n+ q  K) u: B5 g. l                 pkexamine dx240 dy240
               
                local tempminp=r(tomc)
  S1 y$ P# j/ d$ P# q                restore
. {* E* o+ `7 F5 D* `8 J8 F! ?                 }
                 
2 u/ W* _- ]" d; M- k        else {
               
                 di "cannot find minp"
                 
                restore
                 sum dx240
                 local tempminp=r(mean)
                 
                 }
. o8 E' @1 P( j6 h" C  c                }
8 |! {8 l+ l7 b   else {
- G; ~3 U: \( @& ?9 W" E# K6 x) g+ {* i                 
8 D7 ]' A8 t- k- t" t6 j( ]9 N                 keep if ddy240<=r(min)
" r8 ^( q' `+ ]2 }                 
                 sum dx240
                 
$ f/ y6 \, `0 ?; L. b  t                 local tempminp=r(min)
                 restore
                 }
) D# O! S& A7 ^7 Y0 c; G, T% G
/ r6 m" ~& W: ~( {8 _) F; Q
) j1 i# Y* m; [/ G/ w. M
9 |4 T/ }; B5 l3 \6 L2 n. L! I' N   gen minp=`tempminp' in 1
- `; p; r$ l3 K; Q' @   *240 minp
- l( k$ R1 `- B1 U" |
7 c! f* g; G9 m
   
1 ~0 u7 o* W, }/ i; B% v  n0 W   preserve
   
3 V# b  }: a5 v, j   sum v5 [aw =avol]
$ Z9 ?7 c6 |4 t; r/ r   
   keep if dx240>r(mean)
6 ~. n! y3 F+ U   
: y  {1 [) w8 k& m; {/ X! w   sum ddy240
   
4 l" C8 x1 Y/ l0 m  l3 h9 |3 ~: t3 t   if r(max)>0 {
               
                keep if ddy240>=0
               
( a' ~# Q% a" j2 W8 }0 W' H7 Z1 O                count
1 h) r& h% U1 O: n# Y  h               
                if r(N)>2 {
- K/ d+ D2 _: `; G$ J% X                 
) ?- r4 b  r. f5 y3 i  y5 ]                 pkexamine dx240 dy240
               
: |. a' k1 `8 Z( n" k2 O                local tempmaxp=r(tomc)
4 b5 B$ \; K2 S* o- v                restore
                 }
, r+ U0 W4 R0 H+ L, J                 
7 s- y  b9 R, z) L/ P4 t* J& G8 B6 L, ]                 
1 R. x" Y1 H5 r        else {
# G7 [% i: L9 I$ u                restore
                 di "cannot find maxp"
               
% E) u, O9 z- J& ~$ d9 z2 `7 n               
- q+ q% B; |# u6 |! D4 d. O+ F                 sum dx240
                 local tempmaxp=r(mean)
                 
                 }
( n8 |4 h" i+ z3 U7 U                 
) o$ T4 Y1 d$ f6 |4 C: f+ H' _7 n3 [                 }
                 
' p6 U0 N, k# p  A: v) S# Y    else {
+ ?8 P) v! W% @  ]# j1 i; |0 M5 Q8 c       
           keep if ddy240>=r(max)
4 g1 W' I, L- [  t( q& j          
           sum dx240
          
& L  ?& C3 p" t% A; L           local tempmaxp=r(max)
           restore
* u- @; @' {. d* ?- [           }
, k% Z/ }2 d$ O       

  
* X2 r) A) l6 |5 P
5 Z, [0 Z% _! W) {5 }7 ?% d   gen maxp=`tempmaxp' in 1
4 c4 y+ z- Q" t/ y+ U6 J  H9 o: H   
7 n  R1 j% g5 b( l/ `: l% T4 y   *240maxp
   
   
5 @! f7 a. ?8 F  @5 r   *-------------ddy240pminp左侧高峰值-----
   
- D7 r* G$ `* }# h9 Q+ I    preserve
   
( v+ Z5 @" k  d8 D; E9 i  E9 Y% z   keep if dx240<`tempminp'
   
! i$ ^$ O0 K8 V! e* C, U# m   keep if ddy240>0
   
   count
               
                if r(N)>2 {
4 i1 V6 k1 t- f. T# d- x                 
                pkexamine dx240 ddy240
% T0 j; ]2 ]% k               
) q- a# O' B: C/ X2 r& T                local tempminphddy240=r(tomc)
0 r! B# F5 A7 M3 @. u- Z  m                restore
) ]6 l. }! Q- Z% M& y* R1 w                 }
                 
# v+ e  ?# D5 K- F/ y8 o        else {
" u$ J; u0 c- p  h2 A0 e                restore
: f5 B5 L( B4 ?% n: L2 g' ~' G( x                 di "cannot find pminphddy240"
8 G5 O! b8 r0 h# k  c) R                 
0 i- L7 e7 q% g. x% q3 {) M7 Y  i0 s7 P               
                 sum dx240
3 M9 K1 p; ^9 D+ m! k' f                 local tempminphddy240=r(mean)
               
' C$ L) l+ T/ m8 b                 
                 }
0 R* l$ R* u3 `8 M: H2 _0 W                 
   
* X0 g# l, ^7 Z5 \   
1 ~( {! X6 M; A' S( A: S! W   
   
   
' p+ E3 K# p  s. H        gen pminphddy240=`tempminphddy240' in 1
6 @& \" P2 U# i) k9 b4 [; s" A       
4 v6 J) H; j$ U       
9 J& @# A% [2 K        *-------------ddy240pminp右侧低峰值 -----------------
  D! j, z2 U  F0 \+ z" j" ?1 F6 ]    preserve
% l) S, Z; u& |
/ s( R) e9 ~4 [( g' w6 q   keep if dx240>`tempminp'
# M2 _) \! Z/ u+ T1 Q   
0 w& M0 U! ^, R' [' [* B   keep if dx240<`tempmaxp'
- Z, m' p' {1 [2 i# F/ x   
   keep if ddy240<0
# v6 ^- L3 D% Z* t2 J0 P+ @   
   gen temddy240=-ddy240
+ Z) g  N: n: j& v   
' T1 Z7 }* `$ B% Y* U$ j& p    count
               
  P6 l- F" ^& R! q                if r(N)>2 {
0 N6 C3 \# L. [5 k4 {2 ~                 
& u2 o+ q: Q6 ?) A: S! U                 pkexamine dx240 temddy240
/ k1 A& a/ R+ P8 U! S   
( v( v  R; c. X1 p) a( T( _  local tempminplddy240=r(tomc)
5 k# y# @1 M9 I( M1 @7 c, e                restore
" Z* l7 A4 a8 a+ E/ E                 }
1 j* V8 l5 m9 g4 N% N                 
        else {
( Z- S% J! O! Q5 q               
                 di "cannot find pminplddy240"
                restore
                 
# B6 ^/ H' y. a$ S                  sum dx240
                 local tempminplddy240=r(mean)
                 
                 
                 }
   
   
   
$ C* l. n$ L8 }  {. C   
1 x2 T. _. A3 Y% m* L3 L' T   
* O5 b+ a1 [( b8 s0 p' N# ], P8 N        gen pminplddy240=`tempminplddy240' in 1
, h7 v- {" B# g  E( f       
. M* f1 P0 J5 p( _        *-------------ddy240pmaxp左侧高峰值 -------------------
2 b! @$ t- D' s/ m) D4 }" k  x       
         preserve
6 ~* D) ^7 D( P6 N2 a' j% G   
1 H  w0 S) K& ]+ r* S   keep if dx240>`tempminp'
   
   keep if dx240<`tempmaxp'
   
9 K) l/ _/ `: Z  O. |' C. @   keep if ddy240>0
   count
  \; V5 P2 B2 p3 V* i               
                if r(N)>2 {
1 V7 _* Y+ v: b5 b                 
                 pkexamine dx240 ddy240
   
   local tempmaxphddy240=r(tomc)
7 r' e" A: A5 @$ _1 h   
   restore
                 }
" u$ @( j" f3 x7 u# d7 U! a                 
        else {
/ i9 ~+ H- t& S0 d6 J  r  o               
                 di "cannot find pmaxphddy240"
                 restore
8 ~' F- m7 \' \8 R# \  d- f& _               
                 sum dx240
8 j. i( o3 H6 l. q/ C6 y1 m  o; U                 local tempmaxphddy240=r(mean)
/ \) H% v3 \* e* ^3 R& k" K                 
* |5 x/ }8 q# a4 V' {                 
# f! V6 t9 [5 A2 s                 }
   
   
   
- _+ u* O+ v7 Q; W  p; u' s' s   
8 j- ~9 v  _5 q3 Z% \2 C" W   
- ^  ~8 B7 G5 g# X   
        gen pmaxphddy240=`tempmaxphddy240' in 1
- b( w' n, R2 A       
4 f) L& h0 e" s* M# y% y7 o0 k       
        *-------------ddy240pmaxp右侧低峰值 -------------------
       
         preserve
  
   keep if dx240>`tempmaxp'
   
9 g+ t# `7 i$ H: R1 \0 t   gen temddy240=-ddy240
   
7 g$ @1 g# Y5 c' X. B6 q( H2 ]   keep if temddy240>0
; J& n  K  F" _   
   count
               
                if r(N)>2 {
+ j0 P0 h/ m- {: E! O" h                 
                pkexamine dx240 temddy240
   
   local tempmaxplddy240=r(tomc)
   restore
0 V2 N: |# u- q+ A2 e                 }
                 
! u9 r2 m7 |; C6 h        else {
3 f2 E; g" Y# L                restore
                 di "cannot find pmaxlhddy240"
                 
               
                           sum dx240
4 L* u7 d& C' Q& f/ y                 local tempmaxplddy240=r(mean) in 1
* q4 W  F/ g% O3 W+ ^# q               
- {- q" n4 L; A5 ~) J                 
                 }
   
   
   
9 e2 R4 X1 `' l! ]9 w' y& e   
, ~  u- G: [; O8 {   
% f7 r  @0 l: w2 W2 D. j8 S' s* N   
        gen pmaxplddy240=`tempmaxphddy240' in 1
* n7 b" z4 c; Z       
        *-------------------------------------
2 C( z% C. |4 r0 }5 ]6 t   
1 O) h' ^6 o; X& y& X" E   gen price5d=v5 in 1/20
4 }- K$ B; T9 D      
   gen price5yy=_n in 1/20
$ g2 j8 Y9 ~0 v1 y+ P   
   gen price5y=price5yy*(dy240max-dy240min)/20
   
   *price5y &price5d
5 N2 C5 F5 l% Y4 Z   
, E) B* ^, k3 G& h. l1 p4 b   gen priceny=price5y in 1
* x1 O  f8 e$ W. F) s  @5 a* J   
   gen pricen=v5 in 1
7 e$ J/ g0 D7 _. G6 _
0 Z5 S0 c# X- Q$ e  K# L* s     
6 j% q) v2 s6 ?$ Q   *priceny &pricen
   
  k  U& m1 m: }   
2 u" H9 Q  c" p3 p( {1 f   *-------------------------graph------------------
7 _1 @: O. ?+ W! s( z  
  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)

* ---------可选项目----------------------
; {) r' P$ R9 x; V& @*(spike yxis amean,lc(blue))(scatter yxis2 amean, lc(blue) msymbol(+) mlabel(amean))
- h) j) m  U, m$ N/ U3 i4 {* b8 o$ S}
8 U" R. h& a8 @0 R: P" Fend