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

墨轩书童

5 L0 Q( @  m5 F
! u% \( M+ K/ x' E8 j感谢楼主！！

海阔

7 N. g& L5 `# e- S, Q% ?/ G- p/ [8 s
: S! T7 ?0 U0 o$ q1 b按照版主的步骤操作到了最后出现
/ s) k! x$ ?+ o4 ~6 y
( _# e: \) z3 E. _7 ?3 ^. hello 600031
- J1 f4 l- {. f3 q- N1 D8 }avolDo not found
---------------------------------------换成其他代买，也是一样；总是出现这个问题，怎么办呢？奇怪

maodeng

大师，那两条曲线代表什么啊

maodeng

如果我想看最近 一个月，三个月 六个月 的成本，怎么修改程序啊，谢谢大师了

maodeng

我用同样的方法导出数据，发现有的能生成，有的就不行，不知道问题在哪儿，比如300236就生成不了，我的邮箱是19860225gaoyuan@163.com，麻烦大师了

maodeng

. hello 300236
! z; j2 j5 @8 Bstatement out of context
r(119);
' \9 [% X% r* n& a
' a7 C0 B& x& }& |6 Z大师 就是这个

sergio

这贴实在！

awer

maodeng 发表于 2013-1-28 13:13
3 }6 U. N  g! W. p* \$ l4 P. hello 300236
. b) `$ c. Q: Y6 T' Bstatement out of context
+ e- H1 D1 M6 Q% S- o" k7 z& W" Ur(119);

/ p' ^3 H! G" ?( Q* U' l抱歉，最近太忙了

300236 我用下面的程序试过，没问题
. N" q7 `8 _" G1 C
4 a8 F" |6 K$ g3 G你把 D:\stock\test\source 改成你导出的数据存放的目录，再试一试
4 }/ q, \; Z8 G' K% u. S# y* q
, u7 @4 g; n' v* C4 o9 ]& hcapture program drop hello

program hello
) {5 Y( \9 @, T5 x4 `& g0 s" X
( ~1 p' A- y8 ~6 f* ~ quietly {

7 d! g; W+ M0 Y5 P& @clear
! k, o$ Y- f: m: x1 ^4 Y
cd D:\stock\test\source

. f6 V: S) r# j8 o4 minsheet using `1'.txt
% B7 e4 N8 o9 K
2 N: J) |. W+ ]5 ~5 F! o% Q   keep if v5>0
: Y- Q( U; u1 X% K6 @9 w. _   
   gen vt=_n
   
) `- }9 ~; M) i+ ^9 s) k3 Y) q   tsset vt
% S6 d% Z4 x" I" ~% z  v) E
/ u# u; m8 N0 b) m9 S   gen zf=D1.v5
+ t8 A2 q2 W# V$ P8 ?+ X2 z; R- q. d" p
   gsort -v1
3 X& E6 F) P$ R1 S! r% ?" O
   gen avol=v7/v5
   
   keep if avol>0
   
1 t. p) q, m+ u1 Y7 y( ~$ q   drop vt
  
   gen vt=_n
( U: _( t2 ^) y" y7 X" u+ i   
   tsset vt
   
    sum v5 [aw =avol]
# M0 H9 u2 Q. Y$ E0 t% r
   gen amean=r(mean) in 1
! M3 l6 y( }) O3 E   
   * 总加权平均值 amean
   
! S3 m( Y2 j5 [! D   
! C/ V" Y  }* K3 N, a0 s   keep in 1/240

' E# h$ q- g" ^0 [4 ~   sum v5 [aw =avol]
( O* \  v8 L; y" `( k" e2 h
   gen mean240=r(mean) in 1
   
   *240加权平均值 mean240
$ U! o* _) g2 Q' N6 Q   
+ h/ x+ a2 R. x. w" F" m5 J   set obs 500
* Y" ]4 a" X' E1 D, u; x; e   
) O3 h1 j3 `4 Y  a) \   kdensity v5 [aw =avol],gaussian n(500) generate(dx240 dy240) nograph
; U4 M5 D  D6 y" Z1 ]4 {% z, ]7 L   
   *240 kdensity(dx240 dy240)
! S1 S# G8 b3 V8 m   
! ]# ^! X9 b4 V' M$ Y6 E8 p   sum dy240
   
5 J2 m; }( a0 n0 H0 A8 I   gen yxis=r(max) in 1
7 ~# X2 V; M' ~& W( ^4 h' e7 `" T: B- f   
8 Y! Y# X. r  b& t   gen yxis2=yxis/2
   
