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

maodeng

wenty 发表于 2011-10-27 04:17
总结一下使用心得
1）删掉set memory那一行
2）改路径

记事本打开后，怎么用tab

冰火

认真学习

higer

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

6 G4 x, `$ K* {# i: D% b% d* z使用这种专业的技术分析方法容易造成这么一种现象：大家都对导出的光滑漂亮的图像感觉良好，然而在股市中依然无法进退自如。

jway315

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

awer

改了一下，这次运算应该没啥问题了
: o) Y- I) P1 O
' L- T3 U+ t) A2 |, R% O& u
capture program drop hello
- o8 u& F% A- m' Y8 s0 s
0 [$ A5 Z% F' m+ c* H2 A# |* F1 @program hello

quietly {

# ?: J  e- J7 m  d% mclear
7 h) S  N( \" D* R& T
1 l$ p0 n. X& S  l- F2 acd e:\finance\stock\20110101\source

( T  n. |* b' Qinsheet using `1'.txt

7 |. B+ o7 I  u! l1 y7 e: _7 S   keep if v5>0
' P$ V# Q. y. @" j   
   gen vt=_n
   
   tsset vt

   gen zf=D1.v5
" p* s. D, q5 c3 ]
   gsort -v1
( J  y+ \3 G, S
1 l" K3 C$ I! F0 h% p  Y. a" _   gen avol=v7/v5
6 c) r) S+ m, \( q3 b9 G2 X5 _: j   
   keep if avol>0
   
   drop vt
: [) F2 j  h7 i7 i7 v  
   gen vt=_n
  z7 B: \1 O  S! P: E) `   
   tsset vt
$ \' W9 K2 y$ E3 @, W   
4 n4 s8 Q# {1 M8 _    sum v5 [aw =avol]
0 e  Z- B8 H5 ^9 x" ~  c
6 @7 M) e: P  A% L   gen amean=r(mean) in 1
   
: l- Y1 c7 ]7 H4 f  u   * 总加权平均值 amean
   
1 j4 z: \+ K$ L( I+ Q   
$ p( M& F6 w: l   keep in 1/240
+ T+ I- j# m6 a% c$ z
   sum v5 [aw =avol]

; Z( {- Q6 {- N, [% C  c7 ]   gen mean240=r(mean) in 1
   
   *240加权平均值 mean240
8 d% R: ?7 p- X, t: J" J  Y$ h   
   set obs 500
   
0 d+ m' `$ s/ p% y4 t! o) {+ H   kdensity v5 [aw =avol],gaussian n(500) generate(dx240 dy240) nograph
; g: w. I3 O! }6 b3 |. @   
   *240 kdensity(dx240 dy240)
' D- a/ A: `5 q8 W6 B   
   sum dy240
2 y! B2 F3 C* w0 c' j) H   
* u3 `* M3 y% g3 H   gen yxis=r(max) in 1
   
& R  b1 f' m! O% e   gen yxis2=yxis/2
   
   *yxis
   
0 ~# a# H0 q$ b! q# F   gen dy240max=r(max)
   
( t( W  t* Q( F3 T) ]   gen dy240minp=r(min)

    gen vtt=_n
/ x( P" p8 G0 v- P2 c% R* u
6 B9 }% U; |, h  m4 N    tsset vtt
  v- ]% f* C6 h! b+ B- q
2 N; k- g! K# W( Q  ]! X9 ^   gen ddy240=D1.dy240/D1.dx240
   
   sum ddy240
6 C- D7 c% S" t5 C   
. F4 Z  K; K. L5 x   gen ddy240max=r(max)
  }9 L# N# ]' N8 F) [* B   
   gen ddy240min=r(min)
+ Q1 Y- i" v, Q7 I, O- ~. f     
* [0 V0 r- }9 w4 S, B, p   gen ddy240g=ddy240*(dy240max-dy240min)/(ddy240max-ddy240min)
   
   sum ddy240g
   
   gen yxis3=r(min) in 1
   
: M9 t# M1 s7 ]5 E   gen yxis4=0 in 1
  
  *ddy240g
" L  x0 K; Y$ c: f# _; m   
   
/ g! k9 \  e9 h& K7 o( o6 ^  preserve
0 H; y4 k% H# f' C/ s6 M0 d  
  sum v5 [aw =avol]
  
  keep if dx240<r(mean)
  
  sum ddy240
1 W7 e8 ^3 L2 `. R1 j0 h  
7 j+ Q0 n8 x5 ?3 ^+ b; M  if r(min)<0 {
  m) U" G7 Y2 c        
