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

墨轩书童

感谢楼主！！

海阔

按照版主的步骤操作到了最后出现

. hello 600031
avolDo not found
---------------------------------------换成其他代买，也是一样；总是出现这个问题，怎么办呢？奇怪

maodeng

大师，那两条曲线代表什么啊

maodeng

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

maodeng

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

maodeng

. hello 300236
: @# \0 D3 O  u; t1 O0 s: n6 `statement out of context
" u' W- o" i. a* `' zr(119);

4 G, W  f  Z- p$ B大师 就是这个

sergio

这贴实在！

awer

maodeng 发表于 2013-1-28 13:13
. hello 300236
* ?$ n! Y! T: w2 z( H8 Hstatement out of context
r(119);

抱歉，最近太忙了
. l2 U9 W. W: F' m
300236 我用下面的程序试过，没问题

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

! i* i- G$ `1 Zcapture program drop hello

program hello
/ X" t( S+ f; _  T  t/ [2 I- h
9 x# l# ]& Y" o( I1 ], K8 f quietly {

clear

cd D:\stock\test\source

insheet using `1'.txt
- N" X5 F) P, ?# _+ y
   keep if v5>0
   
   gen vt=_n
3 u+ m  O& R4 @8 {6 G; @& v   
; n; y$ ]0 k; ~+ ]: X' n: ~   tsset vt
! A+ P: u4 F" s9 A" h; t! I8 A7 F
* ^  ~9 N( W& [, C/ Z: A: X! c   gen zf=D1.v5

   gsort -v1

6 e. U6 N2 z4 U, p1 k* ?4 s2 z6 L9 ^   gen avol=v7/v5
   
- u: S: @8 \$ l! Z% G: W9 r   keep if avol>0
   
4 ]' O0 ?" v" J5 W0 `6 n8 B   drop vt
  
   gen vt=_n
0 n1 D4 _* N+ x. n' ^1 ]/ u% Z# V   
   tsset vt
5 l7 p& u3 x$ \6 l8 l   
    sum v5 [aw =avol]

2 O; r5 |3 U+ F. k- y6 {   gen amean=r(mean) in 1
   
" T- K+ f$ H4 m$ C0 P! P   * 总加权平均值 amean
   
   
; }) `0 G' N- ^* k- F   keep in 1/240
: r! p" `- Q& h+ Y" B
0 N# Q6 F& p* |4 H$ T; S! r   sum v5 [aw =avol]
0 s2 L4 s/ `; `) r  c
   gen mean240=r(mean) in 1
   
* d  Z! @6 Z% e) R/ L; p) }   *240加权平均值 mean240
0 c7 {8 s. W# T/ E   
   set obs 500
   
, `/ v( D2 T* W+ ^. `( A2 t' e   kdensity v5 [aw =avol],gaussian n(500) generate(dx240 dy240) nograph
) }5 e7 v- A9 k& C% Z9 h   
# V" P7 M& T/ A" r2 m* [   *240 kdensity(dx240 dy240)
   
   sum dy240
   
   gen yxis=r(max) in 1
   
& \* L/ Q: ~! i. O# z  n) Z   gen yxis2=yxis/2
' E( Q6 [% {6 v2 x+ O8 g+ e+ k   
, i! a0 k2 F( l5 C9 e1 J   *yxis
   
   gen dy240max=r(max)
* b& S" S: V+ f+ i) ~   
7 v1 x3 ~' I* \. `" U   gen dy240minp=r(min)

. z; w# i+ }9 ]    gen vtt=_n
: f: P# p% _: N  k+ @/ [4 p
1 c% o6 b" E! E1 j  N, W" ]' H    tsset vtt
5 X/ R9 }' b. ?
   gen ddy240=D1.dy240/D1.dx240
   
  w& H/ R# n; f' y4 _: K   sum ddy240
" L' i) H- X% {  y   
   gen ddy240max=r(max)
   
' L  W+ y- w  Q, ?* d: r5 V   gen ddy240min=r(min)
# C# X( X* D/ M1 e# q# `" c5 h% M     
5 W9 E2 ~9 C, Y   gen ddy240g=ddy240*(dy240max-dy240min)/(ddy240max-ddy240min)
   
