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

maodeng

如果我想看最近 一个月，三个月 六个月 的成本，怎么修改程序啊，谢谢大师了

maodeng

我用同样的方法导出数据，发现有的能生成，有的就不行，不知道问题在哪儿，比如300236就生成不了，我的邮箱是19860225gaoyuan@163.com，麻烦大师了

maodeng

. hello 300236
statement out of context
4 [! x* q" V; R7 B! `+ T+ j. s9 m1 ar(119);
  S3 {2 ^- x4 T' e) d/ _3 Y0 _& X
- g1 m+ l! l6 p9 ~+ J7 S5 L大师 就是这个

sergio

这贴实在！

awer

maodeng 发表于 2013-1-28 13:13
. hello 300236
4 l: R2 y2 K, O/ Ostatement out of context
: ^& T/ p* K& Nr(119);

9 W7 r& D/ P: R; J" B6 L/ r" K8 Q抱歉，最近太忙了
( z# ?( k6 q* B+ |
300236 我用下面的程序试过，没问题

7 ?, L" b+ N; Q/ _你把 D:\stock\test\source 改成你导出的数据存放的目录，再试一试

capture program drop hello

program hello

8 c$ w$ V5 Y  b7 R* m quietly {
: ?5 u4 R. a. [% t) j
- @9 e6 N9 y2 T; qclear

cd D:\stock\test\source
" H7 p. Z( U' _
insheet using `1'.txt

   keep if v5>0
   
   gen vt=_n
   
* E9 F4 ~7 t, T- Q0 T% J   tsset vt
# |' |* c* L9 K6 v" q
   gen zf=D1.v5

) j, a0 X7 @1 `  b0 W$ c) [( }* |5 [  z   gsort -v1
' n. T; Q+ |- S# B" ^
   gen avol=v7/v5
) h! n, V9 y6 K2 A% ~( a  O   
   keep if avol>0
, R; D0 ?4 d# q6 T, ~   
   drop vt
/ L/ C- s' ?. V3 f  
   gen vt=_n
1 A/ B; W7 q7 [   
4 L. I" B! r) _5 W& k2 w' `$ N) J   tsset vt
   
    sum v5 [aw =avol]

+ K+ `" N$ @9 z$ C   gen amean=r(mean) in 1
0 Y" R2 X  D+ n8 l9 Z   
# K- m" G* T4 [; r7 i+ c   * 总加权平均值 amean
5 W  ~3 T; H' u5 B% `   
! N; n, T7 G5 J% u   
   keep in 1/240

6 N) _! U9 s3 @3 ~3 S. a   sum v5 [aw =avol]

   gen mean240=r(mean) in 1
4 ~* [* O/ o0 b' j+ I" ^' X1 V% a   
   *240加权平均值 mean240
   
( A5 d, @, |# G) M* b   set obs 500
   
   kdensity v5 [aw =avol],gaussian n(500) generate(dx240 dy240) nograph
$ b2 z# b1 T- o2 W: A" \   
4 n7 Q: @$ }- D: j& X+ V# v& l' d   *240 kdensity(dx240 dy240)
   
2 l& @/ _, j8 X/ D   sum dy240
   
% c  L% d4 `' e! H* d/ D   gen yxis=r(max) in 1
   
   gen yxis2=yxis/2
   
) J- w/ Y( w' X$ p6 G7 Q2 t( c+ D. K/ ?   *yxis
   
) v: Z: ~" z7 h4 i' ~   gen dy240max=r(max)
   
1 @. |8 f( w! W, m" ^: `) X7 B   gen dy240minp=r(min)
" p; l5 m$ n; B" m7 U2 K! W
5 G9 n+ R' I" A# e+ K9 H    gen vtt=_n

, w0 X( Q# b! k/ h    tsset vtt
  v& N  G8 a: ?
   gen ddy240=D1.dy240/D1.dx240
   
   sum ddy240
, y# n! Q- g* S; S   
   gen ddy240max=r(max)
1 k1 K4 e: T' f/ x# ~0 C1 @   
   gen ddy240min=r(min)
: E1 x/ q' I; [0 j5 g5 D     
8 \; ?3 ?$ E3 V   gen ddy240g=ddy240*(dy240max-dy240min)/(ddy240max-ddy240min)
; J% H$ n4 `8 o/ [6 w  h- H0 \   
5 ^$ c/ D' v1 B( Y6 c   sum ddy240g
   
   gen yxis3=r(min) in 1
! P5 j+ m% d* \   
   gen yxis4=0 in 1
) ~( f5 X4 a8 t9 K1 b( `; i" F" b  
  *ddy240g
   
   
  preserve
  
  sum v5 [aw =avol]
+ Z( Q- C) s- A. t  
  keep if dx240<r(mean)
  
  sum ddy240
  
! Y+ i8 B  k/ J) I6 C9 P1 N& T  if r(min)<0 {
4 ^: r( |5 f  {- Q( a        
                keep if ddy240<=0
               
                count
* |% v' X+ z$ g8 Z               
                if r(N)>2 {
                 
* H, a2 T2 j# e3 o  C; b. w                 pkexamine dx240 dy240
3 p$ q* ]( Z! \% z$ ^7 [7 @* f               
                local tempminp=r(tomc)
4 l* }. z) O5 k3 ]. p$ V                restore
9 p# r1 g$ ]8 _                 }
# [! o9 g( I# T                 
5 m( r/ a/ |: z; _' t2 s  q  V        else {
               
                 di "cannot find minp"
* l- b( I% M, v' |7 @/ l; C                 
- s5 Y6 S1 i" j# W: [! K                restore
' [% Z6 p3 K$ v! m3 L6 K                 sum dx240
' H: ?/ {3 |- p4 Q& t! c: E  V# `8 G                 local tempminp=r(mean)
; ?- O# R  g$ B, k3 i5 T                 
) `3 T& u: k' k$ \( v, v$ [                 }
; w* ^6 ]4 P, b                }
   else {
* S! X' [+ [( A8 P& x/ o4 B' ~                 
: q3 V# O. @3 [: F                 keep if ddy240<=r(min)
                 
9 l! U: u/ A9 q; X8 o                 sum dx240
                 
- G  V  ~, |2 s7 |4 ]7 C2 l                 local tempminp=r(min)
                 restore
                 }

* e  N$ G$ H2 _5 R3 L
3 L7 ~" \/ U' _6 R# H% o4 m
& x& Y3 C# Y& A* L. s   gen minp=`tempminp' in 1
   *240 minp


   
   preserve
   
3 k9 |' ]' d( P' f   sum v5 [aw =avol]
. c6 B- T" v  b# z$ ?2 K8 e% s" P   
   keep if dx240>r(mean)
7 H( ~; {9 }9 _0 M2 ~% Y6 ^( W   
   sum ddy240
6 q; J* a6 m' z+ r* a0 J) Z   
1 s, X! u# u; o& n( ]   if r(max)>0 {
6 u* U3 t4 v  B5 W               
                keep if ddy240>=0
               
                count
; }  C" C8 P" k  \               
9 E. Y0 k; U5 m2 L( b1 I" F. B" w0 D                if r(N)>2 {
" {$ q/ N( R: x: L% U5 i; N                 
                 pkexamine dx240 dy240
1 L. v' D7 @: n               
; X3 p) ?3 O: Y( G                local tempmaxp=r(tomc)
: y2 F3 Z) z* m) H                restore
                 }
  e, [2 U3 v+ V                 
                 
        else {
: F9 ?& F7 Y) v+ [# r% ~$ |/ |! p                restore
* x* d" T3 \1 o, I& R                 di "cannot find maxp"
4 \: H2 O; E* z/ f1 n9 Y( u               
6 G: P" W2 C4 H% ~               
                 sum dx240
+ u" d1 {+ y4 C' i9 K" l  \                 local tempmaxp=r(mean)
                 
                 }
: q  `# V. o) E' R/ u. ~/ C                 
4 z, n8 r- w  a0 Y, U. c" {2 Y! o                 }
4 U2 N1 H' m, f8 T, u                 
/ {- g1 c( G2 D1 N    else {
       
           keep if ddy240>=r(max)
7 T# ?$ w( O; K- B) |0 i          
6 W4 d. \/ e  l& I           sum dx240
          
- K! u2 [1 i' y: x3 ]           local tempmaxp=r(max)
           restore
6 n  W! F* P  O           }
       
: b" t" b# C$ J1 h. l5 k
; ?) k' d6 l4 R" P/ N! R, R  
& n: U  U8 ]5 D, t- R
   gen maxp=`tempmaxp' in 1
( _; w- G5 H/ j% i: L/ H& N   
   *240maxp
  H  g3 \9 ]4 Q8 V2 d, _   
   
   *-------------ddy240pminp左侧高峰值-----
   
: s+ s  A# X/ a. w2 e    preserve
   
3 T$ z+ x: o  j- c4 D/ c% x   keep if dx240<`tempminp'
$ H' P! g% c$ u9 k   
, y4 g  r% c$ a- }" r   keep if ddy240>0
   
: \6 X: X2 T0 m& a3 ]% q7 [, U   count
4 r9 c* `3 Y* y3 q               
# W6 z9 c6 }- a! a& |) ~                if r(N)>2 {
: c2 Z, j/ O$ p, J                 
9 P& a% {+ t. b$ T                pkexamine dx240 ddy240
+ R+ g& p* i! ]% A6 v' d7 c               
& m( w! }3 G$ d+ S                local tempminphddy240=r(tomc)
                restore
- o/ U- o, z$ e2 v4 }; V8 @                 }
                 
$ u# W8 H0 _) q& j$ W4 S        else {
                restore
0 u2 ~) D3 [3 B' m- }- ?# j                 di "cannot find pminphddy240"
                 
               
: ^+ d3 e6 E1 \  s. j; _$ j                 sum dx240
                 local tempminphddy240=r(mean)
: ?& j, r0 O* g$ E! ?' o               
8 J9 X7 U- l" S6 H. X                 
                 }
                 
) H6 {/ o2 Y* J8 Y0 G   
   
   
   
   
        gen pminphddy240=`tempminphddy240' in 1
7 q: r% t" A7 M) B& ]       
       
! |! ^8 `8 x5 `" ~        *-------------ddy240pminp右侧低峰值 -----------------
    preserve

; N/ A+ B5 O; e: z) c8 F& m0 J   keep if dx240>`tempminp'
   
   keep if dx240<`tempmaxp'
8 E7 m4 M/ f# r7 e  T   
7 Q" r5 B( B+ R9 E' O   keep if ddy240<0
   
   gen temddy240=-ddy240
% D% z  V3 s* S* d   
    count
" h8 J+ M6 e; S               
                if r(N)>2 {
                 
+ N# a# [+ L, U. y                 pkexamine dx240 temddy240
! [& G3 K5 g3 h; _) P: s   
  local tempminplddy240=r(tomc)
                restore
6 G4 E, W5 P% J1 o                 }
& p, s% I- [/ K0 @/ \6 H1 }/ [                 
7 o' d  ~4 V+ y        else {
               
                 di "cannot find pminplddy240"
                restore
                 
                  sum dx240
                 local tempminplddy240=r(mean)
0 J& P- U# `- d6 p                 
8 F# @0 ]0 _% A) S) _                 
                 }
, p& x2 U/ _7 \5 N: [   
   
! S" o, I- S; J! U: ]" V   
   
   
) D) s9 b9 w( O; G        gen pminplddy240=`tempminplddy240' in 1
       
3 n2 o( }5 L4 {6 E, J7 P" a# @        *-------------ddy240pmaxp左侧高峰值 -------------------
       
         preserve
   
   keep if dx240>`tempminp'
" {- z9 C" [& Q, R% g! Y2 L   
: E5 e# E" i: A; g   keep if dx240<`tempmaxp'
! Q* f# s! Q% u9 v; q* S* T- r, r   
8 N" K  I! C9 p! y   keep if ddy240>0
   count
+ _4 E8 c& v" c1 X7 @1 k+ N5 x               
                if r(N)>2 {
0 l" ^, \; C, n$ a# t' b0 J" T4 }                 
5 B! a0 |/ K( [6 H0 ^2 D# N                 pkexamine dx240 ddy240
( u4 Q# b. [4 @0 ~# }   
   local tempmaxphddy240=r(tomc)
" n, g, W$ [1 a% {1 @3 ^) B   
   restore
3 _4 M- }' g8 Y* N! w0 e                 }
6 U/ a; A* p% D2 u% \4 j                 
9 \- k. u' J: d4 v3 b8 N2 s* x        else {
; A) H  Q, ?: }               
                 di "cannot find pmaxphddy240"
                 restore
               
2 X3 `1 O2 k! q$ R( R                 sum dx240
                 local tempmaxphddy240=r(mean)
$ M( m/ ~1 H& `- b) N+ p  ^6 c                 
$ n0 l! q8 D8 V4 G' a                 
+ P; J7 x, T* R* K' r7 ^                 }
   
   
! b. i4 i+ O- U3 K% O+ |* w$ F2 J   
( w8 r( G' A& C+ Y   
3 c6 o2 u( w! t3 ]) s; G   
( Z# V8 V8 z1 H. G/ @  ?+ F  ]   
        gen pmaxphddy240=`tempmaxphddy240' in 1
       
       
        *-------------ddy240pmaxp右侧低峰值 -------------------
       
7 p" n3 M& X7 ?$ `0 F/ B% d$ `         preserve
  
   keep if dx240>`tempmaxp'
   
   gen temddy240=-ddy240
   
& B0 [" C0 {, X+ |$ T$ G* B   keep if temddy240>0
   
5 c, d$ F9 C* O( q' L   count
               
                if r(N)>2 {
                 
9 M$ @' ]2 b! u9 E                pkexamine dx240 temddy240
   
1 t: V' Q) E% M$ T4 p  z$ v5 d# |3 M5 e   local tempmaxplddy240=r(tomc)
2 L; R' o6 V" K   restore
                 }
                 
        else {
                restore
                 di "cannot find pmaxlhddy240"
) U4 O; l, v5 X7 I6 i6 i% a                 
               
                           sum dx240
) ~8 N# K) G$ E' _' J2 {: V                 local tempmaxplddy240=r(mean) in 1
               
4 u$ P5 w; I/ }2 o                 
                 }
+ p2 L: T, W! P: ]: e' a# v   
   
* s, D' {$ N% o3 Z0 G% p: b" C   
   
   
   
& P! l$ }3 S7 M' _% {; ?8 J        gen pmaxplddy240=`tempmaxphddy240' in 1
# O5 |$ J3 e* g$ ]/ J2 w) [# E1 D       
        *-------------------------------------
+ v, `+ w4 {) y  _4 S) l& y1 t   
- y$ Q2 D7 u3 a" @   gen price5d=v5 in 1/20
8 C7 a7 `9 x" v      
   gen price5yy=_n in 1/20
   
3 S7 W7 y) x: [7 T7 X$ v   gen price5y=price5yy*(dy240max-dy240min)/20
   
   *price5y &price5d
   
( e) v; I$ }: y3 v1 {. F0 O   gen priceny=price5y in 1
6 o( _  V* D/ E* ?2 H) [/ ^   
   gen pricen=v5 in 1

     
   *priceny &pricen
4 E2 f% z8 Q/ Q0 A- y   
   
, o4 ]' z, f" ]& K   *-------------------------graph------------------
  
  tw(line dy240 ddy240g dx240)(spike yxis mean240, lc(red))(spike yxis maxp, lp(-) lc(black))(spike yxis minp, lc(black)) (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)
- e) i4 F; h, k1 x) L1 J4 _: R9 d
1 N+ ?5 U& I1 [& _" j* ---------可选项目----------------------
- x6 j$ ?- w. O% L*(spike yxis amean,lc(blue))(scatter yxis2 amean, lc(blue) msymbol(+) mlabel(amean))
) q4 C. H( C- p& ^4 O% b+ K}
5 T7 Z  b8 V% `6 R: gend

maodeng

awer 发表于 2013-2-5 14:40
: e# J% s4 w8 J; k9 U' Y抱歉，最近太忙了
$ ^- g( b- `/ V& q& W4 ~8 A# U0 L* H
300236 我用下面的程序试过，没问题

* p! X7 ^; q( @' W6 u哇，谢谢，我仔细看看，这些字母数字代码晕了

喜欢喝冰茶

stata是不是商用的，R里应该有类似的东西吧

喜欢喝冰茶

awer  老虎推荐的一个处理数据的软件，很小，只有10m左右，很好用  发表于 5 小时前

7 T& C. p) {1 k* E) x5 J
4 p/ g+ H5 w" V这东西是开源的？知道有个统计软件叫stata，不过要收钱。

汉水云潇

准备试一试
& j% E3 f* s" U0 n5 b先找Stata软件

汉水云潇

晕呀
; U! E  F! [& X8 E/ \) u论坛复制程序代码
3 O; e9 y) j" @4 `" ?' e: D乱码夹杂在其中
很烦
. @& X, j3 _0 k6 ~3 O$ X& o! ?对了两次
5 Y2 K, ~+ p1 ^9 I以为没问题了
7 f, v8 J! W3 _7 t
结果
! X0 h: N6 }! h9 E4 m执行后
还是出现
! l! Y1 j7 ^6 H2 v) [: _r 199
的错误

汉水云潇

试验成功，下面是三一重工的图

peterbev

先用用试一试，不知道准确不

中秋不圆

噢。。。这个帖子好，比较直观。

bknite

这个帖子值得学习。

jeanny

感谢！！！

Hubert2000

明月回春 发表于 2011-8-13 20:17
- Z: W# k3 V9 n( ^3 c  J% U程序写的不错.
不过有一点要注意.通达信导出的数据,只对价格复权,并没有对成交量复权.所以数据需要手工处理 ...

老哥，也是山寨这篇文章找不到了，不知为何？