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

maodeng

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

maodeng

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

maodeng

. hello 300236
4 F- X# [' I# o4 Ystatement out of context
$ o- w, `( _* V1 f5 Mr(119);

* b. l" p7 o1 y% ?0 h# g' Q9 m  k大师 就是这个

sergio

这贴实在！

awer

maodeng 发表于 2013-1-28 13:13
. hello 300236
8 q9 V2 ^; {3 A( F" a- b6 L8 ]  Z' sstatement out of context
2 M/ Z, A" N+ q/ lr(119);

5 @( u; k8 J- Z2 x& R1 h抱歉，最近太忙了

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

" p# E# A' d, N& ^/ E你把 D:\stock\test\source 改成你导出的数据存放的目录，再试一试

capture program drop hello

program hello

4 `1 h+ c# {) i$ a/ d2 W quietly {

* \# b0 J( w- n# [clear

cd D:\stock\test\source
- a- b; J' U. I
insheet using `1'.txt

   keep if v5>0
- J) E1 s% b2 h2 ?1 U   
- i" U4 S3 k$ v. Y$ @# @3 S2 J   gen vt=_n
  \7 x" |+ U& @5 Y$ p+ S2 x! h, k   
   tsset vt

   gen zf=D1.v5

   gsort -v1
# J, V1 ?+ Z: v/ K
& Z# U' M3 b% e, T6 W1 ]  n$ o   gen avol=v7/v5
) s* i. }6 p1 o' o   
, T2 r  w+ @7 M   keep if avol>0
9 \% e/ O, E2 e; j) q3 v   
   drop vt
  
   gen vt=_n
) w) `2 U# r+ V0 ~   
' q6 m; t# i* C7 O2 E6 S   tsset vt
   
0 N) Z2 W! |0 ?! \( ]. L, C    sum v5 [aw =avol]

' y7 C& w: l1 U% D' |' L3 A   gen amean=r(mean) in 1
   
   * 总加权平均值 amean
   
8 f: o  ?7 Z7 t- y: H( k   
   keep in 1/240
% i( l) E4 \3 l* z1 g- v; B9 d& U
   sum v5 [aw =avol]

   gen mean240=r(mean) in 1
5 T' n/ v1 Z; ?; \" z3 J2 w   
   *240加权平均值 mean240
   
   set obs 500
   
3 T6 t3 p8 M  y8 z9 `   kdensity v5 [aw =avol],gaussian n(500) generate(dx240 dy240) nograph
   
   *240 kdensity(dx240 dy240)
, H  U9 Y8 l6 I7 ?   
  r2 z, N) ?# \' t   sum dy240
   
   gen yxis=r(max) in 1
   
   gen yxis2=yxis/2
   
' f# ~1 a! u: U% ?   *yxis
   
   gen dy240max=r(max)
   
% t/ u3 k5 Y4 A4 C   gen dy240minp=r(min)
; q) T/ H* y' F+ X+ G
9 D" n9 f8 x/ b) |) l    gen vtt=_n
6 X; e1 v- r% a& l
2 V6 h  j+ q( L3 S# t/ x( D5 a* X    tsset vtt

3 M3 C0 z) d2 |* J   gen ddy240=D1.dy240/D1.dx240
( M$ v+ B/ ^& p' c) r   
   sum ddy240
   
   gen ddy240max=r(max)
/ n" }- ?# W' _8 d   
   gen ddy240min=r(min)
     
0 j5 [, D- o% }# U  B8 s% N- h, K1 |   gen ddy240g=ddy240*(dy240max-dy240min)/(ddy240max-ddy240min)
# G2 x" n9 f+ p& d6 h3 v4 O& _   
8 {2 \2 \, E; q  T0 t" o: U   sum ddy240g
   
) V- c+ R5 x# y6 i0 A   gen yxis3=r(min) in 1
   
   gen yxis4=0 in 1
: t, T) {0 ]0 o1 m  
  *ddy240g
+ l$ Z$ y" E4 ^/ I; _* C7 i: G% p   
   
  preserve
  
  sum v5 [aw =avol]
  
" P9 s4 A0 s7 J% o# `  keep if dx240<r(mean)
1 ?0 t/ ]/ k  J  r: n  
% F/ K! n4 Q" }  sum ddy240
  
  if r(min)<0 {
        
                keep if ddy240<=0
               
1 u. a3 ~; x* b9 a                count
               
9 d  Q9 E8 D1 P5 |                if r(N)>2 {
                 
7 a( n' a/ {7 v( [                 pkexamine dx240 dy240
               
                local tempminp=r(tomc)
                restore
. r/ l0 V/ u; X, X: U& ]                 }
                 
7 A: @& H  c* U5 m5 @$ X8 B9 w        else {
               
4 q0 C/ B) e$ A                 di "cannot find minp"
: v3 f: ~) D4 a5 ~. j                 
! O, g+ \3 B. F& m2 u% o- w                restore
9 U& }6 R" t1 Q7 P1 B7 H& x                 sum dx240
" Z) ~/ d6 g1 A. L. ?, C5 |$ |                 local tempminp=r(mean)
                 
                 }
                }
