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

maodeng

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

maodeng

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

maodeng

. hello 300236
' ?. G+ t4 u/ f* ?. zstatement out of context
r(119);

大师 就是这个

sergio

这贴实在！

awer

maodeng 发表于 2013-1-28 13:13
. hello 300236
statement out of context
  \+ Z  s8 H4 s$ O' I$ f) nr(119);

' u) }4 o3 ]$ M$ O. E5 h* a抱歉，最近太忙了

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

你把 D:\stock\test\source 改成你导出的数据存放的目录，再试一试

% ]8 m( V2 ]4 O1 O% Fcapture program drop hello
1 E( \$ J; Z: p  ?  |1 ^0 ?
5 N2 I; d) n. e) F  |: Mprogram hello

quietly {
: @, k: E) `4 H, A9 E$ s
clear
+ t, ], {0 ^! D+ ~
cd D:\stock\test\source
+ x8 K  ^0 m. c2 X
- q8 \- I$ d- z9 z8 ~! v! einsheet using `1'.txt
) U, V- @& X& X7 m  \7 i$ z
   keep if v5>0
6 P' u( s2 u* x7 K   
7 H0 |3 z- m! L9 i   gen vt=_n
   
0 u2 N5 l6 y8 q4 a/ a% `* r, U   tsset vt

: g3 z, f+ b. L* C& ]9 ?   gen zf=D1.v5

   gsort -v1
+ ^( F$ V+ `: y4 [
$ a% @, u% s, u" @) f% D5 h' z   gen avol=v7/v5
   
   keep if avol>0
   
0 }: s: u% g7 b9 E1 z# M  @8 g: P   drop vt
, t% w0 }" R9 n" W1 i" u8 |; \6 E  
   gen vt=_n
   
5 A, D! t0 w+ m( @/ L   tsset vt
   
- d! c8 Y3 O3 L) s. g" U( o    sum v5 [aw =avol]
! y, m/ Z3 A7 l6 y$ h1 C
   gen amean=r(mean) in 1
* f6 |3 H$ m7 E. z3 |4 m  t% i5 g   
   * 总加权平均值 amean
' b+ D$ D! f, b   
3 O6 f! E, d4 X) D6 @5 I   
   keep in 1/240

   sum v5 [aw =avol]

   gen mean240=r(mean) in 1
   
   *240加权平均值 mean240
   
5 G! N4 o) N! v8 S0 j8 G   set obs 500
   
   kdensity v5 [aw =avol],gaussian n(500) generate(dx240 dy240) nograph
   
: w6 Q  j) x5 n' e   *240 kdensity(dx240 dy240)
   
   sum dy240
9 R  \; Y4 S" i   
   gen yxis=r(max) in 1
# o6 C# Q$ x- A& i1 ~# Z4 t   
% Z. Q- ^: M7 u# G   gen yxis2=yxis/2
0 H" L  l9 }$ R) X   
7 o* w0 p+ H+ W+ C   *yxis
% Y, B+ y2 _- ^3 S) P   
   gen dy240max=r(max)
   
   gen dy240minp=r(min)

    gen vtt=_n

    tsset vtt
& Z* j% [/ O* B  O
   gen ddy240=D1.dy240/D1.dx240
2 ?3 o" V# m9 A! O' Q5 Q   
   sum ddy240
3 n8 n* h  y0 p& t1 ]: J/ x7 O  R   
   gen ddy240max=r(max)
9 h/ _! q$ R) M& g   
7 g# e) c' Q7 Y7 P! L   gen ddy240min=r(min)
     
   gen ddy240g=ddy240*(dy240max-dy240min)/(ddy240max-ddy240min)
5 l- t8 b3 O1 G  T+ ?2 D8 Y   
" h: V* X) m/ ^  Y0 v2 m# r" |/ L   sum ddy240g
. ~& Z$ w3 [7 H   
4 |5 \5 }3 k2 c9 A" @   gen yxis3=r(min) in 1
/ E* F; ~4 b; C0 x2 x6 J/ X   
   gen yxis4=0 in 1
