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

maodeng

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

maodeng

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

maodeng

. hello 300236
$ I" a: q5 [% w3 R' Z. }statement out of context
6 M( @, i- h. {# t' C9 K0 or(119);

$ f' z0 @5 j' S" f大师 就是这个

sergio

这贴实在！

awer

maodeng 发表于 2013-1-28 13:13
- D- w0 q+ Z) a& B/ o9 G8 E. hello 300236
statement out of context
r(119);

4 h! v+ x" s% F! b3 e抱歉，最近太忙了
& {& n7 X3 t5 R. i2 Q7 _$ k
300236 我用下面的程序试过，没问题
8 R3 |' ?8 P3 |9 ]9 r) v
2 b: ]4 N/ C% W: |, V你把 D:\stock\test\source 改成你导出的数据存放的目录，再试一试

capture program drop hello
& H2 k) X8 T5 w4 Q9 h
program hello
0 S# W& a; z& Q$ x! n8 P4 q% B, E* H( r
quietly {

clear
( \# g" t3 [* m; O2 V$ {
8 P+ m3 m' W8 Y2 ]& Pcd D:\stock\test\source

insheet using `1'.txt
) [4 h+ F- r  X7 E0 d
   keep if v5>0
   
   gen vt=_n
   
9 a" d3 J: q  ?$ L) p6 [   tsset vt

   gen zf=D1.v5

   gsort -v1

   gen avol=v7/v5
   
   keep if avol>0
2 B7 E' K  ^- A" P" ~' u3 ~   
; Y' ~9 t% j% Y6 u   drop vt
% P3 `* M/ |% \1 B( i0 J$ E  
   gen vt=_n
   
   tsset vt
. e0 @% M1 [' g2 C2 e, o   
- X) ]- b2 y9 R. B5 G) H* h& y    sum v5 [aw =avol]
" q( G) f, o. h! p
   gen amean=r(mean) in 1
   
) F1 x7 C7 ^, K/ E7 Z7 v% H   * 总加权平均值 amean
   
   
   keep in 1/240

   sum v5 [aw =avol]
5 E5 a  W" Q, E1 u* B3 D
   gen mean240=r(mean) in 1
   
: O# [5 F9 H# A( r3 X   *240加权平均值 mean240
   
4 p# C$ R' ~4 _, i( [; Y0 V   set obs 500
5 H2 Z) V- A( b) ~+ v+ C   
  j. ?' j( y6 z/ k   kdensity v5 [aw =avol],gaussian n(500) generate(dx240 dy240) nograph
% W* y) g# `$ T3 e# S   
) n7 N4 c2 P$ |   *240 kdensity(dx240 dy240)
   
. c1 f6 i( Q* y4 j' b9 L" L   sum dy240
  v6 ?2 O6 X% R0 C. j5 m  c0 R   
% Y' k- G! H$ ?' _/ O! y8 i! M   gen yxis=r(max) in 1
" c$ F" X- Q" o- p2 O8 ]' ~, E; p$ I   
2 f1 j% P& S- l0 W7 I   gen yxis2=yxis/2
2 r2 U, u+ q2 M; }' n1 L5 E   
   *yxis
% y7 [6 ~7 Q, ]' t0 z7 Z/ ~/ P   
" t2 r7 Y) F- Q" b! f5 b   gen dy240max=r(max)
  _3 U8 g$ W5 i/ o2 b+ s: m( [   
   gen dy240minp=r(min)

0 }+ X6 d% K! e7 }2 v    gen vtt=_n

