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

maodeng

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

maodeng

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

maodeng

. hello 300236
statement out of context
  m' ^. S+ X6 ~& Rr(119);

7 W' C& U% G, y% l. E) o大师 就是这个

sergio

这贴实在！

awer

maodeng 发表于 2013-1-28 13:13
. hello 300236
statement out of context
9 F* Z3 s4 z: M4 U4 ]r(119);

抱歉，最近太忙了

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

你把 D:\stock\test\source 改成你导出的数据存放的目录，再试一试
1 F/ v/ t% s3 i
capture program drop hello

program hello

- K% o. a4 Q/ k8 z; f4 R quietly {

. ]7 r9 }' P  b% W- R& [1 P0 Tclear

- e% k+ a1 p( xcd D:\stock\test\source

/ b0 A% L  F2 V1 pinsheet using `1'.txt

& q, H" H, g0 Y7 o( B   keep if v5>0
   
   gen vt=_n
; @7 b7 P9 }  M2 B5 P# u   
8 K9 [6 S3 O, d. z) ~   tsset vt
# T4 f$ S: d5 D2 u8 j* T# _2 w
   gen zf=D1.v5
$ n( H- [( Y6 t9 H( w
   gsort -v1
  e0 }! p( u' y9 R- W
. F9 [4 }! F! k+ z   gen avol=v7/v5
   
5 L. r3 ?1 ?% m) O8 ^   keep if avol>0
; V9 d( @) C( Y   
   drop vt
  
4 [- s0 I& ]6 c5 H0 x9 E   gen vt=_n
( ?1 Y6 b9 e( E   
   tsset vt
   
    sum v5 [aw =avol]

   gen amean=r(mean) in 1
! X" Y0 T( d0 S4 D: ^   
) m9 S5 a" a5 Q9 [/ g2 }* A) k+ t   * 总加权平均值 amean
- V1 y2 L6 O, |7 M  h" g) Z7 A   
   
   keep in 1/240

; g$ j( A/ h4 I3 E* }! F3 P   sum v5 [aw =avol]
/ D, y. {) C8 Y" b' A% a
   gen mean240=r(mean) in 1
2 T4 u# X7 Z5 P4 T: w   
   *240加权平均值 mean240
   
   set obs 500
/ O. D  e3 v9 [2 g1 `) {1 Z   
   kdensity v5 [aw =avol],gaussian n(500) generate(dx240 dy240) nograph
   
   *240 kdensity(dx240 dy240)
; L" f! K; k1 l1 _# H   
) G$ F- Z+ L- f1 @4 h   sum dy240
( a, j; r( T+ ]0 {& K   
   gen yxis=r(max) in 1
   
   gen yxis2=yxis/2
9 U/ a) i, B& D6 t   
- E# p  w  q5 ]/ L) d   *yxis
   
% Q4 x) H6 s& [8 B* s. c5 c3 D+ s   gen dy240max=r(max)
   
   gen dy240minp=r(min)

    gen vtt=_n

    tsset vtt

   gen ddy240=D1.dy240/D1.dx240
   
   sum ddy240
- r2 q2 s  X3 m( o# B   
$ G  J3 K. `5 |& S2 T7 }0 w3 A! k   gen ddy240max=r(max)
! G# [2 D' |% |' X" k8 x   
  A( U: J8 @6 e# r, @   gen ddy240min=r(min)
     
   gen ddy240g=ddy240*(dy240max-dy240min)/(ddy240max-ddy240min)
- I* u& ?2 F, B4 D3 W* e; X   
   sum ddy240g
, r5 w2 P. q6 a3 F   
   gen yxis3=r(min) in 1
   
/ m3 M4 U+ i+ K- t! k   gen yxis4=0 in 1
  
  *ddy240g
   
9 o( Z$ O3 Y  A* W7 j; o, }   
  preserve
  
  sum v5 [aw =avol]
$ E1 h+ q3 r8 Y; P4 c- D' ^3 l) c' O1 b  
  keep if dx240<r(mean)
  
4 B+ V3 f5 i2 }5 t" p; M  sum ddy240
  
  if r(min)<0 {
7 w1 R0 N& f8 n5 ?( \        
                keep if ddy240<=0
3 y/ X$ L# H" X* f$ j9 w2 [% I               
( _- t8 \: T" @1 v0 @( m& I                count
               
                if r(N)>2 {
                 
                 pkexamine dx240 dy240
               
$ J  B( k. g! \2 v7 _2 R                local tempminp=r(tomc)
                restore
                 }
! h8 S/ N1 S( o) Q& O                 
; S+ j* ~/ S" G4 P        else {
               
4 R. [3 E9 M2 Q5 M6 k5 C                 di "cannot find minp"
                 
, b6 J7 `# ?6 a, ^# a                restore
! Z/ S( M" z! N/ c- d                 sum dx240
  w' B# ^5 a- l' d5 }9 s7 s( u                 local tempminp=r(mean)
                 
                 }
                }
   else {
8 B1 [" s" |) R) k- q                 
/ G4 j9 a: n7 g- K                 keep if ddy240<=r(min)
" @: I" P. n8 k6 L2 ?5 K( K                 
                 sum dx240
                 
) L' m, ~" U( M# X* x- T) G                 local tempminp=r(min)
                 restore
                 }
