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

maodeng

wenty 发表于 2011-10-27 04:17
  J2 x) n% E; s% p3 [! Z& C总结一下使用心得
1）删掉set memory那一行
2）改路径

: `5 |$ I- }$ P, P' n记事本打开后，怎么用tab

冰火

认真学习

higer

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

jway315

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

awer

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

/ o, |/ l, T+ H8 w& W
* w' B9 h  E" ?3 Scapture program drop hello
; q" {- Q0 n: H
2 M) H  S& C7 T4 P* R/ E) S  Iprogram hello
4 d9 {  v% Q) e9 k( J( a
9 i3 k. D$ ]; l. S quietly {

. e+ U% {- L& a$ e  k6 M3 Kclear
1 V% E2 v3 L4 J' o
cd e:\finance\stock\20110101\source
( B3 n/ A- K2 L0 n/ y7 R# ]& V
& `2 H) {( s* T' R8 T) h- w/ minsheet using `1'.txt

$ P# F, r7 Q5 h9 R   keep if v5>0
2 o" z3 }8 Z- w0 p# t# c; ^   
% H7 `4 D9 D2 v$ R   gen vt=_n
   
   tsset vt
7 c* y( L+ `% X
) O2 g- Y, r) j$ n, s( ~+ h   gen zf=D1.v5
. U) u8 t  j& ~/ ]. f7 Q
   gsort -v1
+ `; F0 E" \' f# @: ]$ A6 K
3 y  S8 C$ b. z5 e4 A   gen avol=v7/v5
6 T5 L) O5 Z7 s6 [+ ^6 |   
   keep if avol>0
4 G) U. L2 z8 g   
   drop vt
  
9 I4 n# k% V' L1 u4 Q   gen vt=_n
- T7 k* z; R' }! R, _1 p   
" Q3 \9 ^" K7 b' W. d   tsset vt
   
    sum v5 [aw =avol]

5 A# f9 `( s9 n; [% E8 O% A# Y   gen amean=r(mean) in 1
   
: O6 ~0 y% f- v5 }" u   * 总加权平均值 amean
7 K# \, I) I! n) F& m7 F1 s" ~   
   
   keep in 1/240

   sum v5 [aw =avol]
, v2 F0 C2 h6 i; Z6 N( u3 P* p4 c  O
: d2 N7 O' G$ ^: X9 u/ A! t   gen mean240=r(mean) in 1
+ s9 }9 w" X* @' o   
. w4 r0 o4 d; O) ]0 r! x   *240加权平均值 mean240
   
   set obs 500
( K" ?" I4 y! [" j   
" O1 R1 `7 K- w7 n$ X   kdensity v5 [aw =avol],gaussian n(500) generate(dx240 dy240) nograph
5 B& Y/ @* k! h& g+ B: V3 K1 _   
   *240 kdensity(dx240 dy240)
   
9 e1 Y1 f3 m9 ^$ P" l  @1 C* R   sum dy240
8 z  e  f8 w+ w. P0 d! `2 _3 n   
   gen yxis=r(max) in 1
# b) I/ a2 Q- X3 v! u% o6 w  U8 p! S6 [   
1 A2 I+ B8 j% L  d  Z/ W. e   gen yxis2=yxis/2
   
   *yxis
   
8 E/ f. Q2 W1 W( o' D$ H7 o   gen dy240max=r(max)
   
8 R$ P1 L- i: k: [7 ]   gen dy240minp=r(min)
1 l& e- E, Q6 K3 {
' S9 O' i) {* X    gen vtt=_n

- ]4 P9 F5 y( J. m) ]# a    tsset vtt

$ g6 R5 X2 J8 X7 Q  ?2 R/ l   gen ddy240=D1.dy240/D1.dx240
   
   sum ddy240
( \8 W0 F* t4 @/ }% G2 s   
$ k' [) y1 W% W& y: V2 \   gen ddy240max=r(max)
   
6 L1 \# L% |$ r4 \, ~3 a   gen ddy240min=r(min)
$ l' U. V# z7 h3 u+ u8 E     
   gen ddy240g=ddy240*(dy240max-dy240min)/(ddy240max-ddy240min)
9 C' F7 x; u& A% S   
: y3 X7 ^: S9 K( O1 u+ U7 _   sum ddy240g
   
   gen yxis3=r(min) in 1
   
   gen yxis4=0 in 1
  
  *ddy240g
" x2 _6 M$ N% E- E8 P* V   
6 f* S% ]7 F# w0 ~7 p6 \6 ^   
# i& l: P: T: L9 Y; K# l  preserve
  
# G' ^' ~) R+ M4 }& {; S( m: f  sum v5 [aw =avol]
  
  keep if dx240<r(mean)
# f+ H- Z) d  E  
/ R1 N$ v0 R- p& Y+ D6 ~1 k  sum ddy240
  
* E8 b  I. |& Q  if r(min)<0 {
: f8 t% S: d! [  x        
                keep if ddy240<=0
               
                count
" h, A5 `- `! w/ X+ D               
                if r(N)>2 {
                 
                 pkexamine dx240 dy240
1 \- a. u% O1 ^. R- V7 Z               
                local tempminp=r(tomc)
                restore
  P% Q- j& F1 S' `$ z$ ]- g! z                 }
                 
