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

maodeng

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

maodeng

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

maodeng

. hello 300236
statement out of context
5 H* I9 Y- \& @r(119);

大师 就是这个

sergio

这贴实在！

awer

maodeng 发表于 2013-1-28 13:13
, }* ?) e5 {6 t, L/ R, C$ A8 L. hello 300236
2 ?9 z/ S- N% Hstatement out of context
0 C9 V" F$ i6 s' V8 \# |r(119);

抱歉，最近太忙了

& n6 ^' `0 ~# ~+ m# t$ H1 g  J300236 我用下面的程序试过，没问题

你把 D:\stock\test\source 改成你导出的数据存放的目录，再试一试
6 k2 y- I' v0 S5 t9 W, W) `/ r
capture program drop hello

$ R  x9 {+ m0 B: rprogram hello

1 U# Z  j3 Q( N& |4 P) @; b0 l quietly {
) S1 f- o  W$ H7 }, X2 W. o& i
clear

1 f: X, l! @# q8 J. a9 p8 t1 G) y  q! Dcd D:\stock\test\source
" F! }/ T2 y: x9 [! D' P0 c
insheet using `1'.txt
( X- Y' e/ C6 L: q. `
  D  v1 O% o! l5 A" }# }   keep if v5>0
2 J9 O4 A0 u# `: ?1 P* b$ Y5 W   
2 L; J1 P: ?( z: S" I% b4 [3 h   gen vt=_n
   
   tsset vt

3 Z. R" `4 V: N, b  U   gen zf=D1.v5

7 y+ Z4 K4 N, j* X( p/ z$ P0 Y$ A   gsort -v1
* }1 `* l' d2 e
   gen avol=v7/v5
   
4 J) W0 g  z) R   keep if avol>0
   
   drop vt
& O+ d2 _8 L8 J+ C% i" z9 j. M& B8 C0 z$ ?  
9 T& U1 |* [; Y5 _7 w3 p' W   gen vt=_n
# T1 ?3 F+ v$ m4 d+ S+ n$ R   
* ?2 i% f6 d5 W% E) p  a   tsset vt
7 t" d5 v% z& I# z* Q8 e+ F! W3 D   
0 S4 [# B4 S/ _0 k  i    sum v5 [aw =avol]
; E/ i1 H: B' S+ m8 g3 j
   gen amean=r(mean) in 1
   
   * 总加权平均值 amean
   
   
9 Y2 j# p! b+ }8 X( f) M   keep in 1/240
# F- N3 m% ?& B$ q1 R& l
   sum v5 [aw =avol]
& h& w5 w8 V  t) b+ w: f
   gen mean240=r(mean) in 1
, C  n7 j4 x; i% R2 B   
1 e7 }' u4 k, J1 w9 N   *240加权平均值 mean240
   
6 n- x0 C3 L$ n* @* r" A$ {* a9 z) b& ]* D   set obs 500
2 ?. M! l( C# O- z0 d   
   kdensity v5 [aw =avol],gaussian n(500) generate(dx240 dy240) nograph
0 R7 Q* J9 b: c   
0 L0 l. \0 T$ a8 w! Y   *240 kdensity(dx240 dy240)
   
   sum dy240
; A9 y: C+ f  X! @& v   
9 q8 L+ u: D) G! e! }7 _, Z   gen yxis=r(max) in 1
   
   gen yxis2=yxis/2
   
   *yxis
: T; T7 _7 M$ ^2 z1 H% A% `$ |   
& @8 P0 k# v( F. ~5 G/ N   gen dy240max=r(max)
   
. }; ^  H, G: f' H   gen dy240minp=r(min)

, \- t! L7 m1 W$ B) X    gen vtt=_n
/ O% I" @1 O- z/ I; T6 K
5 g% O; b: B/ v! X# j1 M+ a    tsset vtt

: V) N7 ?$ W$ y6 M4 x& E5 ~   gen ddy240=D1.dy240/D1.dx240
   
   sum ddy240
% |7 M! O6 F' f3 ^8 _   
# P  k+ x1 N7 K" R# N8 ~5 A   gen ddy240max=r(max)
   
   gen ddy240min=r(min)
     
