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

maodeng

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

maodeng

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

maodeng

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

大师 就是这个

sergio

这贴实在！

awer

maodeng 发表于 2013-1-28 13:13
. hello 300236
4 m6 B* W1 ]# J8 {9 m9 f( A. T. h$ Qstatement out of context
* A$ {- U1 `' ?$ q6 D. V' br(119);

抱歉，最近太忙了

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

: y5 t4 u0 h# c2 E( ]. {  S0 }: J你把 D:\stock\test\source 改成你导出的数据存放的目录，再试一试
3 J+ t! c& A5 }# ^/ d3 Y( y
/ E; [. l8 r% J3 y* y/ Wcapture program drop hello

program hello

/ y  Z3 |6 t/ W2 d9 U/ p6 C quietly {

+ [6 f# a! N' Eclear

cd D:\stock\test\source

) a$ J$ Z- P- k# f& \insheet using `1'.txt
" T% |) i2 u7 Y% B4 b
   keep if v5>0
  j; Y! Y6 u% c; z   
$ X/ u; A# E# q7 K0 w/ H" S6 K   gen vt=_n
# e$ x8 f1 B" Z5 P   
   tsset vt

3 u% v% m2 _* f1 N   gen zf=D1.v5

   gsort -v1
; I9 Q0 [6 j) Q1 R
4 b$ Q! t0 n& k6 P) h1 v   gen avol=v7/v5
8 B8 I9 ^2 h. v3 b" n   
   keep if avol>0
