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

maodeng

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

maodeng

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

maodeng

. hello 300236
statement out of context
r(119);
! ~; d" ?7 t" Z3 s0 c" _
" \( Z9 x5 n- s$ F' g! `3 z+ V大师 就是这个

sergio

这贴实在！

awer

maodeng 发表于 2013-1-28 13:13
. hello 300236
statement out of context
8 K6 _9 f$ B; v  Q$ O9 ]# Xr(119);

抱歉，最近太忙了
# {8 }4 V: ?! i) }9 E) u3 a
3 o# K6 Q% B; I5 @9 _; I+ q300236 我用下面的程序试过，没问题
7 R0 i0 D) x* ~: V% N
你把 D:\stock\test\source 改成你导出的数据存放的目录，再试一试

capture program drop hello

program hello

7 c" b, F9 X& M1 q9 i, ?; T- c/ K quietly {
! @3 |/ ]9 `' Q" f7 a
) K: r9 M. t- Q! D4 \7 W  t" iclear

$ J6 L$ C* u2 g1 ^4 Y0 |cd D:\stock\test\source

insheet using `1'.txt

   keep if v5>0
$ Q6 E4 O, x6 C9 U7 Y4 e   
   gen vt=_n
% O( J4 X+ I# S2 }7 e6 Y   
   tsset vt
  [9 a3 h  E  H/ b  K4 O
% V" ~9 d4 H( B5 T4 u   gen zf=D1.v5
+ k0 i5 s+ h- V
) n8 N3 f' r$ A7 e) e% j  Z/ @   gsort -v1

) I, a) R* e( }0 c* J' }3 E; A   gen avol=v7/v5
   
# ]% w3 q$ r. x$ L9 k   keep if avol>0
   
   drop vt
  
   gen vt=_n
   
   tsset vt
( x3 I: x9 P% ^0 ^) T# {   
7 J6 m6 |+ P& Q; j  W2 L    sum v5 [aw =avol]

. W; V$ d! a# H' v   gen amean=r(mean) in 1
   
% U% H: a7 {4 D+ E8 |" l( u   * 总加权平均值 amean
3 w' f; p0 h2 Y; R' N! G2 e   
1 s4 t7 _+ d. V9 M; _9 W7 j* y   
   keep in 1/240
( @+ ]  f* }* `) S
   sum v5 [aw =avol]

) l, t+ n0 P0 ]% o   gen mean240=r(mean) in 1
- l$ U+ @* G- {5 ~( p# M   
* k6 b. O6 h6 Z9 F- d* K   *240加权平均值 mean240
   
   set obs 500
   
   kdensity v5 [aw =avol],gaussian n(500) generate(dx240 dy240) nograph
% ]8 @& s# ^+ Y3 \. z3 S& j   
   *240 kdensity(dx240 dy240)
& q; s) O: i8 N3 f5 ~$ Y   
   sum dy240
   
4 w, \, E8 \! W  O, p6 l   gen yxis=r(max) in 1
   
& K# K' N' i' y2 J+ ^- n   gen yxis2=yxis/2
( q: E& I  S, p6 B0 ?6 N% Z* R   
   *yxis
   
9 @* {# m# k' ~3 t6 ~2 K0 v   gen dy240max=r(max)
   
$ N! U0 n7 b) `   gen dy240minp=r(min)
% r% z+ S4 e. o# Q0 Q$ Y  h
    gen vtt=_n
7 E/ x1 u6 N3 r5 B: Y
; m% ]7 g" `5 ]& Z  ]; e    tsset vtt
0 u1 k$ M' e& e; }2 O
   gen ddy240=D1.dy240/D1.dx240
   
# F; y  g" T9 q# n/ P   sum ddy240
8 a" ^. o" N  |0 [* C) O   
   gen ddy240max=r(max)
   
   gen ddy240min=r(min)
     
   gen ddy240g=ddy240*(dy240max-dy240min)/(ddy240max-ddy240min)
