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

maodeng

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

maodeng

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

maodeng

. hello 300236
1 }) J9 O" Z8 B. C7 ?: ostatement out of context
5 x4 A, w! K( Yr(119);

大师 就是这个

sergio

这贴实在！

awer

maodeng 发表于 2013-1-28 13:13
# \# S& p. p' O; h5 c4 f. r. hello 300236
' F% c1 n' P, V7 Vstatement out of context
7 G1 i4 {* k) n9 W8 D% C* J, ~6 Rr(119);

抱歉，最近太忙了

+ E3 P) K& R, x300236 我用下面的程序试过，没问题

4 d' B$ J7 E- O/ ~/ i) B  i你把 D:\stock\test\source 改成你导出的数据存放的目录，再试一试
2 G3 h7 H3 N2 G
capture program drop hello

% s9 P( d  ]* g" bprogram hello
# M* q. y/ O. e2 m' ]5 w4 c
quietly {

clear
3 q% \* L7 D+ t, t3 i0 k& u1 d) P& {
8 ?7 a5 F5 D3 ?8 H0 _cd D:\stock\test\source
% F4 k& z' W1 A. P$ |2 \, y" S
6 K0 o, z/ Q) L; finsheet using `1'.txt

! F* {5 q& h6 ~; y0 l0 V3 ~: i   keep if v5>0
   
" Q5 ?$ E8 V( \& L; a0 ^) V   gen vt=_n
   
   tsset vt
, |& O2 s& t' N0 q
  l* `6 y. j" Y   gen zf=D1.v5

   gsort -v1

# \8 l& `5 e: D5 x7 d   gen avol=v7/v5
   
   keep if avol>0
   
, @% N, L! U0 t   drop vt
; V1 I6 l! i6 w! c3 T  g" t  
   gen vt=_n
' x* ]# {- e3 b- l. `" Q& k   
7 a: s% C/ W/ P0 q# G   tsset vt
: ]# u" b' w, V4 P) v   
    sum v5 [aw =avol]

   gen amean=r(mean) in 1
   
   * 总加权平均值 amean
" q7 q( O: n, V3 |# g6 S: \7 V   
% j; Q" A) Y1 O, l! w& h3 l   
   keep in 1/240
; a- X- H' I# Y3 z3 b, T
9 O5 M) _, I: v0 @( f7 @8 K. G   sum v5 [aw =avol]

3 r2 d/ M0 h5 L   gen mean240=r(mean) in 1
   
* |" q% \3 c1 I   *240加权平均值 mean240
) y  a1 D7 u5 c/ T   
$ z7 j# n" w" E; b  D   set obs 500
7 g2 {% B# ~; j/ y, l   
# j8 q' l  B: m; z* M  u9 U+ n   kdensity v5 [aw =avol],gaussian n(500) generate(dx240 dy240) nograph
/ |. ^. f5 z  o' j   
   *240 kdensity(dx240 dy240)
; J3 M% K0 x, R   
  z2 Q* h% Z2 i6 ~& k$ o   sum dy240
4 d9 p: h% |! e5 K6 W# F: q5 L   
   gen yxis=r(max) in 1
   
3 H1 U: j9 V% R; f! m( H   gen yxis2=yxis/2
5 s; n3 ?" D, E   
   *yxis
. d2 W) x& \/ n  x  n2 U! c4 B% ?   
  T7 q9 c! l* K" u* t: m% `1 l   gen dy240max=r(max)
   
% ~  c7 G* q3 A8 I( y   gen dy240minp=r(min)

7 z) b* @5 d$ R' Q$ L( G. I    gen vtt=_n
- _: f( c5 P  V' Z; f& s% I' _
$ f- R) J* ~8 M# C3 r: m5 b    tsset vtt

1 Z( L; L5 `1 j   gen ddy240=D1.dy240/D1.dx240
: U* N6 l5 l9 B% H& m   
   sum ddy240
+ G5 O5 P- ^) y9 J, k: V   
   gen ddy240max=r(max)
   
0 n! i: C5 K8 ~   gen ddy240min=r(min)
     
   gen ddy240g=ddy240*(dy240max-dy240min)/(ddy240max-ddy240min)
   
8 m: `/ z( q0 p   sum ddy240g
0 w5 X3 w3 J. M( K, H8 J( ~. I   
- J  F8 ]7 t# N   gen yxis3=r(min) in 1
! K3 h# o+ f6 @, P- T5 {   
* }; O5 b) K7 ]% x  }5 M% W   gen yxis4=0 in 1
1 ?+ O5 I8 R- r7 B* E  
  *ddy240g
