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

maodeng

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

maodeng

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

maodeng

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

( y4 f2 H# Z/ k. t) m大师 就是这个

sergio

这贴实在！

awer

maodeng 发表于 2013-1-28 13:13
. hello 300236
statement out of context
r(119);

4 ~$ K2 I3 U) B4 J抱歉，最近太忙了

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

你把 D:\stock\test\source 改成你导出的数据存放的目录，再试一试
" V9 v- H: p3 d' H* c
capture program drop hello
! Q* l1 k! f1 F8 i+ }
" _8 T) Q/ O- I/ h! E/ iprogram hello

quietly {
0 _+ ^# Z) e8 [! y1 d* V& g
7 I/ J# ?* q- r1 B3 n- F. d. s* X# Aclear
" Z! ^- T: B% _* [/ R. G% `6 T
cd D:\stock\test\source
9 m1 @( }; G7 r9 ]" {
5 E" d3 ]7 j' @& P, s, m  Yinsheet using `1'.txt

   keep if v5>0
   
   gen vt=_n
   
7 B" S) R* H- E4 H7 f8 ^* l* U   tsset vt

" {' w& v& p0 m5 @" K   gen zf=D1.v5
& a- g9 n9 r; d+ D, {
   gsort -v1
+ I: k# q1 [/ x, \
( {; h9 u& c0 f5 e+ N   gen avol=v7/v5
# W; _! Q3 v- ~: O1 n% P2 `# p   
   keep if avol>0
   
   drop vt
8 P2 X1 @. @8 y8 F4 u/ u6 C2 R; q  
   gen vt=_n
   
# u  [& n$ [! u/ F7 m( F   tsset vt
   
    sum v5 [aw =avol]

   gen amean=r(mean) in 1
   
! j- X( l+ t' f3 Z   * 总加权平均值 amean
0 F. x3 V5 w$ k8 c2 ~: z   
& Y4 ]9 j* C' I1 v   
   keep in 1/240

% V0 V, `4 j' ~$ G& j  R. R4 V! a0 f   sum v5 [aw =avol]

   gen mean240=r(mean) in 1
   
! x( z- a/ e' n! Q( {- X0 J   *240加权平均值 mean240
   
2 V' W% o1 j  k0 l) |! w; ~   set obs 500
6 C4 C# u$ X8 {5 }8 g0 {   
3 ^' e& u6 M2 O4 s. o   kdensity v5 [aw =avol],gaussian n(500) generate(dx240 dy240) nograph
   
   *240 kdensity(dx240 dy240)
: {" c- \; Z5 I9 I* K5 ^2 d: K4 \4 Z   
$ a9 u& I7 q8 [: _   sum dy240
: Z; z- s3 f. ]9 N9 x, j- q   
   gen yxis=r(max) in 1
   
& O, L: P- L1 |7 s; b  z4 {& {   gen yxis2=yxis/2
1 Q+ N/ d! i3 l   
   *yxis
9 |- ?1 H/ I$ _   
   gen dy240max=r(max)
- i  w7 O8 n: }: V$ c8 M7 C" h' W   
   gen dy240minp=r(min)
8 q- a+ E8 y8 n: r
    gen vtt=_n

9 J( y. f5 J% g) I% ^( z    tsset vtt

. Y: }5 w& ]9 ~9 Z% P4 s- f1 B. V6 A' j   gen ddy240=D1.dy240/D1.dx240
   
   sum ddy240
   
, H! g/ v. l. \9 k* |5 a/ x. S+ {   gen ddy240max=r(max)
. A. Y' k/ I6 w; z& \   
1 B' [9 w; m# }3 S' {5 e  c, t   gen ddy240min=r(min)
     
   gen ddy240g=ddy240*(dy240max-dy240min)/(ddy240max-ddy240min)
   
   sum ddy240g
   
   gen yxis3=r(min) in 1
   
   gen yxis4=0 in 1
  
  *ddy240g
   
+ B# X% `3 {1 G   
  preserve
  
% I) p* B( v7 g8 ~( h! X& V, }  sum v5 [aw =avol]
( ?8 Z9 A' I0 f4 i( T3 e8 u  
  keep if dx240<r(mean)
  
  sum ddy240
  
  if r(min)<0 {
# E. J' r' U3 J' m        
: W  H: l* T" q% G                keep if ddy240<=0
+ p: W1 E, y) Y5 S# `- v. N               
                count
6 ?" }0 }% A+ N7 s) [               
+ s. F& y/ P+ j" b; G                if r(N)>2 {
                 
. S# y; Q( I+ X+ z0 h8 l9 P                 pkexamine dx240 dy240
               
: ^/ E, b% B5 {& C) s                local tempminp=r(tomc)
; X. L$ T( Q, R4 R                restore
                 }
                 
        else {
7 j2 w1 d& b/ ?$ R& Y               
( {, K# H, p3 B- p8 c                 di "cannot find minp"
                 
$ K( ]2 X6 e- t. Y- X                restore
                 sum dx240
- A6 P8 G5 n6 {1 Z: h. G                 local tempminp=r(mean)
                 
                 }
                }
   else {
                 
                 keep if ddy240<=r(min)
7 n& F; e6 Q2 E  I" e                 
0 ]4 e7 a/ L$ I3 W0 U( ?                 sum dx240
. I& l. K: m+ @. ^$ t0 k; ?2 _                 
: ?% a$ h  t7 ^9 T                 local tempminp=r(min)
                 restore
4 X1 a1 I! o* ^                 }
+ J/ Z- g9 _1 m( r- z; H6 I
0 {- d3 A% d, ?* F( G  F
" T& T: v) U. k+ A" u1 `7 f# g
   gen minp=`tempminp' in 1
( a+ d+ j/ ^0 h) S   *240 minp


   
   preserve
0 V. {7 y& T- }   
   sum v5 [aw =avol]
   
   keep if dx240>r(mean)
   
, |% T: _5 {0 Z& _   sum ddy240
   
" F7 i' ^' p& }2 G/ x   if r(max)>0 {
/ U" m# {. A( P6 p, M" {2 n+ b3 U               
                keep if ddy240>=0
               
                count
               
+ }" o( Z0 N3 r0 l  d( K                if r(N)>2 {
0 f0 h; I" u' q5 \  L: h& M                 
                 pkexamine dx240 dy240
               
                local tempmaxp=r(tomc)
# v* B$ j3 ]1 W# d6 q                restore
                 }
                 
$ K. r  D4 {7 g+ P5 O7 X% H                 
3 X( @( ^: \& t5 M, D' t( u        else {
! L. K  s, C# B0 ^0 l% i                restore
                 di "cannot find maxp"
               
: |  d* v! c& q! C: D2 M4 }               
                 sum dx240
                 local tempmaxp=r(mean)
( @) y" ?5 @" `9 ]* W; q5 T                 
                 }
                 
                 }
( d$ ~! Q0 s8 ^, G" F4 j( u$ d                 
5 `7 o1 n& s) c; ?5 p! m/ f    else {
( }$ v. u3 p9 O9 J9 L       
           keep if ddy240>=r(max)
          
           sum dx240
" M# ^& I9 ^7 Z6 ^8 n: X! D          
+ y5 v" T5 p) w" H- f           local tempmaxp=r(max)
* u$ |6 r2 F- j2 g           restore
8 t  z. \; o, `6 t0 C* z: l, Z4 y           }
       
! l6 _5 q* \' W- Z' t" t
  
2 x. l+ i- g0 ]& m1 c
$ S% n4 }" A1 \   gen maxp=`tempmaxp' in 1
* ?& n6 s1 o3 \" z9 R   
4 S, Y# l1 w. t  {/ c   *240maxp
   
; h& R+ K; _* i) r   
* E3 ~* L5 S- C$ q( S   *-------------ddy240pminp左侧高峰值-----
   
3 J9 n+ w5 e. F3 ^; A5 B. ~    preserve
   
   keep if dx240<`tempminp'
# D4 `3 Q+ P) k. k   
5 h8 t4 L2 O& A$ ]   keep if ddy240>0
   
   count
' h+ G- O# x# e; T/ Y               
                if r(N)>2 {
3 u$ ^3 [4 Z6 ?- S5 H% Z5 n                 
                pkexamine dx240 ddy240
               
                local tempminphddy240=r(tomc)
                restore
                 }
; H, |8 q1 r& s- ]/ Z                 
        else {
4 S  g% p+ R" `* f                restore
( D  U4 E. o" }  C; }4 k                 di "cannot find pminphddy240"
) P/ V% ]2 p) |- Z8 N                 
8 `+ o/ W/ k' V* k; |! E               
  a; g; K0 w( P                 sum dx240
1 Q$ m7 |3 W1 u+ U/ M                 local tempminphddy240=r(mean)
               
' [! _4 m! c& T1 r                 
                 }
                 
0 h/ b2 Q) u0 D* n3 A   
, {+ N# J3 B  w( a8 u' V, G   
  w+ w: C5 j& {9 U/ L- d( g7 I" z   
8 V( E9 ?, B* L* B$ e   
   
# I, l4 F  i0 \$ ]; [4 O* O) I7 |        gen pminphddy240=`tempminphddy240' in 1
* j/ X! z9 U2 i4 u9 J6 `       
! o" B0 v/ ~/ p6 M       
' V* @9 x/ U% A9 w9 J" z        *-------------ddy240pminp右侧低峰值 -----------------
/ X! j# U: A& |' L    preserve

   keep if dx240>`tempminp'
   
   keep if dx240<`tempmaxp'
  o: K0 s# G( M# p7 g0 W% E+ C: F   
0 L1 q6 i; B7 a5 k+ c# H- Q   keep if ddy240<0
   
- ^5 b+ i' U9 Q# [7 u8 I   gen temddy240=-ddy240
0 x& @/ j& A" O   
' l  ~# {" X" \& g7 ^. `% a    count
               
                if r(N)>2 {
                 
& Y! i" L$ N% R4 B0 T                 pkexamine dx240 temddy240
   
! Y, i- y! D( B2 \; S  local tempminplddy240=r(tomc)
                restore
                 }
: C4 V0 W4 \8 i5 x7 Z                 
" Z, U$ ~* F- a4 d) d3 u  [        else {
               
                 di "cannot find pminplddy240"
                restore
1 E# M, G/ u7 M                 
, B: h3 \0 a* ?6 }% \: b                  sum dx240
                 local tempminplddy240=r(mean)
                 
                 
                 }
/ f+ }0 V" S! c5 @   
   
! F7 g9 C  Q3 P* b2 W   
   
   
( e5 v# [' p/ ~( l        gen pminplddy240=`tempminplddy240' in 1
* W- x+ M. k1 ]$ _- m$ F2 W& h       
. r) X3 d9 w4 d& u) {        *-------------ddy240pmaxp左侧高峰值 -------------------
' V+ O, q8 f2 O+ d       
/ {. g* e) J9 C8 v         preserve
2 D/ E6 ~, b* o   
$ O5 N  J+ o2 W5 _! ^! W   keep if dx240>`tempminp'
, L3 m- ~3 M# [# O0 _4 L   
   keep if dx240<`tempmaxp'
+ ]& c+ ~5 E) B   
* h* W4 W% P) W0 }" H   keep if ddy240>0
   count
               
                if r(N)>2 {
                 
; H) R1 X; R& d9 `                 pkexamine dx240 ddy240
' F" q( r( W8 @2 W6 ]0 P   
   local tempmaxphddy240=r(tomc)
   
   restore
                 }
                 
        else {
               
                 di "cannot find pmaxphddy240"
                 restore
               
                 sum dx240
                 local tempmaxphddy240=r(mean)
                 
                 
* T. A% z, r% w, O                 }
   
   
   
% C/ U+ L3 R& @) t% j# }5 I   
   
   
        gen pmaxphddy240=`tempmaxphddy240' in 1
       
, R( R; S4 A" H; ?- h* E* U       
; M5 o) V' T! X: G: x# H" n1 A* Q        *-------------ddy240pmaxp右侧低峰值 -------------------
& q& w+ ]; c7 C0 u+ a) Q+ U9 k  k5 Q       
1 g! M) N2 X" S0 H         preserve
  
