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

maodeng

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

记事本打开后，怎么用tab

冰火

认真学习

higer

我也曾尝试过使用stata软件对股票平均成本进行分析，但是作为一个非专业的炒股人员，而且日常还有工作的人而言，这种方法似乎显得太专业，而且操作难度极大，不一定实用。
! [! k6 H' L4 w! l- _/ b5 d9 w' }1 d
( f+ Z$ ]# \3 |1 r5 l使用这种专业的技术分析方法容易造成这么一种现象：大家都对导出的光滑漂亮的图像感觉良好，然而在股市中依然无法进退自如。

jway315

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

awer

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

$ |! |8 Z2 o2 _9 T5 ~5 [% I3 h( k
+ j' x$ }# n/ Z& `. x0 |capture program drop hello
- o* t' K/ x5 f3 L5 [. y
program hello

3 p( f* T& `% A; S/ r! t" }% L quietly {
: V5 G2 a& d& }) q
clear
- a5 F1 V2 E$ F) l/ ?
cd e:\finance\stock\20110101\source
2 L, A. ^( H  _8 O( j
insheet using `1'.txt
* `! a! g. E8 {
   keep if v5>0
   
   gen vt=_n
$ |4 x7 B1 Z" N( x- Y0 [3 t' C# {   
! X( j! I4 W. Z) l5 @3 i   tsset vt
: I# }+ P7 }1 m/ S4 ]$ _
) I* p" o1 ~0 f$ L, i& m   gen zf=D1.v5
- s/ {+ H% \! Z5 u+ @
   gsort -v1
$ Z$ g5 N5 f8 o/ X4 H
   gen avol=v7/v5
3 }  m- `5 [" \: Q  Z. M3 h   
   keep if avol>0
   
   drop vt
  
# b' q& g8 o, V& |( ?( R, w3 f   gen vt=_n
   
* C7 g4 O+ F1 M+ p4 \$ n   tsset vt
   
    sum v5 [aw =avol]

   gen amean=r(mean) in 1
   
   * 总加权平均值 amean
; n* g  T3 C% ?   
2 G% }/ t; I$ d3 k4 k" {' e, x* T) X   
3 y& Y! c$ s) y3 e. r8 f0 x+ f% l   keep in 1/240

2 p0 N, [3 X1 T4 l4 t5 m4 a" f   sum v5 [aw =avol]

   gen mean240=r(mean) in 1
   
   *240加权平均值 mean240
# \* e8 p+ [: \+ G   
' e1 c/ }$ M) l- M   set obs 500
   
   kdensity v5 [aw =avol],gaussian n(500) generate(dx240 dy240) nograph
   
   *240 kdensity(dx240 dy240)
   
3 `; `- A4 J, k9 \' q# [) a   sum dy240
0 w; `3 A0 o9 A3 y  C2 v   
; }8 y: W2 V. |% D  f   gen yxis=r(max) in 1
5 x, @; p$ N/ U: b- C. ]8 q   
& z, G$ I& E, z  s  i- C   gen yxis2=yxis/2
" n8 ?. @" \9 G2 b   
" @# ~2 |4 f2 |0 {  U6 ]) {   *yxis
6 l1 ^$ X6 T% ]8 s) M; g# m3 N8 E   
   gen dy240max=r(max)
   
0 o2 W" Z7 l4 E7 W9 r5 [   gen dy240minp=r(min)

    gen vtt=_n
# w3 Z, b5 b* P
# H' f' V% }$ F: q$ V, @" j/ h    tsset vtt
. u) P- C/ K8 n% M' b1 g
   gen ddy240=D1.dy240/D1.dx240
8 c. D# ^- Q6 y8 H" d   
   sum ddy240
3 w6 o: o* l1 }# G7 [   
   gen ddy240max=r(max)
   
   gen ddy240min=r(min)
" N/ O8 E3 R( _" m5 `     
   gen ddy240g=ddy240*(dy240max-dy240min)/(ddy240max-ddy240min)
   
' u+ }; \( [) Z   sum ddy240g
+ b' w. n) e: W   
. f0 D$ C! b& Z: a6 ?   gen yxis3=r(min) in 1
$ F0 r8 T) w3 m9 I( i   
/ x, l: I' @% q# G1 c) P* U   gen yxis4=0 in 1
  
  *ddy240g
2 B) q! X. {! h% g   
6 Z# G4 S' {4 }% T' \* C0 N   
  preserve
8 [$ M5 T* U6 \. M3 p8 F  
  sum v5 [aw =avol]
. O5 k2 z3 Q0 A6 O" B  ]  
  keep if dx240<r(mean)
  
' j5 \1 \" V+ B, ~  U# I5 M  sum ddy240
  
8 U4 N0 e2 F# ~7 a0 z" x  if r(min)<0 {
        
                keep if ddy240<=0
. K; O" N+ v  \) g5 {) _               
                count
               
                if r(N)>2 {
                 
7 D$ w8 W3 U3 I) ~: o  V                 pkexamine dx240 dy240
               
& _$ z1 H' Q/ ~! I; ~                local tempminp=r(tomc)
                restore
+ ~9 F0 l+ J# V: M( l                 }
                 
4 g: n) G  C( L! Q% H        else {
               
                 di "cannot find minp"
7 H& N3 I9 a7 S5 s, d                 
) f1 A, S) ^0 k- S% `! P* r                restore
                 sum dx240
& N8 z) b) W3 h8 O( ?                 local tempminp=r(mean)
                 
2 h& V2 G- \9 ?" w- m* h" m                 }
: m" T5 N$ y0 W8 t" B4 S+ t                }
) D  K# U7 Q: n& C& {* m% y1 z- m   else {
- t9 t, f0 g% o- \. e                 
                 keep if ddy240<=r(min)
                 
                 sum dx240
" y! d# F, B  P) B: F( ^, |+ J                 
                 local tempminp=r(min)
8 ]7 _. W# D+ K. [                 restore
6 V9 Y. K: S. e                 }


# j# S$ T' W% o6 T4 @
  P( S- y/ F# l   gen minp=`tempminp' in 1
   *240 minp


   
   preserve
8 H7 b2 E5 G7 k   
   sum v5 [aw =avol]
/ y+ n! L: Z8 W( @   
; j) O  c$ r+ ]: G$ n+ V   keep if dx240>r(mean)
& Q, o$ i  j6 G/ N   
   sum ddy240
+ {- a+ C" z2 P" ?/ P! w, h+ z- J   
: u( n; r0 t. u5 f# C& F   if r(max)>0 {
               
                keep if ddy240>=0
7 g( J* n1 F+ X               
- Q! [! S4 n( _2 y- {8 R                count
9 a4 Z/ X0 l9 }               
- H( F8 M3 ^: R* e                if r(N)>2 {
                 
& V: z; ~4 i! _                 pkexamine dx240 dy240
* i! E* g  g) }7 N7 p               
                local tempmaxp=r(tomc)
                restore
5 C# _; W! t. B) E                 }
) n/ c1 E( {: k                 
                 
+ q$ {1 d1 F+ J- y" Z( h, D        else {
                restore
                 di "cannot find maxp"
               
               
3 J1 D  S3 s) V! C' L- t. V/ p                 sum dx240
                 local tempmaxp=r(mean)
                 
& _9 y0 g) c8 p                 }
                 
                 }
                 
    else {
, C) M; x) V! d3 m) u8 _% [1 h       
5 i+ e: G3 o7 f/ P/ T6 S+ O) W           keep if ddy240>=r(max)
          
           sum dx240
" q/ R( H+ h* N' ~+ Y0 X          
           local tempmaxp=r(max)
           restore
7 M* S0 S, j  Y8 x/ l9 }           }
" H( }8 V7 g+ a, p! Z) N- e6 y       
( B, Y9 W* }7 N6 C/ K
. s& W1 A5 s8 x6 f' j0 N  
3 B* \4 r5 y7 d& k. N
) B3 @6 c8 a" N4 w7 N# ]" V   gen maxp=`tempmaxp' in 1
1 H) s; z7 F4 M1 H   
   *240maxp
   
5 r! B- a1 O3 ^5 B9 h) `   
   *-------------ddy240pminp左侧高峰值-----
