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

maodeng

wenty 发表于 2011-10-27 04:17
总结一下使用心得
: a* X: ?, v' [) q$ x6 s1）删掉set memory那一行
- L  D+ S, n- B  l& w' _' @2）改路径

3 M9 t# b" i& {# o6 i记事本打开后，怎么用tab

冰火

认真学习

higer

我也曾尝试过使用stata软件对股票平均成本进行分析，但是作为一个非专业的炒股人员，而且日常还有工作的人而言，这种方法似乎显得太专业，而且操作难度极大，不一定实用。
: p/ i+ b. i" G7 a' X# j- X
使用这种专业的技术分析方法容易造成这么一种现象：大家都对导出的光滑漂亮的图像感觉良好，然而在股市中依然无法进退自如。

jway315

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

awer

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


, l5 J4 g& o" d/ Q& h0 Wcapture program drop hello

program hello
5 \; S6 k3 U6 q0 i% y
9 D6 G- I$ O- k3 U4 H' M* q quietly {

clear

cd e:\finance\stock\20110101\source

insheet using `1'.txt

   keep if v5>0
   
   gen vt=_n
   
0 w+ E. ]( \& T$ t2 a   tsset vt
- R9 c+ m" j4 y6 k, g
   gen zf=D1.v5

( @7 b4 J% R: b! g   gsort -v1

   gen avol=v7/v5
   
0 W1 p! D- P) E- S$ A2 R# N( A   keep if avol>0
+ G: ]% `7 O9 n9 O, Y0 Q: Y4 |   
5 o# d( @/ P! {5 {5 P& x! p& M) z   drop vt
  
   gen vt=_n