$ Y1 d% ?8 i. {2 s$ @! U   keep if dx240>`tempmaxp'
  Y& j, C: f# n7 p. t! w; A   
- v3 ~) F8 Y% T6 l) ^% C) g9 Q   gen temddy240=-ddy240
1 |/ h/ y5 t7 K( O# U   
) B: H2 {3 ?! L   keep if temddy240>0
   
   count
$ Q8 L) x; v; g9 G! v% U               
                if r(N)>2 {
                 
; X6 N& P; l* F6 k9 K                pkexamine dx240 temddy240
9 e0 c: P5 h3 h6 R   
   local tempmaxplddy240=r(tomc)
2 b, U2 Q( O0 B; R! l   restore
                 }
7 f% Y; P" T( d$ Y4 m7 S0 A                 
        else {
                restore
                 di "cannot find pmaxlhddy240"
' Y) L' y! F; K3 K2 N                 
               
* o: |5 f! X, Z& i. }                           sum dx240
                 local tempmaxplddy240=r(mean) in 1
               
$ \% j. U3 v5 @3 W* s                 
; e, D0 `7 [% I8 ^2 o                 }
9 b5 G& I+ ^" _% ~   
   
   
   
, T1 k9 K* y7 }7 v% T( s   
* j4 C9 K. w/ R1 m% p  p- j1 W   
: B# J4 @) H7 [        gen pmaxplddy240=`tempmaxphddy240' in 1
       
        *-------------------------------------
$ U5 [9 b4 N# Z9 X   
' |4 N1 N% V/ i) B- \   gen price5d=v5 in 1/20
      
' N2 Q4 n6 J' {, C0 a7 }   gen price5yy=_n in 1/20
   
   gen price5y=price5yy*(dy240max-dy240min)/20
+ m8 \# D' [1 k, |2 V   
   *price5y &price5d
& m* |8 g- a1 J5 t( {   
% {$ i; ~$ @4 n8 @( X" `' \0 _' e   gen priceny=price5y in 1
6 N3 H  O# u9 h  A   
+ C! I: L% S, ^# \' d   gen pricen=v5 in 1
7 s) j1 [, X6 u# s4 {
' y& r9 n. u: h' ]& H7 W; A     
   *priceny &pricen
   
3 d# |; l8 n. l  W- Y, N   
% B; X8 t5 K* n. l* u" S! M/ a   *-------------------------graph------------------
& L2 \8 q4 S! n6 U  ^+ O2 f  
: |7 c9 g" S- _  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)

* ---------可选项目----------------------
" a) D. ], \0 y7 k. x6 }*(spike yxis amean,lc(blue))(scatter yxis2 amean, lc(blue) msymbol(+) mlabel(amean))
}
' L3 w0 y  L/ {: \) k0 L/ Oend

maodeng

awer 发表于 2013-2-5 14:40
抱歉，最近太忙了
  b' S/ r6 d5 q+ @  K4 `
300236 我用下面的程序试过，没问题

哇，谢谢，我仔细看看，这些字母数字代码晕了
4 o# s. J; ]* v5 P+ l6 C+ l

喜欢喝冰茶

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

喜欢喝冰茶

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

: J4 U) E5 B9 h* l. g; G  Q/ F7 }
这东西是开源的？知道有个统计软件叫stata，不过要收钱。

汉水云潇

准备试一试
- G& x) f8 V) @5 J8 y先找Stata软件

汉水云潇

# S( X2 c: S5 J( F7 G2 R
" ^% O) U4 f# s0 x. p: e晕呀
论坛复制程序代码
乱码夹杂在其中
很烦
3 m1 t9 @& C- A: S. o! g  o- w4 n! S1 R对了两次
) ?5 z! ]" U* F4 }; w' v以为没问题了
$ v* x. ]* u0 W
2 P2 ~) s- _8 [结果
执行后
& b) H7 `) g) U3 r1 x还是出现
r 199
. [/ M7 I4 N$ i2 e的错误

汉水云潇

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

peterbev

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

中秋不圆

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

bknite

这个帖子值得学习。

jeanny

感谢！！！

Hubert2000

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

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