2 w* m+ b# J, z' c) g* I) G, l9 j    tsset vtt
0 A7 a  h; h8 [6 ?/ j& _
   gen ddy240=D1.dy240/D1.dx240
! z" B- {) G" D   
3 C& N/ D9 q6 {8 j7 T   sum ddy240
   
7 s# U/ \- `' T7 o0 t( M2 {   gen ddy240max=r(max)
   
. X; P( L* e! ?' J& n3 P; a0 h   gen ddy240min=r(min)
' B: }) p/ W- x     
   gen ddy240g=ddy240*(dy240max-dy240min)/(ddy240max-ddy240min)
6 E" M/ v: a9 t   
   sum ddy240g
   
   gen yxis3=r(min) in 1
   
  X% u0 f$ J4 i# Q  G   gen yxis4=0 in 1
) N; T/ a5 R, _  R& u( j  
  *ddy240g
   
9 _0 Q3 k* f/ G3 D   
  preserve
. S2 f4 J0 a/ O7 Y) n2 ^5 }8 w) G  
5 W, n( b! ~: w- S/ W" ]- W  sum v5 [aw =avol]
  
4 c8 A3 W) `! Q7 \  keep if dx240<r(mean)
  
& @) K7 u4 N  a) d4 u6 \; D* R  sum ddy240
4 M9 {3 t! a+ i3 f" S  
/ s# S6 o( b$ X+ R: ]  if r(min)<0 {
) `, |$ }. M8 o( {        
                keep if ddy240<=0
               
                count
               
* z& M1 C; `2 `                if r(N)>2 {
                 
                 pkexamine dx240 dy240
, E9 f# Z5 Q7 ?2 a! H; r; V               
                local tempminp=r(tomc)
( ]* C! p, K0 x  {- t                restore
                 }
                 
7 e# N# W0 _* p' K7 Q        else {
               
                 di "cannot find minp"
                 
. U4 ~& {; Q+ Y2 S                restore
, X* g( Y* U+ [, ]- P                 sum dx240
7 O5 i; T# ~. S6 c& l                 local tempminp=r(mean)
/ x( ^* ?$ |% F! d9 G                 
                 }
                }
   else {
4 i  m0 R8 W3 d# F                 
                 keep if ddy240<=r(min)
7 _/ h0 S; U0 x: ?+ I* B& M                 
                 sum dx240
                 
" f# S" B% N9 p; z; W6 }2 C+ ~* W                 local tempminp=r(min)
0 V( I( d) n/ V3 a& }                 restore
                 }


4 K0 b5 ]3 L! r
   gen minp=`tempminp' in 1
9 N$ W+ y1 U# |" T   *240 minp

/ p$ p. J" `6 L6 v7 o9 T& }$ _: D
4 Y) Q; b: j; H+ h. F* `   
7 Y0 `3 Z, G* k   preserve
   
   sum v5 [aw =avol]
/ c% [: `9 q& ]6 n   
   keep if dx240>r(mean)
   
# [' x- I8 N1 Q! S; r' W1 g+ ^   sum ddy240
   
" Z! r! e7 F$ H/ x, q) _% K   if r(max)>0 {
3 T$ a! h: f) U! f               
: Q9 q) t/ z5 Q4 G1 n7 d                keep if ddy240>=0
               
                count
               
1 O9 C1 w% V' p5 F% c% N                if r(N)>2 {
                 
/ L$ x4 N# W5 u  K. _& S* e. F5 v                 pkexamine dx240 dy240
               
                local tempmaxp=r(tomc)
                restore
                 }
2 B; W5 U, H1 I! c* Z" q& e                 
; C3 V& b* n, D( j; |$ |5 o                 
        else {
                restore
2 m; m! r2 ?7 w- U; ^                 di "cannot find maxp"
               
8 k# L; G* F0 z0 e$ T: Q( i3 k7 c               
                 sum dx240
+ x  [" d) |  K: C7 Y                 local tempmaxp=r(mean)
8 N( i9 b$ L) H$ c3 e$ X" s4 H                 
                 }
                 
9 _/ s/ y; H7 p                 }
                 
    else {
       
2 F) ^# C$ G4 ?1 _+ |; i           keep if ddy240>=r(max)
          
           sum dx240
          
           local tempmaxp=r(max)
           restore
           }
       

  
4 j& i7 e5 {: y; @% O
   gen maxp=`tempmaxp' in 1
   
3 q% z+ ]+ l3 J0 h/ \: ~   *240maxp
# T+ s* N6 Q$ I" K9 ]   
   
