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

maodeng

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

maodeng

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

maodeng

. hello 300236
statement out of context
r(119);
2 k) p  b+ h' N& @8 e. Y4 X3 L
大师 就是这个

sergio

这贴实在！

awer

maodeng 发表于 2013-1-28 13:13
# I$ O' f; d/ \1 X7 C3 {. hello 300236
* D0 E, J% ]/ q% j. Xstatement out of context
5 ]1 R5 c1 D1 |$ T$ v' t2 l* m9 R6 }r(119);

抱歉，最近太忙了
1 y8 a& Y( Z. _9 n& z+ O
; H+ g$ d" }* A% v+ r6 Y300236 我用下面的程序试过，没问题
, `& L+ {! u* }" R' E9 c) ~+ T  U1 Q
你把 D:\stock\test\source 改成你导出的数据存放的目录，再试一试

& P9 @- G4 b; jcapture program drop hello

program hello

quietly {

' C7 K7 a4 Z1 E" K/ M- @clear

cd D:\stock\test\source

insheet using `1'.txt

0 _; e6 U7 O, U+ n  |2 H   keep if v5>0
   
1 D# W/ k! U6 j+ o   gen vt=_n
) {+ @: x- Q/ @. k( X6 ~! T   
   tsset vt

   gen zf=D1.v5
" w% L( H. o* S2 M3 z
   gsort -v1
/ X8 D  O7 u5 _7 o2 p  ~
8 P1 i4 l% ~* v5 e# V/ m! \   gen avol=v7/v5
   
: J' V9 ~/ |! D; Q- g   keep if avol>0
   
   drop vt
/ N% V- P5 d, Y# J* x3 C$ ^  
   gen vt=_n
' R' u6 E. G. G# h/ b   
- @& b/ h- m) @( d. G   tsset vt
: y/ e$ g( K8 C" _& `   
    sum v5 [aw =avol]
& N" |9 ?; \8 X, x3 ^. C
( m; i+ \; E$ n4 o& d5 q   gen amean=r(mean) in 1
) i" y) N3 ^3 U3 Z, W   
   * 总加权平均值 amean
   
8 Q. {/ r7 ?# H& e! F3 @  I   
& a" G) S  G% ]  H, v: p: |- c; ~: {   keep in 1/240

; c, m; Z+ }& |4 q+ Q6 K  W* L3 z   sum v5 [aw =avol]
- |; d( F2 u+ A; M- O8 e
   gen mean240=r(mean) in 1
; x5 \* r. ~2 T$ C" x   
   *240加权平均值 mean240
   
' g' k& I+ q! i7 d4 s1 z; r' C. k   set obs 500
) E; n3 p  M' d* g6 X& g   
5 }" M- V9 p2 R& }  F   kdensity v5 [aw =avol],gaussian n(500) generate(dx240 dy240) nograph
) Z9 P& A6 @; j! J   
   *240 kdensity(dx240 dy240)
   
1 u& x' D$ ?5 P9 ?/ e1 N1 `3 e  Q   sum dy240
, s! ~- k: J: F: x* \- u. ^   
   gen yxis=r(max) in 1
: l& V4 M. @  _4 C2 ?+ l% P   
   gen yxis2=yxis/2
+ p' y5 _0 M7 P0 b. h$ _   
   *yxis
# ?0 I0 U; W5 u   
   gen dy240max=r(max)
   
   gen dy240minp=r(min)

1 _  h1 L( m3 H# \; V5 J6 R$ o% [    gen vtt=_n

& |* S) w) j. T" e    tsset vtt

   gen ddy240=D1.dy240/D1.dx240
   
   sum ddy240
1 l% H8 W' U1 P7 |6 k   
   gen ddy240max=r(max)
) T& b! G9 B% M   
$ j. K9 ~) p' a. ^6 C0 P% l1 J, j   gen ddy240min=r(min)
' v2 R( @( k0 c0 @2 C5 C) \     
   gen ddy240g=ddy240*(dy240max-dy240min)/(ddy240max-ddy240min)
2 ~9 y( o- b9 O. t) s+ f+ q) d   
   sum ddy240g
   
   gen yxis3=r(min) in 1
   
   gen yxis4=0 in 1
  