/ i/ @4 u0 b; l; O  [( U5 |   
   tsset vt
   
    sum v5 [aw =avol]
! a& y3 A" Q( T
" b# r  O- H4 J- q   gen amean=r(mean) in 1
   
   * 总加权平均值 amean
   
( I# k+ g7 j$ }+ x0 N! O% p   
& B- d" G5 D0 Z! w+ f' ]   keep in 1/240
$ c/ g2 }2 f. }  o4 h
3 f, L  q7 r2 ^; J   sum v5 [aw =avol]

( u, C2 \, k9 U, \6 b8 \   gen mean240=r(mean) in 1
   
   *240加权平均值 mean240
   
- u1 a4 v. E* d8 y+ p   set obs 500
! [0 v+ d$ |9 [6 t; c   
   kdensity v5 [aw =avol],gaussian n(500) generate(dx240 dy240) nograph
   
   *240 kdensity(dx240 dy240)
   
   sum dy240
; S0 [2 k2 [+ S( ]: |   
   gen yxis=r(max) in 1
! W7 j: i1 e- T7 q5 P   
   gen yxis2=yxis/2
. ~  O2 w: ^: E9 k   
   *yxis
   
  _- n0 G% P6 [" j. `4 m   gen dy240max=r(max)
   
9 W. \6 d$ p" n7 h  C   gen dy240minp=r(min)
% O2 S8 S) a/ x, V
    gen vtt=_n

7 d* V; S' N( b# y) Y    tsset vtt
; ^% D! ]- F4 E! [* m6 Q1 P/ a
1 g# T7 r) e  F   gen ddy240=D1.dy240/D1.dx240
   
   sum ddy240
/ e: N) \" z, F; v& K& ~) l( |2 w, s+ M   
$ ~: o. J' E) q: l   gen ddy240max=r(max)
   
6 G5 U0 S- ]! z& H2 R$ O& D( b   gen ddy240min=r(min)
* L# U* t+ O1 F6 E$ s- s7 A1 X9 y     
, ^( H* N3 c/ R7 ^1 |   gen ddy240g=ddy240*(dy240max-dy240min)/(ddy240max-ddy240min)
1 c* k9 o: G  m: n   
5 l# _4 T# F' ^# y4 [6 \) c, g& X   sum ddy240g
. W1 j( X5 Y" S4 z: ?. Y* x   
; O5 ]' K' {0 `3 T1 Q   gen yxis3=r(min) in 1
3 c( {  R4 P# t   
   gen yxis4=0 in 1
" u4 c% w7 M% f: n  
  *ddy240g
   
+ b4 P6 {- g# r1 V4 b/ W; L- R   
( ~% l& E# F* N  preserve
+ n& a7 D3 A& t  u. {  
% S9 b1 w+ j5 a3 f1 p" ~  sum v5 [aw =avol]
& |: M7 h( m/ J4 z  
1 j8 p/ }, E; \, \  keep if dx240<r(mean)
  
  J1 q4 T* M. ~$ r* ?* V  sum ddy240
  
. F* l, U' U# |% A  if r(min)<0 {
8 J& j  V( O2 R        
                keep if ddy240<=0
               
9 [" v! A& c4 x                count
1 e! v% B5 L8 x6 V  p$ L# e& `( T6 f               
9 h2 w3 ]+ j! `' X. ?0 y2 c' j5 y                if r(N)>2 {
                 
* ~& t3 ]' B0 Z, b                 pkexamine dx240 dy240
               
                local tempminp=r(tomc)
/ b& c$ |" m2 P" v* P) Y8 a                restore
+ b' `! L3 r) ~, z6 R5 N6 |                 }
6 _8 f* I' c1 c% }                 
. I5 [+ J. O9 ^7 p0 I5 f! j1 `        else {
- }0 }, N! S& u& t               
7 V0 E7 q/ O' e( K; {( D$ Q                 di "cannot find minp"
                 
* {( p) U" o. y( H                restore
  o+ X' f% R8 W  h1 U                 sum dx240
                 local tempminp=r(mean)
* W6 N4 p- b) y% G! q+ A                 
( C% ~1 O+ A4 b                 }
; L  Y0 C5 i( y' z) V, d                }
   else {
. h* o+ I! S9 h. `" _6 Q                 
# a# V  x% d# P. \8 q& U! ^                 keep if ddy240<=r(min)
2 m4 y: C; C0 n4 O# D                 
                 sum dx240
% a0 f$ @" Z2 K/ c% L% b                 
                 local tempminp=r(min)
                 restore
                 }


" D. m, j4 a5 L/ x% g4 q; e
/ p$ J8 n: ]2 c  T: j, y   gen minp=`tempminp' in 1
- t1 a8 ~2 ^! R% u) w0 |   *240 minp


) }3 R8 b9 q) W2 @7 ?   
   preserve
   
- ~7 {" C( m5 a   sum v5 [aw =avol]
  z) T4 l! @* ?! B9 n   
, u6 N, Y/ d9 K" I2 C: g3 v   keep if dx240>r(mean)
   
4 c! H5 y3 I3 o: `! W" @, \0 k   sum ddy240
% m2 g7 r( I- q; {   
, @, F5 ?6 f$ [: U; j   if r(max)>0 {
8 n  \6 r! v# t8 _) o" b               
( I8 Y& n0 t( o/ W4 |                keep if ddy240>=0
               
- i( ]0 t! a: W: Y, i& H; N  e                count
               
& i- ~" U, J) C  ~, g, Z5 [                if r(N)>2 {
: W" b) I' {: R                 
                 pkexamine dx240 dy240
( g) r1 n, I7 {               
                local tempmaxp=r(tomc)
! E4 _( V0 w% p+ V                restore
                 }
- U. S9 v# {( t, {                 
                 
4 u  X1 N6 l: N) g        else {
                restore
) S4 I* ?4 U5 P                 di "cannot find maxp"
5 j8 f( S; c& v  w               
$ a, c* K* J/ M               
' Z5 Z/ j& K" ]                 sum dx240
" w/ z! b% T4 Y& k4 f; n" o5 c7 p                 local tempmaxp=r(mean)
8 H% Z7 }( k/ Q$ x8 G! ~7 }                 
                 }
                 
                 }
% V4 |6 u  B% R2 ]9 l6 M3 }                 
. ?+ I/ [+ \6 |0 i5 k. z1 G% t  d: j    else {
       
           keep if ddy240>=r(max)
          
           sum dx240
' ?6 d7 L. s3 I# I          
% T, U$ |, y1 V8 I           local tempmaxp=r(max)
0 q% _# C, I* M5 Q7 e7 N6 h           restore
8 q7 {0 F" ~9 \) T8 U, g' e1 ~7 b           }
       
# ?$ Q5 R2 G( d9 c9 @) r
  
: Q  r& x  `0 [/ E' v* R
; |' a) Q# C- O& W& ^   gen maxp=`tempmaxp' in 1
   
   *240maxp
: T# x# v% }/ [; h$ u1 U; _   
   
   *-------------ddy240pminp左侧高峰值-----
: v9 u; l) V' s# x1 S! u3 z- y   
    preserve
6 C% M- h8 E7 f   
; G5 C3 \! \* G8 V  L- r   keep if dx240<`tempminp'
4 S6 L, d8 O# i( z   
/ p  j- y! _' v: v. [! o; E   keep if ddy240>0
( o. E" s( ^/ z# B& q3 f   
   count
# x. W* M) {5 E+ k: K               
                if r(N)>2 {
                 
                pkexamine dx240 ddy240
               
: W! b$ U$ v& q, ?                local tempminphddy240=r(tomc)
                restore
                 }
+ r' i) P5 a' k                 
2 n' s0 B2 _+ C) T: w  L1 m' v+ c5 ^        else {
                restore
, v3 J- I6 A" `5 W2 ?9 }3 i                 di "cannot find pminphddy240"
) l/ n' D9 @- I1 G# r. b' a                 
               
                 sum dx240
6 C$ V" X8 P$ I/ X                 local tempminphddy240=r(mean)
9 X% F) ^/ C& H               
                 
                 }
                 
8 \; `+ v0 s7 l8 X   
) [. f3 \& Q+ W! [   
" Z; f+ v9 N: t8 C, l   
   
# \; o5 Q( Y5 A% ?8 q+ [5 G! P   
) I2 h3 l  _% T; o( ]" v        gen pminphddy240=`tempminphddy240' in 1
$ Q* W$ P$ X0 B# u; l" z  I       
5 ]4 }4 z2 M% e       
! A8 z# t4 M- }7 c        *-------------ddy240pminp右侧低峰值 -----------------
    preserve
! O  y8 R; w8 h, M
' s/ o3 Y8 p: R. n7 w+ |   keep if dx240>`tempminp'
8 K0 q4 F1 M1 B$ G+ E2 F- b. U2 a   
/ z( U: x& }6 V. q% X3 ^- F- q9 Z   keep if dx240<`tempmaxp'
% c: R3 P- D/ N7 L4 ?   
   keep if ddy240<0
# ~7 v  Q; r9 l& O7 u2 Z3 u   
; x3 m- \0 z1 M* l: w- K5 ?   gen temddy240=-ddy240
4 V! {: Z0 t) ?# T   
    count
, r5 A' c" q- Q               
                if r(N)>2 {
  l2 W1 [1 G6 b' C. V* l) j; j                 
                 pkexamine dx240 temddy240
8 |6 \% y* N3 g/ X  a   
) L  w  Q/ l0 y9 V  local tempminplddy240=r(tomc)
                restore
3 r- g. X5 g3 ~1 r6 A! o$ J                 }
9 i7 }3 A5 x* {6 _2 |( J1 y  |                 
        else {
  t4 Z: ?9 b) M, Z/ X               
                 di "cannot find pminplddy240"
* O8 R3 H( f5 g$ a                restore
                 
3 R. T; m0 ?+ X" K+ T                  sum dx240
4 }! L  Y6 W+ z0 G                 local tempminplddy240=r(mean)
7 o8 |' G  b+ C7 \4 X: a+ Y" A                 
' l! `, m& l1 u* V9 \& F7 Y3 N- f9 F# v' t                 
                 }
! F/ ]" L' ^5 U* ^4 O4 W* W( h; R   
) _- N$ y5 s$ [4 _, D   
   
$ k& [% p6 K: O8 D7 {, {, _1 a7 k   
   
. `5 b/ S) L0 M( h- }$ i        gen pminplddy240=`tempminplddy240' in 1
       
        *-------------ddy240pmaxp左侧高峰值 -------------------
       
- ]4 G/ j" V8 P         preserve
   
   keep if dx240>`tempminp'
- M: v# ^# X7 t( h: q   
   keep if dx240<`tempmaxp'
2 U- }" k0 A9 j7 k; o   
   keep if ddy240>0
3 v$ r4 n2 C4 C& k1 c   count
               
! R& w4 z& d7 D: f' C; T                if r(N)>2 {
                 
0 t% }$ o5 h1 d$ o                 pkexamine dx240 ddy240
; g+ a9 ^9 {5 W   
   local tempmaxphddy240=r(tomc)
   
   restore
/ Z! E; e; T; k: |6 J                 }
                 
' O7 x2 O- C5 W0 h        else {
               
' V! ^& k4 T: b* f& K& |! a8 D                 di "cannot find pmaxphddy240"
                 restore
               
6 ?! }; ]) x1 x" e3 X6 ]/ ?                 sum dx240
$ A* x. F8 t' u5 [                 local tempmaxphddy240=r(mean)
                 
) ^: t& W3 V1 t: ^$ W6 T8 Z1 e8 ]( w                 
, l3 ~9 K; `" D' U$ @# \+ c                 }
   
   
   
  N! G3 L& B) B2 z7 M+ D   
+ ~: _. J5 {0 R7 q: X6 Q2 q; @   
3 Q* F! \7 I- [6 L% F   
4 y3 @/ s& D  g+ C, j& q        gen pmaxphddy240=`tempmaxphddy240' in 1
       
       
7 b' \3 `0 {. y- p; i" x        *-------------ddy240pmaxp右侧低峰值 -------------------
       
         preserve
  
9 @1 H& p  ~$ y9 R& m. H% P# n   keep if dx240>`tempmaxp'
  Y3 `; t' U7 A0 I1 q0 [- a   
8 S2 v9 q3 f2 c" l! C. e' m5 b   gen temddy240=-ddy240
1 K. ]  L* i; v" \# K   
4 \' o4 G) u+ n& p; k8 x- `8 {   keep if temddy240>0
. R5 c5 U4 O: J: M% S1 B, ^4 ]  k) b   
' U: R. v9 `/ S1 v3 U8 v   count
               