9 Y8 D  q7 m" P0 {+ _- [" b        else {
9 D8 D. _, {2 K5 D: [. ]% p2 g: W               
                 di "cannot find minp"
                 
                restore
  C4 p, U, K: \% L                 sum dx240
                 local tempminp=r(mean)
' L1 D3 T! M+ J0 O/ B# e                 
1 U5 L  i& P& G% t                 }
                }
   else {
                 
                 keep if ddy240<=r(min)
                 
                 sum dx240
                 
  ~4 w& |9 C' X' O% S: n                 local tempminp=r(min)
/ z$ i; a- g4 T0 L                 restore
                 }

# |8 r6 n; o0 J& \: m
* |6 L5 e1 F" c( w( ^
   gen minp=`tempminp' in 1
4 g8 I- w+ @) G0 R- S+ Q5 g1 ]   *240 minp
0 n: z# _! w9 N6 {; s
- w2 y) \  d- ?& Z5 L6 T
   
   preserve
   
   sum v5 [aw =avol]
   
   keep if dx240>r(mean)
   
3 r4 r5 ?6 C0 m' a. N" W   sum ddy240
   
   if r(max)>0 {
! e5 e/ P+ w7 p7 i               
& n& a* i% n# V$ H8 g8 K% i7 o                keep if ddy240>=0
               
                count
               
' u/ C5 c. }$ L* z1 n& g                if r(N)>2 {
' T7 i/ R$ e& n+ D3 |$ g, Q* I. n                 
5 H) ]: C6 v5 K" h9 e0 ]                 pkexamine dx240 dy240
6 o( F  B0 f3 k+ I               
" o- `1 x! B6 r                local tempmaxp=r(tomc)
+ g- C1 ?# A9 O; e4 S4 A) Z! C                restore
                 }
                 
                 
        else {
                restore
8 i" i4 r; f" o) z1 l) X4 u2 H+ h                 di "cannot find maxp"
               
9 F6 ~% q9 }# c1 c5 |$ h: o! |               
+ R/ E* O+ N+ p* z3 C                 sum dx240
- c4 N' f0 u/ X6 I. }% d5 O0 Q                 local tempmaxp=r(mean)
                 
8 g& b# e# ^8 s  s                 }
9 K* q5 k  E  ]* y                 
& u2 g) z/ [9 l: ]                 }
                 
! s+ W7 B& x2 q! E" y    else {
       
