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

maodeng

wenty 发表于 2011-10-27 04:17
3 `. p. v6 |* k9 V0 u总结一下使用心得
1）删掉set memory那一行
" Z* z% z' f* D$ k2）改路径

记事本打开后，怎么用tab

冰火

认真学习

higer

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

使用这种专业的技术分析方法容易造成这么一种现象：大家都对导出的光滑漂亮的图像感觉良好，然而在股市中依然无法进退自如。

jway315

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

awer

改了一下，这次运算应该没啥问题了
3 C5 c0 w/ a; V$ H
# S, u5 Z6 i& n4 C  u
capture program drop hello

  t' b; O  V. Zprogram hello

quietly {
9 _- B0 G/ _7 z- R* A
- R0 D, x- l3 M0 Vclear
! ]/ f- H- k8 C- G; v" Y) K
: N% ^' }7 ~3 N& {* m5 \9 [" Pcd e:\finance\stock\20110101\source

2 R# i8 \& C6 p3 H/ Ninsheet using `1'.txt
1 a. D" i5 l/ V8 ?
/ F* R3 H8 o: s# }) k   keep if v5>0
4 ]  D; c# p6 d; O) ]- n   
. R) d( _9 j9 F. X' }8 F# T   gen vt=_n
2 `- t5 a/ o+ r2 e' E4 L   
   tsset vt
1 c, }9 {" x" o; b; ?  f
   gen zf=D1.v5

   gsort -v1

9 ~6 f9 K8 u- b9 v! m   gen avol=v7/v5
   
' A/ R. F5 M+ H( Y8 m  w# r6 C   keep if avol>0
   
3 p2 C  R, O9 M3 H  u   drop vt
2 h  Z1 }4 I6 Y* J  
   gen vt=_n
+ i* ^2 [% |5 ?/ v# @  Q, Y5 R& Y" L   
   tsset vt
   
    sum v5 [aw =avol]

; D2 K$ X1 G$ d   gen amean=r(mean) in 1
   
6 ?- U3 ~& P% ?1 W   * 总加权平均值 amean
/ g8 z( k0 s- j3 O0 }* q   
   
   keep in 1/240

% `9 L' W$ M$ {4 ~; I: s# p   sum v5 [aw =avol]

   gen mean240=r(mean) in 1
4 ^6 X) Y  w! x2 b8 c5 g; y1 D   
   *240加权平均值 mean240
+ p' ^6 I5 i" R0 ]   
8 d0 Y- `. C* w   set obs 500
   
   kdensity v5 [aw =avol],gaussian n(500) generate(dx240 dy240) nograph
- f( z& S; ~3 D( I" y* r$ F   
) S; x- N( R, R5 K   *240 kdensity(dx240 dy240)
  S9 [7 o- G; t   
   sum dy240
$ s' {3 w8 b7 }8 d' d8 s   
4 ?6 o& K+ y% p* e. P9 K% z   gen yxis=r(max) in 1
   
1 ^8 O$ g5 l" c# Z   gen yxis2=yxis/2
   
" g9 g; Q9 X3 L: o7 ]) w& a) o   *yxis
3 \& }$ I: D4 k  d   
   gen dy240max=r(max)
( Q, N4 y  C# ~+ l! I0 V   
   gen dy240minp=r(min)
: z- m1 v0 b7 X# U! R6 S
* k  @  x4 V) E& m1 h3 G* T    gen vtt=_n
$ N# f- u3 A! }1 Y( D
    tsset vtt

   gen ddy240=D1.dy240/D1.dx240
  g! S+ e+ W4 g   
/ M  R" |( U9 f' @  Y! t   sum ddy240
   
   gen ddy240max=r(max)
0 y( F/ \, j4 u3 q9 I   
6 X( _2 Q" O% h  _$ d) K' v3 V/ @   gen ddy240min=r(min)
     
8 u, I" q6 B' c+ ]4 Y6 T   gen ddy240g=ddy240*(dy240max-dy240min)/(ddy240max-ddy240min)
/ z+ o% O0 I' b9 n8 R   
5 b$ v/ v) j8 d5 K   sum ddy240g
* q6 B1 u# x6 `  g! Z9 e' \/ Y   
   gen yxis3=r(min) in 1
   
+ H- v7 c! V5 b) ^   gen yxis4=0 in 1
  
