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

maodeng

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

maodeng

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

maodeng

. hello 300236
. i$ ^$ Y- a1 i, i- x3 S- `statement out of context
; E6 `; j- r* M( y2 g" x2 W4 O5 Br(119);
0 O- X. |# U! g+ b  \
大师 就是这个

sergio

这贴实在！

awer

maodeng 发表于 2013-1-28 13:13
8 }2 U, g8 \; @7 N! ]. hello 300236
4 M/ o; M/ r% l( d& Jstatement out of context
5 k& n  |6 u1 Z- rr(119);

抱歉，最近太忙了
# A, n1 r+ T) o8 _
300236 我用下面的程序试过，没问题

你把 D:\stock\test\source 改成你导出的数据存放的目录，再试一试

capture program drop hello

program hello

quietly {
7 Y- F. \8 R4 s8 [4 G$ W$ b
% |0 t1 `1 r2 I, n3 aclear
3 c/ m$ k0 z3 r5 n
cd D:\stock\test\source

insheet using `1'.txt

   keep if v5>0
- s# k7 B( X  a   
! I' c; i" T  p% _' \3 j, f9 T   gen vt=_n
   
/ N" ^# N+ `1 r0 X# H- q0 d  ?   tsset vt

: H! f) _' [7 ?# ]) s  n   gen zf=D1.v5
$ }6 A% y- f& Z
   gsort -v1

   gen avol=v7/v5
   
   keep if avol>0
   
   drop vt
- ~6 E: L: K- R0 t/ T" w! Y+ u  
) K2 ?/ L) L  z' k8 j   gen vt=_n
! K4 f3 D6 z9 W8 w2 }   
. v7 ?' X; t5 n   tsset vt
   
    sum v5 [aw =avol]
. H: l" u6 e) Z+ O2 ^
% q% I5 m2 {; C& X. ?8 _   gen amean=r(mean) in 1
   
. ^. J8 x2 e7 U   * 总加权平均值 amean
   
   
: b) O$ ]; ~3 f# N   keep in 1/240
" W! z, ]0 @8 b4 ?  q0 S% }* \4 j
; t$ R$ {' y; G1 ^   sum v5 [aw =avol]
/ ]: w$ v0 c! L! a, B
3 M2 \: d- W+ j) o8 y   gen mean240=r(mean) in 1
" e2 M1 T8 Y8 S9 g+ [) ?+ E8 S& U7 {- @   
   *240加权平均值 mean240
$ U( p6 t- `3 p& a$ E   
, ?& F' x& k2 J, }' g* E   set obs 500
   
   kdensity v5 [aw =avol],gaussian n(500) generate(dx240 dy240) nograph
   
   *240 kdensity(dx240 dy240)
  f! E" n& Z3 _7 N   
, r7 \, F3 {9 Q- b   sum dy240
' g% F4 b# |4 d/ R& D/ X   
   gen yxis=r(max) in 1
& }7 D5 J& }+ M/ {3 f4 t   
   gen yxis2=yxis/2
% \& o1 [$ m* |   
   *yxis
5 e& h3 g" |6 A" V' C   
( q2 |" ^' P7 Z& u   gen dy240max=r(max)
; w% G' S) S6 @! f- q   
+ `' Q( ]6 |8 a   gen dy240minp=r(min)
6 ]+ E# n9 T9 D  o% [7 ^
" u* e  |& I" O0 \    gen vtt=_n
" ?$ o; W; A, x: h6 u$ Y+ r/ I
    tsset vtt
" ^: C/ i4 Y: ^: e% X' H. a: m
& @6 }, h& B! s5 b2 R: \" S/ a- w" ^   gen ddy240=D1.dy240/D1.dx240
   
   sum ddy240
   
) s. J: O5 j3 E   gen ddy240max=r(max)
   
7 n4 G$ _# n7 I6 n5 \; E   gen ddy240min=r(min)
& x6 s! u2 c* K' ]     
; \6 y9 c4 ?3 T' o; S% Y   gen ddy240g=ddy240*(dy240max-dy240min)/(ddy240max-ddy240min)
   
$ Z7 R6 |9 w5 V% Z7 N, }- O   sum ddy240g
' m# O- ?/ J6 B& _2 ?0 G! H& d9 H   
' I* {, t/ l8 L. t! n0 g+ K   gen yxis3=r(min) in 1
! U: j, q1 c+ u2 \   
   gen yxis4=0 in 1
  
  *ddy240g