( M, x7 |' J1 Q( P$ @! Y9 T# ~   gen ddy240g=ddy240*(dy240max-dy240min)/(ddy240max-ddy240min)
- G* c" m: u3 Z   
+ S+ g% K6 h( _2 `/ p   sum ddy240g
" h- Z; F1 M; n3 K; j6 z   
   gen yxis3=r(min) in 1
   
# j3 F  G) w3 `" j# \/ B   gen yxis4=0 in 1
  
  *ddy240g
% {+ N7 z* u  e" b+ ?3 S   
2 S' ~, E; T' P' d6 k   
  preserve
  
  sum v5 [aw =avol]
  
# }9 v& B/ g7 J  keep if dx240<r(mean)
5 c6 B; z+ n& ^  
  sum ddy240
6 ]" P$ W1 J, `! P. z  
0 G. ]  D5 e; a. _) j1 O7 x/ T  if r(min)<0 {
        
                keep if ddy240<=0
8 U  t6 ?; w- W               
: R% I& Y( J/ P                count
) ]* I/ Q) Y: b+ f& `               
* l: \9 B3 N( @6 ]# G% Y                if r(N)>2 {
' b. s" P( \  R, W. T, u                 
$ Q7 \! w, {4 V% z5 O- u                 pkexamine dx240 dy240
               
                local tempminp=r(tomc)
                restore
+ j  x) V: N8 c                 }
                 
        else {
               
( @$ x' ^0 M7 O; B4 E# q' Z* D& V- H                 di "cannot find minp"
+ z9 O& h% z: X9 Y1 F* A9 `                 
                restore
                 sum dx240
                 local tempminp=r(mean)
( v8 U/ ~* G5 f& O$ E  }$ ?                 
                 }
                }
