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

maodeng

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

maodeng

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

maodeng

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

) H9 I7 e/ L1 }0 F% D. ]: Q大师 就是这个

sergio

这贴实在！

awer

maodeng 发表于 2013-1-28 13:13
0 f( U% H8 L& m: B5 T4 ?6 K. hello 300236
statement out of context
r(119);

抱歉，最近太忙了
( z& y( H0 c8 C7 p7 S1 E7 ]
300236 我用下面的程序试过，没问题

你把 D:\stock\test\source 改成你导出的数据存放的目录，再试一试
7 j* k  k0 Y, ^! Z% t. Y4 `  t3 F& d
* j- u, C" a/ n; J$ Ycapture program drop hello
% Y6 |/ C* t7 L
program hello
. t1 B  Y* g2 s7 j0 n: G
quietly {
/ @1 \# O/ E( M4 I. l2 T
clear
: r. i6 o1 @/ e9 c/ m5 k7 h
cd D:\stock\test\source

+ k0 L4 h0 ^$ ?insheet using `1'.txt

   keep if v5>0
& i3 }8 ~% U) x8 a   
   gen vt=_n
   
/ S0 F1 _  M  e   tsset vt
: z$ \" H7 |5 |* C" `  R% _% E
   gen zf=D1.v5

" G& J# m& M9 T5 a0 O% a  K& W   gsort -v1
- h2 F# ]* F% X" z' V6 T; V* p3 Y1 q
" m* j9 R) `6 b" K, m; E   gen avol=v7/v5
   
$ i+ j8 G: n9 Q% v   keep if avol>0
8 j# ^, g  {! c: G) ]3 G$ e# S6 B   
; h) V" @" C: R8 ?' V3 ^   drop vt
  
9 H; ?8 L5 o) I( X   gen vt=_n
   
   tsset vt
   
4 p, j' O0 X" [& a& J- X: _2 a    sum v5 [aw =avol]
2 f$ D: ^- s0 V# c
   gen amean=r(mean) in 1
   
   * 总加权平均值 amean
9 A& K& C# S/ Q( d; T   
! l  |$ Z% ^1 m% z2 Y( W   
   keep in 1/240
/ d( f; _4 W) h5 O' G2 @% k
   sum v5 [aw =avol]
) j! A, S3 u" A- c- d0 @: @
   gen mean240=r(mean) in 1
8 C# x5 ?# f* O* W, Q/ O- V' j   
" F7 t  t6 y- Y+ M! r   *240加权平均值 mean240
   
3 U) F' \5 a" F; m; P   set obs 500
5 e) {5 H$ h' c8 F. g, }' D   
. M1 l, y/ B- F& E" P5 n   kdensity v5 [aw =avol],gaussian n(500) generate(dx240 dy240) nograph
   
   *240 kdensity(dx240 dy240)
   
5 x$ o$ {  f% S) X: O7 R% |! p   sum dy240
3 I$ R( ~+ ?6 ]! P/ u   
   gen yxis=r(max) in 1
9 V) n; T5 c) {  X   
   gen yxis2=yxis/2
8 s* z+ \) I  J, L1 W   
0 J2 b! N; Z; m, w7 r/ r0 G/ t   *yxis
   
9 A" Q- q$ Y/ {7 y$ q   gen dy240max=r(max)
' h2 o: U+ m/ g+ q/ i) [4 x   
1 e3 l8 R2 }* b1 X# ^2 _' m$ ~   gen dy240minp=r(min)