2 f4 Q7 p9 Y' `4 }. \9 v   
   
  ?  b9 B& V8 F3 u" u6 }& \  preserve
) |8 ]7 [0 K  `  
6 l+ X% v3 h2 ^+ D2 B  sum v5 [aw =avol]
  
  keep if dx240<r(mean)
0 c4 v4 A# l2 \$ O  
( A1 M6 X7 v+ |  sum ddy240
  
) B6 b7 @  g5 R! {7 m  Q6 [) x/ l  [  if r(min)<0 {
$ ~8 Y- E& ^6 [8 r        
7 g" `! q  m+ K& f# i' Q. x6 ]                keep if ddy240<=0
               
                count
               
                if r(N)>2 {
                 
                 pkexamine dx240 dy240
" A1 ?5 ^$ i& Y               
# J+ O- x! i" e                local tempminp=r(tomc)
1 s. I$ z) r' C, f! ~& h                restore
2 Q$ k- q8 t0 P  q5 F                 }
& g* q6 W8 E; t6 L* x8 S% g                 
        else {
9 i; G3 a4 W, w               
                 di "cannot find minp"
! P0 Z5 k7 Y/ w$ |' Z0 U/ R                 
- t/ \. u% }. f7 s: f                restore
. G2 n6 ~) I: w' V2 n& ~9 M6 O                 sum dx240
' }" D$ d' Q6 }! ~. V+ G0 U                 local tempminp=r(mean)
                 
                 }
# ~# s, S" D4 P0 F                }
" z# R0 J9 C1 _: L; R- c   else {
                 
                 keep if ddy240<=r(min)
. ^: {/ `' h: O                 
                 sum dx240
                 
; l/ m. d" i3 P5 A, G' o' C                 local tempminp=r(min)
                 restore
- E/ m+ K( p+ X; |7 k2 G. F                 }
3 {0 I8 p- q6 _" Y$ \3 z) o

: q3 D$ k9 W0 T
8 w2 X' Q* `& q; l; `. y+ K   gen minp=`tempminp' in 1
   *240 minp
) r/ a' Z2 X0 w- B
( b% {: J' @2 I8 ^, x, b  J7 o- |
* I2 t7 l  g% I( x$ V3 B   
" X! \- m$ }/ D) {# n- r   preserve
/ r& Q, ~$ @9 @/ S   
   sum v5 [aw =avol]
5 r1 q4 ?/ O* \8 ~$ K   
  `7 x) R# q' ~2 P% M* M; W* j   keep if dx240>r(mean)
   
! J0 y& r- D% }/ W6 T. y   sum ddy240
   
   if r(max)>0 {
) q5 F' u3 ~0 \/ @5 e. e0 L, }               
                keep if ddy240>=0
               
7 f' U- g9 l6 z5 d9 o/ t' J                count
6 X$ Q7 c- \+ ?$ Q  q               
                if r(N)>2 {
                 
; W# J. _9 b% T+ J# o                 pkexamine dx240 dy240
               
                local tempmaxp=r(tomc)
                restore
                 }
                 
                 
        else {
% }, M# r1 [  s+ k& b: X                restore
                 di "cannot find maxp"
               
' j. w5 u7 B9 X3 P4 u  ]0 J6 Z               
. k5 c- `3 x( K- _9 j! e) ~                 sum dx240
+ b* X5 w' N; o1 O3 X. g4 F                 local tempmaxp=r(mean)
                 
; \0 J6 i) Y. j- M0 w& l                 }
                 
                 }
; I5 C! `" j& U6 R2 w% k2 {4 g& c                 
" E+ L" ^  }1 E. q+ Z4 k% n3 h6 e    else {
8 r- L# Z5 U# c- J( ]       
           keep if ddy240>=r(max)
: F1 u3 Q! @5 y* c          
           sum dx240
          
           local tempmaxp=r(max)
           restore
8 d" D4 \% A6 u' l# ?) w# B           }
& V- S" w) ]" y$ e       

# f# f  B# [3 q2 P4 [6 `  

