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

maodeng

wenty 发表于 2011-10-27 04:17
$ z+ F) Y& N: d总结一下使用心得
& @6 r! g. t1 \8 s( s  V; M' O1）删掉set memory那一行
2）改路径

记事本打开后，怎么用tab

冰火

认真学习

higer

我也曾尝试过使用stata软件对股票平均成本进行分析，但是作为一个非专业的炒股人员，而且日常还有工作的人而言，这种方法似乎显得太专业，而且操作难度极大，不一定实用。
* d% V! s* H; [6 p; R  ^
使用这种专业的技术分析方法容易造成这么一种现象：大家都对导出的光滑漂亮的图像感觉良好，然而在股市中依然无法进退自如。

jway315

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

awer

改了一下，这次运算应该没啥问题了
- e: H% [* K' p6 T) X5 I) j
5 [' J8 a# g5 v+ }2 e- k- L. x& T
capture program drop hello

program hello

quietly {

clear
5 p) U1 C/ m* z# I  }. \; I
cd e:\finance\stock\20110101\source
% z& `9 @) _7 p
insheet using `1'.txt

   keep if v5>0
+ c  \' R* G4 J6 g4 |   
( z6 ^  I2 m* [" d   gen vt=_n
   
   tsset vt
, t( x( s' v# L
9 B& O5 w8 K5 }2 S# U  h! m   gen zf=D1.v5
& d, o7 b  |  _* C  C7 M& F6 r
% c! n( f' e/ C6 N/ y0 N% h6 g$ s   gsort -v1

   gen avol=v7/v5
0 {, I+ l% S. q+ o9 G* {   
   keep if avol>0
   
9 o/ k2 G6 z# d   drop vt
: L& N# Z3 D. p* i0 @: \# p$ T; }  
   gen vt=_n
   
) v2 q3 Z) J, v( R/ `   tsset vt
   
* T6 W6 g' x/ s( f0 D; |    sum v5 [aw =avol]
/ T# t8 I$ @* N  a  b0 w$ g. |/ }1 x
   gen amean=r(mean) in 1
   
   * 总加权平均值 amean
   
. Q; |  t( I: W8 z/ j- E   
5 H& a, X: F; x* O# f+ {( n$ S, P   keep in 1/240
4 P2 p2 c! m; ?% a, g) g
   sum v5 [aw =avol]

   gen mean240=r(mean) in 1
   
   *240加权平均值 mean240
% U. G) p* x) V; Y! h" C   
   set obs 500
   
/ U! F. Q0 B, Y4 W8 w+ L$ f1 Z   kdensity v5 [aw =avol],gaussian n(500) generate(dx240 dy240) nograph
7 ?  w3 T) H2 x+ S9 J   
* U. g/ P7 ^* P. d! R  Q, I5 `; V; x6 B" p   *240 kdensity(dx240 dy240)
   
   sum dy240
" U0 S- y  g8 }( {! h# L   
   gen yxis=r(max) in 1
. z2 \% y( H/ z% R" O# _$ J   
' F! _1 i8 z' w   gen yxis2=yxis/2
4 ^+ n: v% o* g" r   
   *yxis
& N* y& N8 i( `$ z" [   
   gen dy240max=r(max)
   
" ~$ j1 a! E0 n% [' l   gen dy240minp=r(min)
9 F( K6 _4 D+ ~+ L
, Q7 u3 K) `/ u: U' K+ Y7 [9 s% I    gen vtt=_n
6 |. M- g: b/ Q9 Z2 f
3 q( x- Z! n5 o; q    tsset vtt
* H! [) G- w. N. V7 M/ H( x) ?/ S
   gen ddy240=D1.dy240/D1.dx240
- y* W) J2 b' _! ~2 ?  }* `   
   sum ddy240
7 i9 s" \; X& u" `$ Q' b   
& F4 h$ `7 n9 ~1 B4 d   gen ddy240max=r(max)
7 q1 ~9 y$ U4 R: @) D# I   
   gen ddy240min=r(min)
( |( W0 m3 [) `     
   gen ddy240g=ddy240*(dy240max-dy240min)/(ddy240max-ddy240min)
% [. C% n7 ?8 B3 l3 ~   
! g$ V  U# U" c" z   sum ddy240g
$ m& N. m0 j3 y2 A' W8 W# ^   
   gen yxis3=r(min) in 1
/ s3 `2 p4 `+ D. G) i: k% m* G$ \   
   gen yxis4=0 in 1
  
  *ddy240g
, S$ k6 d6 d5 g   
; r7 p4 N4 e- v   
  preserve
  
- v+ P! _# n: P1 S' Q  sum v5 [aw =avol]
' ^4 l3 h6 c& ]0 @$ B  
  keep if dx240<r(mean)
$ X- b& c, j  W% K( j5 z  
  sum ddy240
1 H  x$ y& ?# C( T/ h7 |( a  
  if r(min)<0 {
/ v& e& G2 C0 Q* p3 M& `2 ?* j        
                keep if ddy240<=0
6 \' z# [4 |9 S1 C; l: \5 x               
9 I4 S( Z- L9 I! q& U" n; I* l                count
               
                if r(N)>2 {
                 
                 pkexamine dx240 dy240
  y% @, B+ x- V" n' n               
4 C) [! d9 |0 y9 i                local tempminp=r(tomc)
( p% O+ n3 u! J( b" i. ]                restore
! N# N" }$ a% j" }: `6 k1 ^; Z                 }
                 
$ p  o5 m; E' P1 }        else {
1 ]) d, z6 m( N, {6 u/ o. h               
                 di "cannot find minp"
                 
- |: c/ R+ ]0 P! L, c& A. L6 s                restore
0 H8 B" z4 h* S' h                 sum dx240
- {. I9 }0 Q. _7 N" p8 t' C4 h                 local tempminp=r(mean)
                 
                 }
                }
   else {
                 
                 keep if ddy240<=r(min)
  V5 s( R/ F  _" w- @                 
                 sum dx240
                 
                 local tempminp=r(min)
                 restore
                 }
3 u8 H0 w' L9 F( p  ]/ o
" A1 e; a4 b8 O9 o; `$ O
0 Q* s9 M' a; C& I2 u9 L7 U
1 U: Q8 \) r: Y* @   gen minp=`tempminp' in 1
9 @5 M4 l. g/ }5 y   *240 minp
8 ^$ n& u* N3 I8 |

   
# G) b4 C1 t  x# Z# i- _6 R   preserve
& r5 |! w' D2 v2 U   
   sum v5 [aw =avol]
" }) A/ H4 w! P/ t. n% C5 N   
   keep if dx240>r(mean)
   
' z4 \" {+ C- A5 O1 y' {   sum ddy240
% V! j9 W# j" B- |. Z   
. w% N4 ^5 Q- |$ e2 f' T8 ?   if r(max)>0 {
# p( ]- C; x; U, `, m4 u0 k               
" a  ~* v! s' r% H                keep if ddy240>=0
               
                count
6 I6 a2 U5 B. Q6 P" ]' R               
                if r(N)>2 {
  O3 }1 u& F' W9 Q) f5 d& R3 b4 g                 
8 p3 c; O) O% S; F% Y                 pkexamine dx240 dy240
               
                local tempmaxp=r(tomc)
                restore
                 }
" Z* c( ~3 ]/ p7 z) l# l                 
                 
* c+ p( L6 y: T  D% c5 Z        else {
                restore
                 di "cannot find maxp"
1 U. G" G* X+ U: [9 F# p) ]               
2 r- G( ?* M4 }0 D; ]! g               
  B7 q' y) j( l3 A+ F6 j                 sum dx240
                 local tempmaxp=r(mean)
4 W8 Q; `: C5 t9 Z" d                 
' D* W3 r% g/ _3 k5 b* y) q+ C# f& Q                 }
                 
                 }
                 
2 {8 J7 k1 s+ p; @1 B- D7 ?% H    else {
       
; f$ C. V& {0 h( O, e# v. v           keep if ddy240>=r(max)
$ U9 B0 \1 D: e* j# z2 Y$ ?) e          
           sum dx240
* z' c$ ?) E9 \. l+ h/ O% g          
$ _4 Q; M$ f+ y4 L6 M  |* U$ u           local tempmaxp=r(max)
% s/ L4 e- ?$ y5 V% M* g* Q           restore
           }
8 w  Y: q4 B4 c& [% }1 A. s0 J       
  e' X6 h8 j/ e: y5 X) H
6 s7 d# R& q: c1 q% c$ b  
  w4 x5 O6 S1 I( {4 l
+ ^( V" I2 N! y   gen maxp=`tempmaxp' in 1
   
' h7 ^7 b( o+ n9 \6 O5 Q   *240maxp
   
   
   *-------------ddy240pminp左侧高峰值-----
   
    preserve
   
9 f. r# q: j% N   keep if dx240<`tempminp'
( o  L# E% R4 W/ A$ e   
( h0 T# g* }3 n2 I   keep if ddy240>0
' E& L4 J6 C5 O- Q( m   
, _: m6 x8 P, g" W5 ^   count
               
                if r(N)>2 {
" r. S, p5 d. C( p                 
                pkexamine dx240 ddy240
- ]. Q7 T$ [) W/ W* e- w* k4 f               
5 P. {. l9 H2 d! ~: q; f                local tempminphddy240=r(tomc)
                restore
                 }
                 
        else {
                restore
' n9 s! c% B0 e$ z( Y                 di "cannot find pminphddy240"
( S) q# v. W8 G$ @* ]                 
: H* n" g$ I& D- V; }               
                 sum dx240
. a5 s+ P  z0 o1 B% {/ R8 k                 local tempminphddy240=r(mean)
0 U% i, g/ I# {               
5 O& r7 _% b( V3 [; \                 
                 }
# I9 {' X- y2 ^8 p2 P" b* ~                 
4 k  w7 z6 R- C: \7 G   
) y' ~- C4 p5 I6 v8 X4 j7 V   
   
$ ?/ o1 ^9 G  t' i9 y, x3 y   
   
7 j1 t5 z, }) b5 }4 _& S; R        gen pminphddy240=`tempminphddy240' in 1
) @; l6 z/ w6 u* J: A8 O       
$ p6 V/ v! I2 N, G9 ~       
. c* N% ^+ M. O, o        *-------------ddy240pminp右侧低峰值 -----------------
    preserve

% A& {* ?) |* k   keep if dx240>`tempminp'
   
3 R! j6 B+ }) H7 {  ]; n4 R   keep if dx240<`tempmaxp'
   
   keep if ddy240<0
   
7 q  a  E* _7 \  G% {/ y9 q   gen temddy240=-ddy240
1 Q: v* M% o  J9 n   
    count
" ~/ w- p; H9 h' Q* @               
                if r(N)>2 {
                 
+ |4 ^' a$ M( e, f                 pkexamine dx240 temddy240
   
) i( @5 u+ R; i+ t$ ]  local tempminplddy240=r(tomc)
                restore
                 }
* f3 D1 r7 |! c: \4 d                 
5 x2 s% S9 ^2 ^1 Z6 k% I8 \" j        else {
$ u# b3 ^1 ~% V* n; [$ f2 O. @               
' q$ }. w- {$ l- J8 V                 di "cannot find pminplddy240"
; F# Q) `& u& }$ I                restore
# a; S* n. ^8 U% i4 t4 D1 N. Q                 
8 T2 ~/ g+ G  G6 |, d1 L                  sum dx240
                 local tempminplddy240=r(mean)
                 
                 
1 J5 A% k) y' F1 ^0 S                 }
1 j( D2 _4 e* \& J   
   
1 C2 @" y2 U4 r0 L: ]9 T/ @   
   
- G- n- w9 A/ b0 r   
  ^  B2 |9 `" y* V2 y        gen pminplddy240=`tempminplddy240' in 1
+ Z& `. s: N4 M6 H5 s% @, O1 O       
* a# k. l# C0 [3 E1 F        *-------------ddy240pmaxp左侧高峰值 -------------------
       
% t6 z  R$ v+ J' O         preserve
   
+ J( e& w$ P* ]) Q' m   keep if dx240>`tempminp'
   
   keep if dx240<`tempmaxp'
   
   keep if ddy240>0
   count
, {! ~$ P- U6 \% F- Q& x               
6 j) ]  j1 \$ L" R* \6 l4 O1 f                if r(N)>2 {
4 b) U8 ?8 l% e3 l+ I3 I, I                 
: O9 h$ b  ~) q0 |9 z: _7 a                 pkexamine dx240 ddy240
   
   local tempmaxphddy240=r(tomc)
' Q  a2 s/ H# _/ ~: ?. s   
   restore
5 p8 }, K7 ?% D/ U* u8 a                 }
8 [9 b* O7 X, ~; e: D1 U                 
! w7 k' i3 j! s- ^( E        else {
               
. A9 R0 c, P2 ^( X1 i3 |* y                 di "cannot find pmaxphddy240"
                 restore
               
! J8 y- N- `4 n5 U7 P' s                 sum dx240
                 local tempmaxphddy240=r(mean)
0 o( s% K7 X1 J0 \. ?                 
                 
+ D: z* p1 a) L0 H& J                 }
/ {+ g8 e1 s" g- @9 @' D3 p7 I" D+ Z   
" h5 M0 x0 d+ s% P. l) S2 f   
   
   
   
   
3 K( Z, C1 x3 ^        gen pmaxphddy240=`tempmaxphddy240' in 1
       
       
1 g7 @! P& d, k5 |5 o8 y4 z5 ~        *-------------ddy240pmaxp右侧低峰值 -------------------
' j6 V) k9 H- Y0 p- ?6 u       
+ d0 F4 r" J  n8 e/ Q5 |         preserve
  
6 c6 N1 t) f- P   keep if dx240>`tempmaxp'
   
   gen temddy240=-ddy240
   
   keep if temddy240>0
) o& g6 ~0 r7 Y! `   
5 r5 P3 h. O* U; A! z   count
               
! k3 Y, ~  R5 u0 ^% Q, X8 N9 J, a                if r(N)>2 {
                 
                pkexamine dx240 temddy240
- T7 W7 t  b; U% \! Z& Z& ]   
   local tempmaxplddy240=r(tomc)
   restore
% ^7 H4 Z$ i3 b! s# B; d0 B! d4 `                 }
  o$ I( F: ?' Z5 k3 X9 q                 
  j9 s# w7 U. n        else {
                restore
                 di "cannot find pmaxlhddy240"
                 
# _3 h. ~) |0 P' t5 M               
                           sum dx240
. z/ D; ~2 L% f$ o3 }, Y8 F                 local tempmaxplddy240=r(mean) in 1
               
                 
+ r+ o* n( @6 h1 h2 ^                 }
: n8 g( j) u8 K2 Z   
- F3 [- P6 E1 V% B( O- ?   
   
   
   
4 n) Q: s9 _# J# Y' f2 t% ]2 Q   
        gen pmaxplddy240=`tempmaxphddy240' in 1
       
        *-------------------------------------
* ^$ H  l5 V( r  \' V2 Z   
   gen price5d=v5 in 1/20
" \2 O! L& }# W2 H: d  b! L( {3 O      
   gen price5yy=_n in 1/20
0 j3 ?( {# o) h" Y  h   
; H( @7 F' |7 V5 @0 E, _; G7 K   gen price5y=price5yy*(dy240max-dy240min)/20
   
3 H% {; u! ]* }+ c4 o$ {6 x   *price5y &price5d
   
   gen priceny=price5y in 1
   
   gen pricen=v5 in 1
7 e5 M+ p" f4 z6 ]& O" U, }
8 r# z) m$ c7 [" A1 Y5 m     
   *priceny &pricen
   
   
' P* ?3 S& F  I1 n# p5 Q. ~; x   *-------------------------graph------------------
0 |. q6 t& Y' j' E$ G5 ^  
  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)

! j& h1 ]9 u; b1 Y* ---------可选项目----------------------

* @+ G3 \: T" k. O5 \& _& c}
end

旅途愉快

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

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

awer

旅途愉快 发表于 2012-3-28 23:23
. {9 w: _( f/ J) R" g2 {  O! Q2 J这段代码是重新复制粘贴到hello.do里面的？

对的，就是hello里面的那个

David09

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

. hello 002680
Obs. nos. out of range
" e% [; R% W# l  c5 Kr(198);
+ L7 Y# l. L1 X1 Z, t
7 k$ n3 y; q$ c1 s/ ^+ e# \7 h
( {: _+ Y; M0 {. w9 o+ `' ]我运行后，出现这个，不晓得是什么错误了，请大拿们帮助

qqysc

David09 发表于 2012-7-2 21:35
. hello 002680
; l1 z7 B& J" v) OObs. nos. out of range
. [; }; q' _- ]# F7 J+ N) ]r(198);

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

qqysc

awer 发表于 2012-3-29 11:40
5 N. l; L, t3 ^对的，就是hello里面的那个

1 K3 H5 c' i3 c. [' _8 P复制粘贴有很多乱码，请楼主把上面的代码做成附件以供下载，谢谢！

David09

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

0 @! L+ G/ p+ a& w2 z+ j3 y
我清除过了乱码的啊....

qqysc

David09 发表于 2012-7-3 17:13
4 r) v# _8 p/ x! y% T# W1 a我清除过了乱码的啊....

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

David09

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

$ m  Q& Q0 @8 A8 B+ S- G2 \很奇怪啊，我看stata中间的数据表，都整整齐齐的导入了成交价成交量。但是一运行就蹦出来那个错误

David09

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

碧水芳洲

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

benhuashu

请问，这一个语句：
gen zf=D1.v5
# |* s( r7 ^% n是什意思？主要是=之后的D1.v5，这个"." 是什么意思？谢谢。

海阔

非常感谢；

墨轩书童

4 E0 U6 J3 M1 G7 l' c. L感谢楼主！！

海阔

按照版主的步骤操作到了最后出现
5 V% t, f- c7 f& L" R4 C) {6 z
: {, ~; A2 ?0 V& N: T. hello 600031
avolDo not found
9 ~2 x+ ]0 Q/ `6 Z1 V/ J. X* f+ ~---------------------------------------换成其他代买，也是一样；总是出现这个问题，怎么办呢？奇怪

maodeng

大师，那两条曲线代表什么啊
" D' P$ ~) ~) t5 ?  _$ i% g