* W" a$ C7 i$ j+ u# {, ]9 \, p- h   else {
; S2 U9 W. O# u# W) y/ [                 
* p- _; {8 A4 T4 q0 M                 keep if ddy240<=r(min)
                 
                 sum dx240
9 t0 B! a# _& U1 J                 
                 local tempminp=r(min)
& q* i) P& {1 C9 q3 {                 restore
                 }



   gen minp=`tempminp' in 1
* A: O+ d. [' z( ~' e5 n   *240 minp


: B# x9 k/ }: i* F* w4 G, {   
   preserve
   
" p: C7 B2 w) w! d+ \( J   sum v5 [aw =avol]
4 z7 \0 Q! n) g0 i# m7 X! J   
% R- k0 _& I& b6 ~: i   keep if dx240>r(mean)
, ?% t: a& K6 h8 `( I# \! ~) D- I   
   sum ddy240
( ^- ?! F0 Q4 l" M" R8 ~   
   if r(max)>0 {
$ i5 U4 B9 |" u" A; D1 B               
& a7 {  A( \5 x9 k% k                keep if ddy240>=0
% H  Q: v; t4 g: f& ~               
, L6 m! X- n, {; f( P4 Z+ o                count
) }' [* }, s" u) Z( ~               
4 p. S$ F$ X" V4 A5 s                if r(N)>2 {
                 
                 pkexamine dx240 dy240
/ V0 K2 _/ ]; Z  a0 m; Q& o) k- i, K7 g+ q) o               
                local tempmaxp=r(tomc)
5 V5 p1 k! J0 m" C; A4 _+ Z( g; M                restore
                 }
                 
                 
        else {
0 [7 x% S7 e8 h" Q: U                restore
                 di "cannot find maxp"
( t5 X2 ]. l; e% l) o/ j! R3 X               
               
                 sum dx240
                 local tempmaxp=r(mean)
, ^* q+ _- `9 U% @$ R" b                 
( c7 T# ~1 C4 q% ^4 s; e" v/ @* O                 }
                 
                 }
$ N9 H5 e  h" j# j- B3 ?                 
( |& C+ G, a& ]% f  `    else {
       
) a. u9 {6 V1 ]2 c# G1 f5 q& U           keep if ddy240>=r(max)
) b7 s& |+ ~7 J( F- L" h          
           sum dx240
% G% c4 t6 T! m/ n          
: c! {" P* P1 W5 l: C           local tempmaxp=r(max)
           restore
           }
1 F1 i. V2 x) e2 c- I       
' O& \+ Y3 \' g2 H
  

8 P$ @$ P4 a5 l) P8 n# Z- _   gen maxp=`tempmaxp' in 1
   
2 M9 K5 G3 u9 b: N8 O   *240maxp
   
   
$ K; K3 O. c" D( _/ D   *-------------ddy240pminp左侧高峰值-----
   
    preserve
   
1 K! c. g# Y0 ~; j; G# f   keep if dx240<`tempminp'
   
   keep if ddy240>0
" I1 d$ Q% B- \2 I& X# H   
   count
               
                if r(N)>2 {
( h+ U9 t' c7 X( ^                 
. E5 K7 a7 a2 \$ l; R9 F  X                pkexamine dx240 ddy240
               
# G) `" d3 Z) `8 ]                local tempminphddy240=r(tomc)
                restore
                 }
" V; u6 Z/ T. b" a                 
        else {
                restore
& v' q$ S0 D7 n3 }  i                 di "cannot find pminphddy240"
& K; _" J+ ?  Q; l! ?: v                 
, R, _1 M' ~9 Q               
% c" M* y2 x; r                 sum dx240
                 local tempminphddy240=r(mean)
7 l! a. E3 c6 g* V9 h               
                 
2 M, s5 I4 ^+ ~, r: M- l                 }
9 i( M% r. P7 o. z: A                 
" e. M* K0 |5 L3 I   
) z1 n# A/ {' B& s2 M; L5 }, t   
2 v0 P8 F$ o0 v' G% }   
: ]" G( X5 }$ M3 o& y! q- w. U- k   
2 \2 ?6 d6 ]! S* m0 l% d' f   
        gen pminphddy240=`tempminphddy240' in 1
       
       
        *-------------ddy240pminp右侧低峰值 -----------------
    preserve

7 {7 J( p* p. ^( ]/ ^& H6 {+ M# c3 {   keep if dx240>`tempminp'
7 h+ d+ \6 H# |/ Q6 m" k6 Y' h   
   keep if dx240<`tempmaxp'
   
   keep if ddy240<0
   
   gen temddy240=-ddy240
4 g5 h$ z4 r  k1 q8 r2 b  K1 p   
    count
               
: p' d& L5 z- }4 i: _5 ?                if r(N)>2 {
  P; Z9 y7 v: A4 e6 y) X" w                 
+ S/ b1 T% o  i5 W                 pkexamine dx240 temddy240
   
# u8 ~$ V. E: k) Y/ w6 ?' p  local tempminplddy240=r(tomc)
, M* e' P! E9 j' B0 N# x* P8 [4 w                restore
$ K2 I/ }  x' E* o                 }
+ P1 M' T/ d0 _, L                 
/ p; J2 b7 m% r+ e" b% i        else {
- U0 r' z( U4 N               
0 g; T/ Q1 b, m9 V$ n% ~                 di "cannot find pminplddy240"
                restore
                 
5 K4 n, s) u- m) L6 K" H2 p$ {4 i                  sum dx240
                 local tempminplddy240=r(mean)
                 
' I5 ?) J+ m" Q' Q                 
                 }
+ }0 V% g( S. l   
8 T- I. T2 [$ L   
   
- w6 _/ {; K8 w6 Q# @# p% i   
( H2 M3 @/ m: C; s4 C# d( i$ A, p   
/ g2 ^9 ?; `8 `        gen pminplddy240=`tempminplddy240' in 1
       
        *-------------ddy240pmaxp左侧高峰值 -------------------
       
         preserve
9 t% A# c( d, o6 Q8 i   
   keep if dx240>`tempminp'
   
9 a1 ]1 Y9 n+ r2 t) ^   keep if dx240<`tempmaxp'
   
1 c0 Q; X/ z; i! L8 \' G- J1 A   keep if ddy240>0
. x7 I+ |% \. L0 G   count
               
& M( n9 s+ ]' W6 O! I- u" ~; B5 u                if r(N)>2 {
4 s. Q% o( d2 I* b4 I( ]                 
                 pkexamine dx240 ddy240
   
   local tempmaxphddy240=r(tomc)
   
   restore
                 }
                 
! h; y$ H# C4 d; ?        else {
               
