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

maodeng

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

maodeng

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

maodeng

. hello 300236
: B; O: E. W, \8 e. Ustatement out of context
0 D6 m$ s( s' G! o' G4 f( p' sr(119);

6 D; A) S. |6 D" k9 u大师 就是这个

sergio

这贴实在！

awer

maodeng 发表于 2013-1-28 13:13
7 q: {5 a* w; w$ H$ w4 C0 p. hello 300236
7 D  |/ ^5 o9 g. X- x3 gstatement out of context
' K& J0 z" `  p- I: v# w4 r0 xr(119);

抱歉，最近太忙了
$ l4 V! x1 x& b' l
300236 我用下面的程序试过，没问题

" {3 a. e+ a. K你把 D:\stock\test\source 改成你导出的数据存放的目录，再试一试

$ ~$ s; _, ?5 ?8 X* R% @# \capture program drop hello

, ^6 S# u& T; N( e  d/ v8 Z- O5 oprogram hello

quietly {

4 u' M0 }) }0 @0 A% Xclear
7 h! E4 u( B  L, `) V9 V- S/ E
3 C! \$ O: y  E6 n( V) dcd D:\stock\test\source
3 m1 f7 S! ^+ Y2 N3 y- D3 p" d
. ~$ |2 A2 B& r# E0 s" l: Cinsheet using `1'.txt

) k* V0 J' B/ z- L) T   keep if v5>0
   
- Q. l0 Q0 r3 g+ w" _   gen vt=_n
/ j. ~3 ?, G  ]) K: t$ K   
   tsset vt
; V' O9 t! }7 t+ }4 K) O
   gen zf=D1.v5

   gsort -v1

3 A! p. _2 p4 t) c   gen avol=v7/v5
1 p; z! d: ^2 u* w5 A9 K   
   keep if avol>0
   
( i, p: ^0 G3 J7 S1 U, J! d   drop vt
6 r$ p& t9 W/ w& g5 m* A% E( J  
) w  n5 d% {- U; l$ B/ N' p* U   gen vt=_n
9 h( ^9 ^( s9 V# a' [2 g   
   tsset vt
6 V! Q' i" u% |$ c! t! W  Q4 C6 G   
    sum v5 [aw =avol]
+ f4 @$ `1 ~  d5 z! X, u
   gen amean=r(mean) in 1
* z4 w; u; w/ T) F   
   * 总加权平均值 amean
. F$ J! B5 j! }7 b/ f   
   
   keep in 1/240
1 I  F4 T, j' G2 L; [. W6 p0 w
. X4 [3 N! {9 I( e$ z0 C. X( \   sum v5 [aw =avol]
3 Y6 u* d& q" ^. B
4 ]* o+ ?5 j' m( w# N   gen mean240=r(mean) in 1
4 A2 ?, p% e' Q   
: l' w2 {9 G4 T8 d3 A  z   *240加权平均值 mean240
   
   set obs 500
   
   kdensity v5 [aw =avol],gaussian n(500) generate(dx240 dy240) nograph
   
   *240 kdensity(dx240 dy240)
; @$ `  B- W- o0 E; }. l   
  P/ ]0 r4 _  H, T& t! N7 I9 ?- `   sum dy240
   
   gen yxis=r(max) in 1
/ @- \  T, Y6 k9 q% Y7 @   
   gen yxis2=yxis/2
   
) }: b" D! X: `* Q1 T   *yxis
$ w+ S9 x9 ]1 z   
   gen dy240max=r(max)
* X1 N3 t. L  }   
# `8 N3 Q: s; x. b% g' w' O   gen dy240minp=r(min)
8 |' H) D; b/ ]- Q
    gen vtt=_n

    tsset vtt

   gen ddy240=D1.dy240/D1.dx240
   
   sum ddy240
   
   gen ddy240max=r(max)
   
   gen ddy240min=r(min)
* Z& c2 @5 K, {# n* W4 Y     
2 K3 @' {3 c" s  I% X' S5 Q1 l   gen ddy240g=ddy240*(dy240max-dy240min)/(ddy240max-ddy240min)
& f; p2 x) |+ ?! W5 t/ n5 @* S/ p   
. O* O. P5 {+ z   sum ddy240g
# x7 K" j, q3 d2 }5 ^1 S/ ?   
   gen yxis3=r(min) in 1
& I2 M' E5 d. v5 K$ a   
   gen yxis4=0 in 1
1 h9 r6 S# S: B( Y5 M  
  *ddy240g
0 b& Q9 b" P3 H& f   
; U6 Q8 {/ a' A) J. P9 C' A  R6 M   
8 Y) C6 M8 p; K  preserve
) b5 v$ g0 [1 L% b+ f0 w, b# P  
  sum v5 [aw =avol]
  
- U) ~# W& X% B  keep if dx240<r(mean)
  