( Z6 ~( B3 R% K0 F4 G                keep if ddy240<=0
% s; c) W2 z3 I5 ]  b2 N               
                count
               
                if r(N)>2 {
                 
8 Z0 v% g5 U# C                 pkexamine dx240 dy240
               
                local tempminp=r(tomc)
; j: P0 Z3 W- r: V0 m: R                restore
% D: _" s5 {; [3 j4 L& i2 o& e                 }
                 
- ^- T! a' `8 F% V1 B6 k        else {
  o6 }, W, K7 G3 U$ z0 s2 p2 h               
$ J/ B) ?' b% U. l& V                 di "cannot find minp"
$ L, Z  d/ [+ s* X7 W                 
4 g0 B: M4 n& K                restore
8 D/ ?$ [1 D7 f                 sum dx240
                 local tempminp=r(mean)
+ V' z3 e& r! P$ e' _% S. ^                 
' U6 x" E$ {* F; B: P& P' o7 @: H! V                 }
. K- H* }) n( i/ o% V8 m                }
+ z2 [: Q; s+ R8 z$ l) \; r- ~   else {
                 
                 keep if ddy240<=r(min)
& S$ L/ R0 i8 K3 r7 D                 
. K% h$ V5 L! d+ @2 {                 sum dx240
, f- o- @6 ^$ ^7 Z. [, X; A                 
                 local tempminp=r(min)
                 restore
                 }



: T) \7 W+ m& k- }4 q   gen minp=`tempminp' in 1
   *240 minp
" W' ^. F8 u# ~% z
+ Z3 p" [4 k" v) o1 P. Y5 P5 y4 n
+ X# C0 h* }+ ]; m. A   
   preserve
/ K& [3 u) i+ j2 J. e   
; T: |1 i" w9 Q3 X+ N# d* f  O# d   sum v5 [aw =avol]
  r2 Y: {  B) R8 R   
   keep if dx240>r(mean)
   
. q& y, n' k" |' y: [   sum ddy240
   
8 k8 N1 J  `: d# M0 b* t7 F   if r(max)>0 {
: S$ |: d2 c& |5 L( w               
( g- C3 t2 q! F( K) O                keep if ddy240>=0
               
                count
% s' c: O9 Y9 y  S9 q& X3 Z               
' V: P0 t0 c# v, V; j1 f8 z, K- i                if r(N)>2 {
                 
                 pkexamine dx240 dy240
7 U0 p9 G1 z. \/ q, \* z               
) J( U9 T) @8 L# z7 S                local tempmaxp=r(tomc)
                restore
3 Q7 q* y6 E  w5 Q9 ]$ T5 u                 }
$ m  o1 ^4 M4 U$ {/ x6 Y                 
9 ~5 b# Z: P# c. F- q) F                 
        else {
                restore
                 di "cannot find maxp"
  g; w$ y. v: s- u7 U4 j: i               
               
                 sum dx240
! q7 u- g1 a; V! C$ D5 X8 `                 local tempmaxp=r(mean)
% R) M  {- |6 K7 J+ c9 e  D                 
# [, _7 n* L- {                 }
& ~& F$ R/ P- \* o9 P6 l# u# W, N                 
                 }
- ?6 Z, z- ^. C& `( w. `. O4 A                 
" _$ Y7 @3 `2 W- ^! R+ }* |! p    else {
" M9 q$ M' `5 s! q! b       
$ ~- @( Z$ G/ H8 W- ?- B5 T           keep if ddy240>=r(max)
- t" F9 U- D( H6 g9 E0 s4 K/ d6 m+ a, o          
           sum dx240
          
           local tempmaxp=r(max)
. w3 ]7 `4 x" f* T2 S( i# X! b           restore
           }
       

) `1 R. H0 T2 g5 P- h  

   gen maxp=`tempmaxp' in 1
$ ~/ P1 ^4 T/ t   
   *240maxp
; J5 Y5 H( h. |% D$ k- k   
   
- p( C+ ]: O4 c  a6 U! P   *-------------ddy240pminp左侧高峰值-----
   
    preserve
   
8 T7 h" @- u2 [  R   keep if dx240<`tempminp'
$ K# f+ e4 }% U- m/ n   
   keep if ddy240>0
: v. }: E8 b* m& h" x9 d! `8 H: ~   
   count
               
                if r(N)>2 {
                 
                pkexamine dx240 ddy240
: }# k, c" V8 C3 G" h               
7 F* V. ]+ x) u3 U7 ]' q1 k0 v                local tempminphddy240=r(tomc)
                restore
