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

maodeng

wenty 发表于 2011-10-27 04:17
! [" z; S$ O) s( x) n总结一下使用心得
* M' S8 S. ^6 }- S" `) O1）删掉set memory那一行
, c2 i% M1 ]+ |$ p* Y! q2）改路径

! H+ q" w# D7 l1 D; N( ]( j+ t% T记事本打开后，怎么用tab

冰火

认真学习

higer

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

jway315

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

awer

改了一下，这次运算应该没啥问题了
2 o8 n9 y& X# P6 V3 o" y% d
5 x, w% s4 N  b/ x, p: }& N
/ e& U0 B8 S" U8 L3 ^0 Rcapture program drop hello

program hello

quietly {
6 K0 \" L5 \) J2 B- C8 s
) Q& N) v1 S4 x  Kclear

' J' s9 w+ j& W8 Bcd e:\finance\stock\20110101\source
% R4 |$ ?/ U# a2 I8 R" y/ _
1 y' b/ N8 |& [$ n- s% W: ~' vinsheet using `1'.txt

. u) Z7 W$ c" Y8 X5 Q   keep if v5>0
3 A. L) S, f( o0 w3 H* K6 H   
+ W% a( _: J+ d) g. n   gen vt=_n
: h8 I! p. l% B/ e" ]. Q  u   
! m- M1 O' }) v' N/ B   tsset vt

   gen zf=D1.v5
8 A/ N: B- C5 s4 \/ Q& \
   gsort -v1
( J3 ~  B* M( P
   gen avol=v7/v5
   
7 Z1 o6 w/ b4 _6 [( B. I' T   keep if avol>0
   
0 C0 [5 T+ Z( D: ]0 z* X   drop vt
  
   gen vt=_n
8 [' Y" g' J3 m) |   
4 P/ J6 [0 f  Z) M& `5 E3 R   tsset vt
   
7 c- F1 c( `7 Z+ J, T, F    sum v5 [aw =avol]

& |, {& V" ?5 ~   gen amean=r(mean) in 1
- r1 ~. G8 b4 ^) c$ l8 C4 J   
   * 总加权平均值 amean
   
   
   keep in 1/240
) q) y9 l) G( l& V8 {$ A
   sum v5 [aw =avol]
! Q0 i* m0 X3 r2 q' F
   gen mean240=r(mean) in 1
   
   *240加权平均值 mean240
   
% Y/ y* \4 t) L# P! ~+ U4 p7 `   set obs 500
   
   kdensity v5 [aw =avol],gaussian n(500) generate(dx240 dy240) nograph
   
; H6 g+ x: w* a. `' r   *240 kdensity(dx240 dy240)
   
   sum dy240
   
+ N- |# p% g$ K7 R5 T, b   gen yxis=r(max) in 1
" l1 p) w$ I* p9 M( {) G+ B' i7 F   
   gen yxis2=yxis/2
   
   *yxis
   
" P6 m5 n; R8 P3 h# f, o* ]* L   gen dy240max=r(max)
   
8 g5 k! A" C, t   gen dy240minp=r(min)
0 F  G" w9 R! y  F* O4 _% E
1 p. l& ]1 _' ^5 W    gen vtt=_n
1 X+ U0 ~- C) Z. |( o) p- B
' b7 }3 C/ d: N" S; B    tsset vtt
6 v1 k5 N- X/ l/ y* i3 T
9 z  Q& @( z9 Q$ S+ j   gen ddy240=D1.dy240/D1.dx240
   
& o/ ^) s7 E) Z& e2 `( E2 t1 ]   sum ddy240
   
   gen ddy240max=r(max)
   
   gen ddy240min=r(min)
     
/ P, b* O5 @( E# E; r   gen ddy240g=ddy240*(dy240max-dy240min)/(ddy240max-ddy240min)
& ~2 D; ]- V8 l. e   
6 Y; \- j* U, B   sum ddy240g
) [0 U5 `1 M5 a  V   
+ Q  _- b" s0 ^   gen yxis3=r(min) in 1
$ `6 |" Y, t: l; ?* f7 q5 u4 x   
   gen yxis4=0 in 1
  
  *ddy240g
5 q8 P; D/ \; U7 A% f) U   
   
  preserve
  