, @, d1 U( U- F' s- T  sum ddy240
  
1 A. c1 O  a8 @) M9 Z! N$ L$ r: N  if r(min)<0 {
        
# G4 d: o# s) c                keep if ddy240<=0
               
                count
               
                if r(N)>2 {
                 
) z' _6 Z, {' Y% J                 pkexamine dx240 dy240
               
                local tempminp=r(tomc)
/ Z" [9 J0 A7 z6 P                restore
                 }
7 |+ K9 C, q( J( U: G. V; q                 
        else {
% ]2 m, }; a% `# R: I+ E1 I               
  V) K8 p+ j9 O& s! g0 }3 N8 P                 di "cannot find minp"
4 [' F0 r0 R/ X' Q0 B2 I0 {+ n                 
                restore
1 k7 Q- p/ {  A, u8 o" x* f/ g                 sum dx240
                 local tempminp=r(mean)
! _4 c& n3 |" `) i- ~. \                 
& c5 Z5 s( W; b7 o% |* B/ l! `1 B                 }
" m( s3 x5 _! b5 P. A1 l; B                }
   else {
, u. l- m  }' @8 s: Q( H9 r                 
" N5 a" R( E, {0 \8 G                 keep if ddy240<=r(min)
0 w+ Y. x+ f2 H, g4 i, ?                 
                 sum dx240
                 
/ v1 \: L! |6 W* R                 local tempminp=r(min)
. b8 Q8 k* V* Q; p                 restore
) @" t( |( g9 p# K) k( {  I5 c$ ]: H1 p                 }


" q4 Y8 }3 q6 ]. P2 e
# G% F) }0 P& |/ c( r1 N   gen minp=`tempminp' in 1
   *240 minp

# i& z/ V% f3 c
   
8 _9 p# C6 y' T# o4 k   preserve
   
8 `: K$ p! i' s# E   sum v5 [aw =avol]
   
1 D6 ^6 _$ V8 a& y; K6 B( W/ |   keep if dx240>r(mean)
   
5 p, {  \8 Q& [( }1 D7 K+ E   sum ddy240
) S+ U4 S( `8 p4 ~4 v* @   
9 D, N8 U# O( A+ O; v6 U   if r(max)>0 {
5 N* `: r4 `3 k/ n$ n0 _$ l               
                keep if ddy240>=0
1 `5 V/ B5 w$ ^6 d3 {) I               
                count
1 W, ?8 h' M1 g. `, g! ]* n5 f               
/ I' ?: e: }% u! h8 G2 z: y                if r(N)>2 {
                 
                 pkexamine dx240 dy240
* ?8 R) D# w7 |- P5 J4 Y               
' y1 [7 O: P3 W  w                local tempmaxp=r(tomc)
  @3 Y: L' B- b- ^8 {6 J                restore
                 }
: @$ |' a: A4 k+ W9 Y3 X7 J, {                 
                 
        else {
7 [3 y6 u. S* q7 m                restore
3 ^) S+ S. @, [% D5 K                 di "cannot find maxp"
               
               
                 sum dx240
                 local tempmaxp=r(mean)
' j" T& \( @5 Q; E. w8 v. P9 g( S$ a                 
  z! a9 p- K- v$ M9 Q6 t: X                 }
                 
4 Q) s% I5 L/ k5 G1 R2 ^3 u                 }
1 D& Y+ n6 {- }# n& ^( E/ n0 G                 
- ]) G3 x- d" p9 I5 z% z    else {
7 O2 n* @* f! n: A. i& J       
           keep if ddy240>=r(max)
          
  G9 _; E; A2 H/ f/ @$ V5 M  Z) @           sum dx240
. m( I  o+ V, C" t/ d! _& B          
           local tempmaxp=r(max)
" J  U& O, z+ t: f+ i           restore
           }
       
9 m  t" v( [) S3 _) P- j
0 g" w/ g! O4 _, m1 {- ]) x  

" b, P! C; x. l5 W1 V. R; L+ [   gen maxp=`tempmaxp' in 1
1 u3 i' D3 g* l: O% X   
   *240maxp
   
   
   *-------------ddy240pminp左侧高峰值-----
   
    preserve
   