$ L/ ~6 }$ t' L9 m   gen maxp=`tempmaxp' in 1
, \2 j* s6 @/ u' [   
, a2 A* u8 I' ~   *240maxp
   
   
   *-------------ddy240pminp左侧高峰值-----
   
    preserve
   
   keep if dx240<`tempminp'
9 n# F; Q7 w- v% Y9 g- L0 G2 {   
   keep if ddy240>0
   
   count
               
                if r(N)>2 {
                 
- o1 D4 z; G; o* _7 a                pkexamine dx240 ddy240
0 V9 h! |+ G% d' u3 K' n               
7 _4 B1 V0 p6 C" j1 w# S                local tempminphddy240=r(tomc)
                restore
  J1 w+ R6 V, R                 }
                 
        else {
                restore
( R3 T& N! L7 o2 _. F5 w1 H                 di "cannot find pminphddy240"
$ X) `" H- ^4 O8 n+ T. e7 e1 I                 
               
                 sum dx240
5 U+ E* n2 r6 ]3 v) A% y/ _% G7 z                 local tempminphddy240=r(mean)
               
0 b% D4 G. X& ?  Z9 }* k* ~& w' H9 E                 
9 t% G7 e. b3 L* x' O, y6 @9 q                 }
' b* l5 P" I8 ~( H9 Y% w  _                 
   
