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

maodeng

wenty 发表于 2011-10-27 04:17
总结一下使用心得
1）删掉set memory那一行
2）改路径

记事本打开后，怎么用tab

冰火

认真学习

higer

我也曾尝试过使用stata软件对股票平均成本进行分析，但是作为一个非专业的炒股人员，而且日常还有工作的人而言，这种方法似乎显得太专业，而且操作难度极大，不一定实用。

# d6 C* O  U( }! }6 _' V& G+ d使用这种专业的技术分析方法容易造成这么一种现象：大家都对导出的光滑漂亮的图像感觉良好，然而在股市中依然无法进退自如。

jway315

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

awer

改了一下，这次运算应该没啥问题了
4 o# J* Q9 [# s/ J
( u) D) Q$ j% }; S; y: c
capture program drop hello
3 {' W$ x- ?3 l$ }0 T3 {( P; o: N
* `! h% G/ @& G/ }2 k, Tprogram hello
/ `! N, W  g, U$ C0 u% V
: q; P: i4 U* _- W( T3 f0 z# o7 b6 l  V quietly {

clear

cd e:\finance\stock\20110101\source
  ?' S5 \( Y! c4 w
insheet using `1'.txt

   keep if v5>0
   
" v% Z/ b% e5 S, E8 S   gen vt=_n
6 x) t8 ^; s/ i3 B8 {" B   
0 h& I# q% _& w2 S* Y% r7 |( d" I( o   tsset vt
9 B5 [2 [1 N) D. ]4 |6 |  k
, P! E& S1 Z  n# W   gen zf=D1.v5

   gsort -v1
% e1 d' i' O; ]
$ ]. J. h) `" j% ~4 N) O2 f$ E1 G" J   gen avol=v7/v5
   
) t) c* s. H3 M5 |; X2 d8 s   keep if avol>0
   
0 w, f' `7 ~; I7 a7 ]   drop vt
  
   gen vt=_n
   
   tsset vt
   
$ Y- }6 g8 r3 n7 f    sum v5 [aw =avol]

8 {8 Q1 y! [8 `6 \   gen amean=r(mean) in 1
/ f9 C$ x& T+ Q. y% M. B' u: q   
   * 总加权平均值 amean
' a$ E: B& X) `6 {   
   
   keep in 1/240
0 A2 U- F; K2 W' [$ P4 L; K5 G. G
0 v! y2 Q1 }1 Q$ J: i   sum v5 [aw =avol]

  B" G1 s7 ]$ ^+ M   gen mean240=r(mean) in 1
   
! {! \$ h7 ^" ]& }% s4 Z: q   *240加权平均值 mean240
) x- X# r3 X; k: K) h! m) a1 g0 B- g   
% c% h- y% U# s+ `   set obs 500
  j' o6 u" h3 f% J5 Q; O3 p# [; v   
   kdensity v5 [aw =avol],gaussian n(500) generate(dx240 dy240) nograph
   
   *240 kdensity(dx240 dy240)
, s% U" W' q# X- c, S% Z' b! e   
$ H* }7 w3 `1 U1 A6 J* k   sum dy240
   
   gen yxis=r(max) in 1
# F7 F/ B3 u4 q   
   gen yxis2=yxis/2
2 B( p2 P- V* U   
   *yxis
   
   gen dy240max=r(max)
% b3 Q7 {' y+ H: Q2 N* _7 f   
7 n( T0 M5 r8 j   gen dy240minp=r(min)

    gen vtt=_n

; D0 ^' \. W  s( e3 F    tsset vtt
8 V4 _( X8 Q- _# U) R
8 T, s4 q( R; N7 y4 a( d   gen ddy240=D1.dy240/D1.dx240
, n7 [: u3 D- h7 g3 L( Z5 J! c   
" S1 O( C5 `% C  b   sum ddy240
0 G/ J* q) l: c3 r   
% g3 O1 A7 o3 T6 S  {, c   gen ddy240max=r(max)
   
! Y6 S: N& S, P& x- \! `   gen ddy240min=r(min)
     
7 N- _+ o0 V% w" r- ^   gen ddy240g=ddy240*(dy240max-dy240min)/(ddy240max-ddy240min)
   
   sum ddy240g
, ^- Y$ B4 j- o' Z: r3 z   
   gen yxis3=r(min) in 1
   
   gen yxis4=0 in 1
  
: d$ L* X; H$ z1 H# |  *ddy240g
# O& v7 \  ?8 g, ]   
4 [6 ?! |& r8 l+ G3 n, M* \* A8 f   
9 v! }6 I. L& l! a& I; a; a  preserve
  
  sum v5 [aw =avol]
: i5 ?+ _; z. \: W: Y; H  
  keep if dx240<r(mean)
  
  sum ddy240
, i1 j$ Q  z/ s/ S  
  if r(min)<0 {
+ d* E3 W5 K! L: B' b  \9 a        
) r1 _1 @1 K- N& w& X9 O                keep if ddy240<=0
5 h. J1 m' j5 t* e. k7 D               
                count
+ R" i7 K' \6 g7 N2 t0 ^! v               
! E9 Y8 j: U6 t+ M- m" B' f                if r(N)>2 {
                 
5 ~. |) `% H  _7 ^5 W                 pkexamine dx240 dy240
               
                local tempminp=r(tomc)
                restore
) r% Q+ \0 U1 [                 }
* E' x+ S6 ^" Y! A$ f                 
        else {
               
- c! O3 I8 V, [4 K: e                 di "cannot find minp"
- W( u+ ?  ~% \5 G/ B) M                 
) W8 S4 }0 T: ?; B5 Y: C) }                restore
                 sum dx240