3 Z2 |  v$ Y' U2 q5 x# f5 x  sum v5 [aw =avol]
% H4 `: S; q# S3 h* r, r/ i. ~  
  keep if dx240<r(mean)
! X, r# S* N7 J& I- k7 L2 a  
+ r' p+ ?5 Z, S# H5 }: A  sum ddy240
  
* l# g* {% F3 M) z5 V  if r(min)<0 {
' j% n. H! H" j( M        
                keep if ddy240<=0
               
                count
               
( H  W$ z- [6 z- d( E" v                if r(N)>2 {
                 
                 pkexamine dx240 dy240
) F: U0 ~( k) K3 ?  O" [* r               
1 _7 Z/ f) [. y5 o8 ^% H$ w                local tempminp=r(tomc)
                restore
) b; [& O* A2 }( [# ]' M* J) f                 }
% S0 Y( z) H! K+ H% A                 
        else {
- m8 ~1 V- F7 _               
                 di "cannot find minp"
$ j7 j9 w( u$ V1 k  V" ?  N  g                 
                restore
2 `! D1 o* t+ p4 ~* i8 a* `; x                 sum dx240
7 J5 _9 u' F& m                 local tempminp=r(mean)
                 
- `' p$ s) l) o, h: M2 k                 }
. m) y( f) A" s0 z& ^                }
! l" \0 k4 n& f4 W   else {
7 ^+ I& T7 D( L- l/ ]                 
                 keep if ddy240<=r(min)
5 e( j! t+ o" z( N0 ?* p( ]& U3 k                 
                 sum dx240
2 {$ {+ w- C1 l% D9 N* J                 
                 local tempminp=r(min)
                 restore
! Z6 C  r2 v+ |2 _: r% p- Q                 }


8 {  r0 A4 W0 z
9 `% R, ?  b( W1 @# C   gen minp=`tempminp' in 1
   *240 minp


   
   preserve
   
. M3 N: c' b; R- Q   sum v5 [aw =avol]
% R7 ]$ q, a3 y. P, I8 ~   
+ A  [7 r) P! \" N% j   keep if dx240>r(mean)
   
! y! X% H0 p" ]0 H7 l: n* q   sum ddy240
2 }5 T6 v1 M- }! r# Z# |3 ^* C   
8 h+ B  P$ i" B0 N' Z   if r(max)>0 {
7 {5 x. L; U1 b% Y- {               
  e1 T/ T. P" X9 q7 @8 C; o                keep if ddy240>=0
               
( f  ]: ]4 K7 W' l2 g: J7 P( e, c  p0 s                count
               
                if r(N)>2 {
  q; G; Q! S' S/ r                 
* ]8 S) s  U3 `! N5 e1 U/ B                 pkexamine dx240 dy240
               
                local tempmaxp=r(tomc)
$ k( s1 o" j: V: u" B. e: L                restore
                 }
                 
3 T8 O1 T" r- D9 L+ K5 H! U                 
! P! o: P( J4 w. A6 D        else {
- E: ~+ l0 j8 ?; U; v                restore
+ r- ]% t* u; R0 Y; L                 di "cannot find maxp"
$ ^; m/ v. V7 z: q$ u4 M               
6 ]4 s3 j2 X" o' A- `" e               
                 sum dx240
* n, Y$ q  C* Y3 k                 local tempmaxp=r(mean)
0 J0 Y, \. I+ w& [+ X                 
                 }
                 
+ C  X+ t' R( q% A5 ^& k2 E7 d                 }
                 
' }; H6 @2 D5 o( }+ L, W3 Z* A    else {
       
           keep if ddy240>=r(max)
8 k8 V. a5 W, t$ U- O          
           sum dx240
* j( D& `( }; j3 \! W          
           local tempmaxp=r(max)
           restore
( X3 ^* q9 k! R           }
2 u; z  T. ~5 R  Z6 ]: D+ X' u       
1 T" k4 A* `) y. ?$ n
% r' S" K1 J6 O5 ^/ \7 V; S  

' r; {, j, `6 {/ M: ^' A   gen maxp=`tempmaxp' in 1
   
   *240maxp
   
4 e2 N9 Y1 a* d1 M! q' `/ v   
   *-------------ddy240pminp左侧高峰值-----
  w/ n8 c8 |. u   
9 u% I& C5 S# I- L; y3 x- G    preserve
   
