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

maodeng

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

maodeng

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

maodeng

. hello 300236
9 Q: p8 h' ^' Ustatement out of context
. X5 {. l; g- o6 ]r(119);
& M% x  w1 g% u; x" \
1 Z! w1 C/ s  _2 q2 b) g" J' G大师 就是这个

sergio

这贴实在！

awer

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

抱歉，最近太忙了

300236 我用下面的程序试过，没问题
1 R( J+ l: q& l: M5 G0 V+ ]
4 E* [4 h2 \  g* Q& m你把 D:\stock\test\source 改成你导出的数据存放的目录，再试一试

) q8 ~! L% W7 `9 U/ W- ecapture program drop hello
6 k$ y5 d  B; I/ B
% ^5 G) R; N" w0 _$ V. k7 j8 c' e% Xprogram hello
) W$ I8 l! U' l: ^  q3 U
quietly {

clear

cd D:\stock\test\source
, M8 }' B0 M) s. Y" v+ q& N
insheet using `1'.txt
6 z  D* e" s$ b2 k
8 O5 ]( E: m6 B   keep if v5>0
# s/ h, T, Y* |   
   gen vt=_n
1 ^- i2 k+ [& W1 N8 a1 W  Y, }   
& k% L. y( {8 i1 K5 g- T   tsset vt
" G' ?( U% V; ^# c" ]2 K
   gen zf=D1.v5
. f  v' P9 V  W: M( U
   gsort -v1
6 ?+ r: C* d/ P6 e3 ]1 r1 b6 u
   gen avol=v7/v5
) f2 ]- o( F/ X4 q" }. k+ F, m) b   
  q' k; X  B: E0 \7 M2 m+ l& `   keep if avol>0
   
* P: f5 x) B/ S* l) M4 F  J* M/ E   drop vt
% T# @$ d5 T# d+ N  
   gen vt=_n
, s! j( e" p* ~8 j, g6 K: Z' \$ b   
* Y2 U8 _% u' W2 ?& L* _2 q   tsset vt
   
    sum v5 [aw =avol]
5 b4 ~; G- J) i' n
   gen amean=r(mean) in 1
9 I4 r3 w, X! P' j( q4 A+ R/ _; T   
   * 总加权平均值 amean
. _# C& s9 t5 W1 i2 T   
1 k4 o& ~1 |; w- Z4 A- ^% [: o   
   keep in 1/240
+ \- i$ B; H: ?' M; |, Z
   sum v5 [aw =avol]
! V- k9 d* t. \, m
   gen mean240=r(mean) in 1
  b* @" A" x' g% F   
% B3 h) K2 f2 n2 W5 I; X; e8 O/ p8 }   *240加权平均值 mean240
/ W$ y( g5 e2 [. R, N% I   
   set obs 500
9 X; W$ _. n% Y* I9 H/ g/ V   
   kdensity v5 [aw =avol],gaussian n(500) generate(dx240 dy240) nograph
   
* n2 Q0 g1 Z+ p   *240 kdensity(dx240 dy240)
  T/ s/ T# A% o0 ?   
   sum dy240
. g0 x& B/ J& I: w) N   
) Y: S/ r4 L0 K2 }& [; z   gen yxis=r(max) in 1
5 c2 c) R, r5 z* _, X( G$ ?   
$ H. A* d) z) y' R- v' g: {   gen yxis2=yxis/2
) k, @# G" w' F2 w2 r   
   *yxis