: n/ r( l: I, f   
, z: E, `/ u4 k& U; }' k   sum ddy240g
) }0 e  P% E. ~4 |8 g) v; `" \   
: N0 x! C- [9 `4 A/ d; e   gen yxis3=r(min) in 1
& ]/ _0 }8 c7 k9 U. _$ g   
   gen yxis4=0 in 1
  
7 [* E, c( I; ]* W% N/ K0 O0 \  *ddy240g
   
, Q- M, }. j; w, d5 z, M& d: r   
  preserve
2 y. Y. @: Y" m4 j) P1 i  
  sum v5 [aw =avol]
  
  keep if dx240<r(mean)
6 N! c- v. I& x- k. E2 x  
  sum ddy240
: O0 f/ o7 D9 i( U7 C4 T  L  
% A" V7 u* I2 Z/ k" v& `# }% Y/ k. x( q  if r(min)<0 {
0 [& U7 p- U! F' Z8 w, R8 {        
                keep if ddy240<=0
3 a0 V# i5 H* N! x' |8 A2 E4 Z               
                count
6 j7 T. i& f) R. G( |/ A               
9 A1 I/ g* g8 r- D$ J, P: o- {                if r(N)>2 {
                 
1 x9 D0 z+ z  \3 l* {! \0 h8 U7 k                 pkexamine dx240 dy240
               
6 U) P9 N: j! [/ V. J. B                local tempminp=r(tomc)
                restore
0 p1 y- \2 C6 Y% L' N                 }
                 
        else {
) }$ H2 Y4 w; w* m/ n# e7 u4 L5 f2 P               
3 d9 @7 Q5 h- c# o/ T5 Z                 di "cannot find minp"
                 
' l- }% e6 S4 C2 z                restore
                 sum dx240
& P5 p0 @' s8 `8 f& }: g  l. u                 local tempminp=r(mean)
7 `" S, `0 m5 C+ ^( U9 j* H                 
8 M4 V( I9 a# y5 |/ Z2 Z                 }
                }
: ~2 g. M6 ?5 J3 c   else {
                 
                 keep if ddy240<=r(min)
                 
                 sum dx240
' G! ?5 I' M: X# g9 @                 
                 local tempminp=r(min)
                 restore
                 }
" b) O! w2 Q% v; Y9 k) y

- y/ l4 b* K! `1 T+ g
/ P# \* V, V7 L( q! p   gen minp=`tempminp' in 1
   *240 minp
1 n# M0 p! D$ ^/ ^; C2 y
! T( m  n2 w& F7 |# c
   
   preserve
   
9 p$ ?; j1 m4 O4 x; P- H   sum v5 [aw =avol]
   
   keep if dx240>r(mean)
   
/ [# T1 ]8 C% h- M) m   sum ddy240
# v. z, M: `% Y   
% E. n) ~7 R5 {9 |   if r(max)>0 {
               
                keep if ddy240>=0
4 u& j* V$ V$ i# m               
1 r6 s9 F4 u. `+ M$ m  A1 l/ g1 m! L                count
               
                if r(N)>2 {
                 
                 pkexamine dx240 dy240
  q- `0 M1 Y" g               
                local tempmaxp=r(tomc)
  r* \' w* ~& F) r% j2 u- Q! r. b- r                restore
                 }
" G* z0 R  ^& s- {                 
                 
        else {
                restore
' w  r8 N- a& T( P9 {& I4 u" p                 di "cannot find maxp"
               
  \* x  ]- A9 v) F/ C% v) O7 A) m               
                 sum dx240
                 local tempmaxp=r(mean)
2 J6 [, @# b) N7 g) h' L" ?1 l                 
  \* Q+ S" F8 x( `                 }
' c, Y* T% a, b" S$ d6 |1 ~, F: y                 
; L8 v5 V" k  t  F4 S                 }
                 
( e& H- M5 \8 a" a& ]    else {
4 b% j' P6 R9 w. v$ f: X% D& \       
9 T- S' ?; }' ^, d% {           keep if ddy240>=r(max)
/ [9 e1 X9 f4 U2 ^          
: o9 q3 d) K5 @2 z+ c           sum dx240
          
: s: a  \' Z$ }$ P; H           local tempmaxp=r(max)
           restore
           }
       
0 Z$ H  Z+ Q& u/ I" ?: }7 |
  

- X, J: x# m/ q+ O# D7 E3 L   gen maxp=`tempmaxp' in 1
   
6 c; |; c5 X3 n9 s" V6 Q) l   *240maxp
, s% |: `3 c& P; v8 S5 j  ?   
   
   *-------------ddy240pminp左侧高峰值-----
   
    preserve
   
$ Q$ M8 U( _, U( [3 |( a3 F7 V   keep if dx240<`tempminp'
" z+ {; u9 ?+ L" A+ w   
   keep if ddy240>0
   
$ Q" i$ L" v& ^   count
               
                if r(N)>2 {
. `$ i6 k8 i9 N: u; \% U+ O: e$ i                 
                pkexamine dx240 ddy240
/ J0 O* c; E! j" N* d               
1 `1 ?. [( i5 x1 c, X9 \                local tempminphddy240=r(tomc)
+ C; U6 K# {/ O2 K                restore
  y! k7 G4 O, k$ {                 }
                 
, x: N7 d( x% p0 a8 `        else {
                restore
4 ^2 M" ^: N( C                 di "cannot find pminphddy240"
                 
% o- Y: \8 P9 J               
1 {5 {% j8 V' B2 U                 sum dx240
                 local tempminphddy240=r(mean)
               
                 
                 }
6 f6 y) m8 n, o                 
0 R  g4 C) y/ g" |% U$ l9 X% w6 M   
   
   
' P: X3 i- R% n% ]   
; \; ?# B7 f) h# B* [2 \   
        gen pminphddy240=`tempminphddy240' in 1
! ^) N9 T; f; |8 G4 ~: X       
       
# ~7 W/ i# u  ^9 n# e1 `        *-------------ddy240pminp右侧低峰值 -----------------
    preserve
4 `7 {8 L2 J4 X) H- F0 x) r" t7 C
" i0 S3 h6 G8 ~   keep if dx240>`tempminp'
   
   keep if dx240<`tempmaxp'
   
   keep if ddy240<0
   
7 P8 V( ]8 R& P2 N   gen temddy240=-ddy240
# e( ^' e# l9 k& O- W7 o   
    count
. U; p1 T5 \& z; O               
& m) k& m+ J, X* a                if r(N)>2 {
! P1 r8 O) c. u" Z6 V6 F; ?; R$ o                 
                 pkexamine dx240 temddy240
   
  local tempminplddy240=r(tomc)
                restore
0 O- q4 B  `: `- `, ]                 }
                 
; j* b5 j5 N5 N9 Q* q        else {
               
                 di "cannot find pminplddy240"
2 f4 i5 L4 |0 c4 l, h" ~                restore
                 
                  sum dx240
5 u/ a! R7 D- T3 i8 h                 local tempminplddy240=r(mean)
                 
+ L" U( M# e2 ~3 ]                 
# ~6 p7 ^6 d/ o: W7 z                 }
8 R( Z7 ]* k7 f- o% F   
   
   
   
   
        gen pminplddy240=`tempminplddy240' in 1
       
. U, B1 N, b  A        *-------------ddy240pmaxp左侧高峰值 -------------------
       
         preserve
( U/ ~& c$ y$ ~8 K# @   
   keep if dx240>`tempminp'
8 U% h$ P( A+ `9 {1 K   
   keep if dx240<`tempmaxp'
   
   keep if ddy240>0
   count
               
                if r(N)>2 {
                 
                 pkexamine dx240 ddy240
/ @6 _; N$ W" P$ C+ p   
   local tempmaxphddy240=r(tomc)
7 P( b2 D2 M# o3 N1 x   
   restore
                 }
                 
  X$ F, e' N* d        else {
               
                 di "cannot find pmaxphddy240"
8 \" `2 }+ e$ [5 I, s  M                 restore
               
                 sum dx240
                 local tempmaxphddy240=r(mean)
                 
                 
                 }
   
! R- \7 H# I9 @7 T" @3 Z& S   
( O; W' O# h  x, Q  _) D   
   
   
! B1 \3 m! d- m   
6 {* b) z, d& I0 K/ B( M        gen pmaxphddy240=`tempmaxphddy240' in 1
- N" l" \( L  ^5 D9 B       
, I! P  |5 R5 Y% A8 j, c, Z       
0 U3 j6 B& y$ ^' h$ u        *-------------ddy240pmaxp右侧低峰值 -------------------
       
  _4 }4 c" Z7 f         preserve
  
* L& S( _* x2 B) q/ r+ j# M7 o  s   keep if dx240>`tempmaxp'
   
   gen temddy240=-ddy240
   
$ C/ Z1 @1 ?8 o! {   keep if temddy240>0
" T7 ~# Y1 ^4 {' Q5 r6 f   
" Q2 O4 ~, Q6 d; ]   count
4 p" f  }3 A5 ~+ S* O# |" q9 D               
) w% `5 O  W5 n1 E3 t+ k3 \                if r(N)>2 {
                 
$ }: S: g# l' Q; T2 N, _" e' q+ w                pkexamine dx240 temddy240
7 t  @! u2 e) ~7 `* z. r  R   
6 [! `6 ?" z( T$ f& y8 j4 }" m   local tempmaxplddy240=r(tomc)
! C$ T5 q6 x+ z( [2 K   restore
                 }
# v6 b- B/ m+ X  K                 
        else {
                restore
                 di "cannot find pmaxlhddy240"
                 
               
, v9 W* |1 I6 m                           sum dx240
                 local tempmaxplddy240=r(mean) in 1
( W: |* N& W5 e& Q3 T               
                 
0 N2 t& c/ E9 k                 }
   
8 d7 I/ Z- H7 K$ U  H. K   
   
1 o7 a7 u/ D6 ~& J* b   
3 O7 z" x5 C& [& U, ]- {% |% L( o0 h  o   
+ a9 C* e( n  \6 W" v/ M   
* Z' \/ i( v7 R) I        gen pmaxplddy240=`tempmaxphddy240' in 1
       
        *-------------------------------------
   
   gen price5d=v5 in 1/20
      
   gen price5yy=_n in 1/20
+ |2 @9 z6 W+ n   
   gen price5y=price5yy*(dy240max-dy240min)/20
   
   *price5y &price5d
   
   gen priceny=price5y in 1
   
- V# T. s% \0 [$ d4 G- q7 y# w0 f0 E" C   gen pricen=v5 in 1
% G2 `" ~/ p, z) X* p; ]
     
6 `1 Z" D' |) {- t+ D1 B   *priceny &pricen
# b+ A6 T% B$ C9 d8 p6 D7 S; j   
   
   *-------------------------graph------------------
0 ~9 O- S2 ~$ {* L0 F: E  
  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)

* ---------可选项目----------------------
% f- ~* c9 [+ M! f*(spike yxis amean,lc(blue))(scatter yxis2 amean, lc(blue) msymbol(+) mlabel(amean))
# r% ]/ t# m+ Q- Y$ e3 ^6 _}
4 _7 T2 m5 d1 ~! @end
0 l2 ?) y* H7 H

maodeng

awer 发表于 2013-2-5 14:40
抱歉，最近太忙了

4 c$ a: i7 Z6 @- x300236 我用下面的程序试过，没问题

; p, [, ?* z. l( O' u' @哇，谢谢，我仔细看看，这些字母数字代码晕了

喜欢喝冰茶

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

喜欢喝冰茶

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

. d7 R2 Y/ u0 x' B! R2 g6 h; H5 z这东西是开源的？知道有个统计软件叫stata，不过要收钱。

汉水云潇

准备试一试
3 O# V/ y2 g1 n先找Stata软件

汉水云潇

" K3 r7 H* w& K! b. y8 P3 G& }晕呀
1 v  V1 D7 \8 U4 ?) J论坛复制程序代码
7 k. t* e' W; \乱码夹杂在其中
很烦
5 y3 W4 j) |! d, p0 N/ _对了两次
以为没问题了
+ p$ @) |, s, r" g$ g
结果
1 ?0 r3 q# X* I  Y0 S7 V- g9 ?执行后
还是出现
: N% V$ Y4 v8 Gr 199
的错误

汉水云潇

试验成功，下面是三一重工的图
8 h* j/ C8 G& I4 a
: e; T/ R3 R; x7 `2 y

peterbev

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

中秋不圆

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

bknite

这个帖子值得学习。

jeanny

感谢！！！

Hubert2000

明月回春 发表于 2011-8-13 20:17
程序写的不错.
4 b, ^( u1 L) p: @2 w不过有一点要注意.通达信导出的数据,只对价格复权,并没有对成交量复权.所以数据需要手工处理 ...

( t( w' Z$ p' ]/ u: e! D老哥，也是山寨这篇文章找不到了，不知为何？