/ a2 }& o9 G0 A& [6 C% v1 `9 S0 [   keep if dx240<`tempminp'
   
   keep if ddy240>0
9 f' b  y  q' G7 y4 N( T   
   count
               
                if r(N)>2 {
                 
9 g7 k; ~0 D$ t. y8 g1 m! @                pkexamine dx240 ddy240
7 o5 C" R- H% d4 l6 U               
1 @% D' U  A' q) Q+ ]1 g. E. K                local tempminphddy240=r(tomc)
                restore
8 J* D1 N* H/ p! y$ ~/ V7 i                 }
                 
        else {
0 _& g7 @( U: C                restore
5 p4 |% Y! l0 s, r& Q0 v) u! |                 di "cannot find pminphddy240"
                 
               
                 sum dx240
                 local tempminphddy240=r(mean)
5 O3 S6 h& q& v3 @               
1 Q) i* ^2 X, u  H                 
                 }
                 
. \9 [+ x( x3 _/ Q8 U   
   
* N5 j0 d  M# l$ |# t, b3 s( D- a, r   
   
& F. w0 D2 o9 i  u3 m6 A   
/ n7 ~& O8 B+ Y        gen pminphddy240=`tempminphddy240' in 1
# h0 p! R& R' x- A       
       
, m! C" P  H; I+ k& K2 A        *-------------ddy240pminp右侧低峰值 -----------------
4 c. z4 r3 `  R6 b  n/ y3 ]    preserve

   keep if dx240>`tempminp'
" T, x1 g( c- C# R6 {   
   keep if dx240<`tempmaxp'
; T" t' O: F  G9 W5 J6 S   
4 ?- m% P& }  ^   keep if ddy240<0
. H% Z6 A9 `$ o; B9 W3 I   
9 S# p1 G" p0 ~% ^/ K; P) j! i( X7 P   gen temddy240=-ddy240
' }9 C, O7 I9 [% z. c   
    count
' t  X' b# a/ J# ~  r  _+ @               
/ V) X  L' f: F: b                if r(N)>2 {
' g; I6 S8 g( s                 
* a) X5 e4 T5 h                 pkexamine dx240 temddy240
   
  local tempminplddy240=r(tomc)
                restore
' }6 b8 u2 c' s& U                 }
                 
        else {
4 f0 h2 K/ j( O3 \1 T6 Y               
                 di "cannot find pminplddy240"
                restore
9 s2 q- |  @2 h                 
+ u: h: j' c5 O* W) Q( l- r% k                  sum dx240
                 local tempminplddy240=r(mean)
7 g. ]( [8 i' v% C4 W                 
6 S: f# Y2 I# T+ a, F& Z                 
                 }
   
7 Q+ M  v. w8 {, @  \% P   
   
   
   
2 J) x9 h- \. L; L        gen pminplddy240=`tempminplddy240' in 1
       
        *-------------ddy240pmaxp左侧高峰值 -------------------
       
         preserve
% h$ L- n+ ~3 g   
$ {) C, d6 k5 q* ^4 D4 e; q1 d   keep if dx240>`tempminp'
   
   keep if dx240<`tempmaxp'
   
   keep if ddy240>0
   count
) R7 g/ o& D; t               
                if r(N)>2 {
5 W# a$ M5 D$ d# j0 u% S3 R$ `$ g                 
+ H. Z. e. s: |( P                 pkexamine dx240 ddy240
   
0 T2 ?& V( l7 n* w   local tempmaxphddy240=r(tomc)
0 _: S. |1 N* |4 o2 H   
  \6 j0 V0 v7 X- D$ K   restore
  C# \* u2 E- y% i% J                 }
* _1 x! t  q7 o3 ~/ V+ j: z                 
        else {
               
                 di "cannot find pmaxphddy240"
                 restore
               
                 sum dx240
                 local tempmaxphddy240=r(mean)
' K/ s0 M  L3 _6 n                 
& \  {  |; j- b) N5 n9 `                 
; E5 o6 e  T  Z1 f- w! ]                 }
   
1 q+ V6 l/ C/ B' o" I, r   
5 z* `/ I' C. @+ \   
9 s. b: V* ], H) E; ]$ T: E' u0 Y   
) m' k0 m3 N, d" C( X$ P   
: D6 I* M2 R. w" K/ [   
        gen pmaxphddy240=`tempmaxphddy240' in 1
5 Z% F$ \7 @% G       
       
        *-------------ddy240pmaxp右侧低峰值 -------------------
       
         preserve
  
. e0 ~7 d' a2 N) f+ x7 t   keep if dx240>`tempmaxp'
   
   gen temddy240=-ddy240