& W( {1 H# h* U( d1 K' L   
2 b6 z3 T9 O8 ^# @; m4 o9 w- i   gen dy240max=r(max)
   
   gen dy240minp=r(min)

    gen vtt=_n

    tsset vtt

   gen ddy240=D1.dy240/D1.dx240
   
   sum ddy240
  T: Z! W3 I- ^, Y4 c/ C   
   gen ddy240max=r(max)
   
; e& `$ T3 p' B4 [9 z, I; x   gen ddy240min=r(min)
     
# r9 b2 q1 I! t/ Y* N* E2 }   gen ddy240g=ddy240*(dy240max-dy240min)/(ddy240max-ddy240min)
. s% F0 W0 I+ V   
   sum ddy240g
  z& j8 \, _$ ?$ t) P   
2 m& w' D2 P# o' x6 t, \, n   gen yxis3=r(min) in 1
( p' v6 l2 _( J   
   gen yxis4=0 in 1
  
7 i/ }3 h0 j2 \! S5 ^8 y, r  *ddy240g
  [( ~1 C( x- a   
   
* U: S% \0 M  W7 d5 a3 W# S9 F  preserve
  
  sum v5 [aw =avol]
  
  keep if dx240<r(mean)
  
  sum ddy240
5 v& z. w( o  O- q  
  if r(min)<0 {
  X1 b) [' J4 O9 ]        
7 I/ p1 N7 p+ Y" ]0 \. f                keep if ddy240<=0
2 r3 |* ?- N/ c               
3 c7 [# P( O) ?4 c) q1 _, f                count
               
                if r(N)>2 {
' b. L4 s; t/ F% _                 
                 pkexamine dx240 dy240
               
9 `7 E5 ~; h( Y: ?+ I- V3 A' k. x                local tempminp=r(tomc)
; i, Z: c: F$ R5 p                restore
                 }
                 
3 {4 C# _% `2 {4 k- d        else {
9 H6 E; J6 {/ f! V5 E( k               
) s' U. H6 V0 c                 di "cannot find minp"
                 
$ m5 W! o3 b; X5 O4 o& w                restore
5 S- n0 @0 @" l+ i2 v                 sum dx240
$ F( s5 P2 ?8 h/ i+ f4 q                 local tempminp=r(mean)
5 o; |! y% e" Q                 
                 }
                }
" b, p; \0 U! R, W# {   else {
4 j7 e+ [# T/ k3 L$ T1 `: H                 
% U7 d! h* e. ^4 t7 E                 keep if ddy240<=r(min)
                 
                 sum dx240
                 
                 local tempminp=r(min)
& J4 D: f  J( s0 Q                 restore
                 }

4 r6 o# R+ b) S2 m# h# k
) q( `+ }: G8 E6 [
8 V1 f  |5 ]) i/ s' ]   gen minp=`tempminp' in 1
, c, D2 k  l* c# v: J- m2 T2 G   *240 minp
3 P4 w/ ~% j9 d' \5 b" w7 d
. E+ L( s( m9 a9 x
; ]' K/ O' ?3 _0 e0 \   
   preserve
) F% ?6 l/ r7 s0 S   
5 o: i# N% ^5 q% F  t   sum v5 [aw =avol]
   
   keep if dx240>r(mean)
3 ^/ l! n3 B/ t6 f   
   sum ddy240
5 m! t3 z0 N" W. |   
  p+ B* Z; N) F0 F   if r(max)>0 {
               
                keep if ddy240>=0
$ f6 \8 c. ]9 L( r1 }               
9 q. \+ w2 J& m7 _1 W6 ?                count
               
                if r(N)>2 {
                 
# X- X+ W+ K; X$ y6 s                 pkexamine dx240 dy240
) @- o8 Y. h$ Z0 b3 e               
' g0 I( S9 I0 n& Q" Y" W( X8 e                local tempmaxp=r(tomc)
% ]. Z7 T  k1 S0 f4 L! z5 R                restore
  z. B% y2 A1 ?, |, t/ ]) Q                 }
. p# S$ B* t9 R0 I5 f/ ]  ^                 
                 
- n% ^/ Z$ C$ ?# L5 J; _        else {
                restore
                 di "cannot find maxp"
               
) Y! ]6 o( Y  ~               
# d" b6 f: ~( U1 l( Q                 sum dx240
0 T; q# x' O) [' d                 local tempmaxp=r(mean)
9 j4 k  z& W8 V% L2 t                 
8 s, O2 M6 Q% \! M5 ~7 J                 }
! @2 v" A% K  G& k                 
& A9 `; L) I+ D                 }
                 
    else {
       
5 i) _9 e' I/ }           keep if ddy240>=r(max)
          
           sum dx240
8 T0 z2 V& c& C1 ?5 p. ^          
- o- p- e  ?) F( Q           local tempmaxp=r(max)
0 @4 q  y8 \" C# E3 i' W           restore
% H% G! x& J! r3 N' N# a           }
       

% ?* y2 B# y- J  

   gen maxp=`tempmaxp' in 1
, T6 j  ^  ?' N) o   
   *240maxp
; b% v" [# ~1 L8 c7 N& ^5 }   
   