4 J  S' L- E8 C6 d* L% U9 s: U   
" ~4 l2 v( H/ ^+ ~( F    preserve
   
  M$ }% t" ~2 \' |+ z  ~   keep if dx240<`tempminp'
   
( V; ^( D3 k! G* T! r3 H   keep if ddy240>0
   
   count
               
7 m" B! `+ f6 o+ O: ~. J) T" T. p7 [$ P                if r(N)>2 {
$ i% n/ L/ Q6 X' P6 k/ O                 
                pkexamine dx240 ddy240
               
                local tempminphddy240=r(tomc)
0 F( l4 J5 I+ A; ^3 c                restore
" Z! ^' K" e$ `2 O: ]6 b. F                 }
                 
6 T4 ], e5 j, V        else {
                restore
& B: d, q, N' @                 di "cannot find pminphddy240"
                 
               
                 sum dx240
                 local tempminphddy240=r(mean)
$ K+ i" E% l7 h5 Y& M               
                 
                 }
                 
8 V6 P* w# J' `! n, i   
* ^3 n6 d+ t; ?; n8 ^   
   
   
' c- T% O9 h4 J% c: i7 T+ P   
        gen pminphddy240=`tempminphddy240' in 1
2 H3 q0 g+ A5 S0 f" ~% B$ D. {$ x       
8 b# ]0 F" W9 B6 c% b( R: [& r       
4 t( K. G( ]2 q2 L$ |2 U        *-------------ddy240pminp右侧低峰值 -----------------
$ r+ i  j. _, @' s2 K    preserve

   keep if dx240>`tempminp'
) H! O! d' }. z$ u0 q& K) v   
   keep if dx240<`tempmaxp'
   
   keep if ddy240<0
9 D( F& h# Z5 x   
" \( H: c9 i; G$ i  g5 e' l   gen temddy240=-ddy240
   
    count
               
                if r(N)>2 {
7 s: l* n' r0 {/ i- L                 
* ^3 `3 V+ t4 m  R                 pkexamine dx240 temddy240
   
  local tempminplddy240=r(tomc)
& ^1 w) }) J0 _. o# v8 E- ^1 D: A                restore
2 P, t- y: F4 s                 }
                 