" f7 ]8 A* T+ z3 G" X
, |  H. H0 F3 t7 r  G' d$ P

   gen minp=`tempminp' in 1
, G" p$ F) f* s4 m7 B# e/ u$ e0 \# T   *240 minp


+ V7 C4 F; v9 ]: C# _4 H  E   
4 P) ~' X# s4 ?3 X$ m. G$ d+ D   preserve
   
   sum v5 [aw =avol]
   
$ l5 n, X% W) R/ J* A   keep if dx240>r(mean)
   
  @: x5 @: O/ Y8 {* D; [8 O   sum ddy240
   
   if r(max)>0 {
( t  v$ k: H8 U+ n- `+ W               
" A/ j8 [5 e$ `, r( C                keep if ddy240>=0
# C7 Y& o/ Z' M9 y, Q               
! T5 a# ?, p' z9 O; L6 }5 ~                count
! ~. |, W& c: z1 o' E               
% P7 K* x) W& A/ i& P0 f: g4 r& e                if r(N)>2 {
% p6 v, K7 K* J- V0 |3 H) a) W                 
9 [6 M  B7 x( B$ s8 j                 pkexamine dx240 dy240
               
                local tempmaxp=r(tomc)
7 V# H# T3 }8 t5 \                restore
7 j. p" q8 X. n/ ^9 f) |9 M; Z                 }
7 {/ x8 F4 }& g; v& R2 k4 C. U                 
$ A$ T' V( K% ?+ {! K- z                 
. J: z. ]% ~, ?        else {
+ T1 l1 a" |6 Q+ k, f2 p6 }+ X3 k                restore
                 di "cannot find maxp"
# p! D- @! }7 B$ _+ \* r; G( I               
- U. {8 w' n' o4 g! g               
                 sum dx240
1 I0 K% ?0 y- @. \$ b                 local tempmaxp=r(mean)
                 
  \; c! P& a- f/ Y( z: n; p3 Q  \                 }
: q0 @9 l$ [& ^, g9 V. i2 m                 
                 }
* d9 q, f  T) D, E/ p( D                 
7 e0 R0 E8 N- [2 |* Y7 _4 d    else {
; h! D, A/ U( q# Z) U       
           keep if ddy240>=r(max)
          
           sum dx240
          
           local tempmaxp=r(max)
           restore
+ e- s; j  E" Y( q. \           }
       

# L- X3 S0 e/ s; c5 C7 d" \8 X+ Z- B  

   gen maxp=`tempmaxp' in 1
   
5 M3 F+ ~& |6 n3 C, |9 l! z   *240maxp
* S2 r$ h6 c- \  L0 q: V   
! M! X, a! j0 K1 k5 A   
   *-------------ddy240pminp左侧高峰值-----
   
3 C8 a  ]  b. B# j& x) g- m    preserve
: y% o5 t$ D8 h2 n3 V) a. R   
  ^- c" l* s1 o, Y* q   keep if dx240<`tempminp'
   
  E, s% c3 L, t2 g& }) Z% B   keep if ddy240>0
' {- R+ B0 u$ `2 r2 }" h. k   
   count
7 V. U* @9 P3 y& ^               
8 ~$ g# \) U7 \9 ~( a' K                if r(N)>2 {
4 X: [6 F2 H$ `$ t# |5 _, ~                 
                pkexamine dx240 ddy240
8 z! e, V. `8 L% w# R* [' p- x               
                local tempminphddy240=r(tomc)
                restore
3 F5 g8 K; D8 c$ \" C+ a                 }
& K4 x+ S7 U9 l, O: o4 \' Z4 U                 
        else {
                restore
8 N2 M) D" H  u! F2 G" V                 di "cannot find pminphddy240"
; _9 D9 X4 M  e& D8 y: {# r                 
               
                 sum dx240
                 local tempminphddy240=r(mean)
. C" i) {5 ~2 S% i, ?& r$ d               
                 
1 i8 c% l  E" j" e6 r# D$ [" g+ G                 }
% U2 u' s  J+ z4 h                 
! v1 W8 c! M, |6 Z% E$ l1 B' Q: t   
: L  D. u, D3 v/ z   
   
! r. Y$ U5 p0 J, n7 G" B. v   
   
8 Y- {6 p1 |% [7 W        gen pminphddy240=`tempminphddy240' in 1
       
: r) r0 h9 w. ~7 Y( k$ M' ]       
& X( Z6 m$ f% Z, _. i: t9 I6 `$ o        *-------------ddy240pminp右侧低峰值 -----------------
3 ?/ C! v# P* o$ \) `: K3 ?    preserve
0 I% C7 m0 H; G
   keep if dx240>`tempminp'
0 O, I- q+ h. j* S# Q5 f2 f& O   
   keep if dx240<`tempmaxp'
8 }) Q' U8 E# A" U   
   keep if ddy240<0
   
  ~" B# c7 o: r0 F   gen temddy240=-ddy240
+ M" ^( I$ f( W   
    count
2 N$ D8 r! m0 s( X8 \/ Q               
                if r(N)>2 {
: z7 ^8 F+ I+ Z                 
8 v7 c( w  D! z                 pkexamine dx240 temddy240
   
  local tempminplddy240=r(tomc)
                restore
* ^! n7 y6 y# c+ f0 e0 H: F                 }
& i& _* T9 E& p* M                 
        else {
               
                 di "cannot find pminplddy240"
                restore
                 
                  sum dx240
0 E8 Z/ M7 R8 v: `$ R! L- o                 local tempminplddy240=r(mean)
3 g$ |2 E8 b1 L                 
- F; `8 t  y0 m) a7 {3 X5 M                 
                 }
   
   
3 C, x$ k! e4 T  `- D9 w' F   
   
+ l& u8 w. R9 Q: X   
        gen pminplddy240=`tempminplddy240' in 1
       
        *-------------ddy240pmaxp左侧高峰值 -------------------
       
         preserve
. Z8 \: `0 R- L   
% S' n) R: ^, |0 c/ P   keep if dx240>`tempminp'
, n* Q: B& D8 g7 S   
' T' V. j6 V/ ?& Q   keep if dx240<`tempmaxp'
& i2 p. x3 m* I8 b. R8 o$ S" T) m   
/ X2 {# J4 [: I, [   keep if ddy240>0
+ L3 z9 {8 [6 Z+ k   count
               
                if r(N)>2 {
$ ^. R" Z5 V/ ]9 M6 X' b                 
( @7 z' s. ?1 C" p& Z# c                 pkexamine dx240 ddy240
; ~6 ^9 \5 S5 G2 d   
   local tempmaxphddy240=r(tomc)
- R" U9 d. A; s: ]* G$ \5 f6 U   
; a6 o8 ?4 T1 G2 t# j& y   restore
                 }
                 
        else {
               
4 z8 x0 k8 C1 A* d3 X                 di "cannot find pmaxphddy240"
/ y$ {6 y6 j- c) f; p                 restore
               
                 sum dx240
                 local tempmaxphddy240=r(mean)
                 
/ e1 c/ ?6 q! Q$ z                 
2 B6 |" z  r8 @- ^/ D& N" D) Z/ M                 }
   
   
   
   
   
   
        gen pmaxphddy240=`tempmaxphddy240' in 1
       
# D  k2 d7 q% e/ j. Z       
        *-------------ddy240pmaxp右侧低峰值 -------------------
0 I. p) u4 c. W$ w2 M       
* G2 h, b+ r* N- Q) m, w         preserve
, z2 s* _/ E) ?. _2 B  
& }" {) E6 B; M$ t& Z' \( ~   keep if dx240>`tempmaxp'
. U) S1 y3 r# W   
   gen temddy240=-ddy240
   
   keep if temddy240>0
* A, r' _" a! R$ b9 B, i! \   
# n6 o2 ^& p" B( z7 s, g   count
* m0 ^  t  H* v: ]9 i6 _               
                if r(N)>2 {
7 f- V1 T, M% Y$ Y4 x                 
; ~$ E# o9 M# a: ^( A7 \                pkexamine dx240 temddy240
   
   local tempmaxplddy240=r(tomc)
   restore
" x/ _5 `7 }& n5 }# B# ?                 }
# U2 d0 n/ M9 ~5 Q, w* ?                 
4 |0 E- I3 q/ w: s4 O$ h        else {
                restore
                 di "cannot find pmaxlhddy240"
  c; w  U4 l  |                 
' }3 C! E8 |$ @+ P; `1 F4 [1 L, R               
                           sum dx240
9 \- L, b, c( Z0 ]* j' ]0 M                 local tempmaxplddy240=r(mean) in 1
+ c* i9 P4 ]6 }4 s  U               
' J  F3 O. T6 I                 
5 p" k/ |) A5 A, x2 W                 }
2 u/ c& G% S. a! e' u: k   
   
  {4 Y2 M. b- C/ |" Q' j: O   
   
; ]4 I! X9 G8 v: Q   
   
; B+ i+ }5 i  A9 J  q        gen pmaxplddy240=`tempmaxphddy240' in 1
6 C# _4 f' }$ R$ A; W' V       
        *-------------------------------------
   
9 ]/ F+ `8 A8 e8 `% d# M5 [   gen price5d=v5 in 1/20
! ?9 k5 c& E. q) q$ i      
   gen price5yy=_n in 1/20
   
   gen price5y=price5yy*(dy240max-dy240min)/20
   
1 Y8 d5 f3 W5 [0 ]   *price5y &price5d
   
   gen priceny=price5y in 1
   
' `" D, v  W, P  w- M9 h   gen pricen=v5 in 1
1 o5 L" N/ k, K9 P% a  \
     
+ N6 a% A$ }4 h# G, o* f   *priceny &pricen
   
* H' D7 @# [% |" s+ M/ n2 _7 f   
   *-------------------------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 P' D+ U$ |3 W, b/ o! m
+ K) J: C- k' O0 {' K; S* ---------可选项目----------------------
, b/ r8 ~" U8 G*(spike yxis amean,lc(blue))(scatter yxis2 amean, lc(blue) msymbol(+) mlabel(amean))
}
end

