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

maodeng

wenty 发表于 2011-10-27 04:17
总结一下使用心得
" G; ~) K7 k6 K# k5 x1）删掉set memory那一行
" f. U; S7 p( M2 G2 X% W0 e) ~* x2）改路径

. g4 e" x- ^3 b记事本打开后，怎么用tab

冰火

认真学习

higer

我也曾尝试过使用stata软件对股票平均成本进行分析，但是作为一个非专业的炒股人员，而且日常还有工作的人而言，这种方法似乎显得太专业，而且操作难度极大，不一定实用。
, y. V: C! ^- \: }# m7 |0 v# N
. \1 }# z0 F# x9 A9 _4 c使用这种专业的技术分析方法容易造成这么一种现象：大家都对导出的光滑漂亮的图像感觉良好，然而在股市中依然无法进退自如。

jway315

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

awer

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


9 n; X5 n. h* X0 L4 D/ t# A% scapture program drop hello

7 Z6 U. W. Y" @% `5 xprogram hello

$ l# w+ A; R) x quietly {

clear

cd e:\finance\stock\20110101\source

* g4 y& \5 A" h# _) J  Uinsheet using `1'.txt
4 E2 o1 L2 a7 z( z' Z% l
! L# i3 J$ |# f6 f   keep if v5>0
   
   gen vt=_n
   
   tsset vt
  v3 A$ b7 |1 e( i; I- E8 ?8 r
   gen zf=D1.v5

   gsort -v1
4 I( c- F( ^* ?9 l( r
. `5 t# X% _" @  |$ G$ t   gen avol=v7/v5
3 S/ g; R) p, p4 ~# D/ Q   
' v% `: m% G. F' }- M) K; l   keep if avol>0
   
( J! b8 |  Z( ~, S% b6 ^   drop vt
  
   gen vt=_n
   
   tsset vt
; l6 O/ }7 B1 W   
    sum v5 [aw =avol]
- V6 `- f( T3 `5 w; K
   gen amean=r(mean) in 1
) j9 W' `6 x/ n, Q( g6 n   
   * 总加权平均值 amean
   
5 Z$ u/ ^. O4 E# f+ s' u" k   
/ t7 G# Y" X5 b' Q   keep in 1/240

   sum v5 [aw =avol]

   gen mean240=r(mean) in 1
   
0 P! u6 s8 q% C1 B& d: x   *240加权平均值 mean240
   
8 H+ t9 A! y3 w! e   set obs 500
   
2 M  {: ?8 f0 I2 \4 z7 F   kdensity v5 [aw =avol],gaussian n(500) generate(dx240 dy240) nograph
/ ]" |* y  k& Q6 s3 z/ S, }' K   
   *240 kdensity(dx240 dy240)
   
5 e2 e) r  @# J   sum dy240
: I* C- ~& j% c- R) \- @   
   gen yxis=r(max) in 1
- v+ A- ]8 w# a- j   
   gen yxis2=yxis/2
# r: W' f9 U  E: {1 o   
   *yxis
6 w, Y7 T# |, h7 l1 y, k   
   gen dy240max=r(max)
! v% j7 G+ }. D% ]1 Z   
: ^5 C! A$ l" s* V  c   gen dy240minp=r(min)

    gen vtt=_n
+ z6 w+ P% Z3 N( Q
    tsset vtt
) ]; R/ U9 I$ I
   gen ddy240=D1.dy240/D1.dx240
   
/ c, i, d3 b$ _. F   sum ddy240
   
   gen ddy240max=r(max)
