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

maodeng

wenty 发表于 2011-10-27 04:17
& s! X2 x( G6 L5 \# b2 j. l总结一下使用心得
1 V; t3 t! d3 _1 v, Q$ F8 l: |1）删掉set memory那一行
% f6 v) T& a; z8 m5 \* U2）改路径

1 \0 i* d+ r1 C记事本打开后，怎么用tab

冰火

认真学习

higer

我也曾尝试过使用stata软件对股票平均成本进行分析，但是作为一个非专业的炒股人员，而且日常还有工作的人而言，这种方法似乎显得太专业，而且操作难度极大，不一定实用。
8 o* V( z9 c8 l( E7 [8 b8 a
使用这种专业的技术分析方法容易造成这么一种现象：大家都对导出的光滑漂亮的图像感觉良好，然而在股市中依然无法进退自如。

jway315

谢谢楼主，希望继续给同学们挖宝哈！

awer

改了一下，这次运算应该没啥问题了
0 T8 `! k5 h! ]$ C
' @' M: f& Q/ g4 D* X, }, n9 F
capture program drop hello

4 `* N* w( I9 z# q9 f/ M! oprogram hello
! ?# ]1 n' K# s% p% x/ y& X2 t! M
+ \; L, N7 p+ I) g1 N' p+ E quietly {

clear

cd e:\finance\stock\20110101\source
. G  d: N. a6 ]8 g+ w) D# I2 d$ ~, m
insheet using `1'.txt
: W1 Q# A6 S- G7 F
   keep if v5>0
8 _) H$ l! k$ d9 \" J   
/ k% Y: v8 z& J; N% G, c4 x1 v  w   gen vt=_n
   
! J0 Y# T7 D1 i& e, O8 e   tsset vt

   gen zf=D1.v5
5 k) {8 k8 a6 O- o
+ K4 |! s8 h6 }7 ^# s" }   gsort -v1
4 X) k  L6 j7 I- f1 U! G3 A
6 R! G  \, D, g8 \' I8 `5 h/ l   gen avol=v7/v5
* l6 ]" k0 m6 s) D   
" H- O6 F; Z+ P   keep if avol>0
, g( w2 W6 S' y3 F. G   
   drop vt
, A6 P+ Q; Y+ v$ x) Q/ j1 ^  u0 ]% q/ d5 L  
& N9 m  S) u! @2 `" C4 h! n6 W   gen vt=_n
& }, y; o! Y+ x2 _   
   tsset vt
   
; A: z3 z) j! h% {4 Q- k+ d    sum v5 [aw =avol]

   gen amean=r(mean) in 1
   
   * 总加权平均值 amean
   
: @) Z# e& |3 [- C6 W& x; Q9 }   
   keep in 1/240
; G/ I! V/ g7 _1 M  r1 s
  z' A$ k6 E6 K3 H$ y; V   sum v5 [aw =avol]

   gen mean240=r(mean) in 1
% H8 Y9 X1 V: M7 {9 R( X   
( U3 U+ t4 S6 R+ _& I/ ?( B. p   *240加权平均值 mean240
   
   set obs 500
   
% s- [& b! f( [7 }& a' j   kdensity v5 [aw =avol],gaussian n(500) generate(dx240 dy240) nograph
( j; l1 D/ c9 @9 F8 h* B% ?; l   
   *240 kdensity(dx240 dy240)
6 W; u& [- J: H   
   sum dy240
* G# k) }6 U( t9 U$ ]# [. [   
2 B# ~* H4 M" l$ Y7 a- a9 S& H   gen yxis=r(max) in 1
) ?  N. p: G/ U/ T. E% `3 y6 C   
  X$ }9 c, }8 ~9 o# [. T* K$ T   gen yxis2=yxis/2