3 p1 `. G/ L" l  H: a# x% W   sum ddy240g
+ d" ^1 G" B# F* f' X   
   gen yxis3=r(min) in 1
; w) h$ X$ e# o4 f) E5 a   
. D' M1 b: }& x; y1 _% G! c   gen yxis4=0 in 1
  
  *ddy240g
. _5 F7 R& r( b. p2 N& O   
   
$ w# i3 X6 u+ Q4 I! c# X  preserve
. y; B5 B2 v5 i& Z2 i3 o) ^  
6 j1 o5 s: d0 j, _8 [' g) Z  sum v5 [aw =avol]
8 G/ b; e2 k2 \; ~. O  
  keep if dx240<r(mean)
  
  sum ddy240
  
2 p! K* r: L6 F. A' Z  if r(min)<0 {
        
3 h) S' d! T6 d! }                keep if ddy240<=0
               
                count
               
                if r(N)>2 {
                 
                 pkexamine dx240 dy240
               
3 u. c& ]+ E; u; n3 a& {                local tempminp=r(tomc)
                restore
7 H- f9 R' N1 @8 J7 |5 n0 a                 }
, S: N2 p3 a) R% {  R. V  ]                 
% L7 X8 E& `& q4 |. ]) f        else {
+ d6 J1 B5 y5 ?: Y2 n4 C6 n               
9 ]2 `4 ]" s5 S% {4 a                 di "cannot find minp"
                 
                restore
. B4 j" U+ x: B# O+ u- j- c                 sum dx240
                 local tempminp=r(mean)
                 
                 }
8 ~% V3 l# ], `+ s3 \* u                }
9 M) ?  K9 T' G% R2 \  a   else {
                 
                 keep if ddy240<=r(min)
! `/ v( o; b1 O+ Z                 
% t6 s# J2 j8 s: `$ \" K5 m                 sum dx240
                 
! {3 `/ w% D3 c) b                 local tempminp=r(min)
                 restore
5 G3 ~  a5 e7 g: O+ H2 @$ U                 }
8 F% w4 u- B6 g6 G( B

# k, G! b, V  ^0 h
   gen minp=`tempminp' in 1
! S' `0 X  u  d; E   *240 minp
' V! K% `* E6 f3 x! A

   
* C% P3 z" K2 M; i; t+ H   preserve
   
; K3 F" _3 C3 ]5 _   sum v5 [aw =avol]
, T7 D5 X, l. X! v+ k   
   keep if dx240>r(mean)
   
9 |7 E' b6 n9 f$ e' |   sum ddy240
   
   if r(max)>0 {
2 R6 D: w1 C* t! v% s6 I0 m               
. R  N# ^0 i% r7 S( ?* g! `                keep if ddy240>=0
* s$ T" J) M" E+ `+ d5 d               
# t: g/ w4 u: y. X                count
               
                if r(N)>2 {
                 
                 pkexamine dx240 dy240
' a2 q: O8 b8 m; j               
                local tempmaxp=r(tomc)
' G5 O1 M7 d  h2 I9 ~& `2 r                restore
1 D% Q0 i3 A# ~9 C( O8 `                 }
                 
                 
        else {
) b4 h  V7 V( H+ y. o1 \                restore
                 di "cannot find maxp"
! A( m" L: I' a+ W               
+ ?# P9 H/ ]) q% [               
                 sum dx240
                 local tempmaxp=r(mean)
                 
                 }
                 
                 }
. P$ i" f8 \' J( e. i* D8 o- p                 
    else {
  a% E& q+ X. B  C       
) l. ~( \7 ?. q1 e3 ?7 o) m           keep if ddy240>=r(max)
          
* s8 f6 |( R' A) d3 T( |0 |/ B           sum dx240
          
7 m8 i: |# {( u8 Z; Q  V  _1 G( u           local tempmaxp=r(max)
           restore
. ~& A. C" T, f( L; k6 F9 X. Z           }
       

. E/ B' K8 ]* L0 g, x  

   gen maxp=`tempmaxp' in 1
7 B$ y% c+ q; p, o4 H- r/ A   
9 R3 `! `( }# j- F+ G& q- C   *240maxp
   
6 S! L6 [/ x* ^! F8 G# b   
   *-------------ddy240pminp左侧高峰值-----
- N1 I! z8 l8 O% D. k& p9 j   
3 q' m, n" {# Z3 b1 i8 U, Y/ p$ Y/ ]3 ~    preserve
   
   keep if dx240<`tempminp'
   
- d- \5 S* M2 D   keep if ddy240>0
   
   count
               
                if r(N)>2 {
                 
                pkexamine dx240 ddy240
0 n/ M8 V/ j1 }: R1 p: S9 O' Z               
% x7 ?8 q& w0 Q& W                local tempminphddy240=r(tomc)
3 i& N( K; F0 H( z, D                restore
                 }
) X/ G5 i+ F6 ~1 b                 
        else {
                restore
                 di "cannot find pminphddy240"
, L, F& p, q0 ~6 W. N                 
: R% q3 A& C6 w! i               
                 sum dx240
                 local tempminphddy240=r(mean)
0 v. Q. f' M2 C- t& f* j% e* u               
                 
                 }
                 
   
   
3 I$ U) M6 d2 z5 x' @* H& Q   
: O$ y, D% n5 b. U# T7 @' q   
   
        gen pminphddy240=`tempminphddy240' in 1
       
& e( D% B; V4 ]: i) b! [; X       
        *-------------ddy240pminp右侧低峰值 -----------------
5 H; ]( |7 R- }3 V) l# U0 F    preserve

   keep if dx240>`tempminp'
+ v; z9 f- i' q+ \. j1 l8 k( j   
   keep if dx240<`tempmaxp'
   
. n- z9 z/ i4 z* v7 }( Z" P  d  @4 G   keep if ddy240<0
, @) F2 `! X5 d: L9 J   
   gen temddy240=-ddy240
   
% S) F+ T3 c9 y6 x. N    count
( E, |7 [8 p. K' T& L/ {               
' a' l* w. `# d1 h" X" I" H$ u3 z                if r(N)>2 {
                 
                 pkexamine dx240 temddy240
   
  local tempminplddy240=r(tomc)
; _/ p5 {0 d; R: j                restore
                 }
' n$ v% `( \8 ~) M) |: n                 
) \, ^' @. q; n# Z3 I* \        else {
               
  z- x- }: w& s                 di "cannot find pminplddy240"
                restore