4 `# D5 R) h8 ^: n6 X- ]" g. Y   
* r  w+ Y  W9 s$ n; d   keep if temddy240>0
7 L& J4 d* f4 Y" s, X   
( C, a& U4 B, `* n* H   count
2 h7 Q4 h( a- N# L               
                if r(N)>2 {
$ K% ~' v' o* V, V& ~7 W: D                 
                pkexamine dx240 temddy240
6 N! Z  ?& M! ]3 s! {   
   local tempmaxplddy240=r(tomc)
   restore
                 }
                 
7 K. M% G3 f) b2 V0 `; p        else {
                restore
% C% l& e9 F& O3 a# {                 di "cannot find pmaxlhddy240"
                 
               
                           sum dx240
+ c/ @- ^, ?' m                 local tempmaxplddy240=r(mean) in 1
3 M- b2 i+ h3 c) d" z& k               
                 
- j' o: N9 b. M# [                 }
   
   
   
   
) N. ^$ g6 B: ?- L$ [' Z   
  \- |- O$ L6 }   
        gen pmaxplddy240=`tempmaxphddy240' in 1
# d7 ~* v# d- `6 q: V8 E       
6 v; m( Y3 B; v2 v. E+ k        *-------------------------------------
, F% o7 \& s$ E" C6 n( L   
. }' g2 |( F, Q% H! x" b: U   gen price5d=v5 in 1/20
4 n" a/ ~7 W0 u0 \+ f9 |* E# O9 t      
) f' D2 V" q+ C0 F   gen price5yy=_n in 1/20
: c" k' W5 y" x% m3 ~$ }   
5 x) u. W8 D" x( y5 b$ f/ \   gen price5y=price5yy*(dy240max-dy240min)/20
   
. U3 c+ \, L- [; B   *price5y &price5d
$ c+ C& k% a, \0 ?+ d2 P5 V   
   gen priceny=price5y in 1
% n& I' m- \. N+ F( H5 |   
   gen pricen=v5 in 1

     
7 }& @6 G" @; B& }" @( ]& ?   *priceny &pricen
  s+ L8 I. L# |4 C4 v   
3 t2 s: `* Q- v5 b" q% r   
' z+ h/ \2 A- u2 K2 M( i   *-------------------------graph------------------
  
4 B+ g* _8 c/ V  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)

4 p) K6 Y& x5 t6 y) _; q* Q% \* ---------可选项目----------------------
. E, f# c8 A7 M8 }7 R
}
: c0 k" \$ D3 L9 K% |9 u4 p) Fend

旅途愉快

awer 发表于 2012-3-28 18:06
- s/ f3 C2 P  [/ |改了一下，这次运算应该没啥问题了

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

awer

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

对的，就是hello里面的那个

David09

awer 发表于 2012-3-29 11:40
" m+ U( j  K) X: v* H# T( b5 f对的，就是hello里面的那个

$ u4 u, K: [8 J- Y) v. hello 002680
Obs. nos. out of range
) h2 |# ~" L* z2 `+ A( Fr(198);
' M$ f" Y$ M' B$ f

7 S; Z5 J  @) P5 P" Y1 Q我运行后，出现这个，不晓得是什么错误了，请大拿们帮助

qqysc

David09 发表于 2012-7-2 21:35
" ?, j* j6 M, f& s# d/ `. hello 002680
. c0 V# h3 P7 g5 S& K2 X) JObs. nos. out of range
: v1 }) R- X- u: D9 fr(198);

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

qqysc

awer 发表于 2012-3-29 11:40
% N  P( A; h1 g+ m8 H  c对的，就是hello里面的那个

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

David09

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

我清除过了乱码的啊....

qqysc

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

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

David09

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

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

David09

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

碧水芳洲

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

benhuashu

请问，这一个语句：
& W$ y) R: d$ e; T* ~gen zf=D1.v5
是什意思？主要是=之后的D1.v5，这个"." 是什么意思？谢谢。

海阔

非常感谢；

墨轩书童

感谢楼主！！

海阔

3 }5 I9 E7 k5 V3 T) o# V
3 Y/ a& x/ D$ L% Q) I按照版主的步骤操作到了最后出现
9 {; U/ E+ a( {  y  L( @' E
. b5 {- a3 ^7 @9 ~. hello 600031
avolDo not found
---------------------------------------换成其他代买，也是一样；总是出现这个问题，怎么办呢？奇怪

maodeng

大师，那两条曲线代表什么啊
$ m) x' |' W  o2 \