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

墨轩书童

( I+ a& X* a& ]1 a5 E
感谢楼主！！

海阔

( Q: E6 I, T! f& `
按照版主的步骤操作到了最后出现
& i+ c# v$ v4 L% J6 ^
, T6 h/ Q/ E7 f5 Q7 ]1 z/ i' n- l. hello 600031
avolDo not found
) ^/ c, o: ^5 }& ]" d6 Y2 J& A+ I9 C# P---------------------------------------换成其他代买，也是一样；总是出现这个问题，怎么办呢？奇怪

maodeng

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

maodeng

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

maodeng

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

maodeng

. hello 300236
! X7 ?+ i  h+ _7 Q4 v1 K/ I2 `$ J3 Xstatement out of context
( }8 H' |/ j0 S. |r(119);

' u1 I% O% F. ~0 c  }# T9 ^# m" b$ x大师 就是这个

sergio

这贴实在！

awer

maodeng 发表于 2013-1-28 13:13
! J& |# x% _1 Y0 N! S. hello 300236
5 V( |% S3 c& g+ Wstatement out of context
r(119);

; {7 w& N- ~% V( |抱歉，最近太忙了
, m0 F  H/ j  j3 W! O) t. C% _
) K5 _! A) w. D2 S6 a5 X300236 我用下面的程序试过，没问题
) b0 T/ _/ x( X& s* W* L4 G
3 g2 Z! P; A: R6 [* l/ M! w8 Z你把 D:\stock\test\source 改成你导出的数据存放的目录，再试一试
" C2 b. e( n: i2 z
6 a) g( I. y/ P. r: t) kcapture program drop hello

) b. [1 p! A# e; Z. J$ k- O* Tprogram hello

quietly {
( W, A* l3 M3 k. A
clear

cd D:\stock\test\source

! ]6 o! i' i  _5 v& j. r/ I) C' \insheet using `1'.txt
- w- z) U. d7 y' R$ A! s3 k
6 H7 X# M' o% h+ _7 z   keep if v5>0
( S+ {2 ?4 f$ X5 i) V8 U   
. U9 a9 t& n, C( P   gen vt=_n
   
9 }( k6 [; K9 f4 ?7 y+ z) Z0 q   tsset vt
, H+ l( `$ v) f( P* U
   gen zf=D1.v5
) M0 c' M; X- ^: r+ {
. X& b+ d8 w9 y* Y* k3 Z   gsort -v1
) f9 S" b, F$ j% G
) E! Y0 c# n5 ~  o   gen avol=v7/v5
   
' D. ]+ U8 D/ \  J* m   keep if avol>0
   
   drop vt
  
5 W# p$ x8 ^+ c   gen vt=_n
8 f1 C! Q' L6 z% @   
   tsset vt
. x5 A2 ^5 j6 L. [4 R! a   
, z* s) j2 d6 Z5 h  `% |$ ]% K2 b9 [    sum v5 [aw =avol]
- f- Z9 V$ P- ^, J
   gen amean=r(mean) in 1
   
   * 总加权平均值 amean
+ R; G0 m6 E; m% p0 Q' K+ F" }   
+ |! v+ S6 i* p, y" }+ O   
  O  p1 J  p& o! ?6 U( z5 R/ @+ Y3 R   keep in 1/240
- k" e$ H# m) Z5 G7 g& ?+ U+ }  H
   sum v5 [aw =avol]

. b% S" n( a# U4 D( _4 I   gen mean240=r(mean) in 1
   
& A& v0 ~: P2 V. R# q& s   *240加权平均值 mean240
/ W% M- M; A. Z6 _; e% n1 p   
/ A0 S* D# l, `2 J( F   set obs 500
) f& v7 M" S' P   
   kdensity v5 [aw =avol],gaussian n(500) generate(dx240 dy240) nograph
   
   *240 kdensity(dx240 dy240)
   
$ c8 b: S5 x* ~% k: E, h5 y8 Y/ u. w   sum dy240
   
3 |# M! S' q6 j" |6 A   gen yxis=r(max) in 1
! B- s( n" t" N) t& F2 Z8 B3 J   
   gen yxis2=yxis/2
   
   *yxis
, `. c& `. z) V6 G4 Z* r9 W   
   gen dy240max=r(max)