" p/ J- z% G/ o" ]& k: n                 }
                 
" r! H  a+ U* S& x' N6 P6 _! y        else {
                restore
                 di "cannot find pminphddy240"
                 
               
                 sum dx240
6 q9 R* r2 v8 }7 f" I& E                 local tempminphddy240=r(mean)
8 A  z1 @  x2 z- \               
6 W/ a+ q! Z. P- n; Q8 i; J                 
                 }
8 Z6 d. ~, L: ^  a3 C- B3 d                 
   
   
   
1 P# y2 C  E/ k' ~  x( k   
   
        gen pminphddy240=`tempminphddy240' in 1
       
       
; ^* v8 R0 \* f$ Z1 P6 u$ H        *-------------ddy240pminp右侧低峰值 -----------------
    preserve

   keep if dx240>`tempminp'
   
, R+ |% P; e* C9 r7 d' J+ _4 P   keep if dx240<`tempmaxp'
   
, E+ C3 r4 h$ L8 f   keep if ddy240<0
   
4 d( ~) R- ~% h0 b9 |) t   gen temddy240=-ddy240
+ G& Q3 W% n& Y# E# U8 k4 J0 I   
) {% |3 g$ E, j7 s; W    count
               
( j/ l; v8 {" D                if r(N)>2 {
6 G3 ^  |' r) d' [$ Q/ @4 _1 O                 
                 pkexamine dx240 temddy240
   
  local tempminplddy240=r(tomc)
                restore
8 r. \1 E, R$ x5 K* n8 S* s/ U: ^                 }
                 
        else {
               
                 di "cannot find pminplddy240"
- H' S5 d, @9 i  B" B                restore
                 
                  sum dx240
7 Z& Z  J9 Y* x                 local tempminplddy240=r(mean)
0 v& W& A2 X' t7 f" k                 
                 
- i8 c9 Z2 |7 j                 }
% I" ^0 c0 @$ d  n' v: ]2 y' S   
   
$ T/ K$ D/ E1 J, D4 e% H   
   
   
        gen pminplddy240=`tempminplddy240' in 1
; f) n/ m% }6 E; h& W: i; R       
        *-------------ddy240pmaxp左侧高峰值 -------------------
       
, L5 d( P$ n6 }& }# r$ @         preserve
" f$ R2 B. x8 S% Z( d8 M& O: ]   
- i" q: k9 F3 j8 D% L   keep if dx240>`tempminp'
   
   keep if dx240<`tempmaxp'
   
   keep if ddy240>0
4 z* g( ^+ U, f5 q' K+ X8 O   count
               
                if r(N)>2 {
2 X9 U: b$ g; n                 
/ v+ v5 z, @, m( Y! ]9 v                 pkexamine dx240 ddy240
   
* }) G5 m, |; d8 ?9 m   local tempmaxphddy240=r(tomc)
   
   restore
. O; f9 L" N7 m- S0 Q! f) y                 }
                 
        else {
. `8 k% o+ M5 A               
* h+ G' u- B+ o, G& d% u2 J) h. n! [/ Q                 di "cannot find pmaxphddy240"
) a# q2 m. u! S" d                 restore
/ [& R# h. M  B) K2 ^: d5 V               
2 n0 X& C, A1 X7 v: y1 _                 sum dx240
8 N3 ?! V$ F, L4 `% Z                 local tempmaxphddy240=r(mean)
8 K5 D$ @: v/ q. K$ Z1 l, @4 F" t# i& O                 
& n7 {% G9 `# N& |                 
                 }
   
   
. y2 b9 C2 X+ B4 D   
: L) J9 j4 [8 |   
   
# y$ u5 j! ~7 K9 {& V: [( Z: ~0 x   
, B0 f7 J$ \" \! m7 Z1 t% D7 `. @        gen pmaxphddy240=`tempmaxphddy240' in 1
/ s2 H( b3 i& |7 m. n; I- g       
       
        *-------------ddy240pmaxp右侧低峰值 -------------------
4 Y, B; C2 e) X/ _" m6 _; S( }       
         preserve
% _7 N, u$ f5 @8 Q7 t1 o3 H) f- i  
   keep if dx240>`tempmaxp'
   
   gen temddy240=-ddy240
   
; d" h# F) I9 d3 ]  u. ~" _   keep if temddy240>0
1 H$ Q/ Y2 s. n+ R   
, {4 K4 e& w/ E5 h   count
               
3 F: e6 P- ]1 G$ c1 }+ P3 L( `                if r(N)>2 {
, _4 M* \! _& Q* S7 }" Y) @$ H                 
. p, w& h8 }/ f8 \' d  G4 ~                pkexamine dx240 temddy240
* C" q! B- W+ R: f7 v* d' N   
   local tempmaxplddy240=r(tomc)
   restore
                 }
                 