+ f9 s% E) r' W2 t" L* r% e1 m* x        else {
* D5 V; c6 c! }: f" U               
- U7 E# ^( e7 c4 ?) O6 w5 v                 di "cannot find pminplddy240"
                restore
! k, Q1 C  X2 C" E                 
  Z8 z3 [) q/ c* W1 Y                  sum dx240
0 `$ L& U/ |  @& B: L6 o                 local tempminplddy240=r(mean)
                 
# G$ l/ R2 f8 }: J- }/ R! {                 
                 }
   
   
8 R  b" b  \# J$ p- ]6 P4 C! b( e   
   
   
        gen pminplddy240=`tempminplddy240' in 1
3 X5 F. J$ Y- q# b) }$ P- t       
  H! a  Y! H" d7 J4 z4 P6 P; r        *-------------ddy240pmaxp左侧高峰值 -------------------
- e; f5 ~8 {2 _6 z, Z       
2 E/ z9 q4 L0 k7 S, F# ~9 u' ^         preserve
& g2 _; b' L- h- `1 U" w   
2 w: S$ i9 T' V" z& ]3 m$ b   keep if dx240>`tempminp'
   
   keep if dx240<`tempmaxp'
   
! p' s  F% k" j( d; G   keep if ddy240>0
5 t: q+ B( D; ]2 l' {' ~1 m   count
               
                if r(N)>2 {
                 
& T4 H- B* l' e' ]                 pkexamine dx240 ddy240
   
   local tempmaxphddy240=r(tomc)
   
   restore
                 }
                 
        else {
& I/ x/ `, D7 V2 m( w6 C0 n, x               
                 di "cannot find pmaxphddy240"
8 J3 q$ g. f* I* B% N* ]4 C' N                 restore
               
8 n3 i9 _$ Y( o6 }                 sum dx240
                 local tempmaxphddy240=r(mean)
9 Z# E0 Q  z* Q  U/ c/ y                 
                 
                 }
   
   
   
   
   
   
        gen pmaxphddy240=`tempmaxphddy240' in 1
       
: ~: {. Y3 t: z0 p  m% d. f       
$ ~7 Y( y- V& {2 b' ?$ p# [        *-------------ddy240pmaxp右侧低峰值 -------------------
       
. ?/ i* |# c0 a( |6 ~7 J         preserve
+ k) B1 \; i# ~1 z% d  
   keep if dx240>`tempmaxp'
   
3 u% }4 l% c* o5 O* i   gen temddy240=-ddy240
   
   keep if temddy240>0
4 _1 z3 Q& D0 D+ E; N1 _  T% V; u   
- I, x' s  W  X0 W6 }   count
! _+ B' t+ I/ m: N* M* D7 U               
                if r(N)>2 {
! @+ \3 ?  q3 L* y/ ]8 t3 z                 
                pkexamine dx240 temddy240
( e; I7 v5 Z  u1 K* K   
/ x/ {' [' a% G   local tempmaxplddy240=r(tomc)
   restore
) F* T& d& q3 ]" o                 }
' }. ~+ D/ F" S" U/ S3 h- n& {0 |) D                 
        else {
& }! q6 F2 L1 }; v$ ]5 s- J. G' y; G                restore
                 di "cannot find pmaxlhddy240"
                 