, h$ K* g6 R3 }: D           keep if ddy240>=r(max)
- f) f3 A! G/ t2 X          
, y* P. Y7 H; h- b           sum dx240
; K2 C8 `4 l$ T; Q          
           local tempmaxp=r(max)
           restore
, a3 ]- p: E" A+ S0 p$ b           }
       
8 e' j+ g% y. @. \1 f6 l
( g6 b0 I( P+ k- V1 b3 z, E/ y  

3 N- \* ?3 a2 M+ f( Y. A- L: ]   gen maxp=`tempmaxp' in 1
  [. t$ e7 U. U9 U( c0 {   
! w" P5 }. S, Y$ v   *240maxp
   
$ p: P9 Y6 d1 {' U$ C* s6 \- V   
3 r3 S. \' A& k   *-------------ddy240pminp左侧高峰值-----
   
; v4 |( |9 U- z+ x- U0 o1 X    preserve
1 [7 k2 @" {0 P   
( m1 d) P9 u6 H# d   keep if dx240<`tempminp'
6 d5 C6 V/ P- u& F& \$ a' K   
   keep if ddy240>0
1 N/ I  U! S% p* e$ ]5 C3 Y, o   
   count
* p& s3 H' T/ v  v               
                if r(N)>2 {
                 
$ O+ T. ?7 Y/ Q& w                pkexamine dx240 ddy240
" T' b( `+ f9 w4 B2 n: b               
                local tempminphddy240=r(tomc)
4 i# e4 K8 d( ]0 W, q# Y                restore
  U0 Y6 t9 O9 l$ O( r                 }
. d$ b1 D& g5 V9 i/ V                 
5 i8 c# _( x' X        else {
% t( y* D- C$ r7 T4 w                restore
5 Y9 Q1 r: W8 p7 c                 di "cannot find pminphddy240"
                 
               
                 sum dx240
- \4 M4 ?6 {) }7 S" J( d% b                 local tempminphddy240=r(mean)
0 B$ @# N  W( w+ e9 M# e               
3 m1 b" q  C8 E; p5 }5 v                 
                 }
                 
   
   
   
: e7 h! D% |% d( Q; F7 I6 p# a   
- p! j5 K7 F) D  j! Q   
        gen pminphddy240=`tempminphddy240' in 1
       
       
        *-------------ddy240pminp右侧低峰值 -----------------
' n9 Z5 x- Z' @/ A2 {7 l8 W    preserve
- y9 R) |6 z8 x4 E
) E0 a: C" Q5 j2 {. U. P   keep if dx240>`tempminp'
+ M, d, W9 y& P% l% m   
   keep if dx240<`tempmaxp'
% F$ K5 F6 @  r$ P; ]5 v6 O   
   keep if ddy240<0
6 f% S. R1 {" G% x8 J1 B( r7 b/ N  _8 {   
   gen temddy240=-ddy240
' z$ x0 j+ O' H   
    count
               
' s# A9 U2 n" i* _' R) \; `* ]  j                if r(N)>2 {
                 
                 pkexamine dx240 temddy240
/ P5 ^6 P0 G# E1 \   
  local tempminplddy240=r(tomc)
                restore
6 ?$ S" ]& c* x: M& _                 }
) J1 e3 }3 X9 i                 
        else {
               
                 di "cannot find pminplddy240"
* M$ v! `) V* g9 k2 f                restore
                 
* J: q0 ~. ~$ c. L                  sum dx240
5 Y; T) N+ D5 {- M; |; y* L                 local tempminplddy240=r(mean)
                 
/ q% S% k) O; A( V6 L                 
% d+ C) K9 y* k# T' z# P                 }
- J9 f9 h9 M' R8 j% x6 I   
4 k9 y" C$ ?( f" X* S   
   
2 b; a- R6 k7 Q% d+ E4 [   
. P4 W2 p8 ~  y2 V   
: J5 W7 Z- F! }- Y        gen pminplddy240=`tempminplddy240' in 1
1 h$ j) j4 ?* i$ m2 ~1 d       
        *-------------ddy240pmaxp左侧高峰值 -------------------
. s' u+ O& @# T$ }# S       
! O6 [2 K% ~) M; E% w" F         preserve
% U, |9 N& q; w2 e- ]* d   
   keep if dx240>`tempminp'
   
+ K! E# r  |$ ^& c: {   keep if dx240<`tempmaxp'
   
$ o) ~3 K7 e# N# U9 s; t' R0 \6 N   keep if ddy240>0
7 V3 ~3 O% W' {5 N* \   count
               
1 b# ]4 S5 ^9 @- _. f' }/ W                if r(N)>2 {
7 }2 _% N$ G( S0 T1 c* _2 _                 
4 c; j% Y; Q3 C# d' A                 pkexamine dx240 ddy240
1 P- J/ E3 C  K1 m& ]: i6 p   
2 y2 Y8 f: G+ Z) g3 n+ c: x/ H   local tempmaxphddy240=r(tomc)
& W) f2 y# t" ?) i   
9 n- u2 _9 ^; ?   restore
                 }
  N0 H% A: {* U% z7 Z8 S                 
) p% F" N& }: F+ T& w" F2 m        else {
               
                 di "cannot find pmaxphddy240"
. `( i9 a  p. \' W9 B2 B+ f                 restore
               
                 sum dx240
5 ]$ R2 \, `, @6 \6 ]                 local tempmaxphddy240=r(mean)
                 
                 
- {8 @7 ~3 c' j8 Q+ E                 }
   
8 l) W5 U0 O; V  m! Q1 Y   
   
6 r- G3 ?' p# u% m4 D: r% B, V   
3 W. M% `9 D7 k6 r: [  F' A   
: o. Q8 _, i) [3 r  ~: L   
- ^* C, Y* b. p- Z4 T. O1 [        gen pmaxphddy240=`tempmaxphddy240' in 1
       
       
        *-------------ddy240pmaxp右侧低峰值 -------------------
       
5 q5 V* q( R' T+ K         preserve
  
   keep if dx240>`tempmaxp'
   
   gen temddy240=-ddy240
2 x& M+ o' t- I+ {# Z   
2 z& {; V( E% c7 v+ I   keep if temddy240>0
   
( L0 x0 T$ P2 ~! n   count
               
                if r(N)>2 {
, S/ @; V+ `) w! L; P                 
                pkexamine dx240 temddy240