: V; O% a8 [! Z" x   *-------------ddy240pminp左侧高峰值-----
   
    preserve
   
* ~, o( i0 l; z7 C4 `. k1 T* }: Q3 v6 \   keep if dx240<`tempminp'
( m1 m: n( x- y0 c   
   keep if ddy240>0
7 p' U9 U1 ]+ z# m: G   
1 B. b# ?% v0 |1 B0 @   count
               
                if r(N)>2 {
                 
& J- p4 `1 n2 o% u1 [                pkexamine dx240 ddy240
               
                local tempminphddy240=r(tomc)
, |% ~% ~5 l$ p' W' r1 ?* X* k, T                restore
9 ]+ X% O1 m) R( q4 ?6 D) |: T* i                 }
                 
) T( `* w# ^- W) Q( c  w  ~+ j) O$ E        else {
                restore
1 q$ q: L1 l( A' H/ K' w                 di "cannot find pminphddy240"
                 
               
. M2 n& D! M& ]                 sum dx240
* n! `8 D. F' C2 y                 local tempminphddy240=r(mean)
, X: T) g" \* Z# T, Z/ g3 S               
                 
+ g% j5 T% x" r* i5 \" }                 }
( ]/ m0 r  _- N7 z7 z9 a# T                 
   
   
( K$ C2 t; z4 K# x+ f   
6 t7 h7 u  |2 u8 B0 V/ d0 U   
7 X% y# E/ T/ ?- `% ?9 o. }   
        gen pminphddy240=`tempminphddy240' in 1
4 ^, w5 \  ~" ?( b       
       
4 S1 J/ R$ B. a# F        *-------------ddy240pminp右侧低峰值 -----------------
% i, N$ j  S3 ]: p6 z2 U9 ~! N    preserve

   keep if dx240>`tempminp'
, r3 F$ W6 y- Q. x" ?   
  S1 A3 g' c8 O9 [" o   keep if dx240<`tempmaxp'
" i0 x0 B& b2 O" J: G* v; h   
& T& }, h  l! s1 V( g4 U" a* i   keep if ddy240<0
   
   gen temddy240=-ddy240
# M- D; ^( s4 w  A+ g2 ~4 N# s   
    count
               
                if r(N)>2 {
                 
3 k& A& ?; a+ x, f" P7 ?) i* ]                 pkexamine dx240 temddy240
   
  local tempminplddy240=r(tomc)
                restore
                 }
                 
        else {
               
                 di "cannot find pminplddy240"
                restore
! c* V) d3 p6 H, `2 e4 e                 
                  sum dx240
% F: t" k: F+ R& d( J                 local tempminplddy240=r(mean)
$ i2 R5 b, v2 P1 A2 X) \0 K* o                 
& {- u/ s# c3 i4 x                 
0 G' N* d1 c8 \6 B                 }
3 Z8 \. a( D" ~) }) w' T' v+ E/ n   
0 {  t$ C! [  }& u% J) H   
   
   
" [. z) l  {% X   
        gen pminplddy240=`tempminplddy240' in 1
       
        *-------------ddy240pmaxp左侧高峰值 -------------------
       
         preserve
   
   keep if dx240>`tempminp'
   
5 b: j$ t; i; t5 S5 \   keep if dx240<`tempmaxp'
   
2 y; e5 i( |$ ?2 @& X' ~# O8 z  a# ~   keep if ddy240>0
   count
8 k! Q) w9 g3 B7 f& l1 L+ a               
3 m/ Q- U) o% p                if r(N)>2 {
1 _8 _, U0 l: C                 
; P; K) T) M9 s  D7 `0 S/ w: q. I3 N- }                 pkexamine dx240 ddy240
- E* \" Z# C& m  z; B   
   local tempmaxphddy240=r(tomc)
; @4 k. j$ W6 v% R   
   restore
! u" }- q+ G1 h! {& d7 Q# p( Z/ J                 }
                 
# B7 c( n# _8 ]" _4 W        else {
1 Q! C5 ~1 L5 i- O7 M6 l8 C7 d9 u! Q               
                 di "cannot find pmaxphddy240"
: g( V) w3 s: i7 @  S6 v4 r9 k                 restore
! F4 `, @' M1 I$ L               
8 p2 [/ l; b" d; \! Q5 J                 sum dx240
( B5 l1 j1 _( v9 V                 local tempmaxphddy240=r(mean)
5 g' `8 W; e) B9 b                 
                 
                 }
+ z3 C- C, x- ^: ]/ I+ a; m   
& Y& a3 A3 Q3 ^& `6 t$ l! B   
$ k2 ~/ U$ S: K! v! Z   
   
   
5 }$ a. Y5 `, O! [2 D/ @2 x+ Z, b   
        gen pmaxphddy240=`tempmaxphddy240' in 1
       
       
  L' M7 J0 e6 v" w( J$ a1 s! h        *-------------ddy240pmaxp右侧低峰值 -------------------
- m5 ]. g( V/ k" ^       
* p, P5 H  x9 }1 j         preserve
  
1 g0 ^+ s" V& X7 k   keep if dx240>`tempmaxp'
+ N% N5 G' F& k9 h' R0 ~: q   
5 r. w. \- z1 N3 H! i( j) T   gen temddy240=-ddy240
   
   keep if temddy240>0
+ T7 t9 Z. `+ l+ Y   
   count
               
  G/ Z. w9 z' S) {                if r(N)>2 {
                 
. v7 p9 f" h! h1 o$ u- D6 S3 k$ b                pkexamine dx240 temddy240
   
   local tempmaxplddy240=r(tomc)
   restore
                 }
! |' A) H9 E  P0 M  I. C                 
$ j0 r2 L3 {& w; e8 C3 `& k8 v: k        else {
3 k5 s1 @% G. ]# _6 H                restore
% o% L+ R# R/ y. ?2 p8 |# _# j) h  c                 di "cannot find pmaxlhddy240"
                 
               
7 c& {2 p, w! r. p2 X6 b                           sum dx240
                 local tempmaxplddy240=r(mean) in 1
               
                 
                 }
   
" T8 H3 i& h( e& Q2 @   
3 u  a) t3 Q- m, K5 ?6 m4 w7 x   
   
2 @. X+ r; y# t/ J2 ?! V   
! K. F( x+ p8 K( ~# c5 m% ?: z   
        gen pmaxplddy240=`tempmaxphddy240' in 1
) t% L  w# o7 ?$ |  j, H6 a3 p       
        *-------------------------------------
  S8 r5 W6 Q) h, \9 j, x5 `   
) |; c1 ~, U. t2 R$ ~! s7 j   gen price5d=v5 in 1/20
      
   gen price5yy=_n in 1/20
9 |8 C2 O7 l9 G- H: U& C   
   gen price5y=price5yy*(dy240max-dy240min)/20
   
7 k0 @% x2 [, T, R- `* D   *price5y &price5d
   
3 k' j! ]+ o/ A! N4 b7 g9 G   gen priceny=price5y in 1
( c  P3 {8 P2 ~9 R$ w* D   
   gen pricen=v5 in 1
& ]' Z* Q$ E* W! h7 B* v7 R7 E
3 I- c$ B$ \$ X     
: [$ V( N  J3 z- N0 k- K# x   *priceny &pricen
: f' E0 \0 D; i$ N+ x   
   
; I# o+ A, e- W" N2 p8 D- H! v5 ^   *-------------------------graph------------------
  
1 Z1 Q' C6 }; M" S, k  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)
9 }* H, a  j; J, q4 G
$ c: v$ y# X* S* ---------可选项目----------------------
* |6 B2 e, P4 Z; |*(spike yxis amean,lc(blue))(scatter yxis2 amean, lc(blue) msymbol(+) mlabel(amean))
6 @# P1 X+ f( m6 g' P0 H}
3 z9 i5 t, y2 u+ v/ o* U4 Mend

maodeng

awer 发表于 2013-2-5 14:40
5 ^9 F6 Q/ ]# {3 i4 B5 x抱歉，最近太忙了
4 f$ ?6 P5 w2 w0 p( Q
, B4 ?2 c/ f) r, W+ V300236 我用下面的程序试过，没问题

4 O# `: E% {' N% N( U- @& w哇，谢谢，我仔细看看，这些字母数字代码晕了

喜欢喝冰茶

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

喜欢喝冰茶

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

7 A' [; j8 m# u, s. }8 U0 \$ k
这东西是开源的？知道有个统计软件叫stata，不过要收钱。

汉水云潇

准备试一试
5 q) q6 N/ T4 O先找Stata软件

汉水云潇

7 Q, }7 i6 T0 A# R晕呀
论坛复制程序代码
乱码夹杂在其中
3 D" s( K& U+ S6 b很烦
对了两次
1 n  l" }% e' v3 b4 P以为没问题了

结果
0 ^3 T9 X' }7 g# w执行后
还是出现
+ L* a7 s2 }8 _+ Br 199
的错误

汉水云潇

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

peterbev

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

中秋不圆

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

bknite

这个帖子值得学习。

jeanny

感谢！！！

Hubert2000

明月回春 发表于 2011-8-13 20:17
  k) C0 D9 a+ v程序写的不错.
, _. u2 |6 G9 c+ T9 Z不过有一点要注意.通达信导出的数据,只对价格复权,并没有对成交量复权.所以数据需要手工处理 ...

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