3 a) U9 F: s# i   
/ c2 E* J7 f. r; K1 [5 @" [8 [   gen ddy240min=r(min)
& Y! `* \5 ~* ]     
   gen ddy240g=ddy240*(dy240max-dy240min)/(ddy240max-ddy240min)
   
# R/ X9 x; q  S   sum ddy240g
" L9 _6 F/ l9 y   
$ M/ U) W/ d( M: N+ S   gen yxis3=r(min) in 1
' m# ]% X7 f' ?   
   gen yxis4=0 in 1
' r' S$ m( v7 I2 E6 `4 \( T  
0 W4 |% N& Z- I  *ddy240g
8 s! a3 R, J7 \3 c( q4 {  v   
  q' \: Q; A$ W/ `6 ^4 L   
: r' b; `7 _. J3 Q+ X% g  preserve
  
  sum v5 [aw =avol]
  
  keep if dx240<r(mean)
0 ?, h1 G8 v4 g2 L+ Z/ ]) l7 |  
! H  L8 g7 ~0 G8 F1 O( a  sum ddy240
  
, d: r# t6 o+ q6 t6 L6 J% K  if r(min)<0 {
: a" Z! V, Q7 x" Y' C; Z8 b. T        
9 o+ L+ ~0 T% W. h* N- E- f                keep if ddy240<=0
               
                count
+ A" X/ f. v4 Q  S* ?: Q9 ^& T! }               
                if r(N)>2 {
9 J5 g" E4 _2 k& T9 }2 @  c                 
8 \5 ^, [9 Y% F                 pkexamine dx240 dy240
' t4 K! p  ]) ^6 R               
                local tempminp=r(tomc)
                restore
' q. D+ b! |- q6 r7 R3 R                 }
                 
        else {
% `6 Q, A" l% H( A6 |2 ^               
                 di "cannot find minp"
                 
                restore
3 K+ B3 L1 Q7 d3 ~+ E" N                 sum dx240
                 local tempminp=r(mean)
                 
* y* m6 d% y- A% m                 }
6 A8 r5 V, L6 @5 o  h6 O                }
0 \, e' ]$ P  H8 p4 K   else {
9 a1 v; F  q; m% Q. r5 O                 
# t" J, M, Y0 |# q& `                 keep if ddy240<=r(min)
                 
8 r4 Q. e$ m% k& o0 m8 S; i: R                 sum dx240
6 @3 Y9 H1 M$ G$ Q/ n4 s& V8 {& f                 
+ Q6 e5 b: d  U0 u                 local tempminp=r(min)
* r7 P; r" E. i, Y9 S                 restore
                 }
, Z. n% S" R; z+ }1 X2 _  p) H
8 c+ X* h1 _; G

   gen minp=`tempminp' in 1
   *240 minp
2 B" ~# H) h3 }9 ?" b4 a- L7 C& ^

5 ^6 N$ ^% b; ^1 A: v1 d   
   preserve
   
   sum v5 [aw =avol]
9 u% j  Z; a1 S( ]: M   
   keep if dx240>r(mean)
   
   sum ddy240
   
   if r(max)>0 {
- X2 l- S/ U% U* a3 |; Z               
                keep if ddy240>=0
               
                count
               
                if r(N)>2 {
                 
( @$ ?8 T% ]3 ~7 G- q1 W( S$ [4 w3 H                 pkexamine dx240 dy240
1 F. D6 b" ^' g4 U+ _               
                local tempmaxp=r(tomc)
; I4 X. v3 ~& g0 [* ]4 ~; J9 d  l                restore
- @5 ?) S% P3 U* j                 }
                 
                 
; t  M: }4 _$ H6 f% P9 q        else {
                restore
                 di "cannot find maxp"
               
0 z6 v1 i' K9 L' Q6 u6 _& c               
                 sum dx240
                 local tempmaxp=r(mean)
                 
                 }
, P/ R8 t! I: ^; Z                 
# U6 E1 ?/ m+ i  L                 }
                 
    else {
       
1 D7 ~& x. N1 S  e1 W           keep if ddy240>=r(max)
, Y# R; e- x/ L6 V* O          
3 C) F% B" n& ~' |) V; ~           sum dx240
1 c- h9 \7 y% o; D          
           local tempmaxp=r(max)
           restore
           }
       
( y% F8 s6 C& V
/ e, G0 x4 ?6 y+ J  
6 _0 o5 d3 t0 v
   gen maxp=`tempmaxp' in 1
9 l0 S' S5 U( }8 }& P   
- L9 C5 q6 i6 l8 r0 F/ V   *240maxp
   
2 v. }: {- V8 Q: j4 o   
   *-------------ddy240pminp左侧高峰值-----
0 y0 J* ~1 y. z9 }2 p   
- Y% Y: V* p: b* I. ?" F    preserve
   
   keep if dx240<`tempminp'
   
0 l9 Z& [% r0 a. k7 W9 r   keep if ddy240>0
   
   count
               
                if r(N)>2 {
/ _( V( U" q$ s( y- i6 b! t                 
% {8 F6 F( R. R5 v                pkexamine dx240 ddy240
& G  _0 H+ k+ k9 ^               
                local tempminphddy240=r(tomc)
                restore
7 N+ @" L7 K) F  k. f( ~) ^% e                 }
2 Z+ }+ c8 z5 x! @9 L' @1 u/ a$ |                 
$ |" c: N! Y4 ?; e: f        else {
4 t+ H  B2 @5 b# i                restore
                 di "cannot find pminphddy240"
                 
               
                 sum dx240
                 local tempminphddy240=r(mean)
               
# e/ g% |0 j; Y/ d                 
                 }
2 D- x1 F& _- i5 O( T- w$ v                 
* M, X9 }. O$ d6 X1 K   
   
% h9 `. B# W! x5 p4 F; F1 l   
   
0 `$ ]7 J- k  N7 f2 Z6 u8 e& N   
        gen pminphddy240=`tempminphddy240' in 1
2 F6 P2 M  H; g( @& N  i       
       
* `# |( U# p6 K  v        *-------------ddy240pminp右侧低峰值 -----------------
9 j8 z. x9 _0 _& G0 K: E; R1 u    preserve

! z9 F6 Z9 d1 V1 J) q& z' l/ h4 E   keep if dx240>`tempminp'
   
  e9 y+ i0 {1 a" _0 J& K, Q   keep if dx240<`tempmaxp'
   
4 [1 j. }$ ?: u7 ?   keep if ddy240<0
   
+ p  q# \, P2 W  }% {9 [. m! ~+ U   gen temddy240=-ddy240
   
    count
               
                if r(N)>2 {
                 
                 pkexamine dx240 temddy240
   
  local tempminplddy240=r(tomc)
7 q: y% S. W5 M1 U                restore
                 }
                 
        else {
  M7 |' P) i/ e& z               
                 di "cannot find pminplddy240"
                restore
* g: d; Z! ?$ D2 J+ \) o                 
6 A% V, A4 k" G, `( t+ R                  sum dx240
# O' |# A% \' x8 k2 D( |# x                 local tempminplddy240=r(mean)
' t: t) _, X# L                 
3 f5 C; Q$ D3 V& t                 
                 }
" p9 n0 S6 @1 I7 S" {, [! c   
   
5 o- u$ Q( r% E( F, K! w& H   
5 ~) {: f0 E: c" T3 ]   
- r: [2 q% K. @   
        gen pminplddy240=`tempminplddy240' in 1
. ~! [+ Z, @1 n. F       
        *-------------ddy240pmaxp左侧高峰值 -------------------
       
5 x) g+ W& X" L5 U- i/ B$ v         preserve
' ?8 n" b$ [$ D5 K$ c4 Q& [) R   
   keep if dx240>`tempminp'
9 a1 Z+ s- L. C9 j0 c4 g7 ]; F   
   keep if dx240<`tempmaxp'
   
   keep if ddy240>0
   count
% w  S# j/ }: b( {               
                if r(N)>2 {
                 
                 pkexamine dx240 ddy240
   
- d" M) m, N. i5 g6 N- i; t$ d5 G   local tempmaxphddy240=r(tomc)
   
   restore
- q$ M6 p' l* g                 }
4 p5 F0 G' B* w# f                 
% q6 f* w3 i2 E: x4 F        else {
; s. ~* I" r5 O" {& k, a               
7 ^) `# i. Y; f6 v% E                 di "cannot find pmaxphddy240"
                 restore
0 D' D) Y$ C: S* U  c3 `               
. @0 t5 O: i' s# N                 sum dx240
# @# Q! @5 S( f1 v" O, ]                 local tempmaxphddy240=r(mean)
                 
                 
% E9 s0 U/ E6 w                 }
- l2 e/ l  f: d% Z, ?! L1 r0 a   
9 T) N' J( B9 B7 \3 j+ ~$ |% M   
- [( A1 c, m" ?6 V# [* _& G   
6 j% b- R+ z3 j7 T* n- V0 D   
   
   
        gen pmaxphddy240=`tempmaxphddy240' in 1
       
% E1 E+ w$ f/ t6 @; {       
% ]0 X* X3 N: D- |  x        *-------------ddy240pmaxp右侧低峰值 -------------------
$ f# p0 V5 K) N# B" o$ V       
# R, {. }) X5 }; c* E         preserve
  b8 [9 }, O6 n% S% Q' X  
* y; H) B: r+ A3 \. z   keep if dx240>`tempmaxp'
   
   gen temddy240=-ddy240
   
   keep if temddy240>0
; r3 O# A. g+ h9 p4 V' N   
   count
9 V( c- c( I5 ]' u               
                if r(N)>2 {
5 u1 n' q, ^* c* b; D3 J0 C                 
                pkexamine dx240 temddy240
( o! W1 m, c; h( f& T9 w! V: Y8 v" a   
   local tempmaxplddy240=r(tomc)
8 W* O$ W* i, }   restore
                 }
                 
        else {
                restore
( Y+ ^1 b& c7 U7 g) \: I8 {$ Q                 di "cannot find pmaxlhddy240"
9 N: I2 m* U1 Z8 K- ?9 o! S/ g                 
% T1 w  G! C3 D" ~" F; c: v5 F/ @# M: y               
                           sum dx240
                 local tempmaxplddy240=r(mean) in 1
               
( t, q, l6 T. H0 \( ]5 O  L0 q                 
                 }
   
( e1 y( x  I; w$ S& B   
   
1 W* j8 S$ r  [$ i   
  P) U  Y0 x2 g8 z& S, v, l   
   
) J) c& D9 R. {! M8 D        gen pmaxplddy240=`tempmaxphddy240' in 1
       
$ e& s: f9 B! k% ~. M5 Y        *-------------------------------------
( l- ]  T& J% y   
   gen price5d=v5 in 1/20
      
' G8 B& B* H  x. J3 ?2 `   gen price5yy=_n in 1/20
   
   gen price5y=price5yy*(dy240max-dy240min)/20
- z1 v7 J! c. |) t+ c1 p* G   
( L0 p& M. h2 B   *price5y &price5d
5 e, O- F& P3 Q/ _' t# L   
( ?! R7 [: o" R- T. h   gen priceny=price5y in 1
   
   gen pricen=v5 in 1

. y, l, Y- e1 q     
   *priceny &pricen
   
1 G7 a: L1 |: Z' S$ A$ f   
   *-------------------------graph------------------
  
3 _/ y( x1 [1 ~+ \  a& M  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)
1 b$ v- K; l. N1 D7 U
5 D; E; n- v* X8 y8 r" L* ---------可选项目----------------------
# o7 a- w  a9 c$ z2 M. t
8 z7 k0 R* B/ f# l7 p4 D+ Z' I, E}
end

旅途愉快

awer 发表于 2012-3-28 18:06
8 G! U3 i7 x( Q/ `- y+ k2 q# Z改了一下，这次运算应该没啥问题了

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

awer

旅途愉快 发表于 2012-3-28 23:23
. I# A) ]# H/ A5 L, G# l( M/ G这段代码是重新复制粘贴到hello.do里面的？

. n! t  P5 K% }  I对的，就是hello里面的那个

David09

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

. hello 002680
Obs. nos. out of range
% `; ~: n. D) y% |, C% b  _& F7 Ir(198);
) o! j- G" l$ n/ S( Z) ^+ b

% N/ q2 v. ^- o1 g! Q  ?我运行后，出现这个，不晓得是什么错误了，请大拿们帮助

qqysc

David09 发表于 2012-7-2 21:35
/ n3 K% Q; d9 ~! k/ n5 S4 Y+ ]4 U: j. hello 002680
Obs. nos. out of range
r(198);

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

qqysc

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

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

David09

' j% |5 \% h$ i5 d: V5 e; j

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

- I! t  ^% ?% `" l8 l7 w1 Q1 A我清除过了乱码的啊....

qqysc

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

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

David09

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

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

David09

搞定了，原来是我选了一个次新股的数据，样本数目太少了，引起程序的错误。
7 C% c- Q8 ~' w- V感谢楼主给的好东东。

碧水芳洲

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

benhuashu

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

海阔

非常感谢；

墨轩书童

感谢楼主！！

海阔

* s8 q. j- u( l) E/ I按照版主的步骤操作到了最后出现
& _# [1 P0 p: r4 h% v5 Q. V9 r
& \  m! W* R  Z9 P; @1 T1 R. hello 600031
! t; R. N6 `% G! wavolDo not found
$ v" \' A4 q0 t* {. c+ J- c---------------------------------------换成其他代买，也是一样；总是出现这个问题，怎么办呢？奇怪

maodeng

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