3 K$ `8 h% u7 n3 b- ^; R   
! D$ Y! T0 M% k- @( l; U' q8 Y   *yxis
   
   gen dy240max=r(max)
   
   gen dy240minp=r(min)
/ E. G. J# y4 c) u  S
+ ?5 y7 \# E* m; E8 G2 V4 O1 }# W* }    gen vtt=_n

    tsset vtt

   gen ddy240=D1.dy240/D1.dx240
/ t; j/ A! @# s  s: Z. @) d( c) y6 y5 ?   
* P( a& v8 u2 |+ N$ r   sum ddy240
7 ?, s& r4 z/ c5 t: U2 y   
0 A2 O' ~- o5 M, v3 I$ q0 [; J% g   gen ddy240max=r(max)
   
   gen ddy240min=r(min)
     
   gen ddy240g=ddy240*(dy240max-dy240min)/(ddy240max-ddy240min)
   
   sum ddy240g
8 L3 s! {+ L/ L# d4 s' Q* }+ x   
% W9 \4 }9 \& O8 L   gen yxis3=r(min) in 1
   
   gen yxis4=0 in 1
  
  *ddy240g
   
1 M! ]* f; m. T% {   
  preserve
  
+ I7 J! {, J! Q5 {6 J  sum v5 [aw =avol]
  
- A1 X- i0 ]2 C- _: M0 o6 G  keep if dx240<r(mean)
  
! n0 c2 `8 K/ a) G. ]7 g$ \  sum ddy240
$ f& i& r5 I5 ]* y* x  
  if r(min)<0 {
        
# B3 S9 k8 F) l                keep if ddy240<=0
               
/ U5 T9 K: P/ `& ^( t                count
1 r" H' [& e( p* G0 ]. z               
3 J/ @$ ]1 ~' Q  V) d                if r(N)>2 {
                 
6 z7 x! I& D$ X! Q. a) N                 pkexamine dx240 dy240
               
                local tempminp=r(tomc)
- u9 ?3 z) p0 Z: v, |7 |" T. ]! s* R                restore
                 }
                 
        else {
               
' m5 i" b' N$ a* `0 g* E+ _( \                 di "cannot find minp"
                 
                restore
                 sum dx240
                 local tempminp=r(mean)
6 X  g, \$ v" n/ h0 w) _                 
                 }
2 s0 E: Y7 @) {# x/ y3 i0 S                }
2 j( [. h4 N$ e   else {
                 
& U  |; H1 P0 V9 c* O* J                 keep if ddy240<=r(min)
                 
                 sum dx240
                 
9 N- R) _9 [* n& T% j, |# J4 x+ M3 U6 a                 local tempminp=r(min)
! I& H: I5 f" ?# d2 q& @5 x4 ?                 restore
2 w' a- D% V3 H. x8 b                 }

- n* a: P: _9 v- q# C( K6 u; _! O0 G

$ i! F! P, v8 w9 H! M. D   gen minp=`tempminp' in 1
3 D8 j- z: t, P% x1 J4 G, E9 L   *240 minp
- X3 O4 I2 X1 o# R9 Y

  R( U4 Y5 k, U* ?+ V   
/ [+ f) S# z2 u) R0 s6 h+ |: H   preserve
* e5 t6 L0 E2 E* B9 Q3 y   
   sum v5 [aw =avol]
   
   keep if dx240>r(mean)
+ O7 q# I8 Z. F+ X+ X- }, Y   
  `! x+ Q% m1 H' ]5 K$ B   sum ddy240
9 G8 N6 v2 i. @$ u% W/ \   
   if r(max)>0 {
               
6 `3 o7 l3 `* t4 t                keep if ddy240>=0
5 a2 n! Q6 W3 ^8 \9 l               
/ u8 W9 W8 B" T: Y' y                count
               
$ o) Q6 B! y9 c2 k                if r(N)>2 {
  l) }- U/ r" p$ X: Q1 i                 
3 _, g& w! A- a2 w                 pkexamine dx240 dy240
               
                local tempmaxp=r(tomc)
- s: a4 s- s# B! S# o, A                restore
' X- Z1 X& ^) ?9 D6 t8 x! m                 }
                 
" a: a+ f# z% Y% Q: M; w$ m                 
: v8 ?( z* t0 Y% H3 X6 F: G        else {
- d- l# J( @4 T7 @5 r- J                restore
9 j$ U9 E, b8 N2 }" E                 di "cannot find maxp"
* H1 e4 Q' T7 n; [0 O- B' I               
: u: q, Y5 Y8 A1 S1 e8 X               
                 sum dx240
& A7 t- x+ n/ u- C; v: R& y- q                 local tempmaxp=r(mean)
                 
+ t% S! C& ^2 O8 ^+ W                 }
                 
, N# ?+ ]/ b+ M7 S5 u                 }
                 
& O) q! T! b2 R  \; s    else {
% _0 s- {) p, e' i& K- T       
           keep if ddy240>=r(max)
          
, l( s- W- Z$ I4 q           sum dx240
          
+ L1 j+ I  B' u4 _) E. m# n           local tempmaxp=r(max)
           restore
           }
1 w& _: I5 ~/ \8 n       
' O. ?$ W/ R+ S
  \' _( K* {- l' e* S  

1 H+ }" v% k; Z, G   gen maxp=`tempmaxp' in 1
" j1 m1 ?' M8 O2 f& y1 f. U   
   *240maxp
   
4 a4 z& f& ]# m% b+ G   
   *-------------ddy240pminp左侧高峰值-----
8 |6 i, G4 Z# {) k- y   
    preserve
   
   keep if dx240<`tempminp'
   
/ _! p! \4 H& h! M: Z   keep if ddy240>0
0 K! Y6 U  a8 a( p2 W   
   count
               
& V  c. @& J2 L1 y& w* d                if r(N)>2 {
                 
                pkexamine dx240 ddy240
: C3 M3 }2 B# D               
5 k; o: [/ G! W% k) ?- H0 x                local tempminphddy240=r(tomc)
                restore
1 ~4 k# t! e# h  r& j) @3 Q2 z                 }
                 
        else {
- N8 v4 [+ J) X! x( {                restore
                 di "cannot find pminphddy240"
                 
               
4 F8 k1 {% w4 z2 B& T, W                 sum dx240
                 local tempminphddy240=r(mean)
: K; Z2 c% u4 }2 e% o. g' o9 D3 V               
& `8 E( U4 o  E/ \+ s7 A                 
                 }
3 S) X/ R( k. m/ z7 |                 
   
5 S* p* n" s) X# R# d7 E! }   
. @: C& D" p7 v) D9 [   
   
, G% |$ I# U: v   
        gen pminphddy240=`tempminphddy240' in 1
       
       
7 @9 U5 ?- X1 a7 S% }9 X4 g& T# ?        *-------------ddy240pminp右侧低峰值 -----------------
    preserve

% H7 L# s' i8 _9 b) C3 d% v   keep if dx240>`tempminp'
   
& O( X; B! F# n* i, G- Q. z   keep if dx240<`tempmaxp'
  L4 J1 o, M- X0 D1 J" l   
   keep if ddy240<0
0 `: e2 x6 ~7 Q  M& L7 V   
. |6 S/ `1 z4 Y   gen temddy240=-ddy240
! ~& ]+ m) @0 D8 e% M- |   
    count
               