* ?' n, Z9 I& i# d# Y- O9 W   
   
  preserve
0 ?5 W0 m' Z) E/ Y. v# h  
  sum v5 [aw =avol]
6 u! O# Q% V9 j2 F) ]; L  
& r5 w% \& L( S0 ?! F  keep if dx240<r(mean)
  
  sum ddy240
  
  if r(min)<0 {
" x. q  `8 v! ]- ~) c% `6 e. u8 J        
                keep if ddy240<=0
0 x; p' l) e/ O1 U8 B+ b               
                count
1 J1 C( J) Q% G9 V6 k) z               
7 x- v' Z/ ]' N, o2 N7 G1 n+ e  @5 C; X                if r(N)>2 {
                 
                 pkexamine dx240 dy240
1 a. U4 k; a- g+ w9 J* @               
                local tempminp=r(tomc)
  M  d  o. q4 O                restore
) T. q  p2 q0 h% i                 }
                 
        else {
" B4 [4 t9 D+ a/ c               
) c7 D. e8 s& n# Z6 c                 di "cannot find minp"
3 }; E/ B! p) t5 y  t, J                 
                restore
# B7 ~* z' q; c+ T1 n7 e                 sum dx240
                 local tempminp=r(mean)
                 
6 i) f$ `! T( i  X                 }
! M, ?1 @. \1 u/ h. B6 k' `' l& r                }
( k+ M; l# p) [9 m; Z5 s   else {
                 
                 keep if ddy240<=r(min)
                 
                 sum dx240
                 
                 local tempminp=r(min)
                 restore
                 }

  `8 u' |% p( K* ]% |! `* P

& y. e, S% C, e; V, E   gen minp=`tempminp' in 1
   *240 minp
1 L" u0 V, D0 ?, A; z4 g

5 F3 J" b: R6 I1 M   
1 V+ ~6 R* Z' o   preserve
   
  P) Z! t" H0 Z" @% f   sum v5 [aw =avol]
. q9 G0 z% }; B   
   keep if dx240>r(mean)
. d( k+ r% S. b2 P6 m   
1 \9 g: A! d6 N% ]& z7 O1 B0 j   sum ddy240
: _4 d% n% I. X" B2 Z$ @   
4 p1 k* q$ Q) ]" @( l   if r(max)>0 {
               
                keep if ddy240>=0
7 G+ }; G9 f! A2 G               
! a" J4 e) ^4 Z# E                count
               
                if r(N)>2 {
# K$ O, X7 z1 T9 `3 }                 
5 _& }" x. N" f) G; U) ]                 pkexamine dx240 dy240
               
5 P: x3 V/ n- B9 U. \                local tempmaxp=r(tomc)
                restore
  a+ k1 f5 a$ H# A# T3 f5 p                 }
6 _5 Z6 q4 d* e& H0 j1 n                 
                 
        else {
) |+ t# j* _, [; v- a0 N! H. c                restore
                 di "cannot find maxp"
% N4 n6 x) A" L$ G2 r( z               
/ t/ G8 m$ h& C3 W/ K( a* M4 T3 s               
! E  ]- |6 R# X0 f+ ^8 Z0 E9 i                 sum dx240
                 local tempmaxp=r(mean)
, A% `- M* E4 f) O+ W  f1 I                 
                 }
* g* ?% B' w  u                 
                 }
( V% x, q7 U3 ]7 ]                 
7 X* F/ p! D  p9 L3 N& D0 k0 P    else {
       
           keep if ddy240>=r(max)
          
           sum dx240
          
. ~7 _: J; B7 q: k9 P0 A           local tempmaxp=r(max)
9 m" G7 R8 ?) b2 y. z/ `1 F           restore
           }
       
& M1 e+ _6 T3 ^( l0 E
  

   gen maxp=`tempmaxp' in 1
   
) T/ ?  i8 Z% {# _! z   *240maxp
) s) Y; f- i7 a) X, t   
7 |7 r( s+ y6 N   
   *-------------ddy240pminp左侧高峰值-----
   
4 r- S# |5 ], ]. a5 _8 ~    preserve
   
5 J- m, X8 i6 L+ R8 S, }   keep if dx240<`tempminp'
   
   keep if ddy240>0
   
   count
/ X' F. _/ N2 {% X) {               
# N7 |# U# Q0 V; E                if r(N)>2 {
, L! z5 Z+ b% O' }2 Y& n) p8 i2 j                 
3 H+ e# x/ F3 N' E% A8 J- ~                pkexamine dx240 ddy240
               
+ w! w. e1 |0 S9 ^) K3 S+ w0 k, e# }                local tempminphddy240=r(tomc)
: D9 i( z6 x/ E. ^7 w4 a0 w# ]                restore
                 }