maodeng

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

9 t; [# u$ P& e) q) P300236 我用下面的程序试过，没问题

( ^' `: c8 W. o# R  ^: V哇，谢谢，我仔细看看，这些字母数字代码晕了

喜欢喝冰茶

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

喜欢喝冰茶

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

4 H& k* M! v. d8 `- y9 f8 f
; X( i" i& V$ J  T5 e; u这东西是开源的？知道有个统计软件叫stata，不过要收钱。

汉水云潇

准备试一试
先找Stata软件

汉水云潇

# y# v4 j+ S6 d3 p: X3 z3 n2 S
晕呀
论坛复制程序代码
/ ?  G/ K7 v2 G& F# U+ V乱码夹杂在其中
很烦
  J+ l/ r4 f) I- t对了两次
, R, P. u4 P7 x( p以为没问题了
# r4 h$ K9 a4 ~( o1 u7 P0 D( G
+ V: f: M6 `/ ^. o2 z& i5 B  g结果
$ I$ _5 c6 @) f4 t" G+ W执行后
4 h1 A- E5 w3 ]! r' [: S5 m还是出现
r 199
的错误
( b, j# t* g! t/ P* B7 o2 ~" C0 ]0 H

汉水云潇

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

peterbev

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

中秋不圆

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

bknite

这个帖子值得学习。

jeanny

感谢！！！

Hubert2000

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

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