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

maodeng

wenty 发表于 2011-10-27 04:17
总结一下使用心得
6 g, v! y2 d% ?3 R4 z" n1）删掉set memory那一行
4 J  I2 ^# v  a( n& Q' J2）改路径

记事本打开后，怎么用tab

冰火

认真学习

higer

我也曾尝试过使用stata软件对股票平均成本进行分析，但是作为一个非专业的炒股人员，而且日常还有工作的人而言，这种方法似乎显得太专业，而且操作难度极大，不一定实用。
6 `  w" `3 d5 S) h7 x# L/ x8 n4 S  L
/ ?# r7 f9 T: K, z$ `, w使用这种专业的技术分析方法容易造成这么一种现象：大家都对导出的光滑漂亮的图像感觉良好，然而在股市中依然无法进退自如。

jway315

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

awer

改了一下，这次运算应该没啥问题了
) x9 V0 h5 ~& c' a3 f. P% {
! k' c2 e, O! A4 r4 r) p
capture program drop hello
, Z& X* c; X, U, _: O
program hello

$ v: {! B- h* F, I quietly {
3 |6 v) X" l3 _7 B: v) P3 a
clear

$ p" d! s: p) I$ m; m. t" N4 @cd e:\finance\stock\20110101\source

insheet using `1'.txt
7 W/ v# |. g: X, }6 s. i
* V1 H9 `% U! v   keep if v5>0
) h( ?( r/ D' _  v' l  \; @, w   
% R9 W0 E+ @" [* t, b3 l5 {   gen vt=_n
   
* ^+ q* k5 J2 A8 v2 N   tsset vt

1 V1 U; {- d' k/ {; e   gen zf=D1.v5
$ l; N  p+ }# o0 a, Y
% m9 M+ K$ |& ?9 y. V' s   gsort -v1
. o# V2 t' _. n# p# l" X& l
& a' T  h4 i! @& c   gen avol=v7/v5
! g8 e4 Z  `- f0 d% p1 X   
' T8 A, g8 s% G7 ]   keep if avol>0
   
! b  C* ]* x* a% t# I$ n   drop vt
) ~) B+ }9 z" z1 d) R$ e  
   gen vt=_n
   
   tsset vt
   
    sum v5 [aw =avol]

   gen amean=r(mean) in 1
' D: \0 }* s- l   
0 _9 o7 z  T' w( B' k" k! i   * 总加权平均值 amean
( A+ m- `' ?2 R$ S   
0 J8 u2 H& ~* x$ y' V   
0 u0 V9 z. Y$ ~1 Y   keep in 1/240
# H6 K: X/ a6 x
   sum v5 [aw =avol]

   gen mean240=r(mean) in 1
7 O4 P1 _+ E/ Z. f& S+ x+ `) s, V- x   
   *240加权平均值 mean240
   
: p; I7 I+ I0 c8 a7 ^9 T1 g   set obs 500
' f7 ?7 ^, h7 I" f  Z1 x# m   
   kdensity v5 [aw =avol],gaussian n(500) generate(dx240 dy240) nograph
% A1 q! P+ r2 r   
  t0 J4 V/ N. U* T( M4 \3 T" L   *240 kdensity(dx240 dy240)
& s; S8 b( u) f4 r" a+ S   
: E# D4 S: S% T' r5 w% C7 o   sum dy240
# s& U- z, T# s8 \# w' h- O+ Q! z   
   gen yxis=r(max) in 1
+ J) {( J) O% n  l3 h' @) o   
  O4 v2 f7 E" d; ]$ ~/ z   gen yxis2=yxis/2
; g+ s2 x% E5 D: C   
   *yxis
   
   gen dy240max=r(max)
   
   gen dy240minp=r(min)
5 y4 e2 b3 b( l
3 _) l0 s% s7 W5 Y. Z    gen vtt=_n
1 T  D6 v# |8 B8 i8 w
4 R2 i9 g$ z5 G3 q3 S! ~    tsset vtt

6 j3 N# Y( r. H/ Q   gen ddy240=D1.dy240/D1.dx240
! ~7 i0 \1 g2 T3 c1 v9 x   
   sum ddy240
# \1 r. j& i0 C$ X   
   gen ddy240max=r(max)
  [8 P$ C+ J: f' ]   
# {$ h5 [; k+ J$ ~) }   gen ddy240min=r(min)
. s# h: R/ a3 W, T5 }( C% B. [     
   gen ddy240g=ddy240*(dy240max-dy240min)/(ddy240max-ddy240min)
; j2 w3 P8 a) F   
   sum ddy240g
5 e% P5 [8 b( F) O( W8 H   
   gen yxis3=r(min) in 1
   
   gen yxis4=0 in 1
  
  `- D7 k5 h- D# x. ]  *ddy240g
   
; F! c# x' p7 g   
; x; Z; J$ @8 Z. }! n7 s0 r6 O. o  preserve
5 x, M5 W* D/ f: a* P2 x  
  sum v5 [aw =avol]
  
& Q2 ]7 z! P+ E& v2 u  n  keep if dx240<r(mean)
2 e, t) f! @; R$ U; M# ^: I; Q5 ~  
  sum ddy240