; [4 }" f" v* l/ l        else {
                restore
$ Z$ T; z- Y: D+ I                 di "cannot find pmaxlhddy240"
% Y1 m: m/ @' i0 }5 _                 
               
* i" H7 K9 Z6 C1 h                           sum dx240
                 local tempmaxplddy240=r(mean) in 1
/ Z4 ^: g$ r  \6 D               
                 
                 }
& R: W1 G& |; A& }% m   
9 N3 b" t" A* y9 i  Q   
   
   
   
- `  ]: @6 m# @  \   
- K) N  F0 B; [( r; P. p) I        gen pmaxplddy240=`tempmaxphddy240' in 1
       
        *-------------------------------------
9 c; Y5 v, Z. j% X5 F1 Q3 _   
   gen price5d=v5 in 1/20
1 K/ D- D" B" q. z      
+ q2 v# U1 v  A1 d$ o% @' v   gen price5yy=_n in 1/20
   
  ~+ }4 M7 Z+ X! R, A: ~   gen price5y=price5yy*(dy240max-dy240min)/20
% U. [- b8 ~6 @7 q9 G   
   *price5y &price5d
7 G1 u5 I- L# a7 o; X5 x   
, u0 Q: l+ J9 f8 B7 l$ F  p   gen priceny=price5y in 1
   
% V) P8 g! o1 u# a5 \3 x" [   gen pricen=v5 in 1
7 z8 m! k2 K0 m# I# }6 W9 g
     
   *priceny &pricen
   
" t/ D; _: Z; ~/ u2 @   
   *-------------------------graph------------------
  [! P& W: R$ T8 x  
! \( `1 k; u# e( z7 Z- P( @+ k) E  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)

* ---------可选项目----------------------
) p# F- ^, C4 Z+ e0 ^- a! j3 j
0 \3 z. }( y- v, H6 \}
end

旅途愉快

awer 发表于 2012-3-28 18:06
3 u7 O  [; M5 A- @5 ~1 I- |$ S4 @改了一下，这次运算应该没啥问题了

* p/ |0 M9 ~2 j4 Z" N" d- T4 v这段代码是重新复制粘贴到hello.do里面的？

awer

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

" A  R! F( x3 i3 W9 k# {7 `对的，就是hello里面的那个

David09

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

. hello 002680
- R5 {) F& Q* @& s8 M* HObs. nos. out of range
r(198);
  p0 a% g: F1 O$ M3 I& S

- l8 S3 p" T4 Y我运行后，出现这个，不晓得是什么错误了，请大拿们帮助

qqysc

David09 发表于 2012-7-2 21:35
. hello 002680
1 r1 D& D' \# w' Y+ O3 M+ zObs. nos. out of range
. H7 e9 Y( B+ O- g) G5 e" ^r(198);

, w. g7 A5 f( f: U! ~% }) b; {估计是导出文本格式问题。

qqysc

awer 发表于 2012-3-29 11:40
  H8 ]9 i# n. R( C对的，就是hello里面的那个

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

David09

qqysc 发表于 2012-7-3 12:14
$ |& Q# A% W( i7 G* b, U$ W4 t复制粘贴有很多乱码，请楼主把上面的代码做成附件以供下载，谢谢！

. M; l0 `. H! @4 {
我清除过了乱码的啊....

qqysc

David09 发表于 2012-7-3 17:13
% f0 O1 `- w$ {$ v  s( g3 ]我清除过了乱码的啊....

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

David09

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

  ^# L7 s* H2 k* W: L& p很奇怪啊，我看stata中间的数据表，都整整齐齐的导入了成交价成交量。但是一运行就蹦出来那个错误

David09

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

碧水芳洲

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

benhuashu

请问，这一个语句：
0 g* B( W- g1 t# X; m( cgen zf=D1.v5
是什意思？主要是=之后的D1.v5，这个"." 是什么意思？谢谢。

海阔

非常感谢；

墨轩书童

& a% I; m2 O5 r6 O: J9 Z$ E$ n
感谢楼主！！

海阔

8 _3 d8 x1 J7 L3 ]
按照版主的步骤操作到了最后出现
, W( ^4 f& R0 Z3 c5 e
. hello 600031
avolDo not found
---------------------------------------换成其他代买，也是一样；总是出现这个问题，怎么办呢？奇怪

maodeng

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