( m1 |; o8 ^  x5 W8 n7 e   
   drop vt
  
   gen vt=_n
   
' G1 w: e! k: Q* t8 R1 f   tsset vt
5 G8 p  _- V, `8 c# D   
    sum v5 [aw =avol]
! ?' u" H: k" ~, Q& x1 y& l! |
   gen amean=r(mean) in 1
   
, g- @5 i) M% x$ z" V. p   * 总加权平均值 amean
# m6 y5 F4 K* ~$ g. U2 t5 P- k/ [  `   
" V# O' l5 |, d& p1 }   
   keep in 1/240

" s1 C& H* y( ]0 k0 l. B/ o: I# f, o   sum v5 [aw =avol]

4 ?2 C3 p2 d2 D* p" d7 @3 q7 o   gen mean240=r(mean) in 1
% l: M; m3 U+ b: X: N+ B9 q   
   *240加权平均值 mean240
   
4 p# t. x/ n. \" r+ M  E9 G   set obs 500
- m/ f: O" B! Y   
/ l# K0 @* `7 e0 e5 J4 G0 @   kdensity v5 [aw =avol],gaussian n(500) generate(dx240 dy240) nograph
! T2 T8 d; G* J; z1 V   
8 p8 {! {1 b: j' a2 x7 I1 B   *240 kdensity(dx240 dy240)
   
   sum dy240
   
( D: V7 |- w' m5 h/ \+ [$ t& G1 U   gen yxis=r(max) in 1
* v$ i' W! l9 u( [2 w6 g# Q   
   gen yxis2=yxis/2
, `$ p3 Z. c! D) t   
2 }0 |1 G+ K- i: i. v, D9 O2 R   *yxis
2 m" Z4 l: W" w- Q   
   gen dy240max=r(max)
7 i+ V9 V0 x1 V   
7 w$ b9 a) V0 D( h   gen dy240minp=r(min)

( l! l. O8 z  c    gen vtt=_n

1 W" ~: a$ ?9 v- T/ k    tsset vtt
1 P* n, P: ]  W2 c
   gen ddy240=D1.dy240/D1.dx240
   
% ^* ~& ?8 n( D: X0 i3 C: \   sum ddy240
   
   gen ddy240max=r(max)
   
   gen ddy240min=r(min)
     
( x/ j! n; L* T  p( P% z  D( B   gen ddy240g=ddy240*(dy240max-dy240min)/(ddy240max-ddy240min)
   
/ t/ e* W1 D; C0 i! R; {   sum ddy240g
' K6 a+ ?$ k+ w/ t* W4 T   
; u  j# Q! E, @   gen yxis3=r(min) in 1
   
4 D/ A2 `% {% B1 x+ s  F   gen yxis4=0 in 1
+ F" q, [4 X3 d# r! C  
  f, m, N1 m- s  c- ~+ A  *ddy240g
   
   
- t' T: V& K  |9 `( w! O1 w  preserve
  
  sum v5 [aw =avol]
  
  keep if dx240<r(mean)
  
  sum ddy240
  r" a9 L4 }8 `# v( P3 w" s9 A+ L  
. c& [! `- W" N% N  if r(min)<0 {
        
                keep if ddy240<=0
- m" j0 R, W' u0 f. q               
' e( x) K, I- l. p                count
               
0 U8 Z9 y: s1 w' J+ x: W                if r(N)>2 {
                 
, T6 ]  G6 f  e! h6 l) o6 @                 pkexamine dx240 dy240
               
                local tempminp=r(tomc)
                restore
                 }
9 a; ]0 ^/ j0 T+ A+ d" l& x+ j% X                 
        else {
               
                 di "cannot find minp"
' [' g+ A" c# E7 e" M                 
) A' i1 M2 r  f8 }                restore
                 sum dx240
' ?0 R  v# }# U                 local tempminp=r(mean)
$ M8 k3 ^0 P, B; f6 N                 
                 }
$ W8 N! O* H4 S0 n/ J. ~                }
   else {
) h4 O( o4 e; B% g                 
                 keep if ddy240<=r(min)
% u/ f% h  B% H! B% l& `( w                 
                 sum dx240
$ P5 Q0 e1 x% W  t! G                 
                 local tempminp=r(min)
8 c2 |2 S7 U9 x$ [( B                 restore
                 }



   gen minp=`tempminp' in 1
/ a4 \* r- V. y. \   *240 minp


   
7 Q! |- b) v6 G0 g% R& m   preserve
   
1 D( M- N* A) Y   sum v5 [aw =avol]
   
  _0 P! {4 R, Q   keep if dx240>r(mean)
' Y+ ^9 P6 }  O3 k0 Y   
0 P6 M( L( [  H+ {" ?0 L7 P$ c5 S   sum ddy240
   
   if r(max)>0 {
               
/ T9 Q! n9 e3 O4 P( m                keep if ddy240>=0
               
                count
" \% R7 M. C4 q8 u: J3 Q/ i1 r               
                if r(N)>2 {
. b9 c+ P9 @- d, x" I9 b: m                 
) a0 q3 O' ?/ W+ w1 F4 ]                 pkexamine dx240 dy240
" p$ f3 @/ N7 t; v               
                local tempmaxp=r(tomc)
                restore
! U* e! _) @7 Q$ Z) l                 }
1 J3 J8 i! @. D) O                 
                 
        else {
                restore
3 d+ U. ]4 o' V' U7 r4 g* i# b, ?. f                 di "cannot find maxp"
               
7 N8 e3 `) ]- `6 L3 T9 K: @               
                 sum dx240
                 local tempmaxp=r(mean)
: y" a3 W) E! ~$ C' M+ Q8 P1 f                 
                 }
                 
                 }
* [0 `4 p2 y5 Y                 
* C& n! g& c& E    else {
1 ~1 |9 Y" H) {+ n. O6 s       
" m) y  ~! Z9 ^, B           keep if ddy240>=r(max)
8 n2 j3 @, K4 ]; g! @4 m0 a          
           sum dx240
8 F/ {" H/ A( H$ {$ T7 j& Y          
           local tempmaxp=r(max)
           restore
           }
       
, f% I$ a' q6 p2 h% s
9 s+ `, f. |7 V, V5 A! g  
. U4 |$ A8 u& @% L
   gen maxp=`tempmaxp' in 1
   
   *240maxp
   
   
   *-------------ddy240pminp左侧高峰值-----
" @' B4 K8 M  [7 h  F& }   
    preserve
- e  C6 D. O6 x- P   
   keep if dx240<`tempminp'
! F) Z; M8 c5 J4 w, d+ T   
   keep if ddy240>0
9 s- `" s' Z; ?+ y, [/ X2 m   
9 x- p& R, }' N- ^1 u) `   count
               
                if r(N)>2 {
; s# E3 y) A7 \: J1 l8 H                 
* Y4 _3 o) j* G                pkexamine dx240 ddy240
; B8 G9 N6 N( }               
                local tempminphddy240=r(tomc)
& H4 u1 a$ r& ^5 i3 m! O                restore
; e+ W5 y( z* k, a' r. {" S                 }
; V. l8 k4 u" Z- m* |                 
        else {
" J. v, s, K+ g# ^% ~8 j1 g: e! M                restore
                 di "cannot find pminphddy240"
# X1 ?/ R. c5 }8 {8 T( m                 
               
0 `* S9 @  C: Z3 v3 m                 sum dx240
                 local tempminphddy240=r(mean)
6 [. n7 n* x" s$ O  M- z+ v6 e               
" f) N7 j+ t% e                 
                 }
& c$ H9 G  t. }3 d$ C8 L9 z$ z                 
8 C6 U! A) x9 h( E" O, Z" l7 H   
) ]3 T2 R( n: y+ E1 E9 a& e   
; |. {+ X. j8 C, C   
. Z6 G; t% j" p6 ^   
   
        gen pminphddy240=`tempminphddy240' in 1
+ q  p% v# Q7 P1 Z       
       
        *-------------ddy240pminp右侧低峰值 -----------------
' H% f* [$ V) l& Z0 g% J) f: s    preserve
5 ?9 f: ~; n1 x" L' p$ h  t
   keep if dx240>`tempminp'
9 `+ ]; b  V4 D1 c, r   
   keep if dx240<`tempmaxp'
   
   keep if ddy240<0
' K  h" Y* B% L  A) s! E: `- Z   
   gen temddy240=-ddy240
8 s: M. b" b* M# t   
    count
+ ?. \; a8 F* H# E! p, F               
                if r(N)>2 {
5 x3 l; S1 E" a+ m                 
( Z5 X# C. K( X4 ?% g, J+ B                 pkexamine dx240 temddy240
  o: C" ~0 L) y" P2 y   
: [( J, |  \" ?( R% n  local tempminplddy240=r(tomc)
                restore
+ v8 a# Z' I- f, w/ S                 }
; }0 f( U' e! I                 
        else {
+ Z5 s/ j& E7 G8 J$ K+ `& a( f6 I               
                 di "cannot find pminplddy240"
) {5 g4 W- \5 ?                restore
                 
                  sum dx240
2 q1 N7 H2 \: A7 v7 `2 S" p8 @1 [  T                 local tempminplddy240=r(mean)
                 
                 
                 }
   
   
   
% A) O7 M# \* T9 Y7 L! W$ ?/ k   
8 r3 E/ A7 L2 G, B8 a# z6 k   
        gen pminplddy240=`tempminplddy240' in 1
       
        *-------------ddy240pmaxp左侧高峰值 -------------------
3 _, v' ^5 O+ @9 w       
% M8 b+ o+ I8 W6 r4 v" X         preserve
# I( j  X. h. {: T7 C4 @! e' Y   
   keep if dx240>`tempminp'
   
7 G$ a$ q8 f! j( q) ]' W   keep if dx240<`tempmaxp'
3 e; T! L* r. e. x   
7 }; |) e" o0 C% V& X' b   keep if ddy240>0
   count
               
                if r(N)>2 {
                 
                 pkexamine dx240 ddy240
   
8 A7 E1 I; i: Q; S9 O" k   local tempmaxphddy240=r(tomc)
+ S) S: w5 P+ f4 k+ ]   
   restore
                 }
                 
        else {
6 ~. f/ _/ ^, ?               
                 di "cannot find pmaxphddy240"
                 restore
               
# p  a& w$ i) t  j                 sum dx240
                 local tempmaxphddy240=r(mean)
2 k8 m5 I% K8 W: h' o0 P) P                 
                 
) y6 Q' z- z% }/ o% S                 }
* }! {# z, O0 k% D   
   
/ A$ U4 b" ^3 b2 R. K3 h" z7 N- G   
& Y( H# R: u) A3 b+ j0 q1 G6 E# g7 I   
   
   
        gen pmaxphddy240=`tempmaxphddy240' in 1
- x: |) n7 [/ O  T       
       
        *-------------ddy240pmaxp右侧低峰值 -------------------
       
         preserve
  
   keep if dx240>`tempmaxp'
6 @$ H) X, |$ N6 [3 k# N1 n   
   gen temddy240=-ddy240
   
   keep if temddy240>0
7 ]9 V. K, w2 [# T+ U- e8 f6 Y   
   count
7 P9 Q* @$ K% ]: s! u+ n$ n               
& H1 M/ H# Z& }* h1 Z                if r(N)>2 {
                 
                pkexamine dx240 temddy240
   
9 o% m. i+ d( K$ e3 C7 l( C8 y0 R   local tempmaxplddy240=r(tomc)
   restore
7 `4 y  r  C- S' m                 }
9 h( t8 f! ^. g7 X# t                 
        else {
                restore
                 di "cannot find pmaxlhddy240"
                 
! `$ [" ~3 m7 I               
1 g5 p; V1 _/ z- m; O1 ]5 N. \                           sum dx240
                 local tempmaxplddy240=r(mean) in 1
9 w( M/ q1 ^' y! y  i& Z               
0 C+ m% X6 Q: m5 {" l" q. Z9 _                 
2 ?& z* l9 `, p7 w% ~                 }
+ }6 o1 d$ m: J% {. e   
   
   
1 v& L3 Q" v0 U   
- B! `- Q' y# V2 d$ t   
   
        gen pmaxplddy240=`tempmaxphddy240' in 1
       
+ }7 x8 {; Z! W3 ?( G& v* x        *-------------------------------------
: r+ ]8 k) k7 n: q- g( C- e   
   gen price5d=v5 in 1/20
5 h# h! U" L2 q- \6 [, Q      
   gen price5yy=_n in 1/20
5 s( d8 A7 V( I  }   
   gen price5y=price5yy*(dy240max-dy240min)/20
$ s5 C6 R+ e; x9 z( C   
$ S& l. T9 F# U1 s5 P5 M   *price5y &price5d
   
   gen priceny=price5y in 1
% v  ^/ j7 m+ T. O2 c& @% t* J   
   gen pricen=v5 in 1
8 B* ~) Y% p0 [1 V6 x
: g; w; v2 p% y' r" q, S     
6 k3 A0 v3 y6 I- D$ H2 z: D& b   *priceny &pricen
   
   
   *-------------------------graph------------------
  
6 e; r3 B. ?: ]; ]& n  D  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)
, ]1 K  }- E6 V  A# Q* Y
* ---------可选项目----------------------
*(spike yxis amean,lc(blue))(scatter yxis2 amean, lc(blue) msymbol(+) mlabel(amean))
1 c; l/ a- ?/ m* q# y  u}
end
: i& a" j4 G" r! B! B7 ]8 o

maodeng

awer 发表于 2013-2-5 14:40
! L) J6 j/ h/ D( G: `2 {9 |抱歉，最近太忙了

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

( h% g% B  z: @6 f  k哇，谢谢，我仔细看看，这些字母数字代码晕了

喜欢喝冰茶

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

喜欢喝冰茶

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

3 Z1 p; ~% e1 G
4 ~' e- _' g  w. {这东西是开源的？知道有个统计软件叫stata，不过要收钱。

汉水云潇

准备试一试
" x) M0 V2 h: a0 u& _7 x先找Stata软件

汉水云潇

: p+ h7 p; N# }. x7 [6 q! b: v
# Y* o# E# K" q* W6 J7 K晕呀
论坛复制程序代码
( Q2 p' [0 ]7 s9 F; d8 S' U9 L乱码夹杂在其中
/ R; Q0 {# C0 w8 M8 q很烦
对了两次
以为没问题了
3 [( _- B2 w8 f
) O) x* N& v! U9 J  t, d结果
执行后
还是出现
r 199
2 B% k: ~! z; c8 ^$ L' Y$ y4 a7 {8 C" h的错误

汉水云潇

试验成功，下面是三一重工的图
+ t$ X# x2 G( M9 D: t

peterbev

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

中秋不圆

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

bknite

这个帖子值得学习。

jeanny

感谢！！！

Hubert2000

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

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