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

maodeng

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

maodeng

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

maodeng

. hello 300236
3 X5 u: m/ M4 \$ G" istatement out of context
. F- V6 q1 I& p6 T1 Vr(119);

' ^9 A9 e  e2 l9 i  p4 U大师 就是这个

sergio

这贴实在！

awer

maodeng 发表于 2013-1-28 13:13
' e, G) H5 c0 t$ Z: @% f  B. hello 300236
statement out of context
r(119);

* Q' D/ |. Q9 Q' u2 |4 A& C8 @( u抱歉，最近太忙了

300236 我用下面的程序试过，没问题
. u1 h/ y: S/ v8 g
, |( h. _2 j( H, p1 O你把 D:\stock\test\source 改成你导出的数据存放的目录，再试一试
1 L$ h/ {4 z& w, I9 _* d
" q& Y3 P- u5 f) Q2 N( \& T! _capture program drop hello

8 p$ _2 V" `, K/ \9 z1 H4 uprogram hello

quietly {
/ M) ?& A* C# U: J# T- V9 \$ [/ ~
# T1 s! O2 w0 W+ u- |/ Cclear

cd D:\stock\test\source
- |% E, r/ P3 q
8 X) [' v2 W- ainsheet using `1'.txt

   keep if v5>0
   
   gen vt=_n
7 Z& J" P$ F" X0 Q4 r7 W2 @   
% {) z3 I7 G2 o! b! \! V* ]' k   tsset vt
9 q. Z  ^2 ?7 Z& O2 _& b
! H; p+ A; D0 X* J   gen zf=D1.v5

3 \. H9 O- o6 A" I+ @8 s   gsort -v1
1 I+ U! l" y3 T  y
2 {& T- o% ]  o- e   gen avol=v7/v5
   
; N" ]. R* Y0 q   keep if avol>0
   
   drop vt
  
$ D8 e' l7 i8 g% R( o3 x9 x   gen vt=_n
) J7 X9 G, M3 K( s  }) D  \   
2 C3 a& I7 T- |$ |$ N   tsset vt
   
    sum v5 [aw =avol]
7 m* `9 A2 s# l. j5 A- N
   gen amean=r(mean) in 1
   
+ e2 b0 e' J! B0 F  Y- ]% q: U   * 总加权平均值 amean
' Y, J: F* ?( p8 E   
   
1 a0 A8 n8 J- N/ u9 Z   keep in 1/240
  c" `8 d7 A7 q1 W: W) g
   sum v5 [aw =avol]
+ J' a, g- a6 u# v$ o
   gen mean240=r(mean) in 1
   
   *240加权平均值 mean240
   
7 \9 a4 U$ X5 |7 a   set obs 500
   
# Z1 {9 x" {5 ^" g   kdensity v5 [aw =avol],gaussian n(500) generate(dx240 dy240) nograph
% ^; W" M3 c* r) p: C   
   *240 kdensity(dx240 dy240)
1 i5 y2 m6 D0 |) ]$ t6 D' C   
   sum dy240
   
1 S( |7 |, S3 ^   gen yxis=r(max) in 1
   
   gen yxis2=yxis/2
  i* T3 ]6 n5 S1 \   
% Y* j9 W3 W% j$ `, I   *yxis
   
   gen dy240max=r(max)
   
   gen dy240minp=r(min)
" }7 x* {$ K2 @8 T
% W+ F, H3 ?% A- k8 ^- f+ M    gen vtt=_n

    tsset vtt

   gen ddy240=D1.dy240/D1.dx240
5 j7 T# R0 I$ I6 D3 R! t4 r9 N  B   
8 h# P, ~( X1 d8 E) V   sum ddy240
   
   gen ddy240max=r(max)
   
   gen ddy240min=r(min)
8 R* l/ |, o/ z     
7 T# C: q9 j& B( `   gen ddy240g=ddy240*(dy240max-dy240min)/(ddy240max-ddy240min)
8 J6 y, ?% M, Z& T2 t- q   
   sum ddy240g
$ c$ V3 R0 I8 F/ l) [/ m7 p& [   
   gen yxis3=r(min) in 1
   
. A2 O6 y5 }; Z- x9 Z" q2 g   gen yxis4=0 in 1
- }/ Z+ m2 T+ Q9 n  
  *ddy240g
# c" Q6 u$ D1 i9 I   
   
  preserve
  
0 p* v9 T- ~# P  sum v5 [aw =avol]
  
