爱吱声

标题: 一个用stata 计算股票成本的小程序 [打印本页]

作者: awer    时间: 2011-8-13 18:03
标题: 一个用stata 计算股票成本的小程序
本帖最后由 awer 于 2011-8-13 20:49 编辑 ' N; M  ?- P+ @- V" M/ n

( G1 g, g0 I- ]7 f
   我是从西西河一直追虎大追到爱坛来的,虎大曾用stata 分析过股票的长期成本和一年的短期成本,但是直接用kdensity命令无法得出筹码峰的具体数值,因此我写了一个小程序来计算股票的筹码峰值、长期平均成本和短期平均成本。下面就给出具体实现方法。
计算之前需要有 :
1) 招商证券系统(通达信其它客户端也行,功能都差不多)
2) Stata软件 (我用的版本是Stata/SE 11.0)
然后
1.新建一个文件夹储存股票数据,比如:
e:\stock\source
2. 股票数据导出
a 进入招商证券, 系统-盘后数据下载,更新数据
b快捷键“34”进行数据导出, 在弹出窗口选择
) J+ u$ i- b" J; ~“
高级导出”,
在高级导出窗口选择 “日线导出”,
将导出目录选择为你新建的储存股票数据文件夹的目录(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 0021288 H. R! K; y* `" q* `5 ]3 G
0 V; u: K4 d8 W1 l8 D
(002128
为你想查询的任意一只股票的代码)
几秒钟之后应该会出现下面一张图

& Q; V3 [. T1 d) a7 w
hello.do (5.25 KB, 下载次数: 14073) 其中蓝线为股票上市交易以来的平均成本值(注意!此值没有去掉一级市场发行和国家的0成本)
红线为一年内交易的平均成本值
黑实线和虚线为一年平均成本值附近两个筹码峰的值
黑方框代表的事20内每天股票的收盘价格
红方框代表股票最新收盘价
山寨完毕, 收工!!