4 {$ _6 ~3 a( V3 y9 ]' {                 
        else {
                restore
3 f% S9 o8 \) m! V$ d                 di "cannot find pminphddy240"
                 
               
                 sum dx240
  Q0 B1 ~+ F8 T6 I$ l                 local tempminphddy240=r(mean)
* Z: S! m% f3 N3 p2 ~               
                 
                 }
                 
8 y% D- p2 A3 q; ^! W) c$ b   
+ M- }" |: a  D5 @& n& y   
   
- C+ q- k) T1 h0 y- l0 ^: O( Z   
   
        gen pminphddy240=`tempminphddy240' in 1
       
) M* I( i3 s; [2 B  A       
        *-------------ddy240pminp右侧低峰值 -----------------
! k  U" Q( ^  ^1 R8 y    preserve
- i2 M" p  `* u
   keep if dx240>`tempminp'
   
" Z% {9 m" x2 f# a/ Y   keep if dx240<`tempmaxp'
   
) G+ E# ]9 e& h  |/ y2 x* F   keep if ddy240<0
   
- H2 M. ?) `  n0 P0 q   gen temddy240=-ddy240
  i& Y: _, _4 q& h: p& e   
3 s+ d& j6 {! Z1 {7 Z    count
               
% y; j  C- `. Q# j! f' ~                if r(N)>2 {
                 
- |' d/ ~) l# k+ t/ y$ y' [                 pkexamine dx240 temddy240
   
9 f2 s7 X/ _3 R' t  local tempminplddy240=r(tomc)
, `  C9 k5 u5 e1 L" D  x$ l                restore
                 }
                 
2 m8 S* k# F* l) `! `        else {
+ a5 q. J# C+ Z8 W, [" N6 X% h/ s               
                 di "cannot find pminplddy240"
                restore
; L( M8 T6 w6 H5 {                 
                  sum dx240
2 C4 z% @8 k0 _, |6 ?6 `! r3 x                 local tempminplddy240=r(mean)
( W  m9 h2 F/ F. q8 p( z# p                 
! Y$ x# d  S0 Q! V0 \) |                 
                 }
   
8 M  ]$ g/ J5 L) B% V   
   
1 x4 t4 G& q$ T7 b   
   
1 N' J6 N( W8 D( X- ]2 w        gen pminplddy240=`tempminplddy240' in 1
       
9 C  f. C2 i3 s- b        *-------------ddy240pmaxp左侧高峰值 -------------------
       
8 O& p( j2 ~& y! d+ v! Q( G         preserve
) D- ]3 o, _) l6 z0 P   
   keep if dx240>`tempminp'
   
3 L, T4 f: `3 k. W" }3 p. Z   keep if dx240<`tempmaxp'
; w" C9 X. j: j, d   
   keep if ddy240>0