5 X" \* i( L$ d% n, |1 ~                 local tempminp=r(mean)
                 
                 }
                }
   else {
                 
                 keep if ddy240<=r(min)
, r7 ~( N. ]! S% f& X% c4 K/ u                 
; P# v+ ~  Q$ D% J1 T! q                 sum dx240
' _* E$ w0 I5 n+ k0 ?/ d                 
                 local tempminp=r(min)
! b. a+ W, f: ?( M4 m& _# J) H                 restore
                 }


7 ^7 x7 h( a4 y8 D9 k0 X/ Q
   gen minp=`tempminp' in 1
) ]: ^2 t) W  n6 X$ d# x   *240 minp

" \1 e2 H0 N3 h6 Y6 q. H+ h: ?6 y' o% Q
   
   preserve
' W, A; ]7 \5 C: d. i5 K0 z   
  ]3 P) Y) T+ ?& w; s" I   sum v5 [aw =avol]
   
   keep if dx240>r(mean)
2 f+ e+ w5 [' s  N5 x   
# I6 W1 I& h0 b/ |% a   sum ddy240
   
   if r(max)>0 {
# W- l5 ]) D& L/ I& k' A$ V               
                keep if ddy240>=0
               
                count
               
% T/ N& h9 B* A4 {& w9 y                if r(N)>2 {
                 
                 pkexamine dx240 dy240
               
                local tempmaxp=r(tomc)
$ I/ f) h" T' ?/ S# n0 F                restore
. H' X! N. g1 Q% w& V/ s                 }
  ]. u) C$ n! a# V% q                 
* B4 q0 X9 y2 o% [1 R                 
, w2 s( K+ ]! J9 A/ j8 x% c        else {
7 s( K3 i$ G# `6 Z. R' f  m5 M                restore
& f& [: e$ [* @; g                 di "cannot find maxp"
               
               
2 f4 A+ [- q! U                 sum dx240
+ H0 C' u/ m3 ?  p                 local tempmaxp=r(mean)
7 j$ t% C, d4 u& B0 Z, s; C  S                 
                 }
                 
8 K, W2 V& ?6 D. s" ~- {, M5 Q                 }
                 
    else {
       
  w1 K; \* t9 |: x! X' l5 K' b           keep if ddy240>=r(max)
          
           sum dx240
          
7 s6 ^" E& |0 s6 I" a% I! w2 k           local tempmaxp=r(max)
" H, y- r  b, \: G' o% S           restore
           }
2 b/ ^7 D/ r& _2 z0 U- e       

+ t2 {* b& V& e  i% ?/ S  
1 J6 j" l/ y5 k, P- j. u
   gen maxp=`tempmaxp' in 1
+ P5 m. o: |, k9 }  w0 N; X   
1 ?4 U  e& j/ L; a   *240maxp
* f/ R( L& N7 t$ {7 q/ f" X   
8 n. l/ `* X# y8 L, q& Y8 M   
   *-------------ddy240pminp左侧高峰值-----
   
8 Q$ p& n2 O" l8 z2 u' n    preserve
6 t& Y( v1 i' j& V& f% @   
   keep if dx240<`tempminp'
   
   keep if ddy240>0
/ C" s8 f& k* l. U' F0 n7 Z   
   count
& ]; k9 Q0 d$ P               
% ^! `) V3 Q6 K3 d2 w& @; _! ?( k                if r(N)>2 {
* q* ~: Z" r; _  P# R) x" ]                 
                pkexamine dx240 ddy240
: l0 p8 ^4 F6 p3 Y* t  S0 I               
                local tempminphddy240=r(tomc)
! B8 Y- B! \3 S- a- y6 y                restore
                 }
                 
        else {
                restore
& t- B! d! c% q0 F- C: H                 di "cannot find pminphddy240"
3 }7 E  b  |; R2 |# W                 
               
3 G$ u0 v& D  x. T                 sum dx240
( l: c% y2 T# g7 i. H/ h                 local tempminphddy240=r(mean)
               
7 y- Y+ E% c9 ^8 j% _                 
+ ?" z/ K. B. m# K# J4 Y) z( [                 }
0 v0 N, b4 k& o0 N% y                 
   
) F( l( A3 J6 j   
   
   
   
9 p' P1 _' |  S+ A5 J        gen pminphddy240=`tempminphddy240' in 1
       
       
        *-------------ddy240pminp右侧低峰值 -----------------
    preserve

; [7 c# v, Z  ~- k5 u   keep if dx240>`tempminp'
   
! T" w" |- I+ }7 B0 [* \   keep if dx240<`tempmaxp'
   
# q3 c" R" l! q* l: x   keep if ddy240<0
   
   gen temddy240=-ddy240
0 i& _  D+ f5 b7 o! N4 {   
* T7 q# s7 c+ b    count
               
                if r(N)>2 {
                 
1 u; {6 a% W( i  n                 pkexamine dx240 temddy240
) i! D7 o$ d7 M" r   
  b$ }/ {- i* R# g3 M2 a  local tempminplddy240=r(tomc)
" w+ d! c5 J: O! l3 _                restore
                 }
! t; s, a" H5 V4 y% K5 s. ^                 
' O' e8 G# H$ U4 G: c# O        else {
               
; u5 f: ^) Q' n% {" E                 di "cannot find pminplddy240"
                restore
4 y; i1 Z; q5 W  k, t3 m+ k                 
                  sum dx240
) u4 X7 c" a3 v0 R" N                 local tempminplddy240=r(mean)
7 V& _9 i% K5 k+ Y0 E                 
                 
                 }
# R2 D9 k4 Y2 G/ X" T" s   
3 k9 `0 P, \) k" q+ F- m: Z, r) B   
   
   
   
  \8 K/ d6 V  _4 N        gen pminplddy240=`tempminplddy240' in 1
       
3 d- g9 d6 u& ]        *-------------ddy240pmaxp左侧高峰值 -------------------
7 Y9 j6 Q% ?  j6 O( M       
         preserve
, d- d0 ~5 ?1 Y+ f  s- `; N   
! O: g1 j: T# ]5 L  q0 ?# L: s* t   keep if dx240>`tempminp'
* N) q) ?. G7 K# L8 f* N   
   keep if dx240<`tempmaxp'
   
   keep if ddy240>0
4 f3 m/ z: Z  M. ?# m9 F/ z   count
( k) d# I1 [, W, M& @8 t               
2 y! I7 a  w% C3 l. \3 s( W4 L, v6 K/ G                if r(N)>2 {
                 
                 pkexamine dx240 ddy240
! v# k  h3 }; j0 y3 \0 f% E/ G   
7 y/ e- u& r& `1 S1 y   local tempmaxphddy240=r(tomc)
   
   restore
                 }
1 Y! O9 G& }1 w. D. ^                 
        else {
               
                 di "cannot find pmaxphddy240"
7 u/ [. g' Z1 J' f+ b                 restore
3 g/ J# [5 J' j! O+ I/ K- t               
2 ~5 h- H+ G/ T% y2 q8 M) z& l$ u                 sum dx240
2 M: h3 F1 a7 {                 local tempmaxphddy240=r(mean)
& [& n$ P  [9 e' E: x2 o9 l+ \8 Q  ~) E                 
7 D* t! K6 H+ Y  F1 u0 n                 
' J( g1 j7 K8 F0 g                 }
   
* @0 P, S$ B, D" Q- L   
% n0 V+ h7 T! Q1 b5 [7 Y   
$ O: _: G: J: o   
   
1 n% F% S3 }, o7 n   
) R! B0 [$ R9 c: `. k+ g        gen pmaxphddy240=`tempmaxphddy240' in 1
1 Q( h9 e& K0 v6 o4 L' ~$ L       
       
7 j9 r3 o$ |  x0 }$ q& ~        *-------------ddy240pmaxp右侧低峰值 -------------------
       
         preserve
! V5 S, T. c* k! Z( ?  
   keep if dx240>`tempmaxp'
   
   gen temddy240=-ddy240
   
   keep if temddy240>0
   
7 m% O8 a# R% f1 I3 V   count
               
                if r(N)>2 {
! ^* `( T% _# u                 
                pkexamine dx240 temddy240
1 v0 h) w1 h# n   
# [0 q4 X) z9 m   local tempmaxplddy240=r(tomc)
   restore
' j6 H0 r1 B$ y                 }
                 
        else {
4 F! `2 P; p" A  o+ V                restore
) r8 Z1 s/ w& Q: A1 Y& d                 di "cannot find pmaxlhddy240"
, m& L  z! G; z: G                 
               
                           sum dx240
1 d# V' ^% v* I+ ]                 local tempmaxplddy240=r(mean) in 1
6 b, ~) j! v2 E, @" q$ l1 f               
& H" J7 O  |* U* F                 
( h9 P- F; t. `5 a7 t) t  Q                 }
   
( {$ _4 T- b2 F1 u   
   
. b/ {! ^6 g2 r0 i: y' f   
   
   
7 `1 N% Z- u; E& Z  g        gen pmaxplddy240=`tempmaxphddy240' in 1
       
+ \3 U: l$ \3 w8 F, e* D        *-------------------------------------
+ [6 z# W" y8 B" {   
   gen price5d=v5 in 1/20
9 L  Q( }, b" l! U. c      
5 K" y8 }( `9 e( `( B7 g6 g   gen price5yy=_n in 1/20
   
   gen price5y=price5yy*(dy240max-dy240min)/20
   
" A7 N8 r5 A1 c& Y& G, u   *price5y &price5d
! T& z+ G) K0 ]  q, ^# L" s   
   gen priceny=price5y in 1
. |7 t! W! ~+ r! X9 Z  U   
" i6 d. L, [# [   gen pricen=v5 in 1

     
! S' c& w. ~1 E   *priceny &pricen
, |. ^$ m' L+ N& b7 {- O6 l   
& U6 M# g- c/ G) }, O   
   *-------------------------graph------------------
  
  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)
2 C3 x5 h+ M- v2 I
* ---------可选项目----------------------
2 D" \8 V+ k5 F2 [8 s$ y, o
}
end
7 n$ P) _; a6 u4 f

旅途愉快

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

  ]5 t  ~/ F* u+ R1 S) d8 d这段代码是重新复制粘贴到hello.do里面的？

awer

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

对的，就是hello里面的那个

David09

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

. hello 002680
7 T. Q5 }9 K" qObs. nos. out of range
r(198);

) p8 _. }* t$ o- V" z4 I
8 z% E6 `2 o7 I) N+ b我运行后，出现这个，不晓得是什么错误了，请大拿们帮助

qqysc

David09 发表于 2012-7-2 21:35
  g. o% j, H# ]/ @. H. hello 002680
& ~% R$ M# H+ X; e5 A0 vObs. nos. out of range
r(198);

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

qqysc

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

# E9 V. v7 k- z9 a复制粘贴有很多乱码，请楼主把上面的代码做成附件以供下载，谢谢！

David09

qqysc 发表于 2012-7-3 12:14
  ~' L( ^" |. q7 F9 A0 X# d  d复制粘贴有很多乱码，请楼主把上面的代码做成附件以供下载，谢谢！

( U* b9 Z! j+ D1 k! R
: u( A9 ]- o# I! v* X& \我清除过了乱码的啊....

qqysc

David09 发表于 2012-7-3 17:13
我清除过了乱码的啊....

* d' w7 L4 g* Q& A/ g; n6 j说得不是一回事。你说的问题可以试试在 002680的文本中，把tab键补齐，导出时不要带表头

David09

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

* ~' I! a& I' ^# R- p/ J* Q很奇怪啊，我看stata中间的数据表，都整整齐齐的导入了成交价成交量。但是一运行就蹦出来那个错误

David09

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

碧水芳洲

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

benhuashu

请问，这一个语句：
gen zf=D1.v5
/ m) N3 e, u) \# N! E; X, U( ?, b是什意思？主要是=之后的D1.v5，这个"." 是什么意思？谢谢。

海阔

非常感谢；

墨轩书童

% i2 y7 J; v+ B6 Q+ v9 V$ d: O
' k3 \, q4 j  H' \0 o; B感谢楼主！！

海阔

1 S$ ?' z6 {- s7 {! f3 l7 O
0 \# Y$ r8 A- d3 D! q4 l0 ]按照版主的步骤操作到了最后出现
& ?1 k- t" A# S3 g3 w9 n
5 A) V) P9 p+ N) d. hello 600031
avolDo not found
---------------------------------------换成其他代买，也是一样；总是出现这个问题，怎么办呢？奇怪

maodeng

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