, U7 x! P# }5 j, C; A  
  if r(min)<0 {
6 T9 ]$ s& H0 N' F        
3 E/ J+ E5 _+ s( X% |                keep if ddy240<=0
               
3 W. E" Z) e, d. G! L- \: {                count
               
                if r(N)>2 {
3 \& J1 S/ J+ X6 ~& z- u                 
                 pkexamine dx240 dy240
3 M- Y' z' E5 [8 x6 `               
                local tempminp=r(tomc)
' }; C( Q4 u. Z/ ?8 s                restore
                 }
) x# ?6 R9 S( P, P5 o6 M" j: |2 R+ `                 
        else {
               
* g" o/ d1 X$ E5 ?! |) a                 di "cannot find minp"
                 
                restore
                 sum dx240
                 local tempminp=r(mean)
                 
                 }
                }
   else {
                 
                 keep if ddy240<=r(min)
9 P, B8 F5 E7 }+ g( _) w  R                 
# Y& U# q/ I) q( |+ o/ e                 sum dx240
: U( R$ b) c0 u/ \                 
                 local tempminp=r(min)
. ~/ D8 J1 w: K) B                 restore
                 }
. d! R+ t# O6 B/ o) S+ S

7 k- J' _! E' A
   gen minp=`tempminp' in 1
   *240 minp
3 u: ]5 ]+ R% h- S1 S
  o/ ~/ Z( C  p' ]! ^( I" l% t
   
( {" u* b5 |# j1 k+ X( q   preserve
# G0 S3 h4 Q' F   
7 N* ]- W7 T( [# s  J   sum v5 [aw =avol]
   
: j. j- \* F% ~   keep if dx240>r(mean)
; Y& D0 _* T8 i3 t1 b   
   sum ddy240
+ k" c3 \. U9 e/ T5 `   
   if r(max)>0 {
               
5 t  C/ u, _# `5 [: L! C" G                keep if ddy240>=0
3 u! r* Y& G* h9 o6 }" Z/ {1 C7 b               
: R" _+ R! S  W1 U; V) A9 Z                count
               
" Q; v1 }! b# N9 j5 Z- e                if r(N)>2 {
7 m2 I5 G0 y2 c& y" b                 
8 j% b1 X- p7 H" A! h& I# P0 z                 pkexamine dx240 dy240
, g2 I) }( @6 ]+ Y$ ?& ?               
. T2 K% `( @; ~# {, Q+ a                local tempmaxp=r(tomc)
& ^. W3 {# q7 \. a                restore
3 q' j0 v) O2 q7 {5 Y# ?                 }
; X( F$ |& D. h7 c; E                 
8 G( r. y+ V8 O( N% u* k                 
( p" Q, h9 [, z- P+ \* I6 ~        else {
                restore
1 J/ N9 l5 I% Z( m; s4 M                 di "cannot find maxp"
. p. R1 h7 o8 a4 L6 [9 `; x               
3 i9 c7 T. m) V2 }, E0 [+ G5 [               
8 K! R( @) v- o, u! Q0 C                 sum dx240
                 local tempmaxp=r(mean)
                 
& J5 D, B; v: J5 t5 {                 }
                 
8 g2 h' ~8 B; c, ]9 R  \& q                 }
                 
    else {
! i: T$ r5 L* f0 O2 R       
           keep if ddy240>=r(max)
, w, {- ^. g$ R, ~* I" z& T          
           sum dx240
          
           local tempmaxp=r(max)
           restore
           }
       
( s' ?4 {7 t5 G! c, r
  

; I7 J; o' [% q( }7 y; B   gen maxp=`tempmaxp' in 1
* z) R* n) j8 @) K8 s   
   *240maxp
   
   
   *-------------ddy240pminp左侧高峰值-----
: j/ t! y* C3 a& f* y   
    preserve
   
, ~9 R; h7 P  s# |5 |& `' |6 B   keep if dx240<`tempminp'
   
- k3 c! T* J/ q$ F4 i  C! Q! K   keep if ddy240>0
9 f& y6 D0 s( _2 A) B7 M# e! p   
   count
               
                if r(N)>2 {
2 l" `% g6 |! [% x                 
7 }: M6 }3 A/ g/ b7 ~, e8 H                pkexamine dx240 ddy240
6 G6 R2 H* b* M5 b7 e5 g. S               
6 x* E! X: {/ }3 Y9 J                local tempminphddy240=r(tomc)
2 D, ?6 a4 _/ Z& n: Q. u                restore
                 }