, n  l# t5 t1 B  {$ y9 j. S4 b   count
               
                if r(N)>2 {
8 h" q4 r4 {) K+ W                 
                 pkexamine dx240 ddy240
   
8 V8 b/ ~" c, k0 l' T9 K   local tempmaxphddy240=r(tomc)
5 ^9 [& |2 w8 d( F   
1 ]' p9 s1 j5 D1 A+ s   restore
3 V6 I1 `/ X% J( ~( g                 }
                 
+ u0 F7 Y, S" D. s+ d# X        else {
               
                 di "cannot find pmaxphddy240"
                 restore
5 j7 x& V4 I! j2 F2 f7 o) d: B               
                 sum dx240
                 local tempmaxphddy240=r(mean)
                 
                 
, T5 Y% W# W/ ]4 |                 }
; _- {5 ~7 j7 a; h" V% b; h   
   
   
8 x, T. g- B3 k5 _7 |4 l8 }   
2 x5 Q; ?$ {: L$ l- r' F1 ^4 z   
   
  N+ k3 f' U+ v0 U* Q        gen pmaxphddy240=`tempmaxphddy240' in 1
  P( g  I4 ~4 p8 _" ~; }. i       
       
        *-------------ddy240pmaxp右侧低峰值 -------------------
       
         preserve
/ U$ `5 R; Z% ?% p, u6 ?! k  
+ O# y5 n. g8 r+ K5 t   keep if dx240>`tempmaxp'
7 S. L1 M: _9 m4 J+ T! g  s/ ]; C7 X) Z   
" P1 c8 e% U$ E8 Z& l: E   gen temddy240=-ddy240
1 d9 _1 P2 K0 x# G   
; y9 j6 d- l+ }9 i9 l   keep if temddy240>0
6 [+ R) R5 P: Z2 P, e9 T   
  F  I: o( E) O9 ?& O- e   count
               
                if r(N)>2 {
                 
                pkexamine dx240 temddy240
; X* k3 N; t4 h" X& L   
   local tempmaxplddy240=r(tomc)
3 X9 ?, K4 }3 k7 \, c$ V/ D8 \   restore
                 }
                 
        else {
. _6 O" i! i" R3 x: d( G                restore
5 e# |( p6 w; K. X& o                 di "cannot find pmaxlhddy240"
                 
& H: i& n+ }5 f, ?5 |! `- ~               
                           sum dx240
, [9 q' ?5 N% v" u3 h3 G$ \. [                 local tempmaxplddy240=r(mean) in 1
               
                 
$ Y9 `$ v/ T0 x4 R                 }
   
   
   
# c: E! g9 A( Z- R: T( r5 H/ I   
   
5 Z- x' d; V7 m4 C   
7 C" n* I; y7 s  _        gen pmaxplddy240=`tempmaxphddy240' in 1
       
        *-------------------------------------
   
# e6 u" d, T7 I4 n  r3 p   gen price5d=v5 in 1/20
      
   gen price5yy=_n in 1/20
5 ~& Q: S" P8 r, Q) O. G   
   gen price5y=price5yy*(dy240max-dy240min)/20
( C/ w: }; S1 `% H) C2 }& ?   
   *price5y &price5d
   
0 b4 Z2 E6 Z7 r2 I   gen priceny=price5y in 1
3 w$ O, y+ }" G0 {9 R   
   gen pricen=v5 in 1

6 m& ^) a7 f+ r( L     
   *priceny &pricen
8 v8 Y1 E. ]) [& f" c   
, N  x$ P: [7 W# e) L2 G! n   
2 @7 V4 G% i7 y" b" A   *-------------------------graph------------------
  
+ I$ |( X' n( S  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)

, B2 C$ q  y/ r, ^: M) s* ---------可选项目----------------------
- v) Z# I5 s; l6 f+ E: ?*(spike yxis amean,lc(blue))(scatter yxis2 amean, lc(blue) msymbol(+) mlabel(amean))
3 C! h3 _( t& e* t$ ]- O$ c}
  ~% p8 _& H/ M' x, l: }end

maodeng

awer 发表于 2013-2-5 14:40
# p5 u8 H$ ?. {$ \$ A1 N抱歉，最近太忙了
. w, R: P6 l! u5 d! F2 ^4 r# U& P: _
300236 我用下面的程序试过，没问题

哇，谢谢，我仔细看看，这些字母数字代码晕了
# B) u: L/ d0 e

喜欢喝冰茶

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

喜欢喝冰茶

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

% I, S5 y0 V8 h0 g8 u这东西是开源的？知道有个统计软件叫stata，不过要收钱。

汉水云潇

准备试一试
先找Stata软件

汉水云潇

晕呀
论坛复制程序代码
乱码夹杂在其中
很烦
对了两次
2 Q2 b7 j+ M+ Y$ G) z以为没问题了
$ y% M" x. y3 @
# p6 }" O$ r9 L- d6 X/ }结果
执行后
还是出现
+ u/ q8 z; Z3 d5 F3 i& b- F( Sr 199
的错误

汉水云潇

试验成功，下面是三一重工的图
9 a* C4 ]$ h' k7 F- f$ k
' G- r! s& M' q

peterbev

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

中秋不圆

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

bknite

这个帖子值得学习。

jeanny

感谢！！！

Hubert2000

明月回春 发表于 2011-8-13 20:17
程序写的不错.
5 `1 G3 w- p# Q3 n# _* c# E不过有一点要注意.通达信导出的数据,只对价格复权,并没有对成交量复权.所以数据需要手工处理 ...

! u) b8 l( U6 T9 U- D# l老哥，也是山寨这篇文章找不到了，不知为何？