7 |  ^  R/ p3 g! N* n    gen vtt=_n
6 `! z# d# q( p# E5 Q
# }* x  W$ a( S6 y9 Q    tsset vtt

   gen ddy240=D1.dy240/D1.dx240
   
   sum ddy240
   
0 }# p; Q( T3 c2 E, g% H4 u" ]   gen ddy240max=r(max)
   
. D( x( R' `/ [" s8 s* J9 g   gen ddy240min=r(min)
; V( U# u' u/ M9 i1 O     
7 T! F; l# N, g4 C8 c2 X   gen ddy240g=ddy240*(dy240max-dy240min)/(ddy240max-ddy240min)
  `- J8 d% i' H   
   sum ddy240g
3 b1 y6 F+ h; I5 b' G. s  n   
5 d6 ?2 ~; ]0 p4 W% q/ I9 @) n: _/ s   gen yxis3=r(min) in 1
   
% }$ K1 z( x( G  n, }) S# h   gen yxis4=0 in 1
# r; b  [9 f* G; C0 O1 f  
' y/ ]- c/ a. E4 U" A$ l  *ddy240g
2 R9 v3 Q1 N% P1 L5 K( ?0 O1 \5 g   
  U. Q' O5 }# W! t& y2 |3 F   
  preserve
' _& R5 `# p1 G% }9 ]5 K+ K- A  
( G8 T5 w  G8 J8 Z  sum v5 [aw =avol]
  
  keep if dx240<r(mean)
# R9 o  e, z+ c+ T* V; {  
! q1 j4 k9 u- G2 z  sum ddy240
9 h, Y3 [/ T- a9 p! P$ {  
  if r(min)<0 {
1 e1 c( E6 h1 p; G        
                keep if ddy240<=0
. c: W8 d9 f% e+ q* h               
                count
               
* h6 e5 S* L0 m, W) c                if r(N)>2 {
                 
. O2 w' I( s) _$ ]* Y: k' J' R# p                 pkexamine dx240 dy240
5 k/ [/ e2 l/ i! Q0 A5 w               
0 f7 @2 f8 G3 I+ ~$ _! u                local tempminp=r(tomc)
                restore
                 }
. X& r6 R% ?5 m6 }                 
, B; u9 c7 f) t. W/ Q        else {
               
                 di "cannot find minp"
; |/ {0 l: K- r, d5 S                 
" E/ S* r  `2 i. Q9 F7 c                restore
                 sum dx240
                 local tempminp=r(mean)
                 
                 }
6 D& K3 w" |; J                }
) ^* b/ r! {: H. s   else {
                 
+ w! k8 }9 |8 I                 keep if ddy240<=r(min)
                 
! o  F; d& f4 U8 _0 q, o. C                 sum dx240
                 
0 E$ v# Q2 p( F: D1 s% b                 local tempminp=r(min)
                 restore
! c* @5 Q7 A" J" q2 j7 B$ u                 }



' q6 o/ B- l' H   gen minp=`tempminp' in 1
   *240 minp
4 g  V6 p2 t2 A! A! f
8 x2 i9 I: B9 l2 B& }, E
   
   preserve
( o* r2 w& U$ L3 M   
* h  [" f) ?$ m9 R! V' f/ e4 q/ u   sum v5 [aw =avol]
   
# ?" q+ Z5 S' H) s' [   keep if dx240>r(mean)
   
+ L+ I  t4 q8 h6 ]   sum ddy240
   
0 Z  }, T& w, s3 j' E   if r(max)>0 {
7 R+ O# _5 n8 Y               
( c; H% t4 A9 `# q: S& E                keep if ddy240>=0
               
                count
               
; h# l/ \7 v; R7 |# U" s                if r(N)>2 {
                 
9 i/ z5 g, g/ P                 pkexamine dx240 dy240
3 \  F6 J' H. x; h+ L               
                local tempmaxp=r(tomc)
/ Q) A6 K7 J, f6 J  B                restore
/ v0 ^1 o2 H4 t1 w* S" x! D3 w. C                 }
                 
" Z9 c, U: X" B                 
! O, s1 b% k$ r' L# o. c        else {
& L3 X3 B2 t; o3 @                restore
                 di "cannot find maxp"
               
               
                 sum dx240
                 local tempmaxp=r(mean)
8 q5 C3 y9 x- \" g) s1 _                 
                 }
                 
                 }
6 `2 `2 c" T/ t9 \3 A1 d8 y% e                 
    else {
       
+ J7 o* k9 j( U0 N; o0 }           keep if ddy240>=r(max)
          
! k7 V# X0 I, [5 Z" L           sum dx240
! u) @! n# F; c1 h5 U4 L          
- P& z6 c3 T5 q, Y1 z' \           local tempmaxp=r(max)
           restore
           }
# ^6 s% }& h  d6 v: v" v       
/ Y2 M4 G7 ^; b0 P; K! g7 G
  

   gen maxp=`tempmaxp' in 1
   
, u- a" m. `( B9 t/ t   *240maxp
   
   
' g& B0 y: H" m0 y2 h) A, w! D5 g% {   *-------------ddy240pminp左侧高峰值-----
, d: @5 V4 u9 M! Q2 E% J   
+ W, f- R7 h4 o3 W- [& h  X    preserve
) _' F& ^) K/ g. s5 U1 A  E3 |   
   keep if dx240<`tempminp'
   
   keep if ddy240>0
   
3 r" R& ~9 h$ P   count
               
8 A4 c1 F; A! L! K                if r(N)>2 {
                 
  }( h5 q( b" E+ [) r6 a                pkexamine dx240 ddy240
               
                local tempminphddy240=r(tomc)
                restore
/ f4 x" x: m9 A6 H                 }
                 
0 s( T' _1 ]/ y. C1 V        else {
0 W0 T) y8 Q5 |' x1 g# w' w                restore
9 v7 X3 N8 @9 @0 A) ^3 ~2 u                 di "cannot find pminphddy240"
8 U# i! O% s4 C  H- O. T3 |                 
8 x1 ~6 v! w- l+ }7 w5 e* Q' N, r               
/ K0 b' Y1 Q! g/ p                 sum dx240
                 local tempminphddy240=r(mean)
3 a6 U- G* ]& }+ W4 X$ ]; O7 J* T               
                 
. j7 N$ i3 Z$ w                 }
                 
   
" T$ F8 ]8 u# r6 {4 R   
* E4 H+ H# a# m. D, \   
. J' o& {. I' P- U+ ]   
/ v: p. H$ I* T: ^7 R   
        gen pminphddy240=`tempminphddy240' in 1
       
       
        *-------------ddy240pminp右侧低峰值 -----------------
    preserve
) `* h. w6 _( {  }2 A4 S9 k2 S
   keep if dx240>`tempminp'
   
   keep if dx240<`tempmaxp'
  E9 @0 @4 t% s$ D: P; o# o   
   keep if ddy240<0
" V7 c! G8 Z: C) y- e7 @) h   
   gen temddy240=-ddy240
   
, F$ N1 T3 G, S2 o  ]9 r    count
               
                if r(N)>2 {
  B) X2 L3 g5 }" j8 x& N                 
+ `8 W/ x8 k0 K- }0 R: t                 pkexamine dx240 temddy240
   
  local tempminplddy240=r(tomc)