: m+ b7 R: U9 n7 r& x" O3 T3 {  *ddy240g
, E9 R: s* t1 F7 O! _5 Z" m1 s   
   
; o6 U, O! B) l7 E: @" A* F1 C+ F0 H% P  preserve
  
/ a1 E. S' I; S/ y  sum v5 [aw =avol]
  
  keep if dx240<r(mean)
1 |2 t+ ^5 {# R. ?( r, U9 q  
  sum ddy240
5 A) I8 }4 @) m4 I1 T4 @3 d  
  ]8 j4 S3 B0 z& w- c  if r(min)<0 {
        
                keep if ddy240<=0
" d& n" B$ c  B' P0 D' ^8 I/ N8 L               
                count
               
                if r(N)>2 {
                 
: K- G9 m1 T; _9 W* z2 I                 pkexamine dx240 dy240
               
                local tempminp=r(tomc)
                restore
                 }
. \6 N- q  b3 t/ _! g$ X9 ?) f                 
        else {
% X* D# I2 g8 p# q) N/ ^! T               
                 di "cannot find minp"
                 
& H8 R$ @" ~/ V- e/ ^                restore
2 M8 H% X) p, {5 D) r/ B) d. ]                 sum dx240
                 local tempminp=r(mean)
                 
" b3 l; w% D; Q" r! n                 }
: a- o9 @  Y7 o* W. C- P                }
' x. P. V4 D. n4 K   else {
                 
( m. M: Q! m" @, w                 keep if ddy240<=r(min)
* m  p. r$ E% x) R                 
, t: [+ h% ?4 A2 }; {7 I                 sum dx240
                 
                 local tempminp=r(min)
                 restore
                 }

" S( Z# n* ?* k9 F) {4 N9 }

   gen minp=`tempminp' in 1
0 u$ \: [$ T+ q. Z, L4 E4 E   *240 minp

/ d; Y( L9 Q+ O; P/ ?/ w
' K5 N, t& j- ^2 J2 l   
   preserve
   
( ?4 f, h3 p4 [9 x# O   sum v5 [aw =avol]
   
   keep if dx240>r(mean)
0 }( `* [  c2 {7 x- B   
/ R$ s% ?3 M6 d% d3 U, y   sum ddy240
   
1 `# j( t/ J9 d; ], S* Z! K$ y   if r(max)>0 {
. w& y( u8 M) |" V( D; Z- R               
                keep if ddy240>=0
6 f" I) e4 ~. p1 @               
, Z6 T( R$ E7 N# ?& j6 Q( u2 }                count
               
4 Y8 _0 q0 n2 O  ]/ E- v; a/ }* R                if r(N)>2 {
7 W" P5 O/ C- t) @                 
                 pkexamine dx240 dy240
: I" f( \: m) B! Y               
6 P' L- m" `3 J8 ~8 Q# s                local tempmaxp=r(tomc)
                restore
                 }
                 
                 
( V. M7 v; J6 l  L9 N' `: Y/ x' [  y        else {
                restore
                 di "cannot find maxp"
; H2 c" P  H6 o  c3 l               
               
' Z& @& y0 N- n5 J                 sum dx240
/ K' H$ g% \! s8 A2 j4 i                 local tempmaxp=r(mean)
5 C" L/ V4 T: l+ a                 
1 d4 \0 x) _. _                 }
4 h, L3 ^" v- P' f) {& w                 
                 }
                 
) @& S5 c* p( w9 u9 Z& |    else {
       
           keep if ddy240>=r(max)
" O' M0 O# S* `' F$ l4 E3 I2 J          
           sum dx240
          
+ O: o! \! C( I' R0 m! y1 u+ Z7 N. T           local tempmaxp=r(max)
$ k$ M  n0 F8 s           restore
4 ?4 J: @- f0 p/ g) C2 ^           }
       
1 [& H7 L+ F% L+ R# T; @& `5 ?
  
/ t; I: L2 N/ i
   gen maxp=`tempmaxp' in 1
   
) v- e/ [, W# A% P" ]   *240maxp
   
6 K, [4 E3 R& I6 @& Y5 w   
   *-------------ddy240pminp左侧高峰值-----
8 H* a; u9 n7 |; P; v+ G   
) U+ {2 d- p: b( ]. n7 R  V$ v    preserve
   
   keep if dx240<`tempminp'
1 A4 \9 z  H7 F   
   keep if ddy240>0
$ d. }- D: g2 ]- `$ }: A6 t   
! @/ B" U9 [4 t0 t" Z5 n   count
               
9 M) Y2 \5 ~) S- `% u                if r(N)>2 {
                 
) ~6 F/ w, W: w; K* D* ~                pkexamine dx240 ddy240
               
                local tempminphddy240=r(tomc)