' R. ?! p3 z% d& P                if r(N)>2 {
                 
5 e9 b5 V5 [5 l; M7 N+ Q                pkexamine dx240 temddy240
   
6 l$ E( V. y. ~! I( b8 R+ J& k   local tempmaxplddy240=r(tomc)
- v$ C1 L0 Q+ R! _6 Y4 q   restore
                 }
                 
& ^. G' u; t7 t, L/ Y7 u. z. w        else {
                restore
3 g1 Y& ]! k- i$ s; @                 di "cannot find pmaxlhddy240"
, `  g! r! E5 c& H/ p7 R/ Z                 
               
$ X1 H: ]! k7 H- g/ T; w* F                           sum dx240
                 local tempmaxplddy240=r(mean) in 1
               
                 
9 M. J( A, J) S- U                 }
   
. G1 x4 M" E) f' w! R* E' R" c   
7 B6 e  j9 w9 Y) \4 b# d   
9 ?$ g9 M, X; ~& j   
   
8 \$ ?  Y5 }& X+ J: t3 G   
        gen pmaxplddy240=`tempmaxphddy240' in 1
       
5 H" g! X7 K  U' }( X9 M2 r        *-------------------------------------
   
   gen price5d=v5 in 1/20
6 t) _# ^0 |; s& S+ W9 R      
   gen price5yy=_n in 1/20
   