8 P7 [. H- r! o9 V$ d                 
        else {
                restore
                 di "cannot find pminphddy240"
                 
               
                 sum dx240
% l% h3 j* ]! @4 `7 V                 local tempminphddy240=r(mean)
" R( }8 e6 R$ n               
                 
, o& n! r( t( ?8 i" E0 M                 }
8 z/ Y- K, b0 D7 T% r) I6 {                 
   
   
   
   
, y# m' C. I7 `: t5 R% n& N8 I- c   
        gen pminphddy240=`tempminphddy240' in 1
       
1 L+ \- ]; {2 z( P) s; K' o9 I( n       
( q( l6 ?" Y7 `        *-------------ddy240pminp右侧低峰值 -----------------
# S7 J: L5 M5 ~/ \4 Y    preserve
" Q4 i4 w2 m- O9 |) y7 M
   keep if dx240>`tempminp'
8 ^2 R: o" F) c4 Q( W1 S+ T   
   keep if dx240<`tempmaxp'
" P1 Z% F( }. I0 x. {& U% ?! a# z   
" p" ?" j' u1 Z   keep if ddy240<0
, ^7 L) C1 x0 o0 m   
- x- s! W7 u4 |! x   gen temddy240=-ddy240
   
    count
               
                if r(N)>2 {
                 
+ Z; o9 `" x- q& Z                 pkexamine dx240 temddy240
( |- @9 g$ H1 P  |3 |# Z2 U% @   
& D. J/ E# k+ X- |  local tempminplddy240=r(tomc)
                restore
1 e6 H% I0 W3 F# o/ X) B                 }
                 
! [6 a2 n4 |2 D1 X  ~& p6 b, k. H0 z        else {
% p1 a. S" F2 a               
                 di "cannot find pminplddy240"
; e- j0 |6 T3 H, h, u4 f/ R                restore
) T; R, N- h1 t. X- T- j                 
4 _# Z4 y! K" S8 V) L9 p- t                  sum dx240
                 local tempminplddy240=r(mean)
2 y  u+ }3 p( A: ]- ~3 U* u  i                 
                 
                 }
' ]6 s* r8 j) w) A+ J3 ~   
   
   
1 U$ D3 Y, W5 W2 F4 u   
' m/ B! \. U" F- M: b+ {7 |   
        gen pminplddy240=`tempminplddy240' in 1
% m0 l7 L1 I0 W5 o, D8 W       
        *-------------ddy240pmaxp左侧高峰值 -------------------
       
         preserve
   
   keep if dx240>`tempminp'
5 U$ R- u  G- d   
   keep if dx240<`tempmaxp'
- j. ^' M1 N" R9 Y+ l* R   
   keep if ddy240>0
) w0 G& I/ j& p. T3 d0 W   count
               
                if r(N)>2 {
                 
( N- y# Y1 w1 n- P; Y6 h- B                 pkexamine dx240 ddy240
   
6 I2 I& t; A9 S/ d6 _( C) L0 u8 u0 P   local tempmaxphddy240=r(tomc)
   
   restore
% e  n/ y( a8 I* T+ D                 }
                 
6 i) i/ m5 H& @" ]0 ]* L9 F( A        else {
- B* n. N) P6 N; \0 m: M* _# Z+ J+ O; u               
9 `' w" w: `2 A; f( T% |" V                 di "cannot find pmaxphddy240"
                 restore
               
6 T$ ^0 J- k  f0 e                 sum dx240
' z. M% y# [" S9 @2 v7 ~+ a                 local tempmaxphddy240=r(mean)
: M3 R4 _+ p% r3 @1 j                 
+ a& k& ~. L7 z3 ?# n                 
                 }
   
/ N# @: V* L* X1 a   
& e3 R3 G  ?7 s   
- n, Y3 ~+ G6 J9 w% N8 M   
) `# {) |/ w% M7 ^3 a- |# i   
   
        gen pmaxphddy240=`tempmaxphddy240' in 1
       
       
# L* T& j: q! w/ v/ {        *-------------ddy240pmaxp右侧低峰值 -------------------
       
         preserve
  
4 W! j' q4 B6 g! e   keep if dx240>`tempmaxp'
( X4 L1 j$ t, f& _5 h8 K8 H   
2 ~) r( @8 }# H, q& E   gen temddy240=-ddy240
   
   keep if temddy240>0
   
   count
+ F$ N4 r1 ]7 K3 c" Y' F9 r; d               
                if r(N)>2 {
& @  y+ t+ v) g8 H' j2 ?                 
                pkexamine dx240 temddy240