7 `% C/ `3 S5 n. `* H                restore
0 y  ~, k3 e. d5 T# c) g9 }2 U                 }
6 z- F" W. U. d: P                 
( ], O% y  E; z# z) p% v* W        else {
                restore
# k0 o& X! ^( x* l                 di "cannot find pminphddy240"
3 Y8 H- |( D. l- l( [" R- _7 z                 
( B7 @7 X+ c/ `4 Z$ c+ h0 M               
                 sum dx240
) |4 q6 b+ B4 O$ F( A% e                 local tempminphddy240=r(mean)
' w" l0 C8 c* C. Y               
                 
0 m, N& }* d' P" o# G" P5 h, J                 }
' l$ `9 o  M* B; }* \: x/ B0 ^                 
   
, m1 c" k4 k$ A& \9 W   
/ G: l# Q$ j; }( y   
   
   
        gen pminphddy240=`tempminphddy240' in 1
3 [$ z  U1 p/ k0 L9 x0 K       
       
        *-------------ddy240pminp右侧低峰值 -----------------
2 `7 {& I) c; m/ V1 c$ W    preserve

   keep if dx240>`tempminp'
   
; N8 [4 v/ K! ~# p! i   keep if dx240<`tempmaxp'
   
   keep if ddy240<0
   
   gen temddy240=-ddy240
2 C7 e; p$ r% i: c$ W+ ^1 K   
    count
9 m! f; s& T5 R1 \5 A$ d' e               
( X( M9 z$ c5 V! K% ~) \: H2 j/ l                if r(N)>2 {
, B; \" C6 i- Q- `, \                 
) p$ s: x/ G$ \1 `% [& M) g* X                 pkexamine dx240 temddy240
   
  local tempminplddy240=r(tomc)
                restore
                 }
                 
1 E# A: f/ P; C, j        else {
               
                 di "cannot find pminplddy240"
                restore
                 
                  sum dx240
                 local tempminplddy240=r(mean)
                 
                 
7 Z$ x- w: T2 ^" l* N                 }
% s2 X  E' h9 U$ s& l. s* f   
   
   
1 I9 q/ J7 F# C1 A, p* z- F1 ^   
   
/ W& a- v& ^$ C  v        gen pminplddy240=`tempminplddy240' in 1
       
        *-------------ddy240pmaxp左侧高峰值 -------------------
       
         preserve
& T7 h2 A1 X3 |8 n8 M6 F4 \  R' b   
* I1 P. {2 l" h7 u, h" {   keep if dx240>`tempminp'
/ h  ^0 d* M' D! ?   
   keep if dx240<`tempmaxp'
   
   keep if ddy240>0
* M7 N( l, S2 }+ W. C+ N, T   count
               
                if r(N)>2 {
                 
                 pkexamine dx240 ddy240
: s  \/ t" I7 D9 H) w* k" v   
# d, F, F2 [% H0 J$ p   local tempmaxphddy240=r(tomc)
   
4 t1 }1 D+ d  t; k* m! n, S8 i   restore
" r1 A- O, j4 `4 z, {6 G                 }
; J2 k" n/ w5 \( c8 N+ j! U                 
/ N8 T, N" p$ L/ q        else {
               
                 di "cannot find pmaxphddy240"
                 restore
               
8 m$ _1 s5 `( G% {; O" ~                 sum dx240
& {* ^& {& n' Q. N8 R* Y) R6 G                 local tempmaxphddy240=r(mean)
2 [7 R# W" a- X) @- c                 
0 o3 C( a. l- r3 u9 F/ a1 V                 
                 }
' s% ?2 F" V( q6 l. O   
   
   
   
6 H2 o$ Z3 i5 k5 M   
   
        gen pmaxphddy240=`tempmaxphddy240' in 1
3 l2 C1 c! i2 ]+ W( r% ]3 Z; d! `5 p+ f       
       
3 Q* }6 G) z2 p# p8 P: P9 f        *-------------ddy240pmaxp右侧低峰值 -------------------
. |) {4 H- m$ j$ |7 K' H- b       
  B5 q) s( v. `. q$ x* Z         preserve
  
9 }; Z; [8 B, o6 S+ C   keep if dx240>`tempmaxp'
   
$ ~; p5 m( s4 t: |. ~& p- X: i: P   gen temddy240=-ddy240
   
* A2 ^0 `! o7 w) }9 N1 \   keep if temddy240>0
   
   count
               