% ^" h2 P) u; v: {0 G   gen price5y=price5yy*(dy240max-dy240min)/20
   
   *price5y &price5d
& B1 C4 J* ?' v% H& d& \   
   gen priceny=price5y in 1
   
7 X; H/ Y( ?; O3 s   gen pricen=v5 in 1

/ a0 \( E- v- b0 {+ Z     
# x: ~! G  p) k! d1 f/ _   *priceny &pricen
   
   
   *-------------------------graph------------------
/ z+ W" Q5 ?; I( m/ P9 J  
$ o9 Q' J. }2 e7 W4 L  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)

* ---------可选项目----------------------
- V; H/ P5 K  S% G9 f4 U4 e
}
' ~2 u" T5 D) e4 Y1 a2 s! `end

旅途愉快

awer 发表于 2012-3-28 18:06
) S. n" d) X7 q) Y* n1 Q- c改了一下，这次运算应该没啥问题了

" X; o/ ?$ D9 G, S这段代码是重新复制粘贴到hello.do里面的？

awer

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

对的，就是hello里面的那个

David09

awer 发表于 2012-3-29 11:40
- \1 ~! E1 H5 x# J6 ~5 A5 i对的，就是hello里面的那个

' f! _8 b! R& p8 L. hello 002680
Obs. nos. out of range
r(198);


我运行后，出现这个，不晓得是什么错误了，请大拿们帮助

qqysc

David09 发表于 2012-7-2 21:35
. hello 002680
. j' q% {9 k3 MObs. nos. out of range
r(198);

估计是导出文本格式问题。

qqysc

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

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

David09

3 i& u: ]3 m5 ~0 N* s4 ~* T

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

1 h; l+ ^7 E9 a* t4 P我清除过了乱码的啊....

qqysc

David09 发表于 2012-7-3 17:13
* A" w6 E2 H! n4 V我清除过了乱码的啊....

说得不是一回事。你说的问题可以试试在 002680的文本中，把tab键补齐，导出时不要带表头

David09

qqysc 发表于 2012-7-3 17:56
7 }. P- F" a% E6 ]& ]2 D0 z3 G说得不是一回事。你说的问题可以试试在 002680的文本中，把tab键补齐，导出时不要带表头 ...

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

David09

搞定了，原来是我选了一个次新股的数据，样本数目太少了，引起程序的错误。
感谢楼主给的好东东。

碧水芳洲

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

benhuashu

请问，这一个语句：
gen zf=D1.v5
6 U. G6 ?. v  G& M; z是什意思？主要是=之后的D1.v5，这个"." 是什么意思？谢谢。

海阔

非常感谢；

墨轩书童

$ f1 H+ Y4 s' m5 _/ t+ S
; H% W9 m- ^/ H: ~1 y& J! I7 E感谢楼主！！

海阔

按照版主的步骤操作到了最后出现
" t" h. T7 s, y
; O( n6 @: n9 H  e% z+ M. hello 600031
avolDo not found
---------------------------------------换成其他代买，也是一样；总是出现这个问题，怎么办呢？奇怪

maodeng

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