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

maodeng

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

maodeng

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

maodeng

. hello 300236
statement out of context
r(119);

大师 就是这个

sergio

这贴实在！

awer

maodeng 发表于 2013-1-28 13:13
7 v* @. h+ m, n/ ~) [& \. hello 300236
statement out of context
r(119);

6 O- c! }$ ~! U& B9 S. q抱歉，最近太忙了

& l2 c- ]4 j* Z$ I( i5 }% f$ A300236 我用下面的程序试过，没问题

/ H# r9 x0 V$ y3 U, Y你把 D:\stock\test\source 改成你导出的数据存放的目录，再试一试
/ q# S2 [0 I' R
capture program drop hello
5 q0 q% x. V! D4 B; h* q1 J2 M
program hello

0 o0 @6 S" X3 ~7 q1 o  I quietly {

clear

cd D:\stock\test\source

insheet using `1'.txt

   keep if v5>0
- a6 i7 D4 g9 D& W2 F9 B   
   gen vt=_n
   
# e) S; ^) r# L  m  I! v- [   tsset vt
! W# z6 i4 b1 r- o+ z/ t: J
   gen zf=D1.v5
1 D( ^( f# B- \8 z
; C, v- ?, B/ f7 f   gsort -v1

) G" K& X# ~& `  Q% v: C7 G4 y* v   gen avol=v7/v5
* f* D8 J8 a% A8 s+ C   
   keep if avol>0
   
   drop vt
  
   gen vt=_n
' s; J" n5 j  V! P/ J9 x   
   tsset vt
   
    sum v5 [aw =avol]
) M9 C# P; C- {% I, D0 w8 W6 |
   gen amean=r(mean) in 1
   
   * 总加权平均值 amean