6 u, |% ]1 g# H1 C9 W" T$ \  
1 i2 v7 ]3 K0 m7 z  *ddy240g
   
   
& R( J: m. z- O. K7 M+ W0 r( r  preserve
4 e9 \# r7 s" h0 E  
  sum v5 [aw =avol]
% f4 |! A' M0 D0 q  
* T+ }8 C# e6 [/ W  keep if dx240<r(mean)
4 ~; n' s+ _2 C# z( y! P( ?  
1 ?. u% L- f8 R! X  sum ddy240
  
( L' E4 T4 _9 n: Z, B) T  if r(min)<0 {
  m& w1 j3 }- ^! b4 g        
                keep if ddy240<=0
               
                count
               
/ G7 @  d- {0 Q  n1 A2 S# c                if r(N)>2 {
3 v7 s) L) Z+ _" E* T, W                 
0 Z$ T  F% o% `8 d) A* s  B* C1 Q6 c                 pkexamine dx240 dy240
: r. A7 p! m$ @               
: \( [" N+ G$ ^0 D                local tempminp=r(tomc)
  z; y. b1 @, _4 Z/ t                restore
! G' h* o  U7 K5 N6 Q                 }
% v/ g- o9 Z! Y" h( [7 k2 N                 
) o' b; g. S! v; {& j        else {
               
  D" V% ?3 k; g! J; f7 X                 di "cannot find minp"
                 
- v) r/ @$ s9 A& \: ?4 K                restore
8 @) f$ J: ~- B, k) l4 D                 sum dx240
                 local tempminp=r(mean)
                 
                 }
8 a& f. v3 I* n8 b                }
   else {
                 
. [6 V) A# i4 D9 ~                 keep if ddy240<=r(min)
                 
5 k* f* R' x9 F2 E                 sum dx240
                 
$ ?: h" w7 h. n. ?5 P+ H                 local tempminp=r(min)
                 restore
                 }
% F/ Z$ s( I2 Y$ a. t
7 S2 l; P8 G" x# P! v- d1 F2 {
% `# `$ g6 n) r3 ~/ u+ {- @8 |
8 d' h7 K& H# L1 s- x; E* T' Z   gen minp=`tempminp' in 1
  k( F+ X9 h- c   *240 minp


   
   preserve
   
* D; b( d; l3 n, Q7 i$ u1 C   sum v5 [aw =avol]
   
   keep if dx240>r(mean)
   
! M; X8 f& {- x) g& }" z   sum ddy240
   
   if r(max)>0 {
               
% Y- N- O2 c' ^  Q% G                keep if ddy240>=0
3 d0 \5 Y' O# I# l- g               
# Z, q& ?' h) w, F0 U8 i, |                count
               
                if r(N)>2 {
                 
4 A# s2 O) g+ o, z" i                 pkexamine dx240 dy240
) K: B/ p# Z  A) l$ S               
- ^9 l1 a2 E. d                local tempmaxp=r(tomc)
$ [* N  h( H) G8 C9 m* W2 [                restore
                 }
. J% }6 n7 ~% g& S  a3 c" U                 
                 
        else {
                restore
$ {( w- A* X/ m" Y                 di "cannot find maxp"
               
               
                 sum dx240
" Q6 K, Y+ w2 S                 local tempmaxp=r(mean)
0 n9 w1 a1 c; T7 d7 n% d                 
) ^, d1 u0 A: G7 ~6 d                 }
                 
1 I0 ^. D5 X3 y8 M% [: {  |9 i                 }
% |0 U5 W: H' ^$ @+ a2 a* ]3 u5 W                 
    else {
& w) u1 x2 u7 U3 k2 ?' q       
3 C# f) S) K8 K  v( T, q           keep if ddy240>=r(max)
          
# x5 m/ u9 `: k' s/ o5 E           sum dx240
          
7 v+ ~9 S9 h& B: J) T           local tempmaxp=r(max)
           restore
           }