! {( x6 y1 o9 C   
! Y/ [5 ^* l. m, x& k8 O   gen dy240minp=r(min)
+ C" B8 v  _& U/ j: K
; z0 N5 X. Y# y# {    gen vtt=_n
2 @; j* m/ b  O6 ?8 H$ K. E5 s
    tsset vtt

   gen ddy240=D1.dy240/D1.dx240
" E) m& p# z; Z9 g$ n5 [  o   
0 c$ o5 F4 A& w! y2 k! H  a   sum ddy240
" F0 T3 a9 T1 l   
   gen ddy240max=r(max)
   
   gen ddy240min=r(min)
     
   gen ddy240g=ddy240*(dy240max-dy240min)/(ddy240max-ddy240min)
) ^/ E) F* m1 W3 z   
   sum ddy240g
# ^, p  y& ]4 z8 s3 l# K   
0 u, ]8 K. N: f   gen yxis3=r(min) in 1
   
1 R( D& B0 r. ^  o# z) H* s) z   gen yxis4=0 in 1
3 V; _0 C# l0 e/ v7 P3 H  
  *ddy240g
   
   
  preserve
  
0 J+ J$ U; K& l2 h9 M  sum v5 [aw =avol]
4 F; P" i: y& X. W1 Z0 P* M& Q1 ~/ b  
$ c3 M* b$ R" c: V8 g; S* Y  keep if dx240<r(mean)
4 i- B+ M. a# E, C6 a  
  sum ddy240
7 Z: u% q6 n1 p8 q5 v4 M  
  if r(min)<0 {
, g" D3 f) {; H; I" `, h        
1 l* i) [4 J) Q) k$ G/ z8 R! ?                keep if ddy240<=0
               
5 N+ Q+ ]: z" k+ N- z                count
               
                if r(N)>2 {
                 
                 pkexamine dx240 dy240
               
# }  G" c5 J  m- X( t1 _& c                local tempminp=r(tomc)
                restore
                 }
4 B4 e* J4 ~+ S. w# R7 O2 J                 
4 B7 l% P7 ?1 H, ~$ t9 t5 E' O        else {
               
2 t+ s" i$ d+ a  G                 di "cannot find minp"
                 
                restore
1 s& m' y3 C. W- P4 R% }* i                 sum dx240
/ x" N6 G* P% f+ o. i  l6 K                 local tempminp=r(mean)
                 
/ O. R/ N% H# n) P9 i8 j6 K                 }
                }
   else {
                 
                 keep if ddy240<=r(min)
                 
! G5 A' o0 M4 m                 sum dx240
( j  y, L' E  m# T9 p                 
1 q" O2 v& x+ V# M( ~                 local tempminp=r(min)
                 restore
3 }% @) Z% @' R5 J" ]/ ]: s                 }
3 Z: ?) ]' n- j4 {1 f

: r- ^" W8 k) A, L) Q
   gen minp=`tempminp' in 1
" |& T/ a& c$ l" J6 D# n   *240 minp


   
   preserve
   
   sum v5 [aw =avol]
4 Y7 i; {9 Y( Q   
   keep if dx240>r(mean)
2 w8 g: F5 S, `2 q) u% c   
   sum ddy240
+ W. L/ ^4 X  a7 a/ p   
   if r(max)>0 {
               
                keep if ddy240>=0
. u& c* e3 o* d6 u               
                count
               
; }( O+ e# E4 ^0 y5 E' p                if r(N)>2 {
$ y. v# B2 R/ L8 k' V+ _4 _                 
+ h8 v( w% w! Y; R) R: Z                 pkexamine dx240 dy240
6 L& p9 F  ~: l) q6 \               
6 a- k2 E) d1 H( s1 D: w8 }& |$ p                local tempmaxp=r(tomc)
                restore
: s& |  C" f, ]4 Z' R+ ^* q8 u                 }
                 
                 
        else {
                restore
                 di "cannot find maxp"
               
               
! V) ]4 f- E  O                 sum dx240
+ x3 u) X6 h: ^5 Z, ^                 local tempmaxp=r(mean)
                 
# B9 [1 J7 z1 M% Z0 N' s7 ]                 }
! b+ G& ?5 ?9 b                 
                 }
                 
! V$ j& n5 o8 y6 w! \1 V0 ?    else {
       
           keep if ddy240>=r(max)
          
% t: a* J- I& n- e           sum dx240
' B9 Q& ]' p4 X/ f# K          
           local tempmaxp=r(max)
. S% F% f7 p* H/ Z           restore
( s1 `  ~: H6 v( y           }
       

9 K1 m' z4 s. I* D6 X, x; F  

7 c) y* W, u7 I5 u3 {& v   gen maxp=`tempmaxp' in 1
   
" G/ _  h. h  {3 o& ]' a5 D   *240maxp
) i% a; L0 L5 A% S; f9 }   
   
   *-------------ddy240pminp左侧高峰值-----
   
    preserve
4 C/ S3 K* [& ?/ _' X   
   keep if dx240<`tempminp'
( k# d& z4 V1 A/ a   
   keep if ddy240>0
   
   count
! _% O# b6 K$ _2 I1 x3 c( `               
                if r(N)>2 {
  c) K# f9 Z& f* V* P7 X                 
                pkexamine dx240 ddy240
               
) _% O% p4 [3 l1 n                local tempminphddy240=r(tomc)
4 b4 M0 [7 y+ h( @                restore
                 }
+ K8 |3 M+ Y6 S0 U                 
9 x& f+ V, M' Z+ c( A* B/ f        else {
% I) f* {: ?$ o- K: O                restore
                 di "cannot find pminphddy240"
                 
/ k3 J  ?: [" ~0 ^, v2 U% a               
  L' H6 Z6 F$ l. \6 Q                 sum dx240
1 ?( H$ H1 W; D3 a                 local tempminphddy240=r(mean)
               
                 
                 }
1 i/ U, n0 C6 S                 
7 U' L2 y! f: L/ o% Z   
: q3 o0 x) k: q/ R/ c( c   
   
   
   
        gen pminphddy240=`tempminphddy240' in 1
* M) {2 Y6 B' h3 C+ D/ E       
7 C. o) e- z. m% s; {; I+ ^2 O; T( [       
        *-------------ddy240pminp右侧低峰值 -----------------
: D9 c" ^) ~3 Q1 ]    preserve

" w/ C6 O+ D+ V& c* F3 u) D9 y   keep if dx240>`tempminp'
- _0 {6 j( b0 H* i   
   keep if dx240<`tempmaxp'
( {) s% s, k+ N9 @/ c% m   
   keep if ddy240<0
   
  K1 H& R9 i4 Y" V( b   gen temddy240=-ddy240
- i& o& j5 L" c* j1 R9 ~; B" N   
2 a' N0 v0 h0 F6 k- j    count
# ^5 @- k' x& K0 L  B9 w               
/ a1 A( E0 C& Y) {3 k) R( ^                if r(N)>2 {
                 
                 pkexamine dx240 temddy240
   
4 a- k1 E+ f# e( G4 q4 r. B  local tempminplddy240=r(tomc)
                restore
                 }
& @, L0 R. ^' t' G0 m4 y                 
& w: c0 Y, R, n& Z9 \        else {
3 K4 s5 v1 _+ {- _; F2 ]0 E               
0 S! w: i* w* {- u                 di "cannot find pminplddy240"
; [& [1 a! i$ G- T                restore
                 
                  sum dx240
) P3 H1 j" s3 `7 F7 c3 w& j3 P& O5 I                 local tempminplddy240=r(mean)
                 
                 
2 j3 A+ H6 d0 y! J                 }
, S* t, C+ v, |0 d6 Y# g: k: a   
   
   
   
   
+ }+ `- Y+ y- j/ k        gen pminplddy240=`tempminplddy240' in 1
& Y0 g! S; E" e       
        *-------------ddy240pmaxp左侧高峰值 -------------------
       
         preserve
   
   keep if dx240>`tempminp'
   
9 ~6 K2 ?: j- p* C1 {) K   keep if dx240<`tempmaxp'
   
2 t1 `; V( f  F) k% e+ w: ^   keep if ddy240>0
. e/ O: b$ K, l4 V   count
4 [" M" C" t  g/ z8 j4 a. T               
                if r(N)>2 {
                 
                 pkexamine dx240 ddy240
+ o& q2 y8 A0 |' X% V   
6 r  X! Y5 d% _8 d   local tempmaxphddy240=r(tomc)
   
+ C$ P0 v$ A# B/ f# M* D6 C   restore
                 }
  O. e& G, P" S4 Z                 
        else {
               
                 di "cannot find pmaxphddy240"
! h$ {% H* y6 v- v                 restore
               
0 n! b5 Q9 b+ J  i                 sum dx240
) L  H6 E& [& ~6 N+ X                 local tempmaxphddy240=r(mean)
                 
                 
                 }
   
   
   
5 I" j$ \$ f5 e! u0 T   
! Q% s3 L4 W. [) s" z   
( ^. |$ a: ~4 d; [( W   
        gen pmaxphddy240=`tempmaxphddy240' in 1
* t  J5 G4 B0 E! T- z1 L       
       
" f' l. }: S1 x$ l; Y        *-------------ddy240pmaxp右侧低峰值 -------------------
       
: c+ ]. Y( i' Z, ]4 _+ R8 i         preserve
  
/ s4 A" [& {) J  B, B   keep if dx240>`tempmaxp'
   
   gen temddy240=-ddy240
7 o. D! J( ]9 ]8 V   
7 b. Y$ s# a# J   keep if temddy240>0
   
   count
               
3 t; w6 R5 [. v2 \1 R                if r(N)>2 {
  }# l! W* T) I5 n6 z& ]  j" Y; W# ]; x                 
                pkexamine dx240 temddy240
   
   local tempmaxplddy240=r(tomc)
+ A4 g: o. D4 c9 n3 n7 e   restore
. j# i* U7 D, S- |                 }
                 
7 S% r8 ]1 E# U6 K3 x        else {
                restore
* R' C0 m9 z' d                 di "cannot find pmaxlhddy240"
                 
' N/ F2 e! M( O: l9 D( k               
$ n  f) t. j8 h" a8 u                           sum dx240
5 ^- e3 T  c+ z( e* B                 local tempmaxplddy240=r(mean) in 1
& w1 T  X( @4 D: ~' _7 j! i               
8 B8 }# y! f3 A' n) c% c9 m7 w$ o                 
                 }
1 |) f( k2 r# ~   
$ I3 b0 B9 F: Y  W   
6 \% L6 V% Q+ d2 V/ w   
   
' J. m% j; B& h0 k   
9 _) A  K0 ~6 R7 N   
" s: J; ?5 s2 S. Q/ Z3 d0 @, P        gen pmaxplddy240=`tempmaxphddy240' in 1
       
0 D% v0 l3 u1 E) _# J, j        *-------------------------------------
   
   gen price5d=v5 in 1/20
6 B4 _2 m( Y- W8 w" C; u: e' h      
   gen price5yy=_n in 1/20
   
( B' T: z) ^) D' ]; J- L& i   gen price5y=price5yy*(dy240max-dy240min)/20
   
2 I/ \* w. t" }3 i6 B: r+ ]   *price5y &price5d
% \. `) R: Y  \2 f/ P   
1 `' F) g. G9 ?$ a. Z* P   gen priceny=price5y in 1
' U1 V" b+ q6 E7 E) D   
' I  H7 [! P* I- D   gen pricen=v5 in 1

     
1 i# Z4 w/ C7 E$ F   *priceny &pricen
   
3 q- `' a  U: Y" G3 f   
- L& h9 J1 n: d  D+ m   *-------------------------graph------------------
  
6 y6 u/ c' C$ U) ~6 Q2 n+ {, x7 l  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)
& ^- M  c0 Q3 c  i8 y/ X$ n
: a" s, J7 Z! M" N8 j& u/ a* ---------可选项目----------------------
3 f0 o. k1 J, {6 e- m*(spike yxis amean,lc(blue))(scatter yxis2 amean, lc(blue) msymbol(+) mlabel(amean))
}
end
% S$ Q2 w9 Q$ B* B! G

maodeng

awer 发表于 2013-2-5 14:40
抱歉，最近太忙了
' A  |) ~) L) ]4 p! q) ]3 U$ g$ _
300236 我用下面的程序试过，没问题

哇，谢谢，我仔细看看，这些字母数字代码晕了
" I5 c3 M7 K  b4 w: h( n

喜欢喝冰茶

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

喜欢喝冰茶

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

1 D& _5 b; K) w8 g) ]7 @
5 M& y& j( [. l" d& c9 ]4 p这东西是开源的？知道有个统计软件叫stata，不过要收钱。

汉水云潇

准备试一试
先找Stata软件

汉水云潇

晕呀
! r) |! ]1 e$ N. A, p/ E论坛复制程序代码
乱码夹杂在其中
: n: `0 x" b7 @4 p3 F很烦
' L, Q7 ^0 U! b5 `. t对了两次
以为没问题了
$ h% x( I9 c+ p+ f& F% Q9 X
- \. [; A5 \+ v4 S) _. }$ ?结果
" T& C/ H/ A( v3 B0 \  X8 g执行后
. T# B6 a, y! L! Z! q还是出现
& B- G* F- s. a$ t8 G/ t  A, Tr 199
的错误

汉水云潇

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

& h8 A$ u6 k& k

peterbev

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

中秋不圆

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

bknite

这个帖子值得学习。

jeanny

感谢！！！

Hubert2000

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

( g( k1 s" F. Y4 s' B# q老哥，也是山寨这篇文章找不到了，不知为何？