; c0 r$ ?3 L- E: [, f3 T( h% L3 V' J   
   
   keep in 1/240
& u- ]& F8 E. J. ?  A" L- Y5 e
0 K( k1 L( e& {  i3 X: C/ o   sum v5 [aw =avol]
* h5 ?- U# q* y; t/ z4 r
   gen mean240=r(mean) in 1
   
/ Y6 V( M/ e% g0 r   *240加权平均值 mean240
% U9 N& n, u5 I" a9 X% f   
( Z( \) f# H% n2 C, y, ?   set obs 500
   
   kdensity v5 [aw =avol],gaussian n(500) generate(dx240 dy240) nograph
( l0 d+ z( J: Z% Y- I+ N   
( B4 R) T/ K! p6 C; j7 k  U  a4 ]   *240 kdensity(dx240 dy240)
  R6 U8 _+ W9 W# x' q* q0 |2 w& f   
& z. ~4 {' `1 B2 F   sum dy240
8 {  j  Z; f( c" M   
   gen yxis=r(max) in 1
   
   gen yxis2=yxis/2
+ x1 S' W, t+ s9 W: o% x4 c   
! ~2 {3 b9 F/ p7 ?; G   *yxis
   
, C2 f8 K0 f- t, n% e1 ^9 z3 i   gen dy240max=r(max)
   
   gen dy240minp=r(min)

" o- |) G  `3 B& {  a* S' |    gen vtt=_n

    tsset vtt

   gen ddy240=D1.dy240/D1.dx240
5 r* V0 m7 o- U* {   
+ P2 L, c" u2 R- |; e. _; q) l7 E   sum ddy240
   
! X5 A) E; I% M. B% x1 O   gen ddy240max=r(max)
   
* k4 ]& J3 K1 W* Y# _   gen ddy240min=r(min)
     
   gen ddy240g=ddy240*(dy240max-dy240min)/(ddy240max-ddy240min)
   
   sum ddy240g
   
' N: h% M- W  |' e: B8 k   gen yxis3=r(min) in 1
8 X% I2 p. A/ M. t: V3 g7 K   
# X% X7 ~: I; S, B. P* p   gen yxis4=0 in 1
. `3 X, w& R& p: l  
  *ddy240g
8 s+ v) Z0 u7 A7 P; g- x3 R2 g. W7 Z   
' ]" X4 s! S: E3 f) U. S. S, `   
0 a% `* _; W7 `4 i6 O( J" i  preserve
4 h1 \7 `* d- b, d/ w, Y) n  
/ M: p/ M4 Y4 v6 X" B+ v  sum v5 [aw =avol]
  
6 {1 _( C! s+ B: E; o2 t; \  keep if dx240<r(mean)
9 _- Y3 S0 j" @6 i  
  sum ddy240
  
) x) _& v! R% h; X' d  if r(min)<0 {
        
                keep if ddy240<=0
               
                count
2 B# i: o% ]8 u! n               
; K6 U" R8 f3 c6 R5 N                if r(N)>2 {
                 
+ {+ L, B, m8 n8 B! f2 Q1 `                 pkexamine dx240 dy240
               
                local tempminp=r(tomc)
* Y5 C4 P: B: q8 X' T/ w                restore
" q" C, H5 D& s: F' h& }3 e                 }
* ~0 i% D! m2 Y# i" D                 
+ B9 L: ~1 d1 J8 q7 M3 b* n        else {
  b- l+ k1 `# t9 A7 L3 P! v/ P               
                 di "cannot find minp"
8 o, N" g  X1 M" I) ^  G0 {                 
                restore
% a# C. t& m5 o/ T/ y0 o                 sum dx240
                 local tempminp=r(mean)
                 
                 }
                }
   else {
: |+ X! I2 E# B! H1 C                 
                 keep if ddy240<=r(min)
# C- a) G1 Q/ L8 }4 k8 ~                 
                 sum dx240
                 
                 local tempminp=r(min)
/ l) v: _8 T- S- L% G& a5 T                 restore
: `* x# Q/ Z+ z% F7 H4 Q                 }
$ P; d# G7 W8 k  D1 O& r8 o
* W, s$ P  j- ~$ N& A* w$ ]% o* a# y
; c+ \6 o6 T/ C0 s; \
6 J) G4 F/ I4 F, M8 B   gen minp=`tempminp' in 1
   *240 minp
  i8 l( R+ w$ m, h: h& m

   
   preserve
: p0 r4 [! V: q1 J5 n& x1 Q   
   sum v5 [aw =avol]
   
, \; s5 z: J/ Y* i0 i+ _0 _   keep if dx240>r(mean)
  V- V8 r5 [/ b  f: S( D   
3 Q# C  G4 G# ]   sum ddy240
   
   if r(max)>0 {
( I- G, j3 ]; N5 Q. E7 J- N+ j               
                keep if ddy240>=0
               
                count
6 o# G$ C+ j0 T1 k  a7 O$ ?8 L3 l               
                if r(N)>2 {
                 
                 pkexamine dx240 dy240
               
9 l5 D( c: |# f) ^+ o                local tempmaxp=r(tomc)
& o0 m* U; v* E, z4 ?                restore
+ b% Z& q/ n+ |9 T3 l. s                 }
) q: ^' }( w6 e( l2 w                 
7 h9 z; t) W7 K/ `" @0 l' C, `                 
        else {
2 t9 w4 ~0 ~! k0 `$ c                restore
                 di "cannot find maxp"
/ k7 N( J* p% r/ D9 H7 D5 v               
               
9 R6 j* Z5 d* b" |                 sum dx240
                 local tempmaxp=r(mean)
                 
                 }
# `: N% n" s$ e0 Y                 
5 \/ \# |0 d& _6 b" E/ c! Q1 a7 d; X- K                 }
                 
    else {
8 T1 P0 {: m- D1 p0 s: B; e       
           keep if ddy240>=r(max)
# R1 }2 D  y2 m8 H          
$ g3 G+ t- [' z8 J2 U% E% w! M           sum dx240
8 n2 _# M7 y% }: ^& n2 c          
           local tempmaxp=r(max)
5 x: A+ {7 @  M6 l% ]; ^           restore
           }
       

  

' |( \( u: `% K9 \" F- Z; B   gen maxp=`tempmaxp' in 1
   
. W- ]$ `! I% ]- |/ P- N3 L+ h, Y   *240maxp
   
' D9 K7 y% s" h3 D4 y. l6 ]; t   
$ U6 _" X0 q3 Z1 b) \6 @+ P   *-------------ddy240pminp左侧高峰值-----
8 ^* F& T9 l, r" E. ~! E7 H   
    preserve
/ |6 D/ I8 w7 ]% n" q   
   keep if dx240<`tempminp'
  p+ G( D5 j: D3 X( |   
   keep if ddy240>0
   
   count
3 ]( n. X3 Z6 g% Y               
  J' V6 E( A! I3 I                if r(N)>2 {
5 _3 U/ \' O$ [# m" n                 
                pkexamine dx240 ddy240
4 I1 n# d9 d, L1 X! M0 ]! Y               
# j! f, H# n. _# G                local tempminphddy240=r(tomc)
6 I, i! ?6 h( \2 L                restore
4 G! Y2 ?7 h  Q5 R                 }
0 S' f5 |- m5 v2 }2 k8 @- C                 
; C- }! Q2 V: C' q9 @        else {
                restore
                 di "cannot find pminphddy240"
                 
1 e. E# m+ m! Q( t               
& r. c1 ^, L6 N- G% {                 sum dx240
, j* G' f, _3 m                 local tempminphddy240=r(mean)
+ y6 m8 c1 O- I1 K/ v  j               
1 e/ h, S. b3 N: M" S                 
) u. \2 l# q: e2 K                 }
                 
) V! q/ q" x3 ]! n' J/ g   
) f4 s# R; g, h- P" W! i5 H0 C   
   
   
   
# ]4 X( z. |  O6 i1 O8 F" L( G9 ?        gen pminphddy240=`tempminphddy240' in 1
' C/ S  E6 L0 c       
* h4 Q8 V& E* ~) p2 e) T: |# u/ x       
/ w/ @$ a" R, j; X! k        *-------------ddy240pminp右侧低峰值 -----------------
0 `+ z6 M& W- B+ d$ F    preserve

5 m9 J$ B% N! G8 m+ i   keep if dx240>`tempminp'
. ?# e$ ^3 {6 \' K   
   keep if dx240<`tempmaxp'
   
) t' V! v9 K: r* f( P   keep if ddy240<0
   
7 D6 X0 H& m, o3 w; P   gen temddy240=-ddy240
   
    count
               
# E/ L2 I) @4 Z: o6 z$ ]0 X4 W" b                if r(N)>2 {
                 
                 pkexamine dx240 temddy240
& Q3 d5 K) Z+ y2 [) u, d1 _8 f/ x3 x1 p* e   
/ u5 p7 E, S3 j  local tempminplddy240=r(tomc)
+ |! M! H* x/ G                restore
                 }
: t% t8 y" S- w: k                 
' U- ^& U+ }+ p& Q) x        else {
               
3 C/ Z) p# Z. s( l                 di "cannot find pminplddy240"
+ g* K) Q# @3 M/ m                restore
                 
                  sum dx240
                 local tempminplddy240=r(mean)
                 
                 
( o* x5 k3 [8 Z9 M& M! S: w" m% J                 }
   
   
   
( z2 ?* V/ F" |3 z6 v9 B% c   
   
! i4 l" T" ~) X" c% z3 f        gen pminplddy240=`tempminplddy240' in 1
( d' T7 j5 y- j. ^. q* C       
        *-------------ddy240pmaxp左侧高峰值 -------------------
       
6 z9 w' J4 H+ e7 ^: X         preserve
# y" j& V* ^1 }: [   
   keep if dx240>`tempminp'
1 [6 V& Z1 g' B   
   keep if dx240<`tempmaxp'
1 C% r' I9 E/ C' b   
   keep if ddy240>0
- t! S: \& [& N   count
               
                if r(N)>2 {
" E# v  c) J; o. ?                 
                 pkexamine dx240 ddy240
   
   local tempmaxphddy240=r(tomc)
. S0 e& S/ d2 _" Q   
   restore
9 y) G+ y, S% i# x9 V                 }
7 _  U4 N2 O0 ^                 
        else {
               
                 di "cannot find pmaxphddy240"
                 restore
* V2 b3 }+ `% e               
2 @+ a# g( B- V) c/ S7 G                 sum dx240
                 local tempmaxphddy240=r(mean)
- D: ^* J9 c4 _; o                 
                 
7 S( y4 K9 k) L, O                 }
' K7 [  o+ o* u# B  K' C7 Z0 q   
   
   
   
2 ], y( w* t* i# N7 R( U   
( i0 _+ D( W% s8 @. R   
$ c) _, N+ G+ n3 E6 w& l  R6 K        gen pmaxphddy240=`tempmaxphddy240' in 1
/ c& W0 l4 g# p6 v$ [       
       
: F5 K0 \0 g1 B7 g        *-------------ddy240pmaxp右侧低峰值 -------------------
" W# |* R* t3 S! @9 U  F       
! r  t8 ?& O1 ^3 @         preserve
  W! b3 ]. ]* |9 p% e8 ^  
   keep if dx240>`tempmaxp'
   
7 Z( y6 N$ H2 u0 C/ S; M- n3 Y1 E   gen temddy240=-ddy240
. I! Z# |) n* D; v% ?, H9 l   
1 P4 T3 ]5 u; b2 ^# c0 O   keep if temddy240>0
   
- A3 @1 S4 _9 e   count
2 @( p& S0 ?( {# n$ s               
                if r(N)>2 {
                 
# \7 V& ]% V3 I5 u$ d, L. X                pkexamine dx240 temddy240
) F2 p+ Y. x# _% I" t( S7 {   
7 h  Q# N8 K' x% K3 Z6 {   local tempmaxplddy240=r(tomc)
4 {5 L0 |* K3 a% Z   restore
                 }
                 
' S" P0 _9 ~+ i7 t        else {
. N% d& v# I% a8 P- K                restore
                 di "cannot find pmaxlhddy240"
+ O$ s1 I6 X' b& C/ t, }                 
$ V6 |- e/ I$ ~4 Y, {               
4 c1 [! P* P) B4 P                           sum dx240
                 local tempmaxplddy240=r(mean) in 1
               
                 
. ?7 V6 Y9 u" J/ J1 E                 }
3 P& P# q* M' |4 R% @3 k   
6 `* t, Y7 A, i( B   
  H2 Z7 C6 Z% S* F/ `! U# {; O0 I   
   
0 ~- K- K! Y* }* _) i   
% s* I  K- V% ~   
        gen pmaxplddy240=`tempmaxphddy240' in 1
. m; e# c, c3 m" r) ^       
9 U2 d: F- u% x7 b. p1 c) ?        *-------------------------------------
   
   gen price5d=v5 in 1/20
      
   gen price5yy=_n in 1/20
3 Q* @! U2 {) y7 f7 T" e   
, x. Y& b. q8 D) I   gen price5y=price5yy*(dy240max-dy240min)/20
   
0 V6 S2 f  i8 _" `# }   *price5y &price5d
' `9 K; [; t& b' X9 v; p   
) K" e: r8 h1 H9 S8 f' `   gen priceny=price5y in 1
% `; S) D8 R6 m: x8 f) P. i  U% X  s   
   gen pricen=v5 in 1
8 Z0 W1 V% c0 z6 a
     
   *priceny &pricen
   
3 P9 H8 z9 D, L  s9 l* c   
' x- N7 X9 c* @   *-------------------------graph------------------
$ z0 k# R4 H. Z3 G1 I  
  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)

* ---------可选项目----------------------
: N+ U$ L/ o4 m% f3 V* z" I*(spike yxis amean,lc(blue))(scatter yxis2 amean, lc(blue) msymbol(+) mlabel(amean))
}
7 K$ d  P  b  j" cend
4 g+ ]- @4 o# ]( `( i

maodeng

awer 发表于 2013-2-5 14:40
% i- Q% S3 v& V3 B  o抱歉，最近太忙了

300236 我用下面的程序试过，没问题

哇，谢谢，我仔细看看，这些字母数字代码晕了

喜欢喝冰茶

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

喜欢喝冰茶

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

# L0 t8 J# }+ r5 i, X  B$ @& Q这东西是开源的？知道有个统计软件叫stata，不过要收钱。

汉水云潇

准备试一试
先找Stata软件

汉水云潇

晕呀
! k  L1 O# Q$ _1 D" e论坛复制程序代码
乱码夹杂在其中
很烦
对了两次
" ]& J" J' `2 t1 s" }! Y4 g' r以为没问题了

% L3 D$ l1 a1 L) E7 `8 _7 h结果
( H7 j+ P# K8 j. z/ \; |! f执行后
还是出现
) u% c4 G' U, ^+ o2 i. Rr 199
的错误
& Q2 B6 a* J& I

汉水云潇

试验成功，下面是三一重工的图
4 G4 Q& N9 V! b
3 }( h# |7 x$ o3 j, M$ W, w

peterbev

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

中秋不圆

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

bknite

这个帖子值得学习。

jeanny

感谢！！！

Hubert2000

明月回春 发表于 2011-8-13 20:17
程序写的不错.
7 V3 u2 v  k: M" B$ |不过有一点要注意.通达信导出的数据,只对价格复权,并没有对成交量复权.所以数据需要手工处理 ...

4 Y, @% V0 M6 E1 S# @% n1 x6 w老哥，也是山寨这篇文章找不到了，不知为何？