" J9 N& H; o  a3 V$ u3 T: Y, o   *-------------ddy240pminp左侧高峰值-----
$ V; ~- h; x$ S2 h2 a' I( X- F   
    preserve
   
   keep if dx240<`tempminp'
   
9 L* l# W- }: L" P- S6 Q/ p& G   keep if ddy240>0
7 R: g- p- [  }0 c+ R4 Q   
4 M1 m4 K: A% C4 y. N4 u   count
: Y$ o3 h4 d! C0 F4 Q) Q               
                if r(N)>2 {
: P; ]. y; L3 m* T+ |2 r                 
                pkexamine dx240 ddy240
* i- i: ^& U4 }6 Y               
                local tempminphddy240=r(tomc)
5 j3 X: ~5 |) P                restore
                 }
                 
        else {
                restore
$ L7 Q( C" L- P9 Y" H                 di "cannot find pminphddy240"
                 
3 i6 O, b, w0 J5 b1 p  V& D               
+ g' w9 {: [- n: q; l3 l# ~0 k  A                 sum dx240
0 y# x  [% c1 _9 h8 @' U* P3 Q                 local tempminphddy240=r(mean)
               
                 
1 |  }" ?1 j  l6 ]                 }
                 
& C( Z. \) f. e5 G; P. J   
* r, J7 ]$ O9 k% s7 q- J" {; K# x   
& Y7 \* r' }, L# z/ a   
5 U  ?5 n+ i# J- Y: R6 u   
   
" Y" e( r, b  ~4 G- S) Q. C; `        gen pminphddy240=`tempminphddy240' in 1
       
$ r$ b( M/ u" {/ H* m3 L       
        *-------------ddy240pminp右侧低峰值 -----------------
( L  [  C. C# n% j    preserve
4 O! C6 c  |* e7 B5 N
$ {. ?. I% k) V; `8 k0 f   keep if dx240>`tempminp'
   
   keep if dx240<`tempmaxp'
   
$ H' H* @3 L& i/ s! e; s3 V" X! A   keep if ddy240<0
   