: |% V% F1 Z$ g+ Y- y- @( Q+ m/ V0 a   
" l* _! ?+ _9 B/ R8 q7 P   
/ _# t% p# R8 e, f% k   
9 N- r; {( g: i) {4 [) P   
' _0 e# J/ F4 q& c        gen pminphddy240=`tempminphddy240' in 1
       
       
        *-------------ddy240pminp右侧低峰值 -----------------
    preserve
+ s: v8 \! u9 \$ Y' o  S: \
   keep if dx240>`tempminp'
   
   keep if dx240<`tempmaxp'
2 e: R* ?. g/ h  q  E   
   keep if ddy240<0
   
   gen temddy240=-ddy240
   
    count
8 @# A7 v0 B6 g) }. P) |               
                if r(N)>2 {
- S% K/ X" G+ ^8 u+ Q' b6 `                 
                 pkexamine dx240 temddy240
' a( d1 o& \: V0 w8 n! u6 y! R+ p   
, m  @+ \9 ~( z* Y  local tempminplddy240=r(tomc)
" F0 A: W, H9 [: [                restore
: d+ [' ^( O/ M- J8 D                 }
) j* q2 J3 M& W/ ?                 
        else {
" }1 ?; Q  l/ ^& R" G: k               
                 di "cannot find pminplddy240"
$ i. j# c3 e( \( p                restore
                 
                  sum dx240
) L9 n$ O4 R, Y                 local tempminplddy240=r(mean)
                 
8 v% A  j. K! ~# Q                 
( l9 C; c$ F( \1 x/ g                 }
% u# w" J# B  |! z8 q8 A) \. n   
/ l: h0 F  j/ \0 Z   
   
   
   
        gen pminplddy240=`tempminplddy240' in 1
       
        *-------------ddy240pmaxp左侧高峰值 -------------------
       
$ {1 X0 J0 K- T, }  O+ G         preserve
   
* o% f# z" ^' L; Z9 k  v   keep if dx240>`tempminp'
   
5 {1 p1 ~4 D2 U0 ^/ J   keep if dx240<`tempmaxp'
* L/ s- M7 u2 f' h   
1 @; h6 Z, ]5 Z. h; a   keep if ddy240>0
" w+ I8 E5 F* {, ]3 W   count
               
8 Y+ S: f$ E# d3 u                if r(N)>2 {
, P7 ^/ t1 @4 D+ L9 U0 n  o) T6 Y                 
                 pkexamine dx240 ddy240
   
! ^- D; r0 t' K, u   local tempmaxphddy240=r(tomc)
2 u5 A- r, V/ r5 N8 u   
   restore
0 D$ F% t$ E$ P9 N& H                 }
  S* ?, G" i4 ^7 k                 
        else {
               
                 di "cannot find pmaxphddy240"
                 restore
* t9 m) j( @+ z. w1 e; X7 ?* q               
                 sum dx240
: w- ?) [6 B; G5 R                 local tempmaxphddy240=r(mean)
! Z- X  g5 f' v0 u2 T                 
7 P0 ^3 J0 x5 u# }+ z7 w                 
. z3 G+ `7 M. h, B1 P7 O! O                 }
   
   
   
3 S. N* G5 v( l6 y  O   
   
   
. N) x& \. p' u0 O. U3 t7 \- \        gen pmaxphddy240=`tempmaxphddy240' in 1
       
2 D  ?9 h9 Z1 l) f5 X1 m. R7 s. K       
        *-------------ddy240pmaxp右侧低峰值 -------------------
       
         preserve
1 U$ [4 S& a0 h4 A$ ?  
   keep if dx240>`tempmaxp'
: R% X4 E" U0 z   
2 F% w+ Z* j( v6 p% ?   gen temddy240=-ddy240
   
5 v" A6 d) P: w6 @3 r/ l+ _4 H   keep if temddy240>0
& J/ f% e+ }# L: E$ q" `   
   count
7 H' S2 C/ Y0 r' G, f               
7 g7 Z7 B& e+ p+ s                if r(N)>2 {
                 
) B6 n, u" T! g( x( ~                pkexamine dx240 temddy240
( R: x  x3 @5 ]- Y4 |   
( e2 k4 _  e! P+ D! k4 z# b   local tempmaxplddy240=r(tomc)
1 }& ^" D  C) ]: w6 `   restore
* C* r# r; w0 S6 o4 U                 }
                 
        else {
2 o2 E5 y& J5 a( {$ l' [% {                restore
                 di "cannot find pmaxlhddy240"
, ^! U- A; P: U! _% S, {$ M                 
               
                           sum dx240
                 local tempmaxplddy240=r(mean) in 1
' R) O1 @% k) i- L# ?) m5 d1 M7 p" h5 U               
                 
: z1 }+ ~% j& j' E. Y( T9 R                 }
4 ~: X, v/ [; ?7 s   
& ^) T; j2 g+ J+ p% V0 M* q   
   
. ]2 P+ N( _: G7 C4 H   
   
   
        gen pmaxplddy240=`tempmaxphddy240' in 1
       
        *-------------------------------------
4 ~# U. P+ S( T3 C: ]: a; ^8 _0 a   
4 J7 j" v0 m0 h3 |   gen price5d=v5 in 1/20
" K. N# H  F8 p( t0 U1 z      
  O8 f& u/ `% Z& K  K" ^4 d   gen price5yy=_n in 1/20
) {, g) x" D4 s' S2 |' W   
   gen price5y=price5yy*(dy240max-dy240min)/20
   
   *price5y &price5d
; }" j. S/ q* @, h' I) r: I) V   
: Q. o! L* J5 K; F   gen priceny=price5y in 1
+ t7 A' f5 J2 T9 P! g' U; m   
   gen pricen=v5 in 1

- H5 k' f/ y/ G9 D     
   *priceny &pricen
9 m% C* I4 Y, I4 r5 U$ r" h& z3 w   
   
   *-------------------------graph------------------
6 s. H: Z. |- Y4 m- A3 x  
# A; Q8 T* }1 w& ~  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)

/ W. p! O' P$ `+ g* ---------可选项目----------------------
! _+ P( ^! g3 x6 ?*(spike yxis amean,lc(blue))(scatter yxis2 amean, lc(blue) msymbol(+) mlabel(amean))
}
6 P2 |# u) j9 d( d7 N" x( zend
+ F5 Q* A; ?  Z8 c8 L. q. p

maodeng

awer 发表于 2013-2-5 14:40
抱歉，最近太忙了
6 \1 X( j4 O* D
1 }  I# `+ D. P: P  {300236 我用下面的程序试过，没问题

哇，谢谢，我仔细看看，这些字母数字代码晕了
/ m# F8 r$ G! C( W# j

喜欢喝冰茶

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

喜欢喝冰茶

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

这东西是开源的？知道有个统计软件叫stata，不过要收钱。

汉水云潇

准备试一试
4 l5 z9 Q3 `0 b! a, z1 E  R  o+ w先找Stata软件

汉水云潇

& I6 c5 f7 O' F晕呀
6 h2 \/ m6 q$ R4 }论坛复制程序代码
乱码夹杂在其中
很烦
对了两次
以为没问题了

结果
执行后
7 X: n; F2 w' R% I' J还是出现
8 a! v, c# A$ b! K0 Q, `- ar 199
4 Q' x$ D( U4 d  e5 Z/ a1 [  T: o$ c的错误

汉水云潇

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

! P3 E3 M; h, |% T% w

peterbev

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

中秋不圆

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

bknite

这个帖子值得学习。

jeanny

感谢！！！

Hubert2000

明月回春 发表于 2011-8-13 20:17
程序写的不错.
' v. [; E' I' u1 i3 n不过有一点要注意.通达信导出的数据,只对价格复权,并没有对成交量复权.所以数据需要手工处理 ...

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