/ O; G: U4 G/ p+ e: [1 A4 [                 
) d& l. h6 P& k2 O                  sum dx240
                 local tempminplddy240=r(mean)
2 r- _5 f7 a+ Z# s5 |# p                 
' R& @" Y8 O; C5 T9 ^9 Q/ c                 
                 }
! Z/ g8 B, C# @1 F   
& \- s/ V4 r6 M$ R2 `) y5 k   
( A9 F6 b' a* v* l% }6 n% U) w   
% C9 ~2 O) v1 F/ m0 A   
- E% Q/ h9 E5 d- V) d* x   
; g$ Q; o: I0 C7 O# A* V, D  l        gen pminplddy240=`tempminplddy240' in 1
9 n8 G& r' b. N' b  l$ S) B, O       
7 Z0 D& C  y4 ~# o' T7 D        *-------------ddy240pmaxp左侧高峰值 -------------------
7 x7 G& h9 y5 S. b6 w3 z! W       
, b8 d5 ], n5 |. ?, y         preserve
  J4 q) S. [- W' t   
+ Z  D1 G6 u2 i' F) s3 m   keep if dx240>`tempminp'
4 ~! p6 S; r$ G5 o: h1 P: Q   
   keep if dx240<`tempmaxp'
   
+ `9 }' l; t& |# A5 k, f4 @0 J   keep if ddy240>0
   count