6 V& g! e: _( s8 Y/ {% l# w# C" `3 m  keep if dx240<r(mean)
  
  sum ddy240
  
  r+ x6 @8 q7 `+ _  if r(min)<0 {
        
% n) e5 L0 H0 ^6 x" n                keep if ddy240<=0
2 [' c/ J! o! H3 L! y4 j8 X               
9 r5 d# z+ C- Y! e0 `" r                count
+ d4 M; Z& e1 Q, }/ ?& g               
; q- v. y. K6 h( y  p                if r(N)>2 {
                 
                 pkexamine dx240 dy240
               
9 H# K% R2 L6 g" C2 l                local tempminp=r(tomc)
                restore
                 }
                 
$ v5 V; m4 b0 M        else {
7 X% `0 A) g) w) @* {               
                 di "cannot find minp"
                 
                restore
' K9 p/ e% M6 A$ N$ M, P                 sum dx240
                 local tempminp=r(mean)
: T  A0 U: p+ S$ ~! D                 
                 }
: Y* y6 C5 D) U9 H3 `! H* c                }
3 A4 A, `2 ^3 m1 p   else {
& c+ L( d/ U% H4 \, w                 
                 keep if ddy240<=r(min)
, w1 L7 {- p3 l# M/ p/ E                 
2 ]8 b3 X8 Z$ x5 ?% K/ E                 sum dx240
                 
                 local tempminp=r(min)
; a0 i$ x% f. q- X2 j) }8 R  y                 restore
5 ?# g  s+ z' k& R                 }
, Q# {& g" Q; n+ r! p- r- C
$ k7 X/ \5 g7 c3 J/ i

  C/ i- ~: ~6 a5 o' \7 B0 [   gen minp=`tempminp' in 1
   *240 minp


0 i* M4 e' c* R9 u9 U   
   preserve
   
   sum v5 [aw =avol]
   
   keep if dx240>r(mean)
8 M) |7 j4 J& s" p6 q4 s$ t5 U   
   sum ddy240
   
7 t8 H: h# c) ?7 Q( L) D   if r(max)>0 {
               
$ w+ \- W' g  ?2 b8 V. ]' l                keep if ddy240>=0
               
# F+ V( U$ q6 Y1 B; f9 `7 C                count
, K% i4 _; f( Q+ Z, c- U& I               
, n& U; q; e1 k: C  m0 v1 M                if r(N)>2 {
7 k+ X) J# Q: S6 Q9 {                 
                 pkexamine dx240 dy240
! E; f, Y8 _" ~! g) E               
                local tempmaxp=r(tomc)
" @0 E: S+ z% c. S- Z, n0 S                restore
% Y6 o4 _( I( o                 }
                 
                 
        else {
3 G9 ~% w1 b- O; q                restore
' D& k( s5 _' O( K4 H# b( p, N                 di "cannot find maxp"
               
4 e* D$ e, V' G. S# ]0 @               
$ c; V% e5 D6 e& R: {' c                 sum dx240
. U' G' O) a) j+ i2 r! z                 local tempmaxp=r(mean)
                 
                 }
                 
6 l) Q, d/ c8 i! u+ C7 x                 }
$ g+ g! ^. \% y# t                 
    else {
       
0 ~$ l' s, j4 c& }1 u5 Y           keep if ddy240>=r(max)
4 P2 q0 s; I/ \0 r6 L' f          
           sum dx240
  y; l( K4 e4 R1 b' d          
6 Q9 y+ ~$ M# j5 e, B5 S+ P           local tempmaxp=r(max)
% Y, q! m; k$ d. n           restore
           }
% T1 i+ X; R& u8 p       

* L1 H, d$ A  x) q9 b  

   gen maxp=`tempmaxp' in 1
& u8 S7 C9 O' v+ z$ M   
   *240maxp
7 W0 u$ P3 B9 }2 h& ?  q   
   
1 M- k( Y/ k' T  F8 `; p- l8 ^   *-------------ddy240pminp左侧高峰值-----
   
% U/ m/ ?) e4 N    preserve
   
3 R8 ]( r* u2 D: Z" s0 y. B   keep if dx240<`tempminp'
   
/ ^2 ]' N0 m9 A/ e7 `   keep if ddy240>0
/ p6 O4 ^, G; I! m' o   
, m) W/ B, d, ~! E& z; J   count
               
                if r(N)>2 {
                 
                pkexamine dx240 ddy240
               
                local tempminphddy240=r(tomc)
                restore
                 }