# S+ f' C" j3 [+ [% k& |- ~3 X' o                restore
& N  x+ g% J1 E                 }
4 k" K) L0 x5 Q8 h  l                 
% }8 N, E4 S* D# S9 ~        else {
1 Q! S" K, s( g) l5 z  t               
( d1 J- F  d# \                 di "cannot find pminplddy240"
7 C1 B  j) E2 |$ G/ U2 c                restore
                 
, o8 U: l  a: Y/ l5 D                  sum dx240
                 local tempminplddy240=r(mean)
                 
9 N( s( w6 @  v' H                 
6 n6 T/ T, u( v0 C- N! \                 }
   
# T% ^5 J7 @) I; O+ ^   
# S+ V% k: W: {/ G. ?. S& q# t   
- [+ B( F& E  V7 l   
& c0 @; j! D+ Z) d& n9 V   
! ~, l2 |$ A. S" g% C        gen pminplddy240=`tempminplddy240' in 1
       
        *-------------ddy240pmaxp左侧高峰值 -------------------
       
2 |$ F( u9 z( `, O; ]7 h4 {7 y         preserve
   
   keep if dx240>`tempminp'
# l' p4 M1 n; ?   
   keep if dx240<`tempmaxp'
! x( B( E0 x( K" n! G   
0 u5 W' \' R9 \, I8 Y   keep if ddy240>0
' e( o' O, K7 e" l   count
9 b4 x& ^% S$ p               
7 `) _& R1 q& j5 F  H; I6 n                if r(N)>2 {
                 
                 pkexamine dx240 ddy240
   
   local tempmaxphddy240=r(tomc)
5 y  l- u! a- y  @2 h5 i   
" F& R- j. ], C# w   restore
                 }
                 
( C, n2 u/ n1 _% t        else {
4 A; e3 j' @( _+ ?0 E. I! G2 W               
: K; W: B# [, z' j0 h+ l3 v/ X+ g                 di "cannot find pmaxphddy240"
                 restore
( A& Q# c" Z0 l: X7 I4 K" d" P/ _) i               
                 sum dx240
                 local tempmaxphddy240=r(mean)
  B4 J* [$ a0 e* F3 z$ ~                 
                 
5 o/ p- _+ C# z. h/ ?                 }
   
) x% q- s6 Z9 q, K: |5 E. p   
   
   
! @' f" T5 D2 H: R   
   
. p6 @9 W$ O* g        gen pmaxphddy240=`tempmaxphddy240' in 1
2 C# ~" ?- w+ r- }4 Z       
       
        *-------------ddy240pmaxp右侧低峰值 -------------------
       
         preserve
7 ?+ ?" d" w2 S5 P% H! e, d  
6 ], i- E/ }2 u2 H5 a   keep if dx240>`tempmaxp'
9 n9 W4 f- E6 S9 O5 a   
/ U( I7 l/ {% Z5 w3 _' K6 [$ t   gen temddy240=-ddy240
; c! r0 v# o& I3 I& c2 C   
   keep if temddy240>0
+ [9 q! ~+ E9 p( ~   
! M2 j) ^% {1 x: d8 N* P" {& n2 n0 U   count
) F; j) E" \& t6 y9 l* A  |               
                if r(N)>2 {
; w- Q0 A/ v8 |: a3 j6 M9 g2 L                 
                pkexamine dx240 temddy240
   
+ F; b+ l6 {. @- E   local tempmaxplddy240=r(tomc)
   restore
                 }
% L, V& B( s2 c                 
        else {
! }: e- i8 l3 u, u) L                restore
                 di "cannot find pmaxlhddy240"
5 ^: Z5 p$ h  w! z3 T- p                 
/ g1 Z& p' p# f2 m               
                           sum dx240
                 local tempmaxplddy240=r(mean) in 1
5 s8 E# h' p2 j9 u" z) n               
' j* T, k$ u$ {$ n; W" F                 
) a5 r! ~, O: S: u  w# S8 ?                 }
2 C4 ?  p* f! ^' J- k   
7 I7 v! k8 n- r7 J/ F+ {   
7 E  ^% Y; a8 H# y! d) d" x) c) d9 I   
   
   
4 E5 _  H: N7 L   
        gen pmaxplddy240=`tempmaxphddy240' in 1
       
        *-------------------------------------
   
" \( T6 m/ D7 l' K% K6 B  ?   gen price5d=v5 in 1/20
      
+ C, [: c1 {; F* }! I+ Q- Y   gen price5yy=_n in 1/20
   
" L" V/ b! H# T8 L' s   gen price5y=price5yy*(dy240max-dy240min)/20
# F$ D( h1 Z' M   
" J5 e; Y. |# u/ C& @$ n' ^   *price5y &price5d
6 W6 z2 j" N) F  R   
   gen priceny=price5y in 1
: I- o3 A& B1 j   
( |( [" S( a1 X- J% ~4 |3 R: f: r6 k   gen pricen=v5 in 1

& z  B/ X: R* q% y4 n     
   *priceny &pricen
/ q2 P/ W; R+ K2 Y  S$ h' t8 \* ]* V/ u   
1 ?- U. P% y# l) R' l4 g   
6 V) K# }8 M( n; `   *-------------------------graph------------------
  
  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)
6 x. L5 @9 \. s- T: e
1 k  F$ T6 P! [( Y* ---------可选项目----------------------
: B2 V& Y# d, {" F) @3 l3 [*(spike yxis amean,lc(blue))(scatter yxis2 amean, lc(blue) msymbol(+) mlabel(amean))
* x! l% d  n1 T$ x: J/ ~; T/ K8 N}
end

maodeng

awer 发表于 2013-2-5 14:40
) P- i9 [7 B$ o& W- A抱歉，最近太忙了
9 M) E2 A( ?% A  W
300236 我用下面的程序试过，没问题

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

喜欢喝冰茶

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

喜欢喝冰茶

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

; e6 [/ y1 m& |3 W- e. m) R' g# X这东西是开源的？知道有个统计软件叫stata，不过要收钱。

汉水云潇

准备试一试
先找Stata软件

汉水云潇

0 j1 w- e( S7 r5 [$ \0 O晕呀
论坛复制程序代码
- i$ k  \. N9 {9 U# K乱码夹杂在其中
很烦
* I6 p  I! E0 c, `# n/ Y对了两次
( F9 @5 f0 [, V+ {- q5 p以为没问题了
' B1 c4 s. [2 d0 X& k% n
结果
2 }7 f$ G) u7 F( Z1 a2 X执行后
还是出现
r 199
- W; c, ]' o9 I3 K的错误

汉水云潇

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

peterbev

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

中秋不圆

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

bknite

这个帖子值得学习。

jeanny

感谢！！！

Hubert2000

明月回春 发表于 2011-8-13 20:17
程序写的不错.
1 l' y) K: C/ V1 }& w  o2 K2 D不过有一点要注意.通达信导出的数据,只对价格复权,并没有对成交量复权.所以数据需要手工处理 ...

8 ]& }3 g( y4 F. C* L  V" M. V$ G- }老哥，也是山寨这篇文章找不到了，不知为何？