4 N5 Q& [4 I) s/ J                if r(N)>2 {
* M; s. N( F# {, u. n                 
                 pkexamine dx240 temddy240
   
; R6 S! u9 h  t2 y% O$ k* v% t  local tempminplddy240=r(tomc)
) N% g1 A6 n$ x                restore
                 }
7 m7 n) T* J; b1 I' j" r                 
) x( o2 }6 |$ m- G) q+ x        else {
               
' H# N' H7 M. Z+ M% q5 C8 P+ s. ~                 di "cannot find pminplddy240"
                restore
                 
& B! \2 c) _7 s( ^; P5 y                  sum dx240
                 local tempminplddy240=r(mean)
9 K/ ?( d) z( y* \% y. X. R                 
                 
; c5 U+ a4 ~5 Y2 f                 }
3 T2 b4 Y3 j: z' w0 f# Z   
( r' ^+ j8 u. _, C# r9 i. A0 M   
5 J. `) c+ G7 Q; u   
% G" T) E# C- `7 S  B   
   
        gen pminplddy240=`tempminplddy240' in 1
9 J# M/ J. ^2 o, p! f& S       
        *-------------ddy240pmaxp左侧高峰值 -------------------
2 [8 [/ ], K0 T1 t) D6 @       
3 e) v/ e) m% x' p- j$ u4 f         preserve
   
   keep if dx240>`tempminp'
$ I, y8 `* M  V6 E# V& I' y) q/ T   
1 E( I. M: O4 d) [$ ?/ H+ r# ^7 I   keep if dx240<`tempmaxp'
9 |* g( e& o, R1 \) i1 k# h- K   
   keep if ddy240>0
   count
               
                if r(N)>2 {
5 O+ s1 |6 s1 b  C1 ~7 e                 
                 pkexamine dx240 ddy240
# y/ s' ^7 _) q% W8 T5 c   
   local tempmaxphddy240=r(tomc)
   
   restore
) v/ a: B3 w' x" b                 }
                 
        else {
               
                 di "cannot find pmaxphddy240"
" x$ i* i2 E6 m* Q; N+ ?/ }8 i                 restore
               
( M$ j0 g) w, o. ~                 sum dx240
2 h# Y3 U# l3 C0 Q& S                 local tempmaxphddy240=r(mean)
                 
0 K, i8 |  L: t8 `                 
" T1 [1 r5 R/ g1 I  \, K( @, N( u                 }
: \2 i) @0 |* j  K" u   
   
. q% N# T1 x. J# m8 F1 N/ R   
   
0 g, n3 k, q% v( F6 l+ B   
   
$ m3 c4 g( F) h1 U3 c        gen pmaxphddy240=`tempmaxphddy240' in 1
       
: u' l2 _' R3 g( |       
        *-------------ddy240pmaxp右侧低峰值 -------------------
# g1 G5 v& {9 e       
6 T. k8 C2 |1 c3 u$ _; m         preserve
# i" }1 K. ~0 H0 q5 d6 j  
1 f" K. C7 t4 `# e$ c: C$ y2 }; @   keep if dx240>`tempmaxp'
$ j% @) r. E" o& V   
   gen temddy240=-ddy240
8 ~9 K; }! a5 j, V) [) }: d   
3 u" T7 n: W5 g: K   keep if temddy240>0
   