" o  R9 K$ D1 K/ r" C/ N. a2 {                 
  c1 P% ^* w. O! |& N        else {
                restore
                 di "cannot find pminphddy240"
                 
               
                 sum dx240
                 local tempminphddy240=r(mean)
               
                 
                 }
, S  u- F8 f- d6 h                 
   
0 N3 k3 }+ L7 h& O1 a  Q" E   
   
# j$ j8 K- k( a/ S# c3 n/ n# L   
4 x6 F- V6 e' N! D- r$ v   
        gen pminphddy240=`tempminphddy240' in 1
: }' c0 h8 p; a: X' N1 H! Y       
       
0 y% I; g$ c$ H0 q+ J        *-------------ddy240pminp右侧低峰值 -----------------
    preserve

   keep if dx240>`tempminp'
   
   keep if dx240<`tempmaxp'
   
& D' h8 Q. [- C+ y8 I% X4 ^   keep if ddy240<0
   
; @$ a/ B4 {) Q6 H   gen temddy240=-ddy240
1 `( @8 z% R3 t$ ?   
" t( `1 `, e( _! g. C1 _7 L    count
( ^  X7 F7 u& Z( ?$ p               
- S- A' Y8 f& q5 P                if r(N)>2 {
+ X0 ?$ S" @: _4 K1 O                 
                 pkexamine dx240 temddy240
3 Z0 }8 c6 d4 f9 D   
( K3 H* T$ o* H9 s4 i$ B; E  p2 t) Y  local tempminplddy240=r(tomc)
) Q8 l' y- O" u, N6 B# r6 r                restore
                 }
& h/ Q2 B7 F8 E                 
1 Y; e) k0 W2 k2 ]* i$ M( b7 C        else {
               
                 di "cannot find pminplddy240"
% t. Y1 c, {7 |- r: A' c                restore
# h  D7 n' h( R$ ^& l  u+ e3 U, M                 
                  sum dx240
7 j! T! J+ u; T" _                 local tempminplddy240=r(mean)
' k$ D# ~+ D) @                 
                 
                 }
/ s) s' {3 k% I: @  f   
  o  O1 p6 R* ]5 q, v  g   
) B% o4 a4 {+ x* @   
   
   
3 m) s# u" Z- j% ]% t/ [: @        gen pminplddy240=`tempminplddy240' in 1
       
        *-------------ddy240pmaxp左侧高峰值 -------------------
: r& X% [) v" L( E8 o1 K       
' k# f5 O2 [4 n  T         preserve
! W3 b% K. _4 r; B   
1 B" H6 D# p+ D9 ~4 C+ |' }$ G3 O6 V   keep if dx240>`tempminp'
$ B; t% U4 ]; d   
   keep if dx240<`tempmaxp'
5 Y. `: m, U1 h   
1 b  P% [/ F& H. h* \3 z   keep if ddy240>0
   count
               
                if r(N)>2 {
                 
                 pkexamine dx240 ddy240
3 b7 A" D* r* U* I+ i+ u. T   
   local tempmaxphddy240=r(tomc)
   
; T# [& _* X2 k" P   restore
3 u+ }8 c4 a; z6 g5 ?6 Q                 }
( M0 c/ u! y: L; ~( h. M: M                 
9 y3 B" H; ?' c! f9 Y- |! X) f        else {
& p7 A# v7 Z& V$ m               
. W+ h' s; a/ I- W9 f! q9 \0 i                 di "cannot find pmaxphddy240"
) T6 j! E+ }+ X- h5 v8 f7 _                 restore
2 j# i! Y4 _6 x4 J, {               
  x2 N( ^8 ?9 s( O# B+ U1 m! p                 sum dx240
- |  K, f2 D  j; Q                 local tempmaxphddy240=r(mean)
                 
                 
                 }
- k, z! w3 P* Q( X4 b1 E7 A" V   
) J7 ~% `5 a6 M/ g! s0 }   
   
   
   
   
% G4 V- m: G1 u+ B8 R) k8 u        gen pmaxphddy240=`tempmaxphddy240' in 1
       
7 H/ g2 V; ?1 L       
        *-------------ddy240pmaxp右侧低峰值 -------------------
0 E/ O. ^2 t2 Z' J( n( y% ?       
         preserve
  s; q7 p! g' M5 M$ J+ J7 ~2 ^  
. d! K; o- r) {1 T  c   keep if dx240>`tempmaxp'
% Z( x& y7 u% z6 W8 ]   
   gen temddy240=-ddy240
  |& p4 p. Q  [, s7 y   
: d1 X2 k, p, \- I  s   keep if temddy240>0
) a' S8 `- m, ]* o# S0 C   
   count