+ N0 {  Y% \& Q  G9 c   else {
                 
                 keep if ddy240<=r(min)
                 
$ d! u0 W5 k& Q$ N                 sum dx240
& p* C0 J! Q$ l                 
                 local tempminp=r(min)
                 restore
! d+ l- D  {8 z% A1 ^( f4 |                 }
& j! m4 Q2 @% ^  M
# i4 N/ [! h9 C/ p6 N
! n0 B1 K5 F+ N$ x+ c3 ~7 F2 ?4 S
0 w" T8 A4 c' U. ]   gen minp=`tempminp' in 1
   *240 minp

4 A" q3 e. }4 T
   
& Q2 c( D3 y" \1 F# m' C   preserve
   
   sum v5 [aw =avol]
# T5 |# K8 z% z/ K  E- ^! X3 `$ Z5 ?   
   keep if dx240>r(mean)
* N: G: z3 b- W& {9 k   
4 v, }" l0 c" @+ v3 E   sum ddy240
   
   if r(max)>0 {
! }/ X# c5 k* ^3 `               
3 h, B( M' w! r/ g% M* ]                keep if ddy240>=0
               
" [+ r5 @+ U) P$ C                count
, g. d  ?$ m7 N9 J9 p               
. U- Z# B; K5 v# b5 V8 G! Z5 K                if r(N)>2 {
1 p! V/ @3 ?8 x9 f  Q2 S: B2 S                 
                 pkexamine dx240 dy240
               
                local tempmaxp=r(tomc)
+ S# r/ K: Z5 N6 |                restore
                 }
2 F- H, u5 x1 W                 
                 
        else {
: P5 L. v$ X: p                restore
8 B/ X3 _4 H( ]; Q2 L+ F3 s. d                 di "cannot find maxp"
4 Y$ G+ a& P! u0 K! \  u. V2 O: a               
               
                 sum dx240
4 g& _( ?. U: C$ }                 local tempmaxp=r(mean)
                 
                 }
                 
                 }
                 
/ E5 K  \" ~( D: M9 F    else {
       
           keep if ddy240>=r(max)
          
5 w. i* x  R* b; P$ \! S           sum dx240
          
: U" L" J6 D: v# V' \           local tempmaxp=r(max)
5 A7 W1 L! i2 R( G& L/ B$ f+ P% K7 q- y% f           restore
           }
- E/ U: ~  r5 ^       

" G7 F  q9 n2 @- U0 B' K  
; |2 V' R0 m, ]* _+ h6 a4 P! }
   gen maxp=`tempmaxp' in 1
: Y8 q4 q7 Q' |, s/ l   
   *240maxp
   
   
   *-------------ddy240pminp左侧高峰值-----
   
    preserve
+ F) a' z* T8 k6 O! }   
   keep if dx240<`tempminp'
   
   keep if ddy240>0
   
   count
               
, r7 Q$ l$ K" L                if r(N)>2 {
- Z; v4 S) q8 ]% r2 Z6 I* L                 
                pkexamine dx240 ddy240
1 J9 }. ^- K+ _6 H/ |! U' V               
                local tempminphddy240=r(tomc)