6 T: ]$ O7 e( z& l1 U4 M( a               
                if r(N)>2 {
; Z3 h" z7 S2 n0 n1 J/ V                 
6 `' p+ e) ~2 W/ S1 j- `                 pkexamine dx240 ddy240
1 K$ L7 Q  M2 U. b6 e   
   local tempmaxphddy240=r(tomc)
   
( N1 I, o- Q, S6 a2 A. s   restore
                 }
3 k! B4 y( ?1 D# I6 V0 ~2 k" @- t                 
        else {
               
; K3 s* o+ H7 N  ]6 M: W                 di "cannot find pmaxphddy240"
: a6 N! Y% A/ o' x7 {5 Y( _4 f  u                 restore
  l9 o3 E0 L. ]3 I. E9 [1 C; E: m               
8 u6 m/ X4 T- `( Z. i* V4 y: D4 W' W                 sum dx240
7 Q* y. q7 S8 _$ ?: B$ Q                 local tempmaxphddy240=r(mean)
* p2 O& N$ h. }                 
                 
                 }
   
   
7 b) [0 K. U# [' f/ z7 D   
   
- g* h9 W2 l, ~   
   
        gen pmaxphddy240=`tempmaxphddy240' in 1
       
       
% p' A" e$ s3 J% r% F7 s/ x        *-------------ddy240pmaxp右侧低峰值 -------------------
, F1 D3 U" l1 y  A  n9 ]       
4 Q- C9 u) {8 Y3 r         preserve
  
* }7 k% n) Y6 ]   keep if dx240>`tempmaxp'
   
   gen temddy240=-ddy240
0 ^) Q! x; f  ^) t, W7 m: r   
" i" W5 P: w* i  W9 I# f   keep if temddy240>0
$ x5 N1 O1 r6 q   
5 g# D' {) D# n" O  z4 X/ {- x   count
0 ], U* B1 r! I5 l               
; X( M/ X( O1 _3 u  o+ D/ B                if r(N)>2 {
                 
* d; I. n8 Y  x% e                pkexamine dx240 temddy240
2 z- o4 k. A3 I$ p. i3 _  c) N   
6 u0 e1 _2 ~3 F   local tempmaxplddy240=r(tomc)
/ v6 q  P& y# q   restore
                 }
1 U5 T, T$ N' _$ q. I% j                 
        else {
                restore
" l- O2 `/ w) {: q                 di "cannot find pmaxlhddy240"
                 
/ x' L8 P) A* J- ^% ]+ N               
                           sum dx240
2 n5 B% {+ ^1 }. Y                 local tempmaxplddy240=r(mean) in 1
, S3 O  r0 l; m' t               
/ [  ^7 v$ M( W2 a                 
4 j% H$ j- @" u! `2 G) s, H                 }
   
# w, j& Y, J3 m5 [) ]4 C   
   
   
6 G+ l- N  y. ~# b/ V2 ~0 U   
* a' ?/ Z0 h& s) M( k1 S, e; l   
$ W1 G" m5 M, T+ a: a4 B        gen pmaxplddy240=`tempmaxphddy240' in 1
       
        *-------------------------------------
   
   gen price5d=v5 in 1/20
      
   gen price5yy=_n in 1/20
   
$ _9 T1 H6 V) h$ c8 w# o   gen price5y=price5yy*(dy240max-dy240min)/20
  F1 q" J7 f, d# q1 `1 ^5 N   
( e; w  W3 ]) z9 ~   *price5y &price5d
6 n2 K6 k( B2 ?   
   gen priceny=price5y in 1
7 n8 X1 b6 N. c, {" T$ o& @( C   
   gen pricen=v5 in 1
6 {; g% C+ p7 g- W* p
, r. t/ o# E; Y* A; `( k% b" e+ o     
   *priceny &pricen
   
   
$ s4 C2 A' W7 o% I2 k   *-------------------------graph------------------
  
) t. D# z& k% @$ s  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)

, P2 i0 |+ z" z. c7 U2 K* ---------可选项目----------------------
*(spike yxis amean,lc(blue))(scatter yxis2 amean, lc(blue) msymbol(+) mlabel(amean))
}
2 ]" A% j7 F. y. Z" pend

maodeng

awer 发表于 2013-2-5 14:40
抱歉，最近太忙了
8 A( ]* q/ J5 X) r0 `) Z
300236 我用下面的程序试过，没问题

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

喜欢喝冰茶

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

喜欢喝冰茶

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

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

汉水云潇

准备试一试
先找Stata软件

汉水云潇

) o( V' |2 K7 q& K, w4 M( o) t
. U+ P: C* x) W2 C8 y/ K1 y晕呀
论坛复制程序代码
; K( W- i- A: R* f" k乱码夹杂在其中
  x( d5 R) [4 i( t8 A很烦
对了两次
& H/ l4 c, w7 p' {+ `以为没问题了
1 r4 p( ^8 ]$ ]) k$ j2 A" p
$ U( A, j: f* G! S$ E7 Q: h结果
执行后
还是出现
r 199
的错误
7 ?% {, m, Y$ X6 f' E, H1 e

汉水云潇

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

peterbev

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

中秋不圆

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

bknite

这个帖子值得学习。

jeanny

感谢！！！

Hubert2000

明月回春 发表于 2011-8-13 20:17
6 b4 N' U( e) ?6 b+ a程序写的不错.
; T9 s  t# o6 ?$ c9 m不过有一点要注意.通达信导出的数据,只对价格复权,并没有对成交量复权.所以数据需要手工处理 ...

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