- m' X$ {' x) B                if r(N)>2 {
                 
                pkexamine dx240 temddy240
   
   local tempmaxplddy240=r(tomc)
4 L. ^8 ^( ~/ N% A   restore
                 }
                 
) z1 E) L& t( f: C$ t        else {
+ M; S2 a' y* ^  p  M                restore
                 di "cannot find pmaxlhddy240"
                 
' p) F. U+ j/ H8 W. Y; `               
                           sum dx240
                 local tempmaxplddy240=r(mean) in 1
               
9 G# E- Y- n6 b                 
# g( Y/ R. b6 Q7 C) J( h: U                 }
/ k. d4 N* n+ t. R; A3 C2 x5 n   
) [3 o9 F/ j% s& I5 h, ~& N1 r   
" o5 v& s5 {3 F4 i# v   
: P' M& y7 h6 f/ v4 ]. d! q. q   
# o; u- O+ J5 G( I& B! T   
5 S, s# s6 K+ O3 o. \6 t9 z   
# f% ]: D% x" t" X. g' ^/ c        gen pmaxplddy240=`tempmaxphddy240' in 1
6 g: b: r) h9 M- k& s9 o       
        *-------------------------------------
   
   gen price5d=v5 in 1/20
      
   gen price5yy=_n in 1/20
   
7 G- ?5 g6 ^: C- ]8 J7 S   gen price5y=price5yy*(dy240max-dy240min)/20
   
   *price5y &price5d
0 H6 z$ L# Z: v   
/ f' F6 s$ x  Q9 q   gen priceny=price5y in 1
6 l. U  C, T' @. B   
   gen pricen=v5 in 1
) S$ [3 o& C) p+ d6 K
     
5 V6 r% l, a0 g! D4 ^   *priceny &pricen
5 N- _) Q9 }0 l. B   
   
   *-------------------------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)

4 f2 Z) o4 q5 v4 a* ---------可选项目----------------------
/ D/ [  H' N% K4 Q/ l, C2 ]; E! u*(spike yxis amean,lc(blue))(scatter yxis2 amean, lc(blue) msymbol(+) mlabel(amean))
8 M- V9 v  w6 d' l# ~2 M1 i0 f1 u}
) w, t; [9 s# Q" ?end
3 L* [/ w0 T) p7 ~5 t7 K

maodeng

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

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

2 T3 q) `& R' t7 R1 e哇，谢谢，我仔细看看，这些字母数字代码晕了
6 f' B5 M! r0 Q

喜欢喝冰茶

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

喜欢喝冰茶

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

2 {$ ?; `8 \% I
这东西是开源的？知道有个统计软件叫stata，不过要收钱。

汉水云潇

准备试一试
先找Stata软件

汉水云潇

" _+ E  @! \0 X& G6 T# q
: {: t, w7 }, Y+ |晕呀
论坛复制程序代码
乱码夹杂在其中
( G: R" f+ I/ j. d& ?+ e很烦
) J* {6 z0 x% v- G$ X对了两次
以为没问题了
$ p( y9 A$ `( M6 \) X' ^
结果
执行后
' G) x+ ^# z* G; p, P0 |# \: C, Y: f还是出现
r 199
的错误

汉水云潇

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

, h- d& H! T! H+ M3 D

peterbev

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

中秋不圆

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

bknite

这个帖子值得学习。

jeanny

感谢！！！

Hubert2000

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

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