0 C% V9 n. x; n* l2 b3 T1 o" d7 b                restore
                 }
                 
        else {
                restore
( k: P" j% ~' D; J* R8 Z/ X8 v- x                 di "cannot find pminphddy240"
$ S# }1 L, `. x# F% N2 R6 ^6 ~" E0 l                 
               
# @) b% |0 m4 U# w  |! {* c                 sum dx240
                 local tempminphddy240=r(mean)
               
                 
- A; k1 F. V) {1 h, \9 t; r                 }
                 
   
   
   
   
   
        gen pminphddy240=`tempminphddy240' in 1
       
       
3 |7 j& Z) p1 i) G4 l, L        *-------------ddy240pminp右侧低峰值 -----------------
. ~; x# B9 M* T+ d% q3 z1 x* F    preserve
' w  z5 _" Z* B# L( }+ m" H( B$ |4 w
   keep if dx240>`tempminp'
   
   keep if dx240<`tempmaxp'
; }* F: g7 w( @3 j, Y$ X2 y   
9 \2 H% P# D& u   keep if ddy240<0
   
   gen temddy240=-ddy240
9 q9 K. D9 g( r   
+ X, s/ x8 M+ C( s    count
               
+ S* ]# E! [" l2 ?3 d+ Q2 Y                if r(N)>2 {
. p0 }# O0 U+ C4 @% f- P0 d$ E, ^                 
0 _! g/ b& C6 l1 R                 pkexamine dx240 temddy240
   
  local tempminplddy240=r(tomc)
                restore
                 }
                 
1 L+ K1 Q0 K; ]9 t1 a  [3 K        else {
' O( `( }0 }# i5 {               
                 di "cannot find pminplddy240"
3 a2 W: Q# ~) m5 J* P' N0 I                restore
9 X2 N1 z! J& D7 b0 |                 
                  sum dx240
                 local tempminplddy240=r(mean)
                 
& f' Q( W  X$ F3 |3 m5 Q                 
# w+ n/ t& {+ w+ y                 }
) ^/ y  \1 H" T) H. }8 b   
3 l* }8 q2 b$ P  C+ M   
( l7 @% f+ r6 Y# ]) y   
   
   
5 L# o4 M+ T6 y/ A9 p  X* K  a$ x        gen pminplddy240=`tempminplddy240' in 1
       
        *-------------ddy240pmaxp左侧高峰值 -------------------
2 A- ^3 g# l( U- S7 ]       
3 K3 l3 s( v! I# Q3 t         preserve
# X. \0 X$ t* S   
0 O: U; n! \( ?7 y   keep if dx240>`tempminp'
   
   keep if dx240<`tempmaxp'
2 U$ K4 z7 i, I0 R$ o- g   
   keep if ddy240>0
   count
' \  K- z4 P, z               
                if r(N)>2 {
/ B: G5 C: `/ k3 n  `% l                 
; \4 w6 _0 a" l' i6 o9 r8 U( s0 ?                 pkexamine dx240 ddy240
8 P9 W, S- S5 G7 P) z   
% N6 A8 ]  P- S! h: {! \4 w; w   local tempmaxphddy240=r(tomc)
   
   restore
4 u0 z4 r: O) T  N% U9 ]                 }
1 K5 x1 ^/ N5 C0 H" Q1 |                 
$ L3 H# b) R  C        else {
9 g. c8 a/ M9 s               
0 A2 N$ b3 k  P2 r+ D                 di "cannot find pmaxphddy240"
                 restore
# q$ c7 N1 o; A               
                 sum dx240
2 C8 N5 o, x8 `( T$ g) n+ b& x+ d                 local tempmaxphddy240=r(mean)
9 C$ H# f7 {& V: t+ h' a- |/ J                 
                 
( _* ^/ F/ |( ]6 K                 }
  [, |3 t; q: w2 S; Z4 q8 c   
   
   
   
7 V3 ?( ?. o) r. d( Z   
   
        gen pmaxphddy240=`tempmaxphddy240' in 1
' y6 O& ?4 L  [       
       
        *-------------ddy240pmaxp右侧低峰值 -------------------
       