9 n4 U( D+ Y7 `7 S( Q- W   count
               
7 r( }0 ~+ F- Z                if r(N)>2 {
! _1 Y$ w' c/ N7 C                 
! o& l; D: k( }. g2 I& N5 G2 K                pkexamine dx240 temddy240
   
   local tempmaxplddy240=r(tomc)
7 {1 b5 l( |/ s/ p- z2 m7 N' E   restore
$ C* x) Z  P9 O" x3 X. x7 L! ~; \                 }
# o5 d% a% X1 B+ s6 i) \1 Y                 
        else {
                restore
9 `- @$ d' q# [$ Q                 di "cannot find pmaxlhddy240"
# m! m1 E. L# R& G' \( Y                 
               
                           sum dx240
                 local tempmaxplddy240=r(mean) in 1
+ {, N' d+ g! Y! D5 @1 o$ _# m) G               
                 
                 }
   
   
   
1 ~1 M; j. A' i   
   
   
1 u0 g' \8 s# G8 w& i. s, G3 |        gen pmaxplddy240=`tempmaxphddy240' in 1
       
+ P6 R' z& A' z# K        *-------------------------------------
   
) R* l! _9 R2 Z/ s+ f# u   gen price5d=v5 in 1/20
      
   gen price5yy=_n in 1/20
   