; d8 W8 G' T( J: G$ C( [, _   
   local tempmaxplddy240=r(tomc)
   restore
                 }
$ [5 y# o6 W$ p5 J$ i                 
4 O: a, a3 f. c7 z* u        else {
                restore
                 di "cannot find pmaxlhddy240"
! k* Q9 @9 ]  }; w                 
               
                           sum dx240
4 ?6 \; T5 m3 V* f7 Q$ E0 a6 c: S                 local tempmaxplddy240=r(mean) in 1
               
                 
, t. T5 y5 l% d                 }
( X9 j. F. g% C- I   
) o9 _2 T4 e3 E: o   
8 H8 a: [. l9 t& C3 J0 d   
" Y4 K/ g( X: T; i& S   
   
4 ~# }& D0 P  a% m: ^$ |* a   
* }7 R: y: e" I        gen pmaxplddy240=`tempmaxphddy240' in 1
       
4 u; ~/ N* V& z$ l" J5 c        *-------------------------------------
   
   gen price5d=v5 in 1/20
. k+ |1 P% p1 V2 J: H      
   gen price5yy=_n in 1/20
   
* i# D6 z9 P7 V. Y, s   gen price5y=price5yy*(dy240max-dy240min)/20
   
   *price5y &price5d
! `% w- ]: _" g7 J% ]! ]& R' x2 g4 |   
) q6 D' X( }0 i   gen priceny=price5y in 1
) o; [% Y# |# I# C   
4 h% t. {2 m9 W, n   gen pricen=v5 in 1

! ~2 o, {8 L+ p5 w/ J$ T! A6 H" X     
   *priceny &pricen
! h. U( j+ n' S* B$ Z& @9 d9 o2 v3 N   
5 t& G6 F; [6 \/ ^" H/ c+ \   
   *-------------------------graph------------------
1 n0 c, H2 _* G$ d0 f0 ?* w+ I2 ]  
  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)

* ---------可选项目----------------------
# O. X4 a/ l! h
}
end

旅途愉快

awer 发表于 2012-3-28 18:06
! \8 V! c+ g3 i改了一下，这次运算应该没啥问题了

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

awer

旅途愉快 发表于 2012-3-28 23:23
$ N% P# Z7 |5 p0 \这段代码是重新复制粘贴到hello.do里面的？

对的，就是hello里面的那个

David09

awer 发表于 2012-3-29 11:40
$ B9 _0 h( r4 [( x1 H对的，就是hello里面的那个

: b9 ^  m, `& O! U: H4 l. hello 002680
Obs. nos. out of range
( T3 u2 T5 O/ G- W* T0 or(198);

, T/ X/ ?; g6 `1 c$ @# s4 h8 U6 F
6 M/ v* t( X4 U0 u) Q5 W9 ?( ?- h我运行后，出现这个，不晓得是什么错误了，请大拿们帮助

qqysc

David09 发表于 2012-7-2 21:35
. hello 002680
6 u% K$ b3 R1 B2 w! A% Y; K: ?Obs. nos. out of range
r(198);

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

qqysc

awer 发表于 2012-3-29 11:40
5 G( p1 C5 O% u" z对的，就是hello里面的那个

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

David09

qqysc 发表于 2012-7-3 12:14
3 ~4 J% \2 X1 ?: E; r9 u" I复制粘贴有很多乱码，请楼主把上面的代码做成附件以供下载，谢谢！

6 E  P( N( }* h我清除过了乱码的啊....

qqysc

David09 发表于 2012-7-3 17:13
/ O( N! W3 W3 F; O& M7 K# c9 u/ m我清除过了乱码的啊....

  H9 Z7 t- T# I8 ]5 y' M7 @5 ~说得不是一回事。你说的问题可以试试在 002680的文本中，把tab键补齐，导出时不要带表头

David09

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

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

David09

搞定了，原来是我选了一个次新股的数据，样本数目太少了，引起程序的错误。
6 o/ S& ]7 T! A; q6 j4 S2 s感谢楼主给的好东东。

碧水芳洲

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

benhuashu

请问，这一个语句：
gen zf=D1.v5
3 E( n; A+ l: C! d是什意思？主要是=之后的D1.v5，这个"." 是什么意思？谢谢。

海阔

非常感谢；

墨轩书童

/ A4 U* b: A3 {, ]8 S4 y1 T感谢楼主！！

海阔

6 Z+ {5 @4 H( A( [8 j4 V
按照版主的步骤操作到了最后出现

0 S$ h' ~! [. G6 U& h/ o( O& S. hello 600031
avolDo not found
---------------------------------------换成其他代买，也是一样；总是出现这个问题，怎么办呢？奇怪

maodeng

大师，那两条曲线代表什么啊
' `/ T2 i# \- p: w$ T