- ~- c+ V- M6 L  u8 G1 }+ y               
: y* e8 D9 e% J7 F. E7 H- {# Z9 A, s                if r(N)>2 {
, T; `: c  z# |3 D) z& z$ j$ `                 
& a3 {4 J2 v8 G# O+ X0 A( q8 t( Q                pkexamine dx240 temddy240
- p* P+ @3 ]0 b3 v$ l; D6 S( d   
* V; ?) A7 A6 g0 h5 ]' t   local tempmaxplddy240=r(tomc)
   restore
                 }
/ Y/ p0 n5 W. F$ j1 {) g$ g                 
9 s' z2 z; ]( h+ a! m' E        else {
; a; G; V9 k# x! M                restore
# a1 c* ~/ W' B8 [% Q                 di "cannot find pmaxlhddy240"
                 
  Z4 y9 F" c' C9 k2 Q               
                           sum dx240
* A: P1 k  n- S: l. B9 l                 local tempmaxplddy240=r(mean) in 1
               
                 
                 }
- v5 }3 ~5 ~. [( ~   
   
   
: U( c/ O" c" n. T0 A7 a   
5 f2 e4 I' G: b3 p% y   
: e) M- V+ M, x! _   
$ B: S7 ^& b6 T0 `+ I+ y. O        gen pmaxplddy240=`tempmaxphddy240' in 1
/ U) d) S3 B! j4 w       
        *-------------------------------------
! h* ]  [* J/ y   
5 N+ m2 T2 n5 V, b0 z   gen price5d=v5 in 1/20
7 P2 m2 Z" ]7 c      
) F2 P" ?; h6 T   gen price5yy=_n in 1/20
. n- N  b: {. d' t* j6 |% H' Y$ E   
% t( b. i- W5 V, x- p   gen price5y=price5yy*(dy240max-dy240min)/20
# ]/ m6 X& a7 n% z4 E' o   
5 a. ^' A4 T8 b   *price5y &price5d
1 I$ Y" h0 E  j4 s3 U; S; C1 y   
   gen priceny=price5y in 1
" [+ `8 y/ }+ F  J1 o1 }   
. n, h1 z5 w8 @   gen pricen=v5 in 1

7 \5 h/ {& L  K! G& ]1 k( f: [7 ]     
: i& P* F4 u, c- ^. j   *priceny &pricen
/ d9 F5 i: M" N: E9 k1 k: W/ V* j   
  M' |; k& }* ]& N9 U' O, P   
+ Y7 U6 ^4 v4 P- {   *-------------------------graph------------------
  
% `. w2 B6 E# Z- @9 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)
. u8 }! O) d4 s
8 v  D" a" [; K2 I4 V* ---------可选项目----------------------
8 K' \2 J! d! I; _! V" H* H*(spike yxis amean,lc(blue))(scatter yxis2 amean, lc(blue) msymbol(+) mlabel(amean))
}
end

maodeng

awer 发表于 2013-2-5 14:40
抱歉，最近太忙了
% O! r; u8 {% p8 t2 N7 U
6 a) s1 T- _& N* ]& d; `300236 我用下面的程序试过，没问题

哇，谢谢，我仔细看看，这些字母数字代码晕了
" _6 C- ?) L  n& a" B- U

喜欢喝冰茶

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

喜欢喝冰茶

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

- I' {$ A) j' [
这东西是开源的？知道有个统计软件叫stata，不过要收钱。

汉水云潇

准备试一试
0 H+ T9 R4 }- u先找Stata软件

汉水云潇

2 I8 g5 k3 p5 v3 {
晕呀
论坛复制程序代码
乱码夹杂在其中
很烦
对了两次
以为没问题了
, U, x0 H% Y! E# M
结果
5 l  I3 H8 K) M! N: m' L4 V, y执行后
: a  M( l) f9 P# O还是出现
. m* \) _6 A+ t, c- w! h' Or 199
的错误

汉水云潇

试验成功，下面是三一重工的图
1 L# q1 I/ `7 ?  r, n! c6 K
0 K$ E) D! C' K6 f0 L

peterbev

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

中秋不圆

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

bknite

这个帖子值得学习。

jeanny

感谢！！！

Hubert2000

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

% i0 z5 d! @- z0 A8 h7 X老哥，也是山寨这篇文章找不到了，不知为何？