' g( A9 Z% [! |# Z) e* O& o9 g3 b   gen price5y=price5yy*(dy240max-dy240min)/20
   
2 L4 ]4 q' p7 ?$ v3 h! Q   *price5y &price5d
   
7 U+ L0 P1 B/ t# `5 A   gen priceny=price5y in 1
   
   gen pricen=v5 in 1

& o$ Z5 O! i, H+ |     
   *priceny &pricen
& Q9 w* {1 L( |: l5 ?) @; |1 w9 N, w6 ?   
! b5 L9 v# D6 }2 R   
   *-------------------------graph------------------
, J' B, ?( G. n( M/ _  
/ b: a. W' w2 F+ c% e9 O: E  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)

' x' q6 r% C+ f' v* ---------可选项目----------------------

2 D  k! Q! Q/ _2 a- }) b+ h}
# t" T; R1 j% }; s8 |end

旅途愉快

awer 发表于 2012-3-28 18:06
改了一下，这次运算应该没啥问题了

这段代码是重新复制粘贴到hello.do里面的？

awer

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

对的，就是hello里面的那个

David09

awer 发表于 2012-3-29 11:40
: e4 ]+ A3 m0 F5 P" \; H对的，就是hello里面的那个

. U7 V/ D' w" @# Q2 z. hello 002680
4 M+ }  F2 E. n8 L+ _; }- TObs. nos. out of range
r(198);

+ R- y& O) W. V* Z, d
8 {- I0 A- a5 d' U& }7 M, I- G我运行后，出现这个，不晓得是什么错误了，请大拿们帮助

qqysc

David09 发表于 2012-7-2 21:35
. hello 002680
: B: U3 ~8 z% c( `5 YObs. nos. out of range
r(198);

7 B$ }& h1 w( o! k3 h估计是导出文本格式问题。

qqysc

awer 发表于 2012-3-29 11:40
对的，就是hello里面的那个

复制粘贴有很多乱码，请楼主把上面的代码做成附件以供下载，谢谢！

David09

$ F* a# V" F9 u( t

qqysc 发表于 2012-7-3 12:14
6 w; d8 R; `9 X复制粘贴有很多乱码，请楼主把上面的代码做成附件以供下载，谢谢！

2 \4 O; x$ k% n) g, L3 G( G$ |  G# h我清除过了乱码的啊....

qqysc

David09 发表于 2012-7-3 17:13
0 w5 m5 I  c; y, P+ ?我清除过了乱码的啊....

: ^% Q* B* x: h1 j: H说得不是一回事。你说的问题可以试试在 002680的文本中，把tab键补齐，导出时不要带表头

David09

qqysc 发表于 2012-7-3 17:56
) A" {4 `& \) w! j& j说得不是一回事。你说的问题可以试试在 002680的文本中，把tab键补齐，导出时不要带表头 ...

很奇怪啊，我看stata中间的数据表，都整整齐齐的导入了成交价成交量。但是一运行就蹦出来那个错误

David09

搞定了，原来是我选了一个次新股的数据，样本数目太少了，引起程序的错误。
$ B( _" h0 I3 ^( D6 g感谢楼主给的好东东。

碧水芳洲

怎么没人问曲线是什么 意思呢？ density指标也不是很明白呢，谁来解释一下？

benhuashu

请问，这一个语句：
7 \8 V$ E3 X# P, C& ^gen zf=D1.v5
: p7 s: }7 p9 Z% j2 t是什意思？主要是=之后的D1.v5，这个"." 是什么意思？谢谢。

海阔

非常感谢；

墨轩书童

  z5 h* H1 U, G  h" k, E感谢楼主！！

海阔

. N6 ]% Z- Y) v! \+ J按照版主的步骤操作到了最后出现

. hello 600031
% ]* t8 `0 Z9 j  m5 qavolDo not found
! b6 P) d- O& o---------------------------------------换成其他代买，也是一样；总是出现这个问题，怎么办呢？奇怪

maodeng

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