作者: 明月回春    时间: 2011-8-13 20:17
本帖最后由 明月回春 于 2011-8-13 20:19 编辑 ' [. I' |/ G8 m/ ?8 f
; z  R4 U# k, i6 \) t
程序写的不错.( g* N2 @+ G. `7 L( Y" @! D3 w
不过有一点要注意.通达信导出的数据,只对价格复权,并没有对成交量复权.所以数据需要手工处理.虎大买的是Thomas One Bank.
! U: z1 ^; {4 M/ J还有就是虎大除了用了kde,还有摘要功能,你可以看看这个也是山寨" }2 I4 D0 ^( O, O) H: U

- U6 U7 a) ?9 U
作者: awer    时间: 2011-8-13 20:20
这点我也注意到了,所以我用 成交额/复权后的价格  作为成交量 用于权重计算
作者: wenty    时间: 2011-10-19 14:32
多谢,挺好用的。
- Y. P! C4 Z7 i9 k, o; e2 ]5 {通达信似乎无成交额一栏。导出的txt文件要删掉header,否则stata当数据为string了。
) b5 T# y- y! F5 i* Y) `
作者: 旅途愉快    时间: 2011-10-19 15:40
本帖最后由 旅途愉快 于 2011-10-19 17:44 编辑
& a; T5 `" p0 ?* J6 F7 H$ E* K7 m" U4 ~9 D) c
没太看懂。。。' v/ G* O8 }- i1 N
6 ^. x# i9 A" e) X  R0 h) L. `
露天煤业最近20天的收盘价在20块左右?
作者: 旅途愉快    时间: 2011-10-19 17:47
仔细看了看,发现原来自己理解错了~:)9 D' M7 R0 e9 v: Y& }7 r% m1 C! a* m: z
) p' ~, {9 O' f# I9 V( D
谢谢awer。有问题再请教哈~~
作者: 旅途愉快    时间: 2011-10-19 18:09
还貌似有一个问题啊,为什么我处理2000年以前上市的股票,都会出现:

  1. ! m  e( V4 I1 P5 p6 X; i
  2. hello 000425
    8 v6 f% w! P% y# g8 B
  3. time variable takes on negative values# I) o  a: t4 u% u
  4. r(459);$ r0 {* V* E( W' D, S
  5. % Q# M+ v2 ^) p% k  f- E
  6. 3 y; n2 N2 L: L) r5 J# v8 |( g
  7. . hello 000783
    + @3 b5 K+ a. F( |
  8. time variable takes on negative values
      @" Y* C& S# B& {8 k5 ~
  9. r(459);
    ( K' g8 A8 ?8 h
复制代码
2000年之后的没问题,比如露天煤业。
作者: 冷眼看世界    时间: 2011-10-19 19:25
不错,我也来试用一下
作者: 沉默的前哨    时间: 2011-10-19 20:45
为什么老是提示:type mismatch
作者: 沉默的前哨    时间: 2011-10-19 20:49
能不能帮忙看看:
$ d8 n/ ~( \' K# }% T. _8 M6 f" Q6 v
do hello
2 z" s! t: ^, T0 s3 S9 X' J. f. ~! u6 w; J$ u3 `. t6 }
. , W6 w2 ~" B8 Z1 s2 d
. capture program drop hello
1 e( E' S9 ]. `) E) q9 k& n. U3 K9 }/ @1 v: \
.
8 }% ^4 p, S$ T, i" R  d. set mem 1024m
7 _7 `) N" ]7 Mno; data in memory would be lost. w: z" a/ ]. e, b
r(4);3 q2 E, C8 J# O. J: f8 ^0 K( Q0 B

& t9 A& `7 |- ?& send of do-file
" ~: [1 p7 o4 J9 L: Z# x! T! ?2 w0 m! i6 _7 {7 ]2 A
r(4);
1 K7 v$ z+ r* [8 Z2 l/ T; i5 u8 s4 ], c3 M! T
. hello 600063
! o0 C2 ]6 X& q! |6 c! K2 Qunrecognized command:  hello
5 X( C- Q# C5 `6 s7 q4 h$ Tr(199);: c" F  v+ u8 ?& \

作者: awer    时间: 2011-10-19 22:03
沉默的前哨 发表于 2011-10-19 20:49
# C( o$ P7 H- b能不能帮忙看看:+ ^0 B0 R% x6 l) a8 \4 u
) f& R% E8 B# \3 j
do hello

+ O/ D; G% ^( k8 X9 i" T& {你算得哪个股票?
作者: 阿雷    时间: 2011-10-20 08:52
awer 发表于 2011-10-19 22:03
' [7 M: p* C: P. s3 M) J你算得哪个股票?
; Z" b: k" i( f
那现在露天煤业可是在平均成本以下啊
作者: 沉默的前哨    时间: 2011-10-20 10:30
awer 发表于 2011-10-19 22:03 ) j" w' `. V. d6 O/ k$ ?
你算得哪个股票?
+ f- _1 y! q$ \( H' L$ S
我是用通达信引出的数据,比如600663,在e:\stock\source 里的文件是以SH600663.txt存在,无法执行,后改成600663.txt还是不行,请指正!
作者: 沉默的前哨    时间: 2011-10-20 10:32
awer 发表于 2011-10-19 22:03 , S# \9 r; h7 ?+ P. r
你算得哪个股票?

% ?3 z( C9 L4 S8 j5 ^* b8 l5 t/ kdo hello 后变成8 T& D6 R3 `8 M. }* j, C
......  J4 E  o/ U% E; r5 o6 u

# |( y0 d- z' b# R2 G9 z3 a8 Yend
2 l9 e6 E, W- ~. H" j; f
3 f- {! y9 \# a- i.
4 C# ]/ e9 p1 W& I/ Q- K, j& `end of do-file
5 }0 p8 F. j9 p
作者: awer    时间: 2011-10-20 11:18
沉默的前哨 发表于 2011-10-20 10:32 * T0 H  U: W8 Z' c
do hello 后变成
7 q0 u. O6 n3 |, V5 K" N7 e1 g......

- ?: O1 D' Y* A+ M. e" C& w; c
作者: awer    时间: 2011-10-20 13:32
沉默的前哨 发表于 2011-10-19 20:45 ! [& E6 v2 a1 T$ p9 ]+ `1 {
为什么老是提示:type mismatch

/ z( o; x- d$ S; C0 G" r  M注意 日期格式要选成  yyyy/mm/dd,要不然时间变量格式不对
作者: wenty    时间: 2011-10-20 14:16
"此值没有去掉一级市场发行和国家的0成本" 是指已包括一级市场发行和国家的0成本吗?
  ~, M/ V) F: B& C谢谢
0 _( j2 d. ~+ X; U$ t) k% L+ T
作者: 琴心三叠    时间: 2011-10-20 17:35
下方的棕色曲线是代表什么的
作者: 琴心三叠    时间: 2011-10-20 17:35
下方的棕色曲线是代表什么的
作者: maqiangcn    时间: 2011-10-20 19:13
明月回春 发表于 2011-8-13 20:17 ( }4 G# ]+ X6 {7 T
程序写的不错.
$ |3 |) R4 _0 J" i不过有一点要注意.通达信导出的数据,只对价格复权,并没有对成交量复权.所以数据需要手工处理 ...

7 h" R8 f% y. ~% W7 P" Z你好,这个链接好像不对,能确认一下吗?
作者: 琴心三叠    时间: 2011-10-20 20:03
琴心三叠 发表于 2011-10-20 17:35 % M+ c' X" M1 X
下方的棕色曲线是代表什么的
' Q2 i3 }0 S$ P6 L0 c) r* I" T
两条曲线只是周期的不同吗
作者: 琴心三叠    时间: 2011-10-20 21:37
谢谢答复  还有一个问题  我导出的文件为何只有当天的数据  而不是历史全部数据呢
作者: 偶然    时间: 2011-10-20 22:19
建议楼主给个数据格式模板
作者: awer    时间: 2011-10-21 10:05
琴心三叠 发表于 2011-10-20 21:37 , W, t9 Y; Z; G" `3 Z  N
谢谢答复  还有一个问题  我导出的文件为何只有当天的数据  而不是历史全部数据呢 ...
& d4 ~5 L/ A2 r, y: K- w
有很多问题是数据不全造成的
  K) i: l' O  J. l如果想要全部数据,
% F3 d: E; T8 x( d* m% Z1 W# i系统--盘后数据下载--日线数据--日线和实时行情数据
9 E1 ^' a0 E; w* i( I: S+ a
3 S+ d% T8 J9 R- h把开始时间设成中国股市开盘的那天,下载全部数据
) {. A- u( |/ b5 B8 b' e4 ^# f
  `; L. O; z& U' O下载完全部数据后,在屏幕上使用快捷键 34
* K% }" I& M5 b) |9 U9 e( H$ @
3 w9 b9 Z2 B, p9 P. b( s" L34--高级导出--
& o7 e5 ]& w4 e9 j8 N& W$ `
3 n( X7 P' d2 U* [, I$ ^' _然后选中   日线导出;选定导出目录;文件名设为:XXXXXX.TXT;勾选精确复权;不要生成导出头部;分隔格式为tab();日期格式为YYYMMDD,然后选择添加品种
% l- I" o* N; H  r6 {; X7 [' |
+ B: u$ X/ C5 n1 V0 ^+ y我一般都是把所有数据下载到指定目录下,然后进行hello运算 / g1 i$ |3 r; L

" t% v) J3 i, y7 s
作者: 沉默的前哨    时间: 2011-10-21 12:40
awer 发表于 2011-10-20 11:18
1 o" o( y9 x" ]" P, X3 N) F" g7 F( `
谢了!我自己还是没搞定,我再看看。
作者: 沉默的前哨    时间: 2011-10-22 12:31
搞定了,谢谢!
作者: warbrai    时间: 2011-10-22 16:38
沉默的前哨 发表于 2011-10-21 12:40 8 @7 u6 D! a$ x
谢了!我自己还是没搞定,我再看看。

# [5 \' i. D) f, i# d% D- K别选择 "导出数据头部"
作者: 沉默的前哨    时间: 2011-10-22 23:16
warbrai 发表于 2011-10-22 16:38
" E7 a. s( k* l. ?4 Q* c# c别选择 "导出数据头部"
6 a! c8 R% E, h
是的,谢谢了!现在可以了!
作者: 心静天地宽    时间: 2011-10-23 20:45
老早就想找个这样的软件和方法了,无奈见识有限,一直没有办法,感谢你的无私。
作者: 琴心三叠    时间: 2011-10-23 22:04
抱歉  我还是没明白棕色曲线的意思   黑色线代表筹码分布密度  那棕色线你说代表随价格的变化率   这我不太明白  请说详细些  不好意思  我比较愚钝
作者: 琴心三叠    时间: 2011-10-24 18:45
hello 600085" i- _" O1 n2 d% L8 d3 X, b
unable to change to E:\fenxi
! Q5 `7 p, t% V7 s! R这是什么意思?
作者: 观风望月    时间: 2011-10-24 20:32
已经下载下来用了,非常感谢热心的分享。
9 S2 f0 C: I, a8 ?2 z因为不会使用stata,所以有个小问题想请教一下,我按照你教导的方法已经成功测试了露天煤业,兴业银行等几个股票,都成功出图,但是在计算600188.601999的时候,程序提示“
2 l. Q6 k; _6 a% _. hello 6019998 o9 W1 Z  d* `3 `3 C  e4 e" a8 d
statement out of context" Q: c( Z+ J6 g; g
r(119);
2 V. p$ Z8 _- D8 N( y+ ]; ?! b1 K
/ {) H3 |7 n  b: J) i我查了一下Help的内容,对于这个错误也没有特别具体的解释,不知道是什么原因造成的,因为数据的导出和格式和前面几个股票都是一样的,所以我搞不太清楚到底是怎么回事。$ Q* ^& l8 _: _+ w# W( j% @- A1 H2 F
+ |0 O. M( n2 {
另外,虎大提到过的国家成本和一级市场成本加权可以通过在数据中添加一行交易数据来实现,只不过这笔交易的成交价格和成交量需要自己加权计算,对于不分红的股票来说加权计算还是比较容易实现的
作者: 本拉风    时间: 2011-10-25 15:22
横纵坐标分别表示什么,2条曲线又表示什么
作者: 钱潮滚滚    时间: 2011-10-25 16:08
发个图验证一下和大伙的结果有没有出入./ o. }; A3 e# }7 A) {, N
( M8 g+ N# ~3 c7 _1 J+ L; i5 m  b
收盘价到昨天为止的三一重工:. q  m0 U1 V% P! K1 T. @! N

  V; A' }7 g  N/ ?/ W5 ~& y! _0 m$ @

; X9 W/ i  l. I. s
作者: kevenjames    时间: 2011-10-26 14:05
就楼上这个图,几条线,还有黑点分别都是什么含义啊?
作者: wenty    时间: 2011-10-27 04:17
总结一下使用心得
! D0 s4 t6 \# z( f, [6 _1)删掉set memory那一行7 A4 E0 t0 r9 I* m
2)改路径
; Q4 i% \7 B; O/ C  i9 b3)从软件导出数据时,要先下载,后导出,不选导出头部。
$ V& b( W& I% i4 \4)用notepad看导出的数据每列对齐没,用tab修正。, a5 l/ v( s7 m! R. \7 S* M
希望有帮助。
作者: 沉默的前哨    时间: 2011-10-27 09:09
老祖爱蓝军 发表于 2011-10-26 22:50
' l* e& P7 P7 H" s. _我出现同样问题,请问怎么搞定的?

% g9 g5 e& l( N关键点:! U2 A6 M8 o% k
1、文件名设为:XXXXXX.TXT;7 o1 y" v* O+ t% x! u& o/ I
2、勾选精确复权;
3 h9 }+ j  {3 \7 _' J/ P' y( \3、不要生成导出头部;
$ t8 F1 [5 S& n% o7 c* ]4、日期格式为YYYMMDD;* m/ n& L1 D; v: R
5、如果不行,建议可以关闭stata再打开试试,我认为和操作系统的版本可能也有关;
作者: 范进中举    时间: 2011-11-8 16:10
不错,按照楼主操作完,最后是选择do命定,选hello文件,然后hello 代码就行。6 f9 ]+ Q5 p% h1 R& V+ n% Y+ _
几个问题
$ B( d3 S8 H: |" b+ k2 u
$ Y- A3 Z; y9 E/ Q: U' o, y最边上的蓝色的4.473554是什么意思,红点13.84054是最近20天平均价?灰点15.3842,黄点19.09930是什么意思?
作者: 范进中举    时间: 2011-11-8 16:12
钱潮滚滚 发表于 2011-10-25 16:08 ' ]  Y! b0 M: P4 ?3 q
发个图验证一下和大伙的结果有没有出入.
' q! u8 }" t: h( T5 B  }8 O& P3 h, u% {  o, u7 w- X0 h
收盘价到昨天为止的三一重工:
, F( l  ^# ?! c. z9 A
蓝色线上的8.520041是啥?红点13.85195是啥?灰点15.40402是啥?黄点19.09793是啥?
作者: maodeng    时间: 2011-11-9 21:47

9 m/ h: \' {+ _9 X! m# Q7 K! u- r3 l7 n
7 ]' S) R: G5 u) i3 s# \+ i大牛,小生这厢有礼了,啥也不懂,照着葫芦画瓢,我的数据应该是更新到11月9号的,图形如下,不知道对不对
- ~( V$ B6 F2 X0 p; B( H8 c9 B; r4 ~8 g
$ W9 n* o! k( C6 O0 n1 ]: d' l$ I/ e$ D  z

作者: maodeng    时间: 2011-11-10 19:58
wenty 发表于 2011-10-27 04:17
# V. g) h7 z9 E4 H, ~2 |总结一下使用心得
& b" O5 S/ h0 k% G5 j* }1)删掉set memory那一行  J; U4 I2 F& K, @0 n' y
2)改路径

' L) P3 y: B! n/ \9 w  J, Q记事本打开后,怎么用tab
作者: 冰火    时间: 2011-11-10 22:44
认真学习
作者: higer    时间: 2011-12-9 22:34
我也曾尝试过使用stata软件对股票平均成本进行分析,但是作为一个非专业的炒股人员,而且日常还有工作的人而言,这种方法似乎显得太专业,而且操作难度极大,不一定实用。
$ n9 d$ K! ?5 ~! A5 N1 ~
( w- ~0 C2 Z1 d7 A7 u使用这种专业的技术分析方法容易造成这么一种现象:大家都对导出的光滑漂亮的图像感觉良好,然而在股市中依然无法进退自如。
作者: jway315    时间: 2012-3-24 20:12
谢谢楼主,希望继续给同学们挖宝哈!
作者: awer    时间: 2012-3-28 18:06
改了一下,这次运算应该没啥问题了
  N1 F# n$ s9 ?( f, U: T% n3 ?+ z! `7 S. S1 y2 g. q% h  t

$ e8 K0 \: X8 i1 E9 _9 @, E! Ocapture program drop hello( R9 Q+ _( P% x
$ ~$ l% A% }2 H( x" I
program hello
4 K0 H8 e% R' c
2 q) c& u6 }& ], a0 O& l& p quietly {; T  D& }' w+ B- K' h$ X
$ X& C6 v+ d  v& B# @/ p8 |/ A
clear
5 D+ a2 j" I, X. F6 e# e8 y
3 Q! n7 T, Z6 V* Kcd e:\finance\stock\20110101\source9 \  R! p" @4 B0 `1 S4 ?+ M8 C
" J3 p4 C1 _9 B8 [4 `% J2 [
insheet using `1'.txt
3 h5 Q4 B' s7 t5 r2 ]9 Z: O
3 s% y. r' V, ~$ B' [   keep if v5>0
) u' S& F7 v1 t    : X, o4 B. l+ t0 Z9 J/ N5 g
   gen vt=_n5 _" r& Y- A5 p, G1 s9 V2 C8 }
   
  h) W3 P; k" N! y8 B   tsset vt
% H) e$ A" ~* {$ P* x
& W7 a5 ^6 S/ j6 M- N( K1 D9 L   gen zf=D1.v5# d4 V  q# o: M5 ^8 t

- H" B( B+ K  e4 w   gsort -v15 K: h8 A0 R( I6 n
; v% Z3 }' H8 o
   gen avol=v7/v5
# F- Y* G$ v) Z   
  n" n9 ^8 e# I   keep if avol>00 `* \# F, I0 O. k8 ~
   
' ~2 \9 E( {1 M. A7 j. W# ]; j   drop vt/ |+ g' h7 s7 h, a5 c$ u- L7 X
  9 Z& i7 M) j1 @9 _9 d, Q
   gen vt=_n5 U+ z4 k( D# K- j$ F% z: z  J6 @4 Y
   8 M1 G4 X' x: S3 S. N- k0 w+ {5 h
   tsset vt* A$ C1 l# m* W/ v! d
   
3 _! Z. `  X% |8 f! Z3 S; r' z& |    sum v5 [aw =avol], G# m( B0 e' P( g! I; G7 c& J1 E
4 B- M" Q# H) R  r
   gen amean=r(mean) in 1% P/ {9 C3 f% H" B" X
   * M" o- n/ Y; B( a$ h
   * 总加权平均值 amean
3 q! I9 r3 i. J, O0 V9 \   8 [( Z# F+ T, A' D: J; O' M- S
   1 O5 P/ B- i" o' h; L5 i2 I! s% J9 Y
   keep in 1/240
  z2 Z. l5 R0 r) [3 q1 T- M9 _" X; {7 {  L; h
   sum v5 [aw =avol]
3 C: `% J# G2 e$ p. R9 j; R# h: P, ~8 j  [6 ^% D! P  A
   gen mean240=r(mean) in 16 a4 I- u4 _. u( I
   
5 U0 u3 G/ p+ Y* B" g4 f   *240加权平均值 mean240$ X4 N3 j3 A& U* e, v& Y% a& J
   ' U" e. v" L' l  ?0 r$ `" [
   set obs 500 6 _7 T8 @! s# k6 x& I" H9 m
   ; h! n  \9 r; f6 F+ w& H" R3 q1 }
   kdensity v5 [aw =avol],gaussian n(500) generate(dx240 dy240) nograph5 w# l  k& i  E& G1 W2 A; N
   $ q. I' d' G, w
   *240 kdensity(dx240 dy240)
. Z; |8 _3 T" V. _   
  |' d+ e0 l/ E+ _3 H0 W9 R   sum dy240; ~/ ^7 P  {: |5 G0 o* _/ P
   / U4 L6 n9 T- Q
   gen yxis=r(max) in 1
1 |% Y6 u, I5 W, `   
% m7 r0 B/ z0 g  [/ ?5 \   gen yxis2=yxis/2
7 U3 ~( e$ f, e+ {3 f  P   
0 Q' V9 N' ?% o' X   *yxis( Z+ A  Y; y: ~6 j$ q3 s5 @
   # t. Z  [( H3 B: e2 G9 V9 a
   gen dy240max=r(max). b* e" V1 |# i- I3 `+ R4 b1 ?
   5 n. x: g: q/ l, K$ t6 h) U. U2 P* ]
   gen dy240minp=r(min)
/ |; V. t/ [+ O' E7 ]3 Q2 \5 y' ~ / q$ s7 J* `' ~; L7 D1 z# Z7 W; S, G
    gen vtt=_n
' R) N7 f/ C+ ]  H
# |$ v# G! k, m; e: a    tsset vtt
. G( E5 P1 W0 s4 d1 X
4 m6 u5 s; d) W  c; g   gen ddy240=D1.dy240/D1.dx240
6 U) p9 r+ W/ [  c0 m' n) U   $ A) c9 V0 R& A, y6 s  y: V
   sum ddy240) O4 W: W, L- j7 z' P$ H
   
9 |& s; }# U0 X6 D9 x   gen ddy240max=r(max)7 [2 w1 p1 f4 Q% |: f
   + j% N  H. D6 ^+ u5 c7 _0 D6 M
   gen ddy240min=r(min)
/ y* T* J% V4 ?/ z+ @" w! X* M2 v- s3 e     
$ d4 W7 U& S- d4 z   gen ddy240g=ddy240*(dy240max-dy240min)/(ddy240max-ddy240min)* C: P8 ^4 M8 E: c/ v
   
. g3 w( F4 B* d- q, H* z9 E" B) S   sum ddy240g3 c: {8 j% V/ e7 \4 m: Q' I+ q
   & D$ v( ?& C, J& \& N
   gen yxis3=r(min) in 14 k) u8 M; W5 e9 t" X4 L6 P
   
" y* m0 F' A, F* f" ~2 `   gen yxis4=0 in 1# g$ N  n  y2 ]) Y2 n) r" r: x
  
* F( g& c4 b+ }, q& T. }( o3 Q! o  *ddy240g 4 X* T, K9 P1 N( @
   
0 `+ v! r) S/ @   
2 C+ D6 w% J9 d$ J0 K( x3 L  preserve
& M: i* y# c9 c/ n+ X  , q2 R, h5 a5 C$ `. T
  sum v5 [aw =avol], L: T# o' E9 K( X( E; n
  & v" U$ Y8 H! C" b0 e
  keep if dx240<r(mean). r) a; V/ n8 T, r4 D
  * b6 D2 Z  c( p1 Z" @
  sum ddy240
/ H. k; a( A( l( Q  8 m2 w6 {: B4 O* J. v/ W
  if r(min)<0 {' L: [1 J. a% q5 p' g
        / W0 h. n( R( [, Z0 Z7 }
                keep if ddy240<=0
! x2 s6 g. I6 p               
( u6 n9 i% a, V8 P                count1 ?- v, w7 K& X( h) ]/ ?
               
1 }0 t. g2 O2 ?( Q                if r(N)>2 {& T6 k8 O/ E8 J, z
                 
8 u3 i5 O4 v$ v1 E8 o                 pkexamine dx240 dy240
2 K/ j( l  O. O* ^# G* L               
5 ~0 o% j5 q9 ^                local tempminp=r(tomc)
% }' v+ g* o7 z9 O" `, u1 @                restore
" |3 ]  c9 {% U  Y% G                 }) w- D! O- A3 f; g1 k' I) @
                 * H6 c3 V3 Y6 A( d2 m& w
        else {
3 w+ T; j2 T' \                $ @: m5 B0 v) t+ \! d2 _  S- O9 n
                 di "cannot find minp"0 J  ?0 i5 ?" K" b9 `& r( N2 J7 Y7 ~- e
                 
" B! L3 I7 i$ o5 S                restore * v, h# m. D8 U) @7 z( _
                 sum dx2407 P  S7 o( ~1 D- Q3 [
                 local tempminp=r(mean)
. t8 E- \) q) W$ E                 
+ O9 n) p% R# V# n" B  b8 g                 }
' t4 v9 p. ~9 s0 G                }
+ f7 m0 u" b! o6 c1 E   else {
, Q4 L& Z3 H& w2 j; m/ C                 
; O% q5 O% Z8 }* v! c                 keep if ddy240<=r(min)
: f* T: P; Y7 y) q5 N/ T/ i                 : U# p2 U5 u' V  f
                 sum dx240
0 Z) L/ n( ~0 D& m0 S7 H7 H5 k3 y                 
- ]/ t  F* h. {$ Z9 A                 local tempminp=r(min)
; e- O4 I4 H7 B  Z; ~" @2 d                 restore ; H+ Q( X0 I: B0 {8 E+ B  i. W
                 }
4 P7 o2 T% W% K" A3 Y% v8 m8 p' Q0 L( p3 X& h7 R

% E5 Z5 ~' \0 }/ g( x. Y7 g5 Y2 z( ?- b" x
   gen minp=`tempminp' in 1
" Z- A; `  c- J; I   *240 minp! C$ U2 q2 }3 |1 Q
; e; X7 r4 F$ T6 R5 I7 `; W/ }

5 t# N+ V4 A8 P* U   
9 H* M9 u8 R8 }0 ?   preserve
: s4 B7 p1 V: G9 [7 u) Q   / |/ M/ B, `, o" u) j
   sum v5 [aw =avol]
* j9 U( N. T1 ^7 L2 Z: [   
* o( V2 e3 ~4 b2 T   keep if dx240>r(mean)
, p1 K( l: E9 }   7 t2 x  [" c1 r  v2 q( S! W! z9 {
   sum ddy240
' s; Y6 d6 |+ p7 i+ G9 M   
9 ~0 g7 v9 J% [7 Z   if r(max)>0 {% w  M+ j% H4 U! T3 ~
               
& T" _, k) f# n                keep if ddy240>=0
9 O, J3 K' _0 D) G6 j               
, l5 H& \" d5 J; y                count* d: F: D6 l+ ^( y  k
                5 U+ u0 N1 K+ H* C" B! Z: n: B
                if r(N)>2 {
7 n' Q7 B, U) k  `: e% a5 k  B                 
, [( M6 p# i# E: _/ q" L                 pkexamine dx240 dy240
$ R  P+ m  Y3 y: n$ W! e- a: c               
; Q: c. e) J$ }                local tempmaxp=r(tomc)   Z" t. x, m% g; c/ B; b+ t6 b
                restore ) z& {+ [+ m8 ~" Y5 y# s) N
                 }" p* Z4 S, Z- z9 O& k8 H( K
                 . C3 x1 X, c: q8 t
                 
2 W. ]6 z* R; k7 V2 p) o5 _        else {% Z, D& n) d7 a9 P6 d6 G
                restore
9 U( B0 T1 S, J. t* U% j0 D* c                 di "cannot find maxp"2 S: W: g% _, U# U; }9 J" D5 J& a( n) ]
               
! h% Y/ m  ^3 Y9 X$ c9 M) J                  S% \& I# P; J
                 sum dx240; V1 W/ \% `, K+ r" e8 m
                 local tempmaxp=r(mean)  N: X& ^( y' k; j& w' I
                 4 R& L$ ]) B. W: z, ~2 y! A
                 }* G# y% E) c0 L
                 + v3 U6 P# p5 N4 m. u
                 }
$ I6 _7 T: d6 r+ g                 . I  |$ C( d* W3 l+ K: q
    else {
7 k6 r/ W% _. ~# G  y6 e        $ r, R  [3 G' u& L
           keep if ddy240>=r(max)
- ~6 {$ V" e8 m1 B, ~; k          
$ M" a" m3 v$ S) H           sum dx240/ Q- x' N1 y# z5 H) U* s
          
; Y' v; }7 q& J: }7 V! w           local tempmaxp=r(max) 3 W: m0 r4 @! \+ D# U, T3 N
           restore 4 N; a( |6 t9 l/ t7 c' ?/ k! l
           }; W* x$ H, F* e4 Z& H
        0 X; {0 b: D9 J2 K5 B
! a/ [1 M) |" f
  1 O8 C1 _- P' p3 e# |9 e2 M$ |
6 Y) \8 o1 I0 J3 p9 I; O
   gen maxp=`tempmaxp' in 1
' C+ }5 d- `1 Q   ( F$ D" O7 v! d
   *240maxp
( \& b3 }& X6 C   
2 F$ Z: y- H+ B% w  U. t   
2 o: I( T" R( j- H   *-------------ddy240pminp左侧高峰值-----
3 Z) H( _# S" Z   
* k" x3 f# s& V2 Z5 q0 w    preserve / e8 Z4 _3 Q2 `; C; }- P8 m
   
) T' y! g3 R" o" d( v( J' \   keep if dx240<`tempminp'
* g1 L& n. m/ Z8 @$ K   " ]. }; H, r! T& y
   keep if ddy240>0
: A- t( K3 Z3 \" T& V6 o   
% r1 m$ h( k1 p   count
( v# @( X% W* R4 v$ `) F7 X8 T# N               
+ m; i* j+ E7 s" e! p2 W" [$ g                if r(N)>2 {
& l$ n! \6 p- v7 P                 
, [" F8 s! G* ^3 L8 H                pkexamine dx240 ddy240( [) A( ]$ D$ \2 W* i7 }0 n: Y
               
4 `% N- C+ F* \: Q                local tempminphddy240=r(tomc)
( V3 S% x0 K" e7 v3 c  w5 f                restore
1 W" {- l4 A5 m& Z- S- j. C                 }
' w' ~/ z) j% [. F3 d9 o; H! a5 W                 2 R& Y  M8 i8 S/ ]
        else {
# Q1 A/ I' K- x3 N- e  p                restore / z  l  @6 g: n3 n' ]$ Q% h
                 di "cannot find pminphddy240"
4 g+ Z/ S4 {" |  X6 \8 d                 : C5 y" }; O' Y7 s
                8 {" A* k; \5 ~! G+ E
                 sum dx240; l- O1 w% i: h
                 local tempminphddy240=r(mean)
1 u. N1 I$ u* h" t               
/ [8 u: J7 N2 `# ?# W6 K2 E                 
, e$ e& Z$ s, o6 K7 G                 }
0 m2 x+ _, s2 L* t' J2 U' p- S                 , Z- W  W0 m8 \6 h' M) H) J8 k
     i4 K) l! u* r6 U2 w) `% N
    7 }8 @3 X8 _) \
   
: ]! C, }! ^' S5 x* ]5 Q6 e- |8 ]) q   
4 d  q) {! K! S   7 W6 u4 f: A* R
        gen pminphddy240=`tempminphddy240' in 16 v  \8 e: a% \( L  U% [/ N
       
% X) j9 A$ D% Q# x) N0 y# s+ _0 M  X       
  P% K" \6 z3 U! N        *-------------ddy240pminp右侧低峰值 -----------------
2 l# k% A3 R! O% {% l9 b    preserve
0 V: o. w+ w3 {) h( T' s, j
6 B9 d+ X8 `+ \   keep if dx240>`tempminp'
5 y: H9 y$ C8 `% I7 B, A" g   
) Y& x: ]1 D  b( L4 m   keep if dx240<`tempmaxp'
" X- _) `5 S6 c; k$ a   9 d; K& d! _7 I# r3 a
   keep if ddy240<0
; x1 m( O: z! i; u0 s$ i   " q, b. |6 z0 K) s# O7 _, {3 x
   gen temddy240=-ddy240) `, ^4 }3 i/ y6 j
     C4 \5 ]7 u- B1 I
    count
- [, e. e% H$ O' q: R               
3 s: ]* X- [  q2 x8 D# T# A' D1 V8 g                if r(N)>2 {
3 T7 @1 v$ Y# l6 |2 X                 
& \; {7 P' w6 N5 |2 B6 K3 L; S                 pkexamine dx240 temddy2400 D( }; ]' z( u+ _4 w% N0 n+ V
   * ?7 z" q" J% y0 S1 h& y1 k
  local tempminplddy240=r(tomc)
2 @7 d/ d+ ^9 Z8 f6 _                restore. {+ S' _9 S2 N% K
                 }) o" d% }1 w, T
                 
  H- v! m( ?1 [8 N& U        else {) c- B  `0 D0 k: {4 Z7 s
                / N) u7 S/ `7 Q1 t2 O, H0 t2 D2 C, a
                 di "cannot find pminplddy240"
9 h' I1 {4 l' a0 x                restore
/ s& @8 Q7 \0 d/ K                 : `% b( R, _- Q4 z9 }, x
                  sum dx240
. g. G2 d0 H4 L, }3 Y2 H                 local tempminplddy240=r(mean)
: ?9 N" R1 z+ D6 `2 }                 
; \0 G$ r7 q$ _4 i                 
  t2 z: L. M3 ^: H2 @                 }
$ _  e4 b( W2 X% t8 z1 @   ; L* o# f* [! P9 C4 o$ F
   
. L& f5 s+ e! Z8 b   ! ^% U7 i! S7 n- I! z& ?* T
   ' m/ i% f6 y% ^8 C  h
   / I8 ~" ~- e% N
        gen pminplddy240=`tempminplddy240' in 1
2 l  F- e. v  R6 o       
; |/ ?7 Q( i* e/ R3 i  [7 D& _        *-------------ddy240pmaxp左侧高峰值 -------------------  d' J& H1 V$ j# Z
        * J4 K8 M0 e- M( j
         preserve * g' X3 \. K* r4 ]/ l$ |
   4 A6 x0 f3 `4 G, w7 X5 @
   keep if dx240>`tempminp'1 R" K% s) D) ~! \. x
   
# R3 [6 f6 t9 k: c( D$ R   keep if dx240<`tempmaxp'
8 ?$ R9 M' ^" y. ^# e% X) J7 M, k   
" m* O) w& Z& w% e9 s   keep if ddy240>0$ ]+ f' ^1 x5 m% Z" J- d
   count2 g# z' I' E) l% g1 p: l8 I
               
" k: S1 O1 W. a- O. N8 d0 c1 Z                if r(N)>2 {
- ], P9 V9 U6 P' L$ u* x/ j                 
6 U9 ?1 l  O; c                 pkexamine dx240 ddy2405 X4 `+ c2 S& \5 Y2 i  ]$ z1 U# U" {
   - o3 Q! X3 Z) ?. W5 d2 z
   local tempmaxphddy240=r(tomc) 1 \/ I  M, p7 G- ~  X( Y% r8 @
   
; }" N2 j1 p- A" c, R" ^/ {( p' m   restore
% r. A- ^/ V" O" b$ a' s                 }" ]. X4 x. c# A5 B7 h: m
                 
& u! f- {! d2 a  T2 H9 X4 L        else {+ R& X7 r2 }  N- F
                7 r- a: \3 n3 @
                 di "cannot find pmaxphddy240"+ b* u) w! |0 x: S
                 restore
  N- |, m+ \% I                  e* l, i, x8 k5 X
                 sum dx240
( G6 {  d7 S( p3 {# P( N0 O                 local tempmaxphddy240=r(mean) # g/ \7 |" [$ a( ^; U
                 
5 }4 z6 E; v: R                 / @! f+ n) z* r4 s: ]' ?7 a
                 }# _) R# y! W+ @
   
* Q7 c' g1 l' p8 w) g, `   6 O) O% E& s0 ?% g: W: v) ~
   
* R- C3 T  B" D5 c# }8 S   
2 o! C2 q- J. ?( u   ; B2 d# C) s+ t3 l8 F) G
   0 V% M$ L1 a( J# W/ C
        gen pmaxphddy240=`tempmaxphddy240' in 1& @! P! x* \2 ?4 w$ O3 ~
        ' Y* g$ j. _5 m# X3 |
        * B% V, q4 D. k; p
        *-------------ddy240pmaxp右侧低峰值 -------------------
8 ~) T& y9 x7 J2 a        ( V. r6 t$ A# k
         preserve
# z) [) J1 K+ M- |" p, }: N2 G  4 p/ I: F8 H; t1 i+ b8 H+ m
   keep if dx240>`tempmaxp'
2 {, D5 s  P; F% Y  }   0 A" Y8 X, `) V4 K* S
   gen temddy240=-ddy240; R# e2 V8 f! ~6 p# ~7 N
   
5 n! S# _8 L" q7 j2 d  R7 q   keep if temddy240>02 a0 A' N% w, z2 A( F; ?/ V
   . H2 |0 P: a! j
   count
  F% J4 W" S2 e5 C6 ?                9 G7 N9 E' Z' s* ~! d" ]8 }9 o
                if r(N)>2 {
3 z# _  P3 Z2 d  N                 
4 G) ]* }5 J9 l3 r7 Z5 P$ \; Z( \                pkexamine dx240 temddy240$ F4 n/ N* k, @& X% B
   . p6 L7 |; Y) s: G
   local tempmaxplddy240=r(tomc) / U$ V6 k' x7 ^; u- p' ~
   restore
2 n) }6 [$ j1 Y" y1 }6 n  Q/ B                 }+ M8 X0 a8 I9 K, G4 N0 E, I
                 # ]9 {6 D- G: m/ B+ E: O0 U# v- u
        else {! D3 v% E' M7 G( \0 ]- S8 u
                restore
# [) J+ @; l5 P. T1 Y7 x                 di "cannot find pmaxlhddy240"
' ~  Q1 E  J( Z5 K                 
- a  n. U8 T3 b5 T5 X6 {                8 I# u: \1 u% v$ F2 ?3 {7 ]7 B- M: [
                           sum dx240. Q! h# }4 y9 `; Y8 @* K8 b1 [
                 local tempmaxplddy240=r(mean) in 1
. O4 n' ^  t+ f3 Y7 q               
4 f& x9 |' @( ?3 H0 H0 ^5 A                 
' S1 q+ i* O8 C$ b- B' g                 }6 Y; Y: A8 b7 [% y7 @  j" l
   
( c' L5 W5 I8 s7 M( g: {3 A   $ c; z- L& V% b% S" t3 v* [4 h9 `
    " |4 J" `, X5 N7 c% ^, Y
   % g' t8 d% {% j  W( R7 z. N* {
   
% K0 z5 p' Y  x2 x( m   " [+ W7 h) o9 T6 U
        gen pmaxplddy240=`tempmaxphddy240' in 1# t4 D+ q! H& @8 j0 M
       
# G+ A- |1 ~3 i# i5 r" J1 w        *-------------------------------------1 j  F1 l% ]! j9 R8 Q
   2 T! x6 Y$ e# ^' d
   gen price5d=v5 in 1/207 q0 U, g: b& \- j% y4 t4 q
      * k. o0 ~& o4 m
   gen price5yy=_n in 1/20; m& H5 t4 H+ M( Y* m
   
5 L3 t* A8 o( p& |   gen price5y=price5yy*(dy240max-dy240min)/20
5 s4 S) l( i# g8 h   
3 j- Q1 @" }' y( g   *price5y &price5d
6 \& U! z. s7 a" y: [# w   
1 W  q, V/ @) m# L; z   gen priceny=price5y in 1
" S* V1 q; q$ A3 u   ! r, a0 V0 B8 j
   gen pricen=v5 in 1
- z4 q1 J( F0 F9 L
& L2 G( f. C6 ?- ]# g  N, n     + ]5 I3 _/ y2 F: m/ H/ Z
   *priceny &pricen
0 B* R7 n! j& Q: u2 k# A) h8 S* e   
% ^/ b* A" V, A2 `; |# v9 H- _1 G   2 p5 [* l$ g  x* C; P% U$ @5 D9 c
   *-------------------------graph------------------5 Y7 D9 o( C& Z+ D" z( ?
  
+ ?+ D& m) X( P6 Y" E6 ?( V  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)" s, b, h! f7 O- V+ h7 t* _+ U
) a- C- U& i8 R: @1 q. n" [8 p
* ---------可选项目----------------------
* W2 [- V9 v) j6 A
: h* u2 l4 m1 P5 o* t  Q}
( T% ^8 d* f  g" R2 s& V7 Bend
+ Y$ _1 ?2 _; v; c4 k3 @$ n
作者: 旅途愉快    时间: 2012-3-28 23:23
awer 发表于 2012-3-28 18:06
* i1 T. ?0 _/ K2 J. m改了一下,这次运算应该没啥问题了

/ s( u6 Q- k  P, }8 @8 B, J这段代码是重新复制粘贴到hello.do里面的?
作者: awer    时间: 2012-3-29 11:40
旅途愉快 发表于 2012-3-28 23:23 - u/ X! m% \' F4 w6 ]$ W1 m# ^
这段代码是重新复制粘贴到hello.do里面的?
$ b8 m8 D7 q: m8 G" F* d2 ^8 g
对的,就是hello里面的那个
作者: David09    时间: 2012-7-2 21:35
awer 发表于 2012-3-29 11:40
8 q+ e3 Q# H1 @* G. b对的,就是hello里面的那个

; M# z9 O* N& I' Y. hello 002680
6 J: A( p" W$ \( i9 L9 `Obs. nos. out of range
, K& Q) h( H% F. {3 c% or(198);6 V/ v9 p& p1 G) ~" m3 t+ m

5 L; Z% `. r/ @- v( x' c1 o# n
* ^- u) x' f9 M  V( o5 Q/ G我运行后,出现这个,不晓得是什么错误了,请大拿们帮助
作者: qqysc    时间: 2012-7-3 12:12
David09 发表于 2012-7-2 21:35
( i$ T% ?' ^7 L- |/ F% b9 x, z( }, G. hello 002680
2 ^  L  i$ i0 C+ K# z/ @2 f# FObs. nos. out of range' g/ O- |) |$ M6 j( n
r(198);
8 j' a4 u( R* I
估计是导出文本格式问题。
作者: qqysc    时间: 2012-7-3 12:14
awer 发表于 2012-3-29 11:40 5 l$ o/ v6 i$ Q% `
对的,就是hello里面的那个

! A' |/ Y* w& F! _5 x( b9 R+ J& |复制粘贴有很多乱码,请楼主把上面的代码做成附件以供下载,谢谢!
作者: David09    时间: 2012-7-3 17:13
本帖最后由 David09 于 2012-7-3 17:28 编辑
  W. W( A9 z+ B
qqysc 发表于 2012-7-3 12:14
+ n5 ^) d7 ~4 b7 E% I复制粘贴有很多乱码,请楼主把上面的代码做成附件以供下载,谢谢!

% i1 g: a0 a% h. J
( z8 ]3 i7 x0 p我清除过了乱码的啊....
作者: qqysc    时间: 2012-7-3 17:56
David09 发表于 2012-7-3 17:13
* s" y% {2 p6 l, r我清除过了乱码的啊....
0 k7 k  G- l( j, f" R7 d6 y- w0 H$ V
说得不是一回事。你说的问题可以试试在 002680的文本中,把tab键补齐,导出时不要带表头
作者: David09    时间: 2012-7-3 19:45
qqysc 发表于 2012-7-3 17:56 6 e+ J9 V6 G% v' t' [& G# a, ?4 E. Y
说得不是一回事。你说的问题可以试试在 002680的文本中,把tab键补齐,导出时不要带表头 ...

) Y2 L6 J9 M4 ?! m很奇怪啊,我看stata中间的数据表,都整整齐齐的导入了成交价成交量。但是一运行就蹦出来那个错误
作者: David09    时间: 2012-7-4 14:44
搞定了,原来是我选了一个次新股的数据,样本数目太少了,引起程序的错误。
2 ?6 Z8 U" e' f+ ^7 o8 q感谢楼主给的好东东。
作者: 碧水芳洲    时间: 2012-7-13 13:40
怎么没人问曲线是什么 意思呢? density指标也不是很明白呢,谁来解释一下?
作者: benhuashu    时间: 2012-7-20 15:16
请问,这一个语句:+ q" l! r( x- J+ [
gen zf=D1.v5% ^; o$ I/ b" _4 a5 S
是什意思?主要是=之后的D1.v5,这个"." 是什么意思?谢谢。
作者: 海阔    时间: 2012-11-5 18:37
非常感谢;
作者: 墨轩书童    时间: 2012-11-7 11:46
本帖最后由 墨轩书童 于 2012-11-7 14:37 编辑 & E8 W! a! {. E7 U0 T

+ M, Q! i. _% _. V+ Z感谢楼主!!
作者: 海阔    时间: 2013-1-20 17:03
本帖最后由 海阔 于 2013-1-20 17:08 编辑
+ O4 |" }# l( c3 n. q# T- j% u! s
" M3 K9 d$ T) G7 J按照版主的步骤操作到了最后出现+ Q0 B" C: y- F5 h  ?# a- R' n3 ~
$ |# R3 H5 D' a! z- J& U7 P
. hello 600031
- F0 k" D- o/ [5 ~* \avolDo not found8 U5 T5 l) x- {8 o% ~
---------------------------------------换成其他代买,也是一样;总是出现这个问题,怎么办呢?奇怪
作者: maodeng    时间: 2013-1-25 13:15
大师,那两条曲线代表什么啊  k7 t) U/ K7 _% n

作者: maodeng    时间: 2013-1-25 13:15
如果我想看最近 一个月,三个月 六个月 的成本,怎么修改程序啊,谢谢大师了
作者: maodeng    时间: 2013-1-25 13:18
我用同样的方法导出数据,发现有的能生成,有的就不行,不知道问题在哪儿,比如300236就生成不了,我的邮箱是19860225gaoyuan@163.com,麻烦大师了
作者: maodeng    时间: 2013-1-28 13:13
. hello 300236
3 h" O" p1 X; R& M( vstatement out of context, W) J. d9 T( U5 h0 L. s' N
r(119);0 j5 @, D# f+ ]1 U: t$ g- M2 E# `

5 V' u/ q: Z0 g6 H; p大师 就是这个
作者: sergio    时间: 2013-1-28 21:11
这贴实在!
作者: awer    时间: 2013-2-5 14:40
maodeng 发表于 2013-1-28 13:13
- p% S) \4 T  m& f; c! L* D; R. hello 300236
: [7 x" p5 c7 k  bstatement out of context2 ~- z! h$ a+ i4 D# a8 k
r(119);
# g( M. k: l* i( T) L" z0 y- q$ m
抱歉,最近太忙了1 ]; I9 ?4 I3 n2 {( |7 ]0 R

  H# F" s  |8 W# x7 X* c300236 我用下面的程序试过,没问题; ?6 ?- v$ L0 T, e6 C' K" Q" }
/ c) t; O7 p. }8 N8 t1 [
你把 D:\stock\test\source 改成你导出的数据存放的目录,再试一试
) p" @: ~' j& a- v0 H  W9 |$ b
. k! d, F3 `/ F. gcapture program drop hello
/ w9 f5 U9 |, g0 [/ A2 e0 d  D
1 t# g9 G$ W- k2 uprogram hello
& r2 O5 _9 B( a7 t8 Z
9 P3 ~. R+ }! O& L( O quietly {
7 M2 t/ g: A2 I( o$ M$ ]8 N5 ]! W0 Y+ @$ }5 x7 ?, y4 N1 J
clear
# X: o( j2 d( }( \+ i! |% y7 {7 l
! S2 O0 B. o2 Y$ ucd D:\stock\test\source
1 e7 |, ?: Z+ a  R- K( J! M+ @- v" w, Z
insheet using `1'.txt
. {2 \  v, S9 ~- f; B, a2 m4 A4 _* T5 G0 R- O. y
   keep if v5>0- Z9 m0 b# q+ B# O2 m6 y
    ' j  [: r+ w  W( V5 C7 R
   gen vt=_n
; L  x/ ]$ l0 F' f7 U   ; n* i9 V3 G) ]: u2 Q8 N
   tsset vt
: K. h; j' j' ]
! {8 [  Z) ^8 l" [' ^   gen zf=D1.v5
. Y  r" r$ l7 V0 q1 r3 s* h! |4 F4 ]- m$ |# f! S! ?
   gsort -v15 t7 u9 d. m$ M7 a( G3 r; L

  v) f6 y* n, b# \9 U3 h! w8 i   gen avol=v7/v5
. Q" V, W/ N7 g' B6 F( \   
7 V. Q1 v: t; r1 ]. t  u2 I   keep if avol>0: n$ ?) O0 L. y( ~: Q4 b( m! u2 \
   
( {& [8 h; ?+ S2 Y   drop vt/ b- h: y" `" v' g
  
# E8 R% k1 M- X8 [   gen vt=_n
' z9 r% H, d# t# m; f3 r% s   
7 U9 B+ y  f- K$ a( ?/ |) n0 d   tsset vt
% a1 _. u% z& b3 F   
$ k' v$ B& P6 v; z3 {! O7 C. s    sum v5 [aw =avol]& _" _% x0 C, a$ U& P
# t# X7 W3 l4 P& l$ }
   gen amean=r(mean) in 1
% D/ g* U' e2 J( W7 l   ) J6 |/ q' }: ]8 P1 u
   * 总加权平均值 amean* y9 @8 y- P7 E5 }+ ^
   4 z0 ~2 x4 _) j3 r3 y9 X
   + ~. r) N7 }2 p4 o2 @
   keep in 1/2408 Q8 A/ Z, X, }4 o1 j* V

7 M6 _9 k# i' Y% T! `6 Q' X5 F   sum v5 [aw =avol]0 Y+ }! j- y+ H+ X) H, l3 c

9 D% _' T* b2 ^5 _" l0 _# ^   gen mean240=r(mean) in 1& D: k9 U4 J3 M' w2 m) j
   * @1 v9 E3 F( F- |( l
   *240加权平均值 mean240
. a' R$ @" o5 c" M1 ]; @1 q   
8 B6 M% G) B# W# [   set obs 500
# Q" f& |( |* ^) c  S. B" [) N7 i+ N9 Q   # R2 L7 J5 g5 l7 Q
   kdensity v5 [aw =avol],gaussian n(500) generate(dx240 dy240) nograph* C( \4 F, L. T- E
   5 s4 Q! N5 O) \7 M% y4 l
   *240 kdensity(dx240 dy240)
' U" h: O- o& c. L   6 D, u1 q3 h. Z: l9 b+ |
   sum dy2404 Q0 m- U7 R$ i% ^, [
   % L7 L) I4 P. ~! c+ X3 ]
   gen yxis=r(max) in 1
7 P, n+ q3 t5 X   
% G) @4 m) `/ ?. K( O   gen yxis2=yxis/2) D: B2 E/ V& Q; [! _1 e1 Q$ T2 Z
   & q. f; T" e2 e4 i) B
   *yxis+ d5 @4 v& ]1 Q/ c3 g' ?: d
   
# C6 T* T2 K# o   gen dy240max=r(max)
) r" \: f7 V* \$ G) [6 e   
8 z) ^) [. Y8 s  Z. P! }   gen dy240minp=r(min)' |; F6 r5 x; I8 \( n
5 q2 G/ E5 \9 q# i9 v# ~
    gen vtt=_n1 k7 O* G8 Y! q* B0 j
) J: T: R8 h0 p9 B
    tsset vtt
- V( L; E) l  I: ^% e+ W+ J: j" d  s* K
   gen ddy240=D1.dy240/D1.dx240, \' Z/ Z! A. O. E8 B& ~
   & |% n  n1 k& E3 W1 `) u6 o+ `
   sum ddy240
$ y# \; Y! ~8 C, m9 f& ^   & z, [5 Y* w' {& K. P+ K4 Q
   gen ddy240max=r(max): o. P0 E& n/ x: }
   6 j; Y& c8 Y' ^# z* K/ @. t9 I; T/ y* o
   gen ddy240min=r(min)# H2 u: l4 C5 v9 r4 a8 h2 V
     
. W8 l: I4 V" M! t3 ]   gen ddy240g=ddy240*(dy240max-dy240min)/(ddy240max-ddy240min)
, W7 Q0 M$ ], l0 t. |0 U" L$ S+ N6 K   
9 E( H) k6 X/ c; B% [2 ?   sum ddy240g
/ S' m6 C  k8 R7 e+ y   5 o' k$ y* H  q" J9 U1 l7 _, B* l
   gen yxis3=r(min) in 1
( Q  x1 _" F; Z4 ?   , s/ |% u) X: i. c" T/ y
   gen yxis4=0 in 1
2 T: L$ v, h, G) r6 l  ( u8 O! t% T) ~1 Z1 g6 }( z$ x
  *ddy240g
7 {7 Z! M% e7 r! F9 c  x+ `/ X   - M+ C% G3 W3 F
   . i5 W% k& u% D
  preserve
" f- f( a$ g) w9 g+ m% e  ) _% u2 f9 i# n2 }8 n3 g, X
  sum v5 [aw =avol]
' q3 x: V5 M' F. ?8 x& h* P  , g; A6 u1 H$ e. e0 F4 K! L
  keep if dx240<r(mean)
: S* F. V  p/ R  
3 K( F* v7 K2 i: G6 u  sum ddy240
1 ^2 Q' D$ ~. I  ( ]% ?' [; O+ Q( W* Y% s
  if r(min)<0 {
/ s6 f8 B6 t  a- |6 U3 X! t        3 b( y/ h$ e3 }1 N( u  w  X0 W
                keep if ddy240<=0% R' ~& }( |  S4 ~& ?# ?0 R8 L$ o
                / g( L9 W! t  v4 X2 r
                count
, M% @0 G' |' }" y7 [/ @                ( }* h4 h1 G6 r+ J& {% l
                if r(N)>2 {
$ }7 c# e! e. r5 f, [& @" A                 & W2 m$ H! v, b/ |
                 pkexamine dx240 dy240; s* u0 z1 }& H6 q+ @: b
                , j+ s$ H) ~9 j% V/ O' r. w
                local tempminp=r(tomc)) \' ^" D9 P2 x4 _" k$ A2 @
                restore ; i" Z: B0 [4 Q/ j$ j7 W
                 }
8 x2 V8 ^  L( S                 " w1 \, ?& n6 i3 l; S$ o) ^, p) n$ Q
        else {, U5 s' a/ s: O$ |2 U
               
0 y4 m( e$ @6 [1 U: M3 |* w                 di "cannot find minp"
0 R; d& X  t" v9 O% e; u                 ; j+ B* u! b$ E& A, e' h
                restore
$ [# Z" @* b! c3 b                 sum dx2408 n) C% S6 g' \$ ], q- I
                 local tempminp=r(mean)
/ q. }0 S( S( B6 Z) [6 K: C                 
( g8 K1 A) H9 n- \                 }- }8 r6 b& n: z, O
                }- s8 e- m' l) S3 ~
   else {
& y: _& q8 P6 ]+ a                 
* ~# v5 E$ L: Z5 q7 W# c2 b                 keep if ddy240<=r(min)
, F( ?' P2 d; ^- N: V                 
; q3 J5 z* h) N2 K6 N- I                 sum dx240
8 v5 }% A5 _! r7 O, a                 ) f3 D8 Z- g2 n
                 local tempminp=r(min)4 e# u8 u* B! {7 i: w1 k; F0 s
                 restore
: ?' J7 v; E" J* ~                 }
3 C5 [/ P# m7 K; k4 v& W, ?' c: h* k" N. I3 R, w) g8 V

3 t. c5 o! B0 e! q* L3 Y: H0 ~
( I5 Z7 n3 E, [. N' z9 J) V! H   gen minp=`tempminp' in 1* Z7 e: N* v$ `4 u4 ?$ I; v
   *240 minp
# s. P6 M1 |1 w+ {! E8 q- x- Q: l6 G% x) z4 x+ x: b$ m( A

/ h/ k  D) u, A' {6 t3 ]$ L   
$ B8 U1 r3 P# Q! |0 B. Y& F   preserve $ J1 h" o1 B9 R6 g4 H
   
0 g3 p( o! ^: M$ Q   sum v5 [aw =avol]
# }, S" {3 ?4 V   
4 ~* h+ w, h0 j0 I8 p   keep if dx240>r(mean); F" j; U7 e" X$ |& k/ y  T1 z$ J
   
0 r3 J9 C, [3 F   sum ddy2406 A9 }" E' k% N; i5 f* U+ N
   - [6 ?$ K# D  [1 G* l6 t
   if r(max)>0 {0 [6 O* Q6 E, W3 z
               
$ P0 ?1 J4 s. U; D2 L& v                keep if ddy240>=04 R3 F$ _2 ^! ^6 J
               
+ i/ [) F* c0 ^$ v& D7 _) q                count
4 z; P5 b: h* N) r                - B! u+ N( ^: |5 n) |6 Y
                if r(N)>2 {' u& Y& f* T" w% H/ ]+ v
                 - C7 |/ O( Y( M
                 pkexamine dx240 dy240) v7 B+ e/ X2 t8 V' D* G+ [
               
" e/ q$ k% |7 O- w8 W                local tempmaxp=r(tomc)
% j: H: x( ^' p                restore
: ?8 S* n2 d! A                 }
: a" u3 `9 q4 T) ^                 8 J2 L- g% x* ?- p( Z4 k, @- Z
                 
! m7 j6 w/ Q1 e; Z        else {# a5 }5 O2 T" e) H/ _
                restore
1 O& j. n0 b% z- P/ M1 T1 k- k                 di "cannot find maxp"7 o, c2 P  W8 ~* q$ p+ D
               
4 J" f. H0 k0 T  h' N. J                  n; g0 S) P9 P7 @* m0 n0 \
                 sum dx240( n+ n2 [; t, p1 }
                 local tempmaxp=r(mean)7 ?* a4 {1 O! |4 U
                 
0 l6 ?4 i1 h9 c/ P4 Z                 }
4 |8 U; g0 C( }, h4 x                 5 R" X* B# k$ d$ e
                 }, k- B& `" U7 r* D  I
                 / D3 N' Q9 b' c+ ~# L/ |
    else {
0 T* m" e* r5 j3 D# a* ^        3 A1 T: E' q3 I: L
           keep if ddy240>=r(max)
' ^6 ^4 x  U$ F9 t' i7 {! l! D! v           9 k: Y  c) ]1 g, P+ S" H0 S! z
           sum dx240( J# ^: x& B9 i2 @# V
             \7 o; z( S* Q/ r( D# y" c
           local tempmaxp=r(max)
; `9 Q/ }% F2 l, N7 {, E( x  {' S2 A           restore 8 {, U4 f- X) S" c3 p  \5 O, s
           }
/ n: D5 ]- S9 N& {       
. |1 K/ H: c2 j1 f1 ?# Y- R' c) O) k; @- ]
  
; f/ s2 z2 ?& W7 X8 G+ B! g* ?1 [4 R9 U* D
   gen maxp=`tempmaxp' in 1
8 o8 D* M0 A( G% G  q! a   
$ [, N' `( [8 Z5 T* _   *240maxp
. n/ O4 R6 U: ~" t! J   / w$ b0 n3 O. O' V5 \
   ' [$ N+ U8 A3 Z) J& i5 j
   *-------------ddy240pminp左侧高峰值-----& I; I9 d) [' X
   - ^: I* P( Y% G3 i
    preserve * m; l$ r2 Y9 m9 C) r- v
   & U* {& \: V4 u9 G4 V6 r
   keep if dx240<`tempminp'
8 r+ Q4 d% z/ [6 ?   3 f7 C: \2 V6 F* B6 B
   keep if ddy240>0
! }2 J/ M( C, X. ]% E, y   ; P3 }* ?: V- G" U2 n, e- Y
   count
7 B3 X# Y, z8 u5 n5 p9 z  u7 v0 W                9 |8 t6 r" P8 x8 w% d5 T
                if r(N)>2 {
4 \6 Q! R2 Y" ?2 V9 ^$ Q* a4 a                 
9 E; A: X* h/ E2 t. N                pkexamine dx240 ddy240
; Y( V+ d9 F+ C6 N4 ^: b. h* F               
9 Z$ ^. {/ F. B6 @8 k/ ^' n                local tempminphddy240=r(tomc) * ?, z$ w: i- X2 Y
                restore
5 E; U1 ]! ^* d0 T6 W/ E3 x                 }' C( N: o& U0 |5 `/ K7 v' h
                 
! }! Y: X3 g6 r; A4 I0 s8 \. `        else {/ w0 k8 t" n  i! A0 ~8 z6 [
                restore $ n" k, w* b& i) h4 U
                 di "cannot find pminphddy240"
0 W/ x7 i5 P1 ], J7 k) K# T) `                 : s  |4 _) e3 A% b* K. @' ?) G
                ' \! z4 S# T5 v6 t2 A: ?% L5 b4 Q
                 sum dx240
: z; |1 |$ u, M# B                 local tempminphddy240=r(mean) - t- A& R5 i: |$ i
               
5 r. q7 Y; L2 P" A                 , ?+ Q( Z7 C4 ^! G
                 }5 B1 ]  q4 b0 M8 n  v/ j5 P/ ?
                 
9 O) E  x* v1 A+ M& i+ a% M   ; S# l# |# ~/ @8 W, b
   
9 }9 x; p& n( C7 D" L$ Z8 T4 Y; O   
& V) t" L2 ?+ w   
. ^  P2 h* v- t0 k; m   / _9 L9 o' d$ P- ]
        gen pminphddy240=`tempminphddy240' in 1$ G' x5 M  N# D5 K& J6 T$ v" N1 w
       
2 ~0 ?. S) V& d       
: }- [( O# e1 U        *-------------ddy240pminp右侧低峰值 -----------------9 o" b1 \, l( l" u% Q, `
    preserve
3 A/ w$ x. W( [: R # i. n( k% M$ b: t* |+ @
   keep if dx240>`tempminp'. G6 H" d+ H2 U8 h4 }4 ^, S9 h
   
# w% Q  E, y& ^/ K3 w# [" Q5 m   keep if dx240<`tempmaxp'
' N& Y2 P, K+ w1 c; `   6 r4 x  z& T4 a8 e( j
   keep if ddy240<0. g% D( G3 ?7 Q( R* K
   3 u5 T4 M. o! X3 r# c. L
   gen temddy240=-ddy240
! p& F! A9 g8 c5 x   
' m0 }' M/ Z4 o9 j    count
% f+ `4 b1 U4 V2 M               
& w* j  c! `+ s                if r(N)>2 {5 m" j! u' m  u! k8 e
                 
; D& p) D" L6 Z4 l                 pkexamine dx240 temddy240
$ ~. c& y* u% Y3 o' P8 d   
2 R7 \( d2 I9 Y* h! i% Y  local tempminplddy240=r(tomc); |2 ~4 j2 H: _2 b# W- a
                restore
/ B( {7 Q1 W2 L6 n                 }
7 _6 {% D2 X, z) b! K& L& y$ m                 
" J* `% r. S/ D+ F7 d        else {
/ \$ x2 @" e$ [( L: p# n" e               
' U: u. r$ i3 B" [                 di "cannot find pminplddy240"' D' v" X+ y; K/ o+ m9 r) _0 F
                restore, [1 R, M& H5 E2 P
                 - }, M/ Z/ l0 w) n  @/ W
                  sum dx240# \' x* k* l/ E
                 local tempminplddy240=r(mean)
2 P9 U: s& {$ n& _                   ]  Y2 a, R2 {) M% k( X
                 5 `# c. l0 c& S
                 }, i# J7 R& n- h" v! K
   * F) y! ]2 W4 @
   
& L/ t) j- x; T' ]$ P$ v' @   
$ v& S$ }! L: n4 N+ B   5 J  ]# Q! c, T4 M
   4 y) H5 Z! B; I0 t, `2 v
        gen pminplddy240=`tempminplddy240' in 1+ ?& j: A# g  ^# w8 d0 x
       
& G) B  H! ]0 n! E2 D        *-------------ddy240pmaxp左侧高峰值 -------------------+ v% f2 c$ y( X3 e' d" v
       
! w7 t2 N, X3 O1 F( |. J         preserve 0 i' B3 x! y- Q7 U. C
   * X& ?8 N- I; F* x# D( a3 V
   keep if dx240>`tempminp') C8 z" M! N' Y* z
   
5 M2 k. B' P7 g) U, L   keep if dx240<`tempmaxp'
: a4 ?9 \& N3 c* j1 x1 z   
2 L! u0 @0 F* B   keep if ddy240>02 K0 w0 \* B8 D& A, Z  o
   count
& P8 l% r) B& P! T; x& k                * \! b0 a8 n1 |5 D# O
                if r(N)>2 {* n" [" f& [# J- }" }
                 
3 ?( S% i- v4 }2 b* @' ^                 pkexamine dx240 ddy240& d' r! a1 D8 v% |% E
   
" a0 P. i; B: n% y   local tempmaxphddy240=r(tomc) 5 t9 J4 ?4 U, A& z
   $ ]; r$ a# {' V% Z  s  u" [
   restore3 R* R, q) ^5 k& S) O) y( F/ M
                 }/ b! A0 [' d5 s8 o
                 1 H8 n  E8 j9 V9 j! w. z- f
        else {" A# c" h- X9 d% O. z% `- ?  h
               
) ]8 K! Q4 l+ K                 di "cannot find pmaxphddy240"/ g3 O3 l3 }' d6 ?# g; [
                 restore
# q( @+ Z4 d7 ~6 I. J2 D               
* r) V1 C- W* I8 w: {/ v                 sum dx2406 u0 c' A! O) E% w" h+ B
                 local tempmaxphddy240=r(mean)
2 w4 w' Q0 s* ]. Q" }! s                 0 b( [4 A+ ~( C* u- p" Z1 r) K) w. J8 z
                 # m# T7 y% P. X- \
                 }$ g- L$ n. `# F( ]
   6 h1 g8 @( ]! ]" r4 M! ]
   * w3 U; q+ n! {8 |+ P( J& B
    ' y0 r% n* |" E  O$ P( s2 \2 P
   
8 k6 I0 R/ G" I6 g   
* N  S# D. q- K! h$ Q   
9 y, A% N% [/ F/ v* f6 M1 I& s        gen pmaxphddy240=`tempmaxphddy240' in 1
) j3 j& J  f8 Z0 Q& g* N, O5 z        ( \( d7 f) m6 U5 p
       
. h( g3 z6 A3 |* }        *-------------ddy240pmaxp右侧低峰值 -------------------2 N! U9 D( }  d1 T
       
  N2 x/ b; F4 j' d* A! |4 h( z         preserve ) u( h% k5 Y* r0 U2 c6 {: V
  3 A) }; {. D/ w' U: l4 O, }
   keep if dx240>`tempmaxp'
) ?9 p! L2 k, Z/ U$ O   
- I( A# q. |/ w8 K   gen temddy240=-ddy240
5 m& M$ D  J. X/ j& S. a/ D   
. M* \1 D+ v1 H% j( J' \/ O   keep if temddy240>0
, Q) A: ^$ d* F/ p   0 P3 Z' o" G( P! Q
   count
& l- P$ n( H3 c- b, U                " h; q# p8 F% ]2 H
                if r(N)>2 {
* }$ S7 ?& y0 I" L                 * t  O, _" G; ~5 ^. N+ v7 a$ k
                pkexamine dx240 temddy240
% L" ^: ]2 j" i5 [) A) o+ n   
1 }! r9 `3 I0 {& e   local tempmaxplddy240=r(tomc) $ ]' w7 T4 {- b1 e# J6 C
   restore- C! v' J; x$ _# K2 Z( W$ B
                 }4 \1 B) Z- ]* X* d6 a! p1 y, c
                 1 z2 C) H) ]: ^
        else {! m7 r2 j" O, j' @
                restore
  o+ }! y0 m1 ^) h; \                 di "cannot find pmaxlhddy240"
% E& J7 h, U: E0 H1 C' D                 
! }! I3 A/ u6 A! y8 g                  O. Y) ^1 S- m) A4 I8 v
                           sum dx2400 R* Z: V6 ^  p" c" k* p
                 local tempmaxplddy240=r(mean) in 1* _* B, c( C6 o, W) z1 T! ?' v+ A
                + w' i( E+ m4 F: D/ H; j4 D* {
                 5 @+ T& A! E! n- Z/ K  J5 w2 U3 u) z
                 }
6 c. y. `/ F% i   
. i( [) p( s4 R) ]   - ^% J4 w! H6 m- b
    1 c1 T8 E* ^: b2 q
   
* [7 L- w* ~4 x& [; X   9 z. A% n. p! Z2 X: v4 p- f
   
2 A& u  F3 y/ K! ]2 \+ K, i$ Q, ~6 [        gen pmaxplddy240=`tempmaxphddy240' in 17 e7 Z; r( ~# m6 d: d
       
1 F4 y5 ]  }# C* R3 z, q6 p% b        *-------------------------------------7 ]9 L0 Y5 a/ g: p
   , Z% ]: o/ r4 m
   gen price5d=v5 in 1/20
1 E  ?& R* [$ ^- k      . V# e2 A# e6 A2 f
   gen price5yy=_n in 1/20
( H- g) T4 U, M. d3 M  B- w   
! n7 {( x/ d& y" K8 R   gen price5y=price5yy*(dy240max-dy240min)/20% d3 e$ G* ]- l4 D( C
   ) y8 `$ C( V2 R9 ?! N
   *price5y &price5d
2 s& t/ Z9 Z- Q/ a) F- r4 m   3 X5 F/ M& n( a! u5 B+ t7 X
   gen priceny=price5y in 1
) v5 ^+ z0 G5 A9 [   
+ P* \! V( Q# M" F! ^6 T: B   gen pricen=v5 in 11 x7 \/ _7 b, v
5 Z( b9 g' H4 {: l. {
     + n  A/ B3 ^3 N. y/ W5 ?3 d
   *priceny &pricen- i9 T( u8 b- v/ a) _1 N
   , p6 A% e2 f0 w6 G1 D+ w: D
   
  `4 f$ d1 U4 R% _   *-------------------------graph------------------) L. Y; V: ~/ e
  , J3 Z! i" q6 \6 c! ?
  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)  x! s: p  g' Y0 O

; U' `$ f* b3 ?5 ^4 g! k* ---------可选项目----------------------
1 ~- i3 Y/ L3 E*(spike yxis amean,lc(blue))(scatter yxis2 amean, lc(blue) msymbol(+) mlabel(amean))
( _8 e) F/ F5 w. ]* h}9 l! C' Y3 j9 ^
end( h$ \% r% T. \! Y" U  f

作者: maodeng    时间: 2013-2-7 00:03
awer 发表于 2013-2-5 14:40
/ o% ]) g8 R# p: \  D# B: C抱歉,最近太忙了
* l& @9 k( C: e
& O& h, @$ T7 i7 C3 b. p, t% C300236 我用下面的程序试过,没问题
7 Q# t" t  e& X
哇,谢谢,我仔细看看,这些字母数字代码晕了- H" a( j4 v  G$ Z3 g

作者: 喜欢喝冰茶    时间: 2013-2-7 04:05
stata是不是商用的,R里应该有类似的东西吧
作者: 喜欢喝冰茶    时间: 2013-2-8 07:16
awer  老虎推荐的一个处理数据的软件,很小,只有10m左右,很好用  发表于 5 小时前
5 E& P5 w  ^# y% K# Q5 d4 }. q

3 O3 G5 H8 I0 W5 A) X这东西是开源的?知道有个统计软件叫stata,不过要收钱。
作者: 汉水云潇    时间: 2013-4-17 17:14
准备试一试( R& ^: p& O& F+ }' Q/ s: k. ]
先找Stata软件
作者: 汉水云潇    时间: 2013-4-17 22:21
本帖最后由 汉水云潇 于 2013-4-17 22:27 编辑
& d' u* `+ C  F0 S& J5 N6 t
2 o0 g% v/ h5 n7 u晕呀
3 T: J. a6 R. _; c2 H* S% F论坛复制程序代码
! K5 i8 U2 t; ^* C. h/ J乱码夹杂在其中
) g4 \8 V% {9 z/ p7 \0 g很烦
/ |5 q. F2 K& K. l/ E对了两次
/ |( D; f; x) Z/ c6 @9 W& S以为没问题了
3 w  l# _6 l7 P4 C+ W, y) z' i% B8 A7 M" m0 J$ i3 b% g
结果
( q% G; p5 q& P& W6 y执行后3 E  X  F5 \8 ~) T7 [& G) x, G
还是出现
2 @4 R4 b* C, F: H, X3 \& t% mr 199
* v9 h6 x3 @  q9 l* W的错误( E: v' y4 F& [7 _6 `/ t+ _! l( ]* Q

作者: 汉水云潇    时间: 2013-4-18 07:50
试验成功,下面是三一重工的图; o7 u) F0 M3 |9 K9 T

0 Z- Y5 K0 |- \- s+ ~
作者: peterbev    时间: 2013-4-26 12:46
先用用试一试,不知道准确不
作者: 中秋不圆    时间: 2016-3-12 15:46
噢。。。这个帖子好,比较直观。
作者: bknite    时间: 2016-3-21 23:51
这个帖子值得学习。
作者: jeanny    时间: 2016-3-22 10:41
感谢!!!
作者: Hubert2000    时间: 2016-3-23 10:23
明月回春 发表于 2011-8-13 20:17
1 u: g/ b8 H5 G" Q# D程序写的不错.1 }1 i" S" U/ X9 h6 X: X
不过有一点要注意.通达信导出的数据,只对价格复权,并没有对成交量复权.所以数据需要手工处理 ...
: _  l! p6 Z9 H1 i( k& e; P8 H2 T
老哥,也是山寨这篇文章找不到了,不知为何?
作者: axiutime    时间: 2016-8-3 13:28
提示: 该帖被管理员或版主屏蔽




欢迎光临 爱吱声 (http://www.aswetalk.net/bbs/) Powered by Discuz! X3.2