! o3 k; ~& q6 D, w4 a+ S7 C/ F   *yxis
* V$ p2 I: j9 o0 k& u   
   gen dy240max=r(max)
4 l# x0 g+ N; l/ m+ D   
+ I  N) m# T2 y   gen dy240minp=r(min)

    gen vtt=_n
- m7 a: ~9 g! m3 A9 V/ G, A
' X! j5 Q3 Y6 T8 n; U    tsset vtt
0 `  O" ^) \$ p0 i. `* \
   gen ddy240=D1.dy240/D1.dx240
+ v& x' Z! q# c3 Z& l   
   sum ddy240
   
   gen ddy240max=r(max)
   
5 I' }9 h# i7 S( q5 i2 W   gen ddy240min=r(min)
$ n3 _: b' v6 m% a4 K  f     
   gen ddy240g=ddy240*(dy240max-dy240min)/(ddy240max-ddy240min)
6 a3 G8 [1 W. a$ K2 k- C- m   
: S; i7 s8 }, ~% w  L4 W$ T. W   sum ddy240g
   
& z% K  d( t4 }. h5 Y0 m9 V   gen yxis3=r(min) in 1
   
/ N6 {: i* }' o6 g. ^$ \   gen yxis4=0 in 1
1 B& H" v# f% c  y, ?  
  *ddy240g
1 w. I( n3 Q) P; C2 B; _   
" K2 v, |. R  q  C7 [8 G/ X2 u% E   
  preserve
7 J5 d" O, W/ u; H/ E! j  
# k) p; G2 }* s8 c, i( B4 P  sum v5 [aw =avol]
9 ^8 T" N" ~" u( M# }# o) b  
; `( f3 i( e( b  s7 _  keep if dx240<r(mean)
$ a& G% H( {0 T# Q1 d$ `  
  sum ddy240
3 _" P8 E4 h) P5 e+ ]1 J  
  if r(min)<0 {
        
# h/ Q6 [. i  v7 H                keep if ddy240<=0
               
                count
! z, C! }) [* h) d" X% l7 t; P               
                if r(N)>2 {
                 
                 pkexamine dx240 dy240
0 q/ @# k2 F! i               
% h0 k0 D6 C3 ^" _6 A                local tempminp=r(tomc)
" ?% B5 H4 e7 H7 N                restore
0 [, g) _, l3 ]                 }
  B+ m- X  I% ]7 B3 r' u5 w                 
) S( W( E# P1 [' J$ A        else {
               
                 di "cannot find minp"
9 K/ e, b  e4 |) `( Z                 
+ Y( w5 X; Q& _8 |                restore
                 sum dx240
                 local tempminp=r(mean)
                 
                 }
7 C+ _; G( S6 D7 C  a" K! C4 J                }
   else {
                 
: ~) p; x5 T8 m* D7 ?  M# [+ Y0 _                 keep if ddy240<=r(min)
                 
6 |6 C6 l0 T3 H" X7 U# N                 sum dx240
+ E: E% F* K# V                 
9 n6 A+ O, g+ n' j  ]  B: K% Y: L7 m                 local tempminp=r(min)
                 restore
                 }
3 r! e! \2 Z2 P, X

3 H. T. `5 B/ L. h
0 {3 z9 u6 V+ ~( P' W' c8 T   gen minp=`tempminp' in 1
   *240 minp
- r3 }) N1 I& Y" e; Q
( ]! I( L- i0 F% c# Y9 Q
   
( y( o9 R: I. _4 `6 x$ l4 _! l& v   preserve
' X$ I4 @8 V5 N) o8 d3 w7 Q* B   
   sum v5 [aw =avol]
3 I- J! u% R3 }6 D% u$ H! i   
/ E% w# ^' Z* K* ~' c8 F0 M   keep if dx240>r(mean)
   
- Q2 m* f; _. l3 W   sum ddy240
- o! d; Z: g$ a) U) Q- R  j   
   if r(max)>0 {
               
: j# _6 i2 D7 L9 l9 b" a6 q                keep if ddy240>=0
               
/ v/ m1 }$ K( v' r! e. k- e/ a                count
               
7 w: I1 Y- p+ Y, Z                if r(N)>2 {
                 
2 I: B; r+ ^' W  V2 d) t                 pkexamine dx240 dy240
               
                local tempmaxp=r(tomc)
7 F- r7 K8 Z6 Y8 y7 A                restore
                 }
) L! U4 t2 W- v3 t  ]' G4 N6 u                 
0 w; f# W; r+ S) z4 G6 x                 
        else {
                restore
) v6 b6 W% y! r                 di "cannot find maxp"
               