3 S( h7 W0 \' o1 E) _5 K               
                           sum dx240
3 g4 F: z2 }) q8 ]; l5 U1 k3 b                 local tempmaxplddy240=r(mean) in 1
  F$ V0 f7 N* X+ \$ T% Z               
7 Q# o; T2 ]* {9 e% H4 Y                 
  p/ W8 ]+ c1 i2 ?8 i& |                 }
   
   
   
, o; s; G/ l9 a  R9 q2 D" R   
4 v  U1 r7 A# \5 N: O   
   
        gen pmaxplddy240=`tempmaxphddy240' in 1
. l/ X4 }( L% M# z$ L) E1 L; w; U9 ?0 s       
        *-------------------------------------
   
   gen price5d=v5 in 1/20
, D% N- J% g1 W) ?9 W3 c1 h2 U      
   gen price5yy=_n in 1/20
8 L) o3 T1 R) ?+ s) R& p- i4 ?+ X   
9 X2 A% |: g- J- U+ X$ g$ a3 W$ A   gen price5y=price5yy*(dy240max-dy240min)/20
) ~! M% \/ ], z- l   
   *price5y &price5d
   
  M1 \) n* B2 @( @) U8 b   gen priceny=price5y in 1
   
  d4 q5 T* r! e% n6 L3 q   gen pricen=v5 in 1

     
- i5 u9 _) j  ]8 f   *priceny &pricen
5 _6 y  l& r0 ?5 h9 w   
   
   *-------------------------graph------------------
6 j& o6 S: `& o* a* H  
8 u* ~$ J% G/ I* q1 J4 D$ H  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)

* ---------可选项目----------------------
9 w$ S2 I% z& h. [" G- L
}
end

旅途愉快

awer 发表于 2012-3-28 18:06
6 v! i) P9 ]; g- C9 _5 g: I改了一下，这次运算应该没啥问题了

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

awer

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

对的，就是hello里面的那个

David09

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

8 O' _" w$ Z+ z+ l. n8 O. hello 002680
Obs. nos. out of range
6 l! b- E: I2 K9 x" a: b* qr(198);
* f- n) q: J# b6 {* @4 T  L. m. t' S

" i; y6 t! y5 W  w8 h& I+ Q! p我运行后，出现这个，不晓得是什么错误了，请大拿们帮助

qqysc

David09 发表于 2012-7-2 21:35
0 [: `* H, E, B. j  U* q. hello 002680
) p) L$ a" }3 e* LObs. nos. out of range
r(198);

- e' c* j3 Q  V& q估计是导出文本格式问题。

qqysc

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

# K/ v3 e8 I0 {5 u+ Z; r( a复制粘贴有很多乱码，请楼主把上面的代码做成附件以供下载，谢谢！

David09

qqysc 发表于 2012-7-3 12:14
9 \, s8 n: m6 |5 ?# v复制粘贴有很多乱码，请楼主把上面的代码做成附件以供下载，谢谢！

$ [% A8 Y& n8 V" \8 q# t
/ X; r. c9 `+ x8 G, n$ i# L我清除过了乱码的啊....

qqysc

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

! }/ I7 T; L- T9 j# e& P说得不是一回事。你说的问题可以试试在 002680的文本中，把tab键补齐，导出时不要带表头

David09

qqysc 发表于 2012-7-3 17:56
# h, i) k* o# J- Q8 Q, N说得不是一回事。你说的问题可以试试在 002680的文本中，把tab键补齐，导出时不要带表头 ...

" z. ]& |4 v0 E% Z很奇怪啊，我看stata中间的数据表，都整整齐齐的导入了成交价成交量。但是一运行就蹦出来那个错误

David09

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

碧水芳洲

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

benhuashu

请问，这一个语句：
gen zf=D1.v5
是什意思？主要是=之后的D1.v5，这个"." 是什么意思？谢谢。

海阔

非常感谢；

墨轩书童

6 }, U& J5 D" V; B9 Z; q感谢楼主！！

海阔

' F8 z! C0 \. U- ^
  w# y% E  w- b' D! C' c按照版主的步骤操作到了最后出现
% P$ N6 E7 n8 n- N4 ~* _) \! M5 m
. hello 600031
* J/ H) O; B* h+ J8 l: u3 UavolDo not found
---------------------------------------换成其他代买，也是一样；总是出现这个问题，怎么办呢？奇怪

maodeng

大师，那两条曲线代表什么啊
! }0 V* b' i4 U