1 t& k; g8 I( Q6 Y6 h3 U       
: Q! H- U5 {  S" U. t$ S% X4 ]
  

: r0 U) o# V  h7 x   gen maxp=`tempmaxp' in 1
) J( H* U& c' G% c, i0 o   
   *240maxp
2 d. [6 D' F) l' w9 k   
7 |9 h0 O5 l" Z7 Z- {! I5 N   
   *-------------ddy240pminp左侧高峰值-----
   
    preserve
   
   keep if dx240<`tempminp'
   
   keep if ddy240>0
# \/ b/ z6 o5 n+ v/ j- F( |, E   
% ]) i8 T, \5 C3 X6 j; u   count
               
                if r(N)>2 {
                 
                pkexamine dx240 ddy240
4 U/ A# K3 L& Z9 N0 k               
                local tempminphddy240=r(tomc)
                restore
7 t+ Q) ?* F+ r5 _/ z0 u8 m                 }
# K: V3 y  Y9 x7 q3 m; [8 I                 
        else {
$ ~7 R* U# ]: }                restore
                 di "cannot find pminphddy240"
1 O3 ?( a* M" P                 
+ r- H2 |4 W9 l% ]               
                 sum dx240
                 local tempminphddy240=r(mean)
5 q2 R$ O. l0 L/ K6 j. ^               
                 
5 h9 }/ B$ E8 y" S0 Y: ]                 }
! N4 }9 T; M1 r+ X# M                 
   
4 c+ T. l* z, j4 O# V' S# k' v   
   
3 X- v% k2 e' Y1 D0 ~# b   
' s, ?" H- P, @  {   
8 S0 s3 U9 r( E% s" F* \1 B        gen pminphddy240=`tempminphddy240' in 1
       
       
) w5 }2 v2 _0 E' `" M        *-------------ddy240pminp右侧低峰值 -----------------
    preserve
5 n, Q8 H, V& C' I8 C, _
  e2 k1 S7 E4 u' L  {   keep if dx240>`tempminp'
   
   keep if dx240<`tempmaxp'
8 E5 z1 I, Y3 o  v+ J! m, t/ n   
   keep if ddy240<0
6 D: N$ B3 E( X. f8 \! g   
   gen temddy240=-ddy240
7 M* ^, c, t! g: A   
    count
               
                if r(N)>2 {
; \( v/ B8 g+ l: H" p" W                 
                 pkexamine dx240 temddy240
   
  local tempminplddy240=r(tomc)
                restore
                 }
5 W' Z9 N$ k' ~0 ~1 Q" D3 p; B/ n                 
        else {
2 D/ }1 p# |6 ], Z8 j: b2 `               
                 di "cannot find pminplddy240"
                restore
                 
& R; d: u' q) s4 q) ~5 r- t                  sum dx240
                 local tempminplddy240=r(mean)
                 
7 B0 u- V6 Z# g% M0 s                 
0 Z& Q- h0 [- @$ T6 Y9 V) @5 Y) y                 }
; \" q' ~1 \! n3 t$ Z   
   
   
/ L) z* R' y7 u9 h% ?; [! Z   
: r& c2 t# n9 S   
) t  q; A* A; B  B- X        gen pminplddy240=`tempminplddy240' in 1
       
        *-------------ddy240pmaxp左侧高峰值 -------------------
! \( @1 z( C5 ?5 a       
         preserve
   
' s# q0 V& s3 \. D: b& w# ^5 F1 Y   keep if dx240>`tempminp'
, c# X8 d# }$ v/ b: j   
   keep if dx240<`tempmaxp'
   
   keep if ddy240>0
* z7 I& C* \9 f) a) s( c9 @4 W7 M   count
               
1 b: l& R, s2 |4 N+ M                if r(N)>2 {
% K* r# m6 n% P, E! A7 p8 Q! b2 Z                 
8 N' k( u9 G. |. R0 z4 ^                 pkexamine dx240 ddy240
   
   local tempmaxphddy240=r(tomc)
   
+ d2 @6 Y( `0 a. e4 C. g( y% b7 N$ w   restore
0 Y, T" u6 j/ z5 G6 k* n                 }
                 
2 o1 Z5 `% B: p5 |0 C" e9 s5 T        else {
               
) X+ W4 z0 v0 A' F# ]                 di "cannot find pmaxphddy240"
                 restore
               
                 sum dx240
- H  l, g- b* u& C+ y4 a" q                 local tempmaxphddy240=r(mean)
                 
) {: q& o: \/ n) n& @5 n2 l; U                 
* Z1 J* u  G% z/ O8 }" J2 Q% t                 }
   
6 e2 d# w  e" d" Q; o, c1 [   
   
: s. o% d7 x- ^5 X, }3 e  ^   
   
   
/ k5 x$ M9 e4 L8 o  x* b* E        gen pmaxphddy240=`tempmaxphddy240' in 1
       
       
        *-------------ddy240pmaxp右侧低峰值 -------------------
       
         preserve
* K7 P9 z; Y) k9 z' `  
5 B) h7 a, H9 C" P6 [8 }   keep if dx240>`tempmaxp'
- Q6 L- B+ d* g+ \$ R   
   gen temddy240=-ddy240
) C3 ?. {) D+ O   
: \8 T5 Q. r: D4 o& p   keep if temddy240>0
   
   count