/ d5 a9 f& y( A* Y. w# j% m   keep if dx240<`tempminp'
   
   keep if ddy240>0
   
/ U+ b+ {1 l6 |! m, e* @   count
  c8 ^' G9 t  x# C               
                if r(N)>2 {
                 
! k3 f1 U. Y5 @3 p4 r- I& c/ D0 N* x                pkexamine dx240 ddy240
7 b' A8 L( E* L! s: {2 {! V6 r               
                local tempminphddy240=r(tomc)
' I5 W2 l5 H) p                restore
                 }
                 
        else {
                restore
                 di "cannot find pminphddy240"
                 
               
/ p6 H% ?/ d' }( M* D                 sum dx240
* l: X: S8 `1 s: w                 local tempminphddy240=r(mean)
9 _& q4 X$ z  G8 }* ?( [               
                 
                 }
                 
6 ]! `* S% I) Y3 \+ K+ g8 t8 E3 w4 ^! k# \8 y   
% G' L2 W6 x2 P. V! w) }) P! }; F   
1 v' R$ U2 Q9 J6 |6 `2 {/ W, u   
* w, ^  c' }) H: \7 H- s# B. ]   
, G  U* }! e" W4 g: ?   
        gen pminphddy240=`tempminphddy240' in 1
: {" x& ]2 A" N$ |; I       
! H) n- b* [7 {- ]% g# _- U       
9 y; c5 F/ j& x: z- y8 ^7 Z) L        *-------------ddy240pminp右侧低峰值 -----------------
    preserve
% @% U' b. Y. ~9 w! v! O7 Z) [3 f; ?: t, z
   keep if dx240>`tempminp'
   
, k, R/ a% n5 @0 X  O+ {0 r8 i6 r   keep if dx240<`tempmaxp'
   
. x/ u. |( f1 f- C6 }" S   keep if ddy240<0
8 P1 T3 ], R. ?. T6 v   
6 z( A; e# P+ r/ |& s6 |0 f% ~+ j   gen temddy240=-ddy240
6 {5 O8 j6 N: W, n# }9 x# z% Y7 z   
) g- O3 Q5 M+ z- V8 A5 v4 T: `! t    count
               
                if r(N)>2 {
                 
" [& i7 k! m4 \" p% I% H0 m( L                 pkexamine dx240 temddy240
   
  T% O, P6 U) y% g/ S/ w, l0 w  local tempminplddy240=r(tomc)
                restore
8 R% E% Q! d' t3 v                 }
                 
        else {
               
+ D$ t4 G5 [( i* \3 Y                 di "cannot find pminplddy240"
                restore
                 
+ @# ]) U  P- b4 X0 s& l  ~. n  s, O                  sum dx240
                 local tempminplddy240=r(mean)
                 
+ D& l; k& j: y/ h* U! L( d) ~2 D                 
0 x0 Y6 O- u& T2 `& M6 l1 P. k                 }
5 i# F; y  w4 d# |+ x, ~   
   
   
: j- }1 L1 B5 ~. \   
1 u* H& y& D' o' ~! `0 e& u6 L   
) B$ \+ i6 a8 d0 M) U" v        gen pminplddy240=`tempminplddy240' in 1
; m8 b# N$ V% U" E6 m9 Z       
  E' X6 ~6 N7 V! l3 C        *-------------ddy240pmaxp左侧高峰值 -------------------
. ]( u% O8 S5 y4 w       
3 @; S6 g  h) E9 o$ n9 o! u         preserve
   
/ }7 K# N2 ^$ l# J* \- a4 ~   keep if dx240>`tempminp'
" ]2 F& o+ s  n" P, x+ `   
   keep if dx240<`tempmaxp'
   
0 l* H$ u' U9 d) o0 ]* J   keep if ddy240>0
1 Y, m; l3 V2 p( w/ O' m   count
               
" P) B7 I% j2 ?                if r(N)>2 {
                 
                 pkexamine dx240 ddy240
8 X, I6 {4 ^$ ~5 v1 f   
   local tempmaxphddy240=r(tomc)
" X7 b6 R- D8 D4 g   
* M" A. H) _( L4 @   restore
- p% a" q: V  E                 }
                 
        else {
               
+ F2 I3 y* Q  t                 di "cannot find pmaxphddy240"
  \% w0 X( z% c0 ?) ?( c! ^                 restore
               
                 sum dx240
# A6 y" ~$ g2 `; M3 Z                 local tempmaxphddy240=r(mean)
) }2 v- z! J6 C2 O+ u+ J0 o. k7 Y                 
                 
* y& Z9 s$ f& _' X/ p3 L                 }
0 C7 ~- E; S$ U% G   
4 k0 w' c7 M9 i* }. ~   
: m5 ~  c* Y" j1 B+ l/ t6 x   
   
* P6 Y; B3 `6 X- f  r* o   
   
+ R8 l8 n/ c) k6 K        gen pmaxphddy240=`tempmaxphddy240' in 1
' m0 p% ^: n6 f7 J       
4 X# e/ C% j1 j5 x/ k7 ], e- Q9 b       
        *-------------ddy240pmaxp右侧低峰值 -------------------
2 a7 A) d% O/ G0 u       
& q. T" A3 I: B2 Z- h0 y! q8 ?         preserve
3 B6 [4 {: Y; c$ ~7 F$ i% o  
   keep if dx240>`tempmaxp'
8 a0 y% O. B5 l$ p2 C   
+ w" s+ z6 b0 B! Z5 @   gen temddy240=-ddy240
   
   keep if temddy240>0
   
   count
& o8 h0 ?. _) H7 o7 z. p. |               
- Y' N( \7 f3 q# u                if r(N)>2 {
                 
                pkexamine dx240 temddy240
   
   local tempmaxplddy240=r(tomc)
   restore
8 R+ A! s. H/ t% ^2 [                 }
( c+ ~& A, M+ U( H5 a0 Q. a$ V                 
        else {
  y' e. }( ?. g, z& p                restore
                 di "cannot find pmaxlhddy240"
                 
               
. }+ `* `  \! l1 j5 l3 X- y                           sum dx240
" \$ B5 O* y4 v8 S& F& ?                 local tempmaxplddy240=r(mean) in 1
               
                 
                 }
   
   
   
   
5 W8 `' x% ?4 P7 t! R# v6 B   
" A: f) ]# }! V; P2 P. E   
0 U' H4 b/ A, E# ~0 N, K5 C        gen pmaxplddy240=`tempmaxphddy240' in 1
- ]" U. u+ D, }3 W/ X: f       
        *-------------------------------------
   
; l6 S, P1 k+ T# M) E0 I8 L) y7 U   gen price5d=v5 in 1/20
      
8 J8 F0 ]/ T, I   gen price5yy=_n in 1/20
" v; M( d( P: J" r0 n7 d   
   gen price5y=price5yy*(dy240max-dy240min)/20
+ u6 j. j% `$ [1 H- a- v   
   *price5y &price5d
   
, k8 |. p* F! u8 C0 a   gen priceny=price5y in 1
   
   gen pricen=v5 in 1

0 K" ~7 B0 F- r, A( j- ]$ d     
   *priceny &pricen
7 v$ D, w3 k/ ]: F1 ?   
/ h5 M% d3 c8 ]% w   
) [) e3 l  t  r+ x* m   *-------------------------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)

* ---------可选项目----------------------
- _2 U$ m/ [: @2 s* T' M4 Y*(spike yxis amean,lc(blue))(scatter yxis2 amean, lc(blue) msymbol(+) mlabel(amean))
# {* A: T: b/ p1 v4 \# S  D' y}
end
8 c6 z  M# @7 }3 \

maodeng

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

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

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

喜欢喝冰茶

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

喜欢喝冰茶

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

- `! ]  v7 u- f! t
这东西是开源的？知道有个统计软件叫stata，不过要收钱。

汉水云潇

准备试一试
先找Stata软件

汉水云潇

4 G9 m& [9 c! t0 U; C% \2 v
晕呀
论坛复制程序代码
乱码夹杂在其中
# I+ e! ?3 F5 Z$ S6 w6 Z很烦
对了两次
以为没问题了
6 M7 K/ o* d* Q4 V- H+ k7 V4 F
结果
, P2 x/ f. c- e% ]& V+ N执行后
还是出现
% h. {- H( Q0 Y5 I+ q  k; zr 199
5 i# `6 O0 h6 x5 Y) m的错误
0 i/ d) m6 h* M

汉水云潇

试验成功，下面是三一重工的图
! u; x8 @2 N5 r& R8 X
9 b# ?& P1 z7 U: K

peterbev

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

中秋不圆

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

bknite

这个帖子值得学习。

jeanny

感谢！！！

Hubert2000

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

4 M: ]2 q4 Q1 O) x7 c老哥，也是山寨这篇文章找不到了，不知为何？