7 y4 b6 ~7 B$ f8 Y5 K$ e# o  *ddy240g
   
   
( u' U' F/ i9 J. J. }2 F  preserve
+ R0 a6 r% |- {" i  
) C, [& J# Z6 o6 i" O  sum v5 [aw =avol]
  
  keep if dx240<r(mean)
  
0 c9 w2 Y" y! M2 _  sum ddy240
" D) n+ i  |# ?3 Z8 H  
' o6 i- O+ |9 Y. g3 c& K7 m) I2 X9 \5 e  if r(min)<0 {
        
. G3 D- v7 l$ {6 R( C3 j                keep if ddy240<=0
               
                count
8 P  F- j. Z% c, ]               
                if r(N)>2 {
                 
* ?- R! k9 O3 d                 pkexamine dx240 dy240
               
* v$ U" o9 I. q% m                local tempminp=r(tomc)
& G5 G! _: h, H( i, Q2 N                restore
                 }
  @/ U( W& ~  {# p# v- p                 
        else {
' T' V! U: f& O$ z               
7 L3 S5 v% @) H                 di "cannot find minp"
                 
( {- z- _. v* l5 R% a, y                restore
                 sum dx240
+ ^& U" ~; r3 A- {                 local tempminp=r(mean)
9 g7 b/ ~4 p2 {' B  m                 
4 B$ a& U3 {4 S/ K7 h$ e) D! `                 }
4 b0 B: j1 l" |* M! h                }
8 |6 c. f- s6 Y2 ~   else {
7 X! ?& ]2 w' e( J. b" |3 c+ v                 
                 keep if ddy240<=r(min)
                 
4 `- }( H6 n* D' D                 sum dx240
6 j; g5 G, D  I! }0 m9 V                 
' X  t2 v) X/ ^( r$ \                 local tempminp=r(min)
                 restore
3 H4 X( w) u: }4 b' c, z/ P                 }
# V5 |( B  w9 q" x1 z

' V1 Q7 l9 K. Y" a+ }
5 Y/ S6 N5 j, M: A" |8 `, n   gen minp=`tempminp' in 1
9 d+ z1 e6 U) t4 t% O, p5 t, j& G4 _   *240 minp
$ H* d, _. k0 P1 {7 S$ i: {$ u! ~* ]
: w$ s: B; M3 x0 W2 V: U- {' v5 B
   
   preserve
$ C% ?$ d5 |  W) _/ ~+ D   
8 ~+ h6 d  P2 W+ s0 j   sum v5 [aw =avol]
   
   keep if dx240>r(mean)
   
% G# E' t9 [" `! c- ]   sum ddy240
   
   if r(max)>0 {
- K9 }" X/ M. w; T1 `               
                keep if ddy240>=0
4 [+ i) G$ T, ]# E               
- y4 Q! p4 D" N6 k; {  W                count
6 j2 Z0 x5 T0 y) b0 s& a+ N               
                if r(N)>2 {
% v$ {( \6 O& H; b" O                 
                 pkexamine dx240 dy240
: m1 X# o1 V0 L8 b) J$ c1 f" X$ C' _               
# D; L, N1 W, a                local tempmaxp=r(tomc)
                restore
                 }
                 
                 
6 P6 i3 k9 I& a2 @0 J/ E        else {
                restore
                 di "cannot find maxp"
( W. j; ?) N# |+ m9 H0 I               
               
                 sum dx240
4 t8 ^1 r9 c+ k  I/ ^9 `% Z4 H8 ~                 local tempmaxp=r(mean)
; O9 ]" H8 E4 d% Z" Y                 
' ?0 L9 _5 h2 o; O9 \                 }
                 
1 A3 B! G! F4 E! {                 }
                 
7 }/ B3 x& B8 x% y( i/ Z# s    else {
       
           keep if ddy240>=r(max)
2 y: p0 N( n' a& Y          
           sum dx240
          
/ c  W5 H5 ]- x% \           local tempmaxp=r(max)
           restore
           }
       

# o: i1 s) Q5 D8 q9 P  

   gen maxp=`tempmaxp' in 1
   
   *240maxp
% {3 J( X# w, X: t6 p- S  y   
* d* ?3 b  h' w7 d! L* J   
+ F, o% S' y# `& \2 e# r   *-------------ddy240pminp左侧高峰值-----
1 Y) i3 Y1 s- c: I   
    preserve
# g1 w7 z8 V: }4 D2 G& T   
   keep if dx240<`tempminp'
   
- B; _" c$ a, g0 i' Q   keep if ddy240>0
" g$ d, R9 j( S  C   
0 u1 }; \" C3 C3 H, J8 @. o0 C   count
               
                if r(N)>2 {
                 
                pkexamine dx240 ddy240
' r, g3 G5 N2 W: `6 ?& ~               
                local tempminphddy240=r(tomc)
                restore
                 }
, z# H8 _  U! N0 i                 
' J: Q& p! j* q  X* H) w2 f        else {
                restore
1 z4 S! m$ G5 j/ m                 di "cannot find pminphddy240"
                 
+ Q1 @" Z) Y+ {2 ]# h               
                 sum dx240
' o  J! h2 U. L1 [( G                 local tempminphddy240=r(mean)
               
                 
% |/ v8 z2 M  q4 V2 [                 }
                 
& K5 r! Q/ _! S# O7 e: t   
   
   
   
% P- K- g1 G# k( ^+ l   
" n+ a( K' ]6 Z' Y        gen pminphddy240=`tempminphddy240' in 1
8 B+ ^; ~, s$ I# l9 r3 ^5 d* s       
) {# `# c3 O! [& @/ y4 X       
        *-------------ddy240pminp右侧低峰值 -----------------
    preserve

& b, B  Q, x5 g/ f3 Y8 A, z   keep if dx240>`tempminp'
   
& [5 m% q; D3 ~2 G: u( W   keep if dx240<`tempmaxp'
( J& y! \) E6 u" G   
* b' I: R& |$ w( u) _   keep if ddy240<0
   
  R, L3 {- [! d* o( ]( {& i# V4 f9 D   gen temddy240=-ddy240
   
. T/ R# F: e/ T; T0 O4 h/ r4 _, O) f! q    count
/ H' Q2 E" `- B$ p( u               
: X# c" I9 ~7 l* C/ B) ]& l                if r(N)>2 {
                 
                 pkexamine dx240 temddy240
% d& i) q" t& n7 |* ]9 k   
  local tempminplddy240=r(tomc)
0 d* w' p' f4 K% Y! m' G7 D3 y! A                restore
                 }
3 n3 `% U, W3 N8 l! @                 
( i: B, J3 l- n4 n6 ~* ?        else {
               
                 di "cannot find pminplddy240"
# ^& i) i; S  D                restore
) Q% }$ S0 a* e6 ], u! Y$ b( s7 R                 
9 m0 p" s4 d+ p0 K$ E                  sum dx240
- C  F4 t& B* X* U2 r                 local tempminplddy240=r(mean)
                 
                 
                 }
   
   
4 x3 Y: B3 e! S0 y9 T% x   
   
3 b2 |. J1 F0 `9 S+ T   
        gen pminplddy240=`tempminplddy240' in 1
  q) j0 v3 m" J  ~0 a       
        *-------------ddy240pmaxp左侧高峰值 -------------------
7 |! p4 U5 I8 N+ h( x# @       
+ x- D  x" ~5 x0 U         preserve
+ e6 k7 [4 u' O# T6 H/ }   
" `* s) \3 e6 E& D9 ^   keep if dx240>`tempminp'
( b& J5 d8 h; O( `3 H+ j* g   
% Q3 m; ^/ S8 q   keep if dx240<`tempmaxp'
  ~3 h4 f$ `1 ]1 z5 @   
   keep if ddy240>0
   count
% c0 N0 l2 {, V, I5 ?! [               
                if r(N)>2 {
0 e3 L* E( f4 a! @" L3 K; K                 
                 pkexamine dx240 ddy240
   
+ }" c  T# e! O2 n9 S   local tempmaxphddy240=r(tomc)
) z& O( @- H  e1 [0 U7 K8 ]: |1 z0 y   
   restore
' g  y6 a* |( a; P$ W" L6 e" }                 }
, A/ W- L- b0 ^7 @& \' G. a: i                 
        else {
               
                 di "cannot find pmaxphddy240"
) r) Q( v, ^* X. v                 restore
               
                 sum dx240
1 B/ p0 \0 R! S5 F, v                 local tempmaxphddy240=r(mean)
7 a/ d1 {, _3 O                 
                 
9 R$ v( O2 M" K                 }
! Y4 @1 I- C; A  n$ N% ~   
   
   
/ ^! W4 R# G! `   
( N; G+ M: c, K6 ]: n. y% T' T  s   
& P: w: P9 w' Y' P$ O   
        gen pmaxphddy240=`tempmaxphddy240' in 1
0 K3 Y. i& M0 l( a5 a( q, S, g& A       
       
3 I' ]4 O- o% e        *-------------ddy240pmaxp右侧低峰值 -------------------
       
         preserve
4 b/ y$ i4 {  [' I  d  
0 f6 b/ I. s' E/ z   keep if dx240>`tempmaxp'
   
: A( E# F; B  D7 a7 g   gen temddy240=-ddy240
   
   keep if temddy240>0
   
   count
               
                if r(N)>2 {
1 g1 w, Q, n( z) _# W* h1 Z' [                 
5 x4 O$ o" V( Y$ h/ A* f9 W                pkexamine dx240 temddy240
   
   local tempmaxplddy240=r(tomc)
   restore
2 U0 P0 @, O8 u: Q5 n/ t+ f# L                 }
                 
        else {
                restore
                 di "cannot find pmaxlhddy240"
                 
               
( A" ~% [- H" N# q7 x3 W                           sum dx240
                 local tempmaxplddy240=r(mean) in 1
* L& g% j+ l, ^- m; J               
) i0 u, K" N3 H) m% y. X                 
" q; Q& [; C2 {% h                 }
   
   
   
   
' k0 t# t# W% y3 d0 D' h" J. t   
( m# X, [, ~6 N8 K, ]   
/ d8 k, Y5 |" S$ v        gen pmaxplddy240=`tempmaxphddy240' in 1
% K* k! T4 I& q. Q: R* ^( h1 D- }: |       
0 o% _/ v) \" T9 v        *-------------------------------------
   
   gen price5d=v5 in 1/20
      
   gen price5yy=_n in 1/20
$ r% ?/ _% i0 ^: w, G( b6 s. k+ @   
( t9 @& A$ X( {   gen price5y=price5yy*(dy240max-dy240min)/20
4 B5 |' V% E* B0 T7 n& E   
   *price5y &price5d
# e0 ]0 x0 m$ d! y4 m7 G   
: q$ d6 E/ N% N   gen priceny=price5y in 1
' P4 l- Z' }6 X' ?( o   
/ k) {) B8 F7 I7 s0 \' t" D   gen pricen=v5 in 1
- T8 e9 r! J8 U$ k. P$ S
     
   *priceny &pricen
   
4 M, }* g' `0 `   
   *-------------------------graph------------------
  
" }, g- P! g- Z$ y) s; w" @  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)
: ?/ H. B9 ~0 h/ P! b
' F+ H& M+ S) B) M1 N* ---------可选项目----------------------

7 T$ R" ^! N) g6 p3 T}
end

旅途愉快

awer 发表于 2012-3-28 18:06
* q) X# R# }5 p改了一下，这次运算应该没啥问题了

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

awer

旅途愉快 发表于 2012-3-28 23:23
; o* t5 ^- k+ I- i4 f这段代码是重新复制粘贴到hello.do里面的？

对的，就是hello里面的那个

David09

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

. hello 002680
Obs. nos. out of range
r(198);
" }7 Y7 w- |5 M2 i5 G! G0 E

6 I1 y, {% P1 f: V我运行后，出现这个，不晓得是什么错误了，请大拿们帮助

qqysc

David09 发表于 2012-7-2 21:35
. hello 002680
& a9 K; n6 M. FObs. nos. out of range
$ }  t1 W- @3 Er(198);

+ f8 X: s: }! F$ e- k4 Z5 S估计是导出文本格式问题。

qqysc

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

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

David09

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

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

qqysc

David09 发表于 2012-7-3 17:13
5 o7 \" c" O0 l% c9 `2 ?我清除过了乱码的啊....

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

David09

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

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

David09

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

碧水芳洲

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

benhuashu

请问，这一个语句：
- _, ?+ r: q% sgen zf=D1.v5
, }% P. G) F. \$ v7 X( {是什意思？主要是=之后的D1.v5，这个"." 是什么意思？谢谢。

海阔

非常感谢；

墨轩书童

% y7 f4 g: O0 C* D/ |% r6 V' A感谢楼主！！

海阔

7 ]: E9 W! ~, n0 E按照版主的步骤操作到了最后出现

. hello 600031
avolDo not found
---------------------------------------换成其他代买，也是一样；总是出现这个问题，怎么办呢？奇怪

maodeng

大师，那两条曲线代表什么啊
0 U/ z+ N2 p6 n