4 l7 v: F. w+ l: r               
6 B* H) o# @$ J, e6 W8 [                 sum dx240
. F: X  e6 w2 v7 t+ H2 g, a5 U                 local tempmaxp=r(mean)
5 |- j2 h2 y$ l& h- ^9 S4 W                 
7 R6 k- }/ n7 z( C$ W+ l1 f) S% r( \# ?                 }
                 
                 }
% {! ]$ n5 g8 k! Z6 E! c                 
    else {
& q6 I, y* g* U" Z1 P$ O       
; C7 J* e2 Q$ u) [           keep if ddy240>=r(max)
# [% N& n$ E, @2 v* {/ S5 O9 ~          
           sum dx240
5 x, d1 J, V3 c# m  C+ |* Y5 ]' f          
           local tempmaxp=r(max)
           restore
           }
  A9 D8 M6 D5 q2 j) j       

  
* z. i8 d; z" q5 Z4 C
4 E+ T, D1 ^" Y0 {) A) K* v   gen maxp=`tempmaxp' in 1
; n, j: ?# ?/ C   
& k5 H. `* }- t* [7 g   *240maxp
   
3 _' Q) D* `  ~( a- B' }9 ]   
) e% {' ]; H8 R# J( J   *-------------ddy240pminp左侧高峰值-----
   
' _/ R5 [9 j' Q4 b4 K7 L    preserve
+ A% w" t! w1 [( W7 [$ ]   
# N5 c2 ?  y# F1 v9 W; d1 e9 |   keep if dx240<`tempminp'
4 r, K. r4 M# q   
" ^( t" a; P1 D4 o+ I9 Y5 P   keep if ddy240>0
3 D  l; G8 H. @, D! j9 h) k3 o1 J. Z   
/ i0 g& R3 u4 L1 v5 m   count
$ l; C- y  n  y! {/ |% i7 F# @* w               
                if r(N)>2 {
, a  W) Y  p+ O, h                 
8 }& Y4 G- U9 z5 n5 S9 L) S                pkexamine dx240 ddy240
/ H5 g  s; ^  g5 z; Q               
4 Z5 Q  ?' E& X/ x) l                local tempminphddy240=r(tomc)
: r3 z+ I6 S2 Y  ]! y/ G                restore
                 }
                 
& I5 e# H/ \# \* m$ {        else {
                restore
                 di "cannot find pminphddy240"
                 
/ y0 ]0 L5 ^$ O- I7 R) Q               
# w, g$ W# q4 i9 W$ |8 O" p) a                 sum dx240
                 local tempminphddy240=r(mean)
               
                 
1 p0 C( i0 f/ R) T                 }
                 
9 b/ k4 \, A% C& V! L   
3 ?; a3 j% n7 u% w2 a  d. i/ g   
   
: C1 O9 ~% H/ p* \4 ?- M: t' s2 r+ i   
, v0 z# D1 J, E) r/ I   
        gen pminphddy240=`tempminphddy240' in 1
" [" J) v0 n, e( V: N+ A% ~       
3 {$ I2 |4 m# n+ d) d8 i# U       
6 l! G/ w/ v9 `' g        *-------------ddy240pminp右侧低峰值 -----------------
    preserve

   keep if dx240>`tempminp'
   
; [8 _( `, @$ V) `0 ~5 p   keep if dx240<`tempmaxp'
4 q+ u* p& g# q# G/ |& m, ]- h) ?( D   
   keep if ddy240<0
   
   gen temddy240=-ddy240
   
    count