+ o- C' W9 r% o  J   
4 d; c+ }8 }8 H( E   local tempmaxplddy240=r(tomc)
   restore
5 P. y; K0 T1 _: ?0 [( ~# H                 }
                 
        else {
  Q4 j/ A4 {+ j+ l1 m$ v* w                restore
7 L) B( E. ^$ y                 di "cannot find pmaxlhddy240"
# X; n% B4 T" W& C                 
( [$ V3 ?) w2 ]2 o5 Q5 k               
% c4 ?6 w0 A0 [  r1 [3 Q5 I                           sum dx240
! q2 B; c) M+ R4 z8 P) c                 local tempmaxplddy240=r(mean) in 1
               
                 
                 }
   
   
   
   
   
) U- o* w) N$ q   
        gen pmaxplddy240=`tempmaxphddy240' in 1
3 X+ F% l$ U* c& i       
  Q) A9 v6 `9 O        *-------------------------------------
8 z4 j* m4 U# m1 m9 `0 @   
2 B8 Q) i8 N# ?5 K) g7 a& H9 Q8 Q   gen price5d=v5 in 1/20
      
2 C  D) o; i% y! ^, n  p6 S  s   gen price5yy=_n in 1/20
   
# ~' l7 z! U: R4 t5 E7 V+ M   gen price5y=price5yy*(dy240max-dy240min)/20
   
  _# b0 z0 @( J7 h9 \' ~2 v, d   *price5y &price5d
0 ]! S  `6 `2 X# o" m, o+ m   
& H6 {) r+ v8 G/ r; Q3 A   gen priceny=price5y in 1
# ]- Y* }% S+ e2 |: f. ^   
   gen pricen=v5 in 1
) Z3 `  L' ~. O+ ?- }
2 o* V2 d3 d; W: i     
) a( Q3 k3 v: o. Z8 L0 k4 f4 I; K* o   *priceny &pricen
7 d) e% q2 J( X/ f  w' L   
   
   *-------------------------graph------------------
# S9 L  S4 \1 n  
! l0 H2 m. T0 x  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)

* ---------可选项目----------------------
/ {6 j8 I' s# N9 h
. \/ l+ v$ b8 H* o}
2 v3 N7 H! f$ F4 _; Pend
  {5 A# I: Q6 }( x) R' X& r

旅途愉快

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

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

awer

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

对的，就是hello里面的那个

David09

awer 发表于 2012-3-29 11:40
2 s; g! N: X7 W1 V2 q对的，就是hello里面的那个

3 ~5 y: a# u' l/ D4 R* Z$ o. hello 002680
; v+ r! F( N: M8 G; UObs. nos. out of range
r(198);

( z7 ]' k) K8 ~( U( j7 y" d
6 r" [" ?- q5 F) C# Y( P$ i我运行后，出现这个，不晓得是什么错误了，请大拿们帮助

qqysc

David09 发表于 2012-7-2 21:35
. hello 002680
' Y( g8 b' ]) i. l% M- sObs. nos. out of range
r(198);

$ ^. o0 K- k% p0 p7 o7 T# o- Z估计是导出文本格式问题。

qqysc

awer 发表于 2012-3-29 11:40
) ^7 C; a1 S+ V6 v8 k, ~对的，就是hello里面的那个

" I6 l) X* G: \) M复制粘贴有很多乱码，请楼主把上面的代码做成附件以供下载，谢谢！

David09

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

6 f' `4 a: T  t, B% w6 `我清除过了乱码的啊....

qqysc

David09 发表于 2012-7-3 17:13
8 |, M: K5 s  B/ u" N我清除过了乱码的啊....

; G6 j) b, s2 X# d) s' u说得不是一回事。你说的问题可以试试在 002680的文本中，把tab键补齐，导出时不要带表头

David09

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

. H& G7 {# d. m* P很奇怪啊，我看stata中间的数据表，都整整齐齐的导入了成交价成交量。但是一运行就蹦出来那个错误

David09

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

碧水芳洲

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

benhuashu

请问，这一个语句：
gen zf=D1.v5
* e$ `2 c. J3 Q' L5 U" b是什意思？主要是=之后的D1.v5，这个"." 是什么意思？谢谢。

海阔

非常感谢；

墨轩书童

感谢楼主！！

海阔

按照版主的步骤操作到了最后出现
/ G; G3 u! x1 K* Q2 m4 j0 a
& J$ G$ j6 A7 g. hello 600031
9 I" v, r, v& M( d/ k& J: navolDo not found
---------------------------------------换成其他代买，也是一样；总是出现这个问题，怎么办呢？奇怪

maodeng

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