1 y# {! h) J7 r  n+ c) p. G: \         preserve
  
: J# y3 d. ^9 J( y5 y0 P6 Y: t   keep if dx240>`tempmaxp'
" N! [* w: i( g7 k$ d   
   gen temddy240=-ddy240
   
   keep if temddy240>0
) T, R9 K- H3 Y1 r) @2 }) b   
' H6 n/ l& A, l) D   count
' Z( l! \8 d% a( y" W               
                if r(N)>2 {
                 
! n% h: m' [4 ]) p. ^, ^                pkexamine dx240 temddy240
  T/ b: Y8 U6 H: F' N   
   local tempmaxplddy240=r(tomc)
. ]: A) ^& `& p$ ~& t   restore
                 }
* x. G' o/ z3 A, ]: D7 i, K                 
        else {
                restore
. R- d4 P1 o6 T: n) q                 di "cannot find pmaxlhddy240"
& c& n  s% s3 q* x+ H! i4 g& o                 
+ l( E7 G, P/ u. ^# N               
                           sum dx240
                 local tempmaxplddy240=r(mean) in 1
. I4 S" o9 G; y8 _3 ]2 M; a               
( i0 v1 N) v2 b% ^2 M$ \+ Q                 
                 }
; V6 F: _' B) ~% q+ E3 U+ S  N   
+ N- X& ~; R! g0 v# `" R   
   
* C/ W! q4 n1 t& U   
3 E; Q. Y: i+ O  J) \' p# n   
   
        gen pmaxplddy240=`tempmaxphddy240' in 1
       
        *-------------------------------------
   
   gen price5d=v5 in 1/20
$ z5 Y6 U( W+ T- e3 q# `      
- Z2 d1 F% j: g' G" N$ v' Q   gen price5yy=_n in 1/20
   
2 s# [0 J8 W4 [8 B! n! t   gen price5y=price5yy*(dy240max-dy240min)/20
; U1 q( E* G- d2 d5 V/ {  n   
! o8 `! s9 A) S4 A   *price5y &price5d
" k. |0 M6 u8 Q2 D+ u   
4 u+ _2 Z* L2 ]   gen priceny=price5y in 1
9 C9 p' ?' I2 z# F( @: j- s   
- ?- t' M  f) ?8 }' k  h  h   gen pricen=v5 in 1

     
   *priceny &pricen
2 k: [8 G* _9 v0 ~   
   
   *-------------------------graph------------------
7 U/ A$ e8 z8 _3 X5 w7 J* d1 ]9 ~  
" m- Y" d+ |3 ?. ?; A  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)

* ---------可选项目----------------------
+ j4 l( t- S0 B8 @*(spike yxis amean,lc(blue))(scatter yxis2 amean, lc(blue) msymbol(+) mlabel(amean))
}
end

maodeng

awer 发表于 2013-2-5 14:40
3 P- ?( q: ?9 `5 F9 }  i% C! Y抱歉，最近太忙了

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

哇，谢谢，我仔细看看，这些字母数字代码晕了
, r; v- w1 e. I2 ^% F1 c: r

喜欢喝冰茶

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

喜欢喝冰茶

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

这东西是开源的？知道有个统计软件叫stata，不过要收钱。

汉水云潇

准备试一试
. g0 Z; S2 T0 M* `7 B, q先找Stata软件

汉水云潇

' b7 R0 s# I6 _晕呀
论坛复制程序代码
4 _/ q* L9 B4 w1 \乱码夹杂在其中
; k2 U) J0 K0 s8 N+ t0 z很烦
对了两次
+ q8 N* K# B8 X' p! f0 l" s3 C以为没问题了

结果
执行后
! W; d; \0 h, ?还是出现
5 N& W2 i' e7 X# fr 199
) `; M2 Q) J4 L) V( [' q的错误

汉水云潇

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

, Y+ F" M" F* v7 T4 d3 v: K( O2 p

peterbev

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

中秋不圆

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

bknite

这个帖子值得学习。

jeanny

感谢！！！

Hubert2000

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

- I! b1 K/ L% n' S1 X, H老哥，也是山寨这篇文章找不到了，不知为何？