! [+ \' R0 B! {/ o% j2 O               
3 G4 ~0 t* v3 x                if r(N)>2 {
                 
                 pkexamine dx240 temddy240
   
  local tempminplddy240=r(tomc)
                restore
                 }
                 
% |7 }$ v. Y6 M8 H) b: d5 |% [4 K) y        else {
               
! Q" u" x  F3 q, {2 i                 di "cannot find pminplddy240"
                restore
2 [8 K, a/ [3 S6 ?                 
                  sum dx240
, S& `2 H: f0 P( i% ]0 e0 y                 local tempminplddy240=r(mean)
# T& F+ x0 V/ ^  M6 I$ E- A                 
                 
; \! d0 K* P2 C% R; a# z* O                 }
   
9 P5 A1 f/ P! x. ]' u   
   
   
1 Q: a" f4 O8 y4 [' G' O   
/ r+ I% p# Y9 ]5 Q        gen pminplddy240=`tempminplddy240' in 1
( ?; m- t$ d! E/ x, C; ]4 _# L       
        *-------------ddy240pmaxp左侧高峰值 -------------------
       
% F0 r* X" H- r/ \& \$ M# f         preserve
: k" ~" m! `* ^   
' d4 u6 j' o  d6 t; Z. c   keep if dx240>`tempminp'
   
) [+ @, s" p* ?. {* c' P# H   keep if dx240<`tempmaxp'
   
   keep if ddy240>0
+ C, t  C3 S  G   count
; `' _) t6 Y8 z! N5 I$ l               
0 r% Z' ]5 O' k) ~2 f                if r(N)>2 {
! v0 T( z% N1 |2 C: }* }                 
" V1 a1 x$ t+ X0 q' \1 {! Q                 pkexamine dx240 ddy240
0 g: p/ j+ u1 D& w" t   
   local tempmaxphddy240=r(tomc)
; c4 H+ o" \! x% S8 [( C   
   restore
                 }
0 [* b0 K: U/ r$ [! p                 
. q4 Z1 I. }- B+ Z% x. V        else {
# f4 o" y4 ?7 y( F1 T3 t2 ?* Q               
4 c5 p% k7 p) [, G( `, i; n" R! K                 di "cannot find pmaxphddy240"
2 O1 Y( R+ u6 ~9 m                 restore
               
: Y* \. Y, j& k; A; o                 sum dx240
                 local tempmaxphddy240=r(mean)
% r, g$ Z. c$ ~                 
                 
5 ~7 y1 U5 w' G  F# K% t                 }
   
4 X% q: X  z/ e" T6 h/ \7 C2 `   
   
5 W8 m, g) \) M, P2 u, T4 x% T" D# [   
   
  L  h9 U% T8 A# \5 C0 K  d   
        gen pmaxphddy240=`tempmaxphddy240' in 1
       
4 E! H( ]3 k! R       
' B" |" W  T; X; J  B        *-------------ddy240pmaxp右侧低峰值 -------------------
% N/ n5 }+ u3 Y& Q( N; U  Y       
         preserve
- y4 P( z0 t) S8 n  
   keep if dx240>`tempmaxp'
9 r. Q4 X: P6 {: [% }   
; `+ _: P  C1 h" C/ g' {$ `/ `   gen temddy240=-ddy240
: N- |# ]" R" O: e% R: _   
+ w& {! o" y+ j1 @- N0 d. B+ K4 n   keep if temddy240>0
& y! h6 N6 J4 @9 a8 u* G6 O/ P' X   
, F  e0 W: k0 i2 Z( j+ c   count
               
                if r(N)>2 {
9 o1 P) N) Z% l* D2 |1 E" w! |                 
                pkexamine dx240 temddy240
7 X" {; S  Z* Z; E, ^7 P   
: W& Q2 r8 v- N  I- P" t   local tempmaxplddy240=r(tomc)
   restore
; ^  B+ K1 Y! D. V  |: n( f0 s                 }
                 
. D0 \% E4 C- x2 _4 m        else {
                restore
                 di "cannot find pmaxlhddy240"
                 
2 x* A! S0 z. M9 l. T' i               
' r* r6 x  ]4 j. B: u6 j! `! c3 z                           sum dx240
                 local tempmaxplddy240=r(mean) in 1
  |; Q/ |3 e$ C% m8 _               
8 }( X; [1 d. ~' H% G4 k, r                 
! X- \3 O6 h- D" s! t                 }
   
   
   
   
6 I( K& F" T; q  v   
" d, ^7 ~3 N5 m+ p   
- W+ n$ ~; ?! K1 ~: _7 E3 a        gen pmaxplddy240=`tempmaxphddy240' in 1
( S6 ]' z2 w" E2 V. W- I% D- a$ j9 j       
        *-------------------------------------
/ |7 ]2 u6 N: ?6 ?   
   gen price5d=v5 in 1/20
      
   gen price5yy=_n in 1/20
   
4 t  X* @, p9 i; T% z, E  N   gen price5y=price5yy*(dy240max-dy240min)/20
! t) L! u9 ]4 X, M+ f% j" U   
3 L, _+ z; H, G   *price5y &price5d
0 L8 H0 Q) t4 h7 y5 G7 z   
, l/ X9 c8 \9 H3 ~   gen priceny=price5y in 1
/ c) \2 H! a$ ^   
   gen pricen=v5 in 1
. ~5 `2 Z. E/ R( @; P
     
   *priceny &pricen
0 k/ D2 q( C4 M2 s   
   
   *-------------------------graph------------------
- a$ g! ]5 l/ H& _. u9 f  
  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$ q: u  X* a7 t& W6 U$ f$ B% {* ---------可选项目----------------------
8 c- t7 c6 q5 n% S3 n*(spike yxis amean,lc(blue))(scatter yxis2 amean, lc(blue) msymbol(+) mlabel(amean))
}
end

maodeng

awer 发表于 2013-2-5 14:40
抱歉，最近太忙了

300236 我用下面的程序试过，没问题

$ s5 w% Z" R: X哇，谢谢，我仔细看看，这些字母数字代码晕了
. r# D: x) j, ?7 v. \4 T7 B  o/ H

喜欢喝冰茶

stata是不是商用的，R里应该有类似的东西吧

喜欢喝冰茶

awer  老虎推荐的一个处理数据的软件，很小，只有10m左右，很好用  发表于 5 小时前

这东西是开源的？知道有个统计软件叫stata，不过要收钱。

汉水云潇

准备试一试
# P5 F+ H8 l. |. d+ {先找Stata软件

汉水云潇

8 K4 K0 j8 }1 f6 x0 j9 J( O
! I5 Y/ h1 Q5 V$ c5 f) p6 M晕呀
论坛复制程序代码
乱码夹杂在其中
很烦
对了两次
0 }$ f4 B- e! g以为没问题了

- a, B$ f$ _& Z4 ~4 Z( X  K6 H8 Q( U结果
" x" I3 C* F9 O, s; \& e! y1 B执行后
  y! e) C5 Y3 c  `还是出现
r 199
5 e; _! V& r! g  |* }: A) }的错误
5 E* `  V* w% L, j* J. I$ J

汉水云潇

试验成功，下面是三一重工的图
+ e2 b3 i  g: M2 d( S
0 o0 a# |, y* p( D

peterbev

先用用试一试，不知道准确不

中秋不圆

噢。。。这个帖子好，比较直观。

bknite

这个帖子值得学习。

jeanny

感谢！！！

Hubert2000

明月回春 发表于 2011-8-13 20:17
. Q+ v, W$ ^: \  H$ E程序写的不错.
1 ~( o) B8 s- y: G0 V不过有一点要注意.通达信导出的数据,只对价格复权,并没有对成交量复权.所以数据需要手工处理 ...

4 u3 ?  x' }4 N0 W4 d/ E0 m老哥，也是山寨这篇文章找不到了，不知为何？