: g  }& O% I5 o0 S3 K; ^: {/ [9 M   gen temddy240=-ddy240
   
    count
               
                if r(N)>2 {
                 
3 G7 H4 y, t" c* T0 E0 S                 pkexamine dx240 temddy240
   
' a0 s$ x& w# x  local tempminplddy240=r(tomc)
; Z6 X, N; R7 u0 `8 l( H                restore
, @) q: l" @' {7 b! R5 v                 }
3 G6 g2 D# C" W) S) A; l- u% K                 
        else {
               
- u" F0 d+ {2 A$ m                 di "cannot find pminplddy240"
                restore
1 u/ X3 X, X0 D2 }# v$ p, G, a8 P                 
- L4 T3 \8 U* W& M% d4 v9 A                  sum dx240
! O/ R/ b' E, L8 Q                 local tempminplddy240=r(mean)
4 z! v$ `& v! K4 A* f1 m( [/ _* s                 
7 e+ M$ P6 R! G                 
                 }
   
0 c% e4 V- ?) D2 Y. Z   
6 [6 S8 T* ^% V   
# x( f+ e$ n3 i8 @2 c' b' P! C   
   
) Z" ]8 v7 j2 F8 {! @0 h. z1 i2 l        gen pminplddy240=`tempminplddy240' in 1
' N! M3 _. t9 N       
& j7 z) c, A( k, ~! w        *-------------ddy240pmaxp左侧高峰值 -------------------
       
         preserve
  \1 n- |2 V8 N* O6 e8 T3 b5 |   
3 d/ b# I* W$ u$ o   keep if dx240>`tempminp'
* F' N' @( E' |) T+ X   
   keep if dx240<`tempmaxp'
( P" t8 Z  y- k* {   
   keep if ddy240>0
   count
! {: S6 K& b& n: }9 ~5 d  D               
9 _; e7 {: {( d& _1 H0 P) {3 \                if r(N)>2 {
                 
                 pkexamine dx240 ddy240
   
   local tempmaxphddy240=r(tomc)
   
$ C8 E) G. h9 j! m: q6 d   restore
# I5 X/ G$ a: o7 ]- Z! W                 }
4 m1 F1 b8 q, D4 m, D$ E' Z" }" I* |                 
        else {
               
                 di "cannot find pmaxphddy240"
3 E2 O% M5 R# u2 k: y+ F                 restore
               
( E) C+ K- [/ J$ p5 U                 sum dx240
/ [6 N; [3 U) A, p- u                 local tempmaxphddy240=r(mean)
                 
                 
! X( o; y1 b; r6 h; b1 F5 t                 }
3 \" v: K& {/ |' V   
   
$ n4 w& a) s% P% J5 J   
   
   
9 U* n0 o$ ]9 r# n$ s* d7 ?   
        gen pmaxphddy240=`tempmaxphddy240' in 1
( h7 C6 x/ C, N" z1 v' O* b- ]       
( A, F6 s" Q3 ]       
        *-------------ddy240pmaxp右侧低峰值 -------------------
       
         preserve
( h+ ^- V0 Y2 r0 Y1 x& m  
   keep if dx240>`tempmaxp'
   
: V: \2 j5 a4 T, @   gen temddy240=-ddy240
' F7 p4 n4 B& E: Y, r) E   
   keep if temddy240>0
   
- @: e/ n6 U" w7 {/ {+ E0 Y7 \: w   count
               
                if r(N)>2 {
                 
' S1 ]  h# ~  C- n* N                pkexamine dx240 temddy240
, i2 B3 G" S- H2 d0 G   
5 L! u6 B- g- ?2 z4 F   local tempmaxplddy240=r(tomc)
) {/ g2 J( f3 n' c5 o+ b0 p. A   restore
                 }
4 _. D$ `6 u* Z; q9 S                 
        else {
                restore
  o: N9 _) g$ y4 J                 di "cannot find pmaxlhddy240"
                 
' }  m5 ]8 P7 l0 N               
                           sum dx240
                 local tempmaxplddy240=r(mean) in 1
               
                 
: O$ z3 B0 J2 f0 d0 h9 X! Q                 }
   
3 q4 t' ~5 ~6 H6 M& ]1 K4 K3 ^   
   
3 Q! e: F) X$ S% v; S   
   
   
5 \* T# @: {7 o; i        gen pmaxplddy240=`tempmaxphddy240' in 1
       
        *-------------------------------------
   
   gen price5d=v5 in 1/20
. E! }/ O- r  T3 K" \* g! B) H0 z      
8 J# x! d  r1 ^4 d% C, I   gen price5yy=_n in 1/20
   
: F$ {2 D: L) n2 L6 M8 z   gen price5y=price5yy*(dy240max-dy240min)/20
) W  ^* c+ Y- Q7 X   
* t$ I- w2 o& W2 i   *price5y &price5d
   
   gen priceny=price5y in 1
, N7 h8 Q# D, I   
   gen pricen=v5 in 1
1 g+ S: h0 G+ g* X0 P: o7 j9 g
     
   *priceny &pricen
8 A' [. `, K" M+ w6 |. ^   
   
* Y" D- j9 ~6 U  y( e' I; |   *-------------------------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)
) h% a. S' q! i% h) |( o
* ---------可选项目----------------------
0 `9 c2 }; S* }, J% N4 C*(spike yxis amean,lc(blue))(scatter yxis2 amean, lc(blue) msymbol(+) mlabel(amean))
}
) R( z! y  W3 N% U2 h% Uend
& Q, y. Y- {1 x+ a- Z  N* \) {/ D

maodeng

awer 发表于 2013-2-5 14:40
抱歉，最近太忙了
' X3 G! P7 E, D* z0 U: ]7 e
% z( T/ d* }, S9 N9 i300236 我用下面的程序试过，没问题

哇，谢谢，我仔细看看，这些字母数字代码晕了
4 B: Q  {6 e% Q5 u/ H! X

喜欢喝冰茶

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

喜欢喝冰茶

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

# v1 K( [0 {" H' x这东西是开源的？知道有个统计软件叫stata，不过要收钱。

汉水云潇

准备试一试
先找Stata软件

汉水云潇

$ t; s- ]' D% R# _, ?. Y2 K5 N
晕呀
论坛复制程序代码
乱码夹杂在其中
很烦
对了两次
1 F' N/ R/ W3 W& l. Z5 T- i# C以为没问题了

结果
- {3 D+ H: K: U& o1 p9 v7 R; T7 v6 z执行后
还是出现
$ g) ~3 `! @& u- ?% j( _" ?5 p* O/ ~r 199
的错误

汉水云潇

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

peterbev

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

中秋不圆

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

bknite

这个帖子值得学习。

jeanny

感谢！！！

Hubert2000

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

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