9 O3 t( \4 C* R6 V  [5 ^9 C                 di "cannot find pmaxphddy240"
                 restore
               
( h) u1 X# S' U) L                 sum dx240
) u- a% A, R& I  q4 h                 local tempmaxphddy240=r(mean)
                 
' T+ l+ \8 g% N' _. ?2 \* [                 
5 v0 Y$ X, h1 w# g( b  O                 }
8 ]9 m2 P1 ?$ f  O+ G5 u   
   
   
4 s0 N& t9 w  O+ P1 d   
   
   
        gen pmaxphddy240=`tempmaxphddy240' in 1
; U8 D! ~5 ]. w+ N7 i       
       
- v! [6 B" n9 _( J- k% L0 g5 Q        *-------------ddy240pmaxp右侧低峰值 -------------------
9 y: D# F. |6 o. e       
         preserve
  
   keep if dx240>`tempmaxp'
* `; U$ E" c4 b4 a   
& L1 C, f5 J& R$ F' q$ b$ X   gen temddy240=-ddy240
+ J7 v6 Y( F$ o  K/ k! v3 i. |" p   
  ^" G8 D: J: ^+ {& B- [   keep if temddy240>0
! s, t2 [3 K8 ~( e3 G   
   count
9 I2 S( H9 [& ~" W4 c               
                if r(N)>2 {
                 
                pkexamine dx240 temddy240
   
: ?4 y% e1 l( A, J& ?6 x   local tempmaxplddy240=r(tomc)
   restore
/ V: T$ @0 s0 E7 P                 }
                 
        else {
                restore
                 di "cannot find pmaxlhddy240"
                 
               
                           sum dx240
; P1 e9 C- ^. G% W                 local tempmaxplddy240=r(mean) in 1
               
                 
                 }
6 h# a; \  X# N/ ^, @5 Y' ^0 l   
, c) u$ \2 p% O   
2 s& _6 q' N& P) }, K  J- ?3 h   
0 o0 \) J; e+ O4 y* Q0 \   
   
/ y+ q$ K* ~: ?, Y% u; A   
) M0 ]2 ~( ~1 b$ i" {6 q        gen pmaxplddy240=`tempmaxphddy240' in 1
       
& l9 y6 d5 J" P: x' ~+ H# w5 _! l5 |        *-------------------------------------
   
+ i$ L' D2 q5 Q5 E- F6 P( g8 \   gen price5d=v5 in 1/20
* K5 D2 @* i: {$ v0 P0 ]" y      
+ D, ^# F1 \5 m+ X( ^1 R   gen price5yy=_n in 1/20
   
   gen price5y=price5yy*(dy240max-dy240min)/20
   
   *price5y &price5d
   
   gen priceny=price5y in 1
   
/ E- \8 V1 l+ I  S' _& G   gen pricen=v5 in 1

. ?1 b. u! D# {& S# D) `     
   *priceny &pricen
   
4 O- c/ N# K3 p- k   
   *-------------------------graph------------------
! P0 R+ O& {! @6 E; d* l0 H  
  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)
' T0 A% H+ x0 J+ {$ b$ b9 @
; `" v5 q* [1 v" Z5 z* ---------可选项目----------------------
' R/ \: W0 T7 g; x7 b) v1 d- ?*(spike yxis amean,lc(blue))(scatter yxis2 amean, lc(blue) msymbol(+) mlabel(amean))
5 C3 C, ]1 s' _- |+ z: K}
! }3 ]( G5 ~6 f! E$ cend
; w1 T9 ^3 F  d8 x" X

maodeng

awer 发表于 2013-2-5 14:40
! y/ X+ V' Z% K# c抱歉，最近太忙了

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

0 r5 Y+ L/ z5 H5 R# g, ?5 r6 g4 X哇，谢谢，我仔细看看，这些字母数字代码晕了
, k7 B  @" v5 E

喜欢喝冰茶

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

喜欢喝冰茶

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

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

汉水云潇

准备试一试
6 v4 i/ Y* j+ t) Y6 K先找Stata软件

汉水云潇

# w1 S1 r: i; x晕呀
; h# [7 H3 g; g, W: r$ |论坛复制程序代码
乱码夹杂在其中
( O) L8 Y2 T4 J" ^2 [0 p很烦
9 l& k) G4 t2 L; z对了两次
- F4 @! P. h6 @8 k$ P& N以为没问题了
3 |2 y! y+ m# T
8 I, y' y) Q1 {结果
执行后
. o! R) G# X' e  i还是出现
r 199
的错误
9 y3 p8 P8 p% ^- z+ X1 N6 Q

汉水云潇

试验成功，下面是三一重工的图
: w8 _8 e5 Q" Z/ o! ^- i- g, Y
( i& a: F4 |& Z

peterbev

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

中秋不圆

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

bknite

这个帖子值得学习。

jeanny

感谢！！！

Hubert2000

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

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