+ s1 C6 }: Y9 w# x4 L; U  ]: U% p               
3 Y6 b+ X2 S( Q' m                if r(N)>2 {
3 m* ]; j' C1 \  o0 f8 A1 K                 
. M7 X: l* }) l' K: x                pkexamine dx240 temddy240
   
0 x+ L  q  Y& x, G  B   local tempmaxplddy240=r(tomc)
   restore
3 p& H4 b5 g( v+ T' M( F( ]  b5 t                 }
                 
4 f7 w$ y0 X1 C! b. ^" P0 o        else {
) c4 W/ p; q. k* a% O- k: b7 [                restore
! H  d: O8 S3 h8 l$ z  L. n/ @                 di "cannot find pmaxlhddy240"
7 n  Y; ^5 V" A# H! B- q" r8 A                 
               
                           sum dx240
9 n. T5 X$ O5 m1 j- N! K. ^: E' A# t                 local tempmaxplddy240=r(mean) in 1
4 {: ]  j# b9 o" ^# K               
                 
. B0 ?! z& ?" V7 c' N                 }
: O' u6 c6 o# f/ ~/ l   
# b; S- E- y2 @% p   
   
   
2 H5 q% g4 A% h0 C& i" t% u# S+ P   
' p) @/ }2 Y, q' t* }   
        gen pmaxplddy240=`tempmaxphddy240' in 1
! n5 l! e# N0 J+ E# l       
9 K* l1 I9 K) f: g  E; Y& @        *-------------------------------------
5 `  C) o  E$ N$ a   
   gen price5d=v5 in 1/20
      
   gen price5yy=_n in 1/20
   
   gen price5y=price5yy*(dy240max-dy240min)/20
) M6 C; V1 ~" I) {5 u  ^: M   
% K  J0 v$ {/ c$ e/ r, e' [4 S   *price5y &price5d
% O. z4 V3 x" Y$ f- c   
5 ^5 q: l, J7 M2 R. d   gen priceny=price5y in 1
8 C1 |9 E: a+ w) E: Z   
  G* T: {  X% \( U2 r4 i" m   gen pricen=v5 in 1

, t, y9 H$ K0 \+ k1 ~     
+ l3 ?2 U: n* T+ ?, b4 c' _1 T   *priceny &pricen
; }8 q- l. M/ P1 c4 l8 X- E   
: N2 _4 k& y2 e" z9 G   
   *-------------------------graph------------------
5 c6 h" ^4 W7 q) U- Q: L  
, h7 d: O+ x7 r" |4 A' X3 _  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)

/ Q) X) V6 O9 w9 N1 |* ---------可选项目----------------------
*(spike yxis amean,lc(blue))(scatter yxis2 amean, lc(blue) msymbol(+) mlabel(amean))
}
end
5 W4 D. `5 \; ?) _! O

maodeng

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

2 t& N$ b9 M! D+ ]& R8 }300236 我用下面的程序试过，没问题

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

喜欢喝冰茶

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

喜欢喝冰茶

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

" V. [( I" p( c; d5 f" P# k
3 A+ n8 |, p8 @+ ?& ^这东西是开源的？知道有个统计软件叫stata，不过要收钱。

汉水云潇

准备试一试
先找Stata软件

汉水云潇

晕呀
论坛复制程序代码
; L* Z' [0 u1 j. R0 S. D乱码夹杂在其中
) G. p$ z7 ?1 N很烦
# q3 `$ w5 T$ u% K. @对了两次
以为没问题了
  E) }; Q% d* Z9 O2 W# L
结果
: X) [9 C% q7 n$ c1 @执行后
还是出现
r 199
6 U- `7 [/ D/ b7 Y的错误
: m; |- y' a" j2 l

汉水云潇

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

peterbev

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

中秋不圆

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

bknite

这个帖子值得学习。

jeanny

感谢！！！

Hubert2000

明月回春 发表于 2011-8-13 20:17
2 X% C, {3 s+ |4 J- |, Q3 b程序写的不错.
0 F/ A( J1 D' F1 W+ ~1 ]8 _7 J. [不过有一点要注意.通达信导出的数据,只对价格复权,并没有对成交量复权.所以数据需要手工处理 ...

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