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

maodeng

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

maodeng

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

maodeng

. hello 300236
statement out of context
r(119);

$ k" ~3 e6 E$ m" F; j; I- Q. b大师 就是这个

sergio

这贴实在！

awer

maodeng 发表于 2013-1-28 13:13
. hello 300236
/ \1 \- L7 @' o9 ^- H7 |statement out of context
2 e& r! |& v, w+ e* v! wr(119);

抱歉，最近太忙了
! j/ q) z& B4 D' U7 r$ F# T
300236 我用下面的程序试过，没问题
* M% I  c% L6 o: O) \+ `
  L+ P! E4 }3 W; ^& U' ?你把 D:\stock\test\source 改成你导出的数据存放的目录，再试一试

  d# P  t! m8 u/ w8 x, Tcapture program drop hello

# D& E2 @( G# k" m# Rprogram hello
6 g' Q8 p/ v* {9 s' R
& t7 S- `. A, P# b8 O: ] quietly {

( o; S0 w: \6 \, l# f8 dclear
1 a6 x, \- G9 R$ t
. e2 j% v2 E4 L2 K" D0 Wcd D:\stock\test\source
9 E4 y" j7 x% ~7 ^4 Z0 N5 E; l5 P
' }$ q1 h& ]3 w0 p% o% H9 cinsheet using `1'.txt
4 @( H9 D& E  _: k5 F  v& ^
- _9 y/ o/ V4 {, x' }   keep if v5>0
9 U2 g# P0 g' x0 B   
$ u( J0 W) A4 B+ Z! s7 O/ C7 D   gen vt=_n
* Z/ ]# z7 `$ ^8 ]9 `   
   tsset vt
! e" Z# n# n9 X# i
   gen zf=D1.v5
$ I# }, @# H% K9 ^' s) R
! H' A: h6 [2 D0 J' w; G! B   gsort -v1

3 K. `! u- n- t3 e" _7 y   gen avol=v7/v5
1 p7 t' H( Y5 X" Q/ W: G# |   
   keep if avol>0
   
9 T$ R8 ], Y4 D& k; Y  z- T   drop vt
* u2 h, C4 c; N9 u  
   gen vt=_n
' E# u& w! p/ B. I3 `3 W" c   
   tsset vt
   
8 N2 B/ ]5 q8 R: N    sum v5 [aw =avol]
: l. g9 d) D' o) \2 M6 L% O9 Q
   gen amean=r(mean) in 1
9 u- N, h" U1 Z3 I' Q- l   
/ \. k+ P/ }; d8 |   * 总加权平均值 amean
   
   
   keep in 1/240

   sum v5 [aw =avol]

   gen mean240=r(mean) in 1
" t0 h2 N" q3 ~$ g5 w7 m   
/ L8 ~. u' R- c- O) ]   *240加权平均值 mean240
% n) i4 H  ]1 F0 ?" K. i6 ^   
( e% l6 p! O3 b  k  q8 v4 D   set obs 500
$ b# U2 |( O5 a) `* Z' p8 I- |0 t   
   kdensity v5 [aw =avol],gaussian n(500) generate(dx240 dy240) nograph
8 A  d2 c( f9 O* G) l   
( i! Y* Q- D# L" R  \   *240 kdensity(dx240 dy240)
3 `: t4 n+ e0 d% T# x   
9 N2 s: v6 p* }* Z; o" G& y   sum dy240
   
   gen yxis=r(max) in 1
- o/ I3 Q4 O/ a( ]! s) v( O7 t: k   
   gen yxis2=yxis/2
! I# H5 z! u5 }% t3 Y+ {   
   *yxis
   
   gen dy240max=r(max)
) @  x9 Z0 h4 |6 W   
   gen dy240minp=r(min)

& x1 T) `2 B% y$ M! l    gen vtt=_n
8 |* S+ |$ `1 s9 p
    tsset vtt
1 B9 S# G7 b6 S" ^9 h
   gen ddy240=D1.dy240/D1.dx240
& e( Y1 P* u5 i, I' J2 X! @! }   
' B& D7 q' e7 N% Y   sum ddy240
3 h9 K3 A8 ]1 @$ J   
   gen ddy240max=r(max)
   
7 u1 {# T# b3 E% u   gen ddy240min=r(min)
     
   gen ddy240g=ddy240*(dy240max-dy240min)/(ddy240max-ddy240min)
/ j2 e. c: j0 |& C" N, c   
   sum ddy240g
/ ?+ ~( {+ g  w/ m! t   
: v7 Z: {* v1 V7 p2 k   gen yxis3=r(min) in 1
0 Z/ E3 [- r. Z$ q   
   gen yxis4=0 in 1
  
  *ddy240g
! B. r- u. I5 G# e3 R: _   
' M, V: F" M; T& `   
  preserve
  
8 e9 B1 [+ x0 }/ I- k  N  }3 n9 L  sum v5 [aw =avol]
0 E2 E# C& J2 a/ @  
  h  K6 V& ?; T& ?! h  keep if dx240<r(mean)
6 h! }$ O2 U: F3 s1 J  t3 F  
  sum ddy240
  
! \" r+ [9 v' v/ n: z  if r(min)<0 {
        
                keep if ddy240<=0
               
                count
: r" [8 E: d( d- Q: C% i               
                if r(N)>2 {
. ^3 d- l1 V2 g2 q" P3 n                 
                 pkexamine dx240 dy240
               
' m6 V7 l/ i; M- l                local tempminp=r(tomc)
! s9 g) C/ J3 L1 j# `! z) L: U                restore
                 }
                 
        else {
               
                 di "cannot find minp"
                 
4 F* h+ ?: \# {  L' ~( }- d# m& k6 M                restore
& j% F% ^5 w! ~7 }' n) |: n                 sum dx240
, J8 o: b4 _8 s7 ^                 local tempminp=r(mean)
% Q* Z' G) T6 t+ y5 C" r                 
                 }
                }
   else {
% V3 V+ S. k) ^7 K- A                 
                 keep if ddy240<=r(min)
                 
4 {; z8 e! k$ ]( V- J& v                 sum dx240
                 
1 G5 C: i8 m  U8 h' |                 local tempminp=r(min)
7 T  c9 A, V" O                 restore
8 ?* Y3 Y5 M# B5 D$ ~5 T, Q                 }



2 q: w; U- x0 y) u2 ?   gen minp=`tempminp' in 1
   *240 minp

" X# n' y' g' p; u3 \) I
& f5 t8 X6 L: a$ @- a6 v# c- \   
- Z, n+ s" s+ l8 H3 u2 k( S3 L1 x   preserve
  T6 Y% `. R0 K- x. O* H+ a   
   sum v5 [aw =avol]
   
4 D* [7 q* S; k. H; p   keep if dx240>r(mean)
, k9 G% F2 z4 ~3 B   
' Y2 U4 k* l9 ^   sum ddy240
   
   if r(max)>0 {
* j3 \) h0 g3 c/ O# D  D8 T; ?               
4 g$ \1 u6 V3 K$ J  Z4 c9 o4 T                keep if ddy240>=0
               
/ _4 V  h( m2 b$ b$ i( J                count
/ Z5 a5 ?) K  ~3 n7 {               
0 S9 c* q! Y% {+ T; E! K                if r(N)>2 {
                 
                 pkexamine dx240 dy240
               
                local tempmaxp=r(tomc)
                restore
* P' {) d; @; x4 L3 f& E7 W                 }
                 
! ^- N8 D8 o$ M                 
        else {
/ t' M; k8 z; R% u, `9 e, l# |$ n, M                restore
                 di "cannot find maxp"
1 V; z+ G4 _4 F- ]+ e  {# f, G  b               
               
                 sum dx240
                 local tempmaxp=r(mean)
                 
                 }
8 {, J4 w1 @$ F9 N                 
                 }
8 B% I" U7 g$ V. Z& c                 
7 H& A* \1 l" T1 @' [    else {
- M! u9 z6 @8 x9 y0 }5 }9 q. J1 {  a       
           keep if ddy240>=r(max)
5 ?6 w. A2 ^- d* z' m% J          
9 r# m" S+ k8 E- x           sum dx240
2 O; z) K. j9 g' n  U; T: J5 a5 L4 l          
           local tempmaxp=r(max)
           restore
           }
% y& ~5 X) h( [: t' ^  P       

  
6 g7 C5 q! v- a# g, b' W
1 g5 u. c' F0 }8 f3 I   gen maxp=`tempmaxp' in 1
   
# h: t( o( }) S7 ~1 i; m/ O- F   *240maxp
' ?+ l: @! t1 L4 R) t   
   
9 @5 Q, j8 ^0 X1 T9 h$ q8 d   *-------------ddy240pminp左侧高峰值-----
   
    preserve
( T# M% @4 ~1 U- Q# t   
   keep if dx240<`tempminp'
   
( m/ \* E0 Y; ^9 H; l% ~8 c# ]6 E   keep if ddy240>0
4 }* Y& e  x4 H# s   
- c; S0 b. z$ P: W' e   count
               
' i/ a4 F8 A8 e+ w0 F) i& z$ G                if r(N)>2 {
% T4 F: x' S" }2 S# ^( F+ V                 
                pkexamine dx240 ddy240
               
                local tempminphddy240=r(tomc)
                restore
                 }
                 
        else {
1 {3 K0 g; }0 Z+ a                restore
                 di "cannot find pminphddy240"
                 
               
- e4 {5 {  Z: x6 p$ Q                 sum dx240
                 local tempminphddy240=r(mean)
! \& t: |* B' x: g               
) {" ~$ u7 m( B$ s0 [' ?                 
                 }
7 a0 }7 i" \, L                 
   
1 i$ }( ^7 D4 Z, z. N2 I. I   
" Q9 v0 j% m) k   
   
   
9 M: J8 W) l# k# G4 |# V2 U1 }        gen pminphddy240=`tempminphddy240' in 1
       
       
9 [$ ^! b" j+ V: B4 q' u, A        *-------------ddy240pminp右侧低峰值 -----------------
    preserve
7 l( M: L, w1 t' h
   keep if dx240>`tempminp'
% N2 |5 f. s& Y4 \   
   keep if dx240<`tempmaxp'
   
   keep if ddy240<0
/ \  O7 z( K  i( `' @7 l& ^   
   gen temddy240=-ddy240
. o1 r% |" e  H# i$ d# i   
    count
9 v; r- \2 N- d" z               
                if r(N)>2 {
                 
                 pkexamine dx240 temddy240
   
1 [: h( b5 K, y  local tempminplddy240=r(tomc)
6 o+ `5 d0 Q8 w& h- x/ r                restore
                 }
                 
$ |. t6 Y+ u: h0 e        else {
               
                 di "cannot find pminplddy240"
) ]4 J( I# M2 d. |3 A/ z& Y                restore
5 L. v3 K5 b) T7 [9 r8 d8 Z& G1 u                 
                  sum dx240
                 local tempminplddy240=r(mean)
                 
, X: v0 ^! o5 U( h3 |) x                 
                 }
! R% j8 ^- }$ \   
   
   
$ k# |9 E, I2 P0 W   
% A' B/ Y: e7 Y1 j, C   
        gen pminplddy240=`tempminplddy240' in 1
       
+ l) Z1 Y0 q, s# A+ d" Y% J# y5 G' h        *-------------ddy240pmaxp左侧高峰值 -------------------
       
; y: y* K7 ~2 b6 S8 J; c6 ]         preserve
0 v* _/ N0 C  O- _   
   keep if dx240>`tempminp'
   
   keep if dx240<`tempmaxp'
- @$ q) t5 H, ~7 A/ \   
; s; p" N) j5 y6 m" p. x* S2 f   keep if ddy240>0
   count
               
* O# U2 z6 D" g                if r(N)>2 {
; ~% m: H  ~# i                 
                 pkexamine dx240 ddy240
$ l0 D# H# \& ^6 e/ W1 i+ |   
   local tempmaxphddy240=r(tomc)
7 _. y0 q2 k6 I' w   
   restore
                 }
                 
        else {
               
- I$ X! q$ W) P( f) C8 O" n6 q                 di "cannot find pmaxphddy240"
0 ~0 l4 ?$ ~+ Q                 restore
' `9 B2 a: ]; c# _; ?+ w1 o3 j               
( x' U- \/ e7 N0 t. X                 sum dx240
                 local tempmaxphddy240=r(mean)
                 
, _- K- I- q* I7 z                 
                 }
5 j7 _2 S- G- L( k6 t   
   
: {: y: K' M6 D* a$ g   
0 b* y0 O' ^, _6 S8 K4 W   
$ S. X2 `* E# F( c  G   
   
        gen pmaxphddy240=`tempmaxphddy240' in 1
  f8 P! L$ g  H4 F( E8 ~8 r       
       
( r8 B3 n9 Y, j        *-------------ddy240pmaxp右侧低峰值 -------------------
       
         preserve
  
   keep if dx240>`tempmaxp'
! p- t+ t- M8 L( g2 X3 |   
   gen temddy240=-ddy240
   
   keep if temddy240>0
   
6 P$ a& K' R/ g4 B1 b   count
% M/ V1 c7 D; ~% b  m. \4 X) W               
3 }. f3 B- Q( e$ l% O) z5 n2 [                if r(N)>2 {
4 q) V* K; \2 o6 _( ]                 
                pkexamine dx240 temddy240
6 C+ Z2 ^0 [# v6 {   
1 N5 ?+ ]0 K/ _5 y   local tempmaxplddy240=r(tomc)
   restore
# g7 [' c- A) m                 }
5 m  K6 Q% o  i- r6 U% y                 
        else {
2 ]: ?6 w) y3 X6 y6 P- S                restore
                 di "cannot find pmaxlhddy240"
                 
5 a" ]. R) u6 r) V4 i5 {               
" o$ W' E  l. H+ ^                           sum dx240
0 K& m. }, I3 m                 local tempmaxplddy240=r(mean) in 1
0 Y6 k: j0 b# t( i/ h- i               
                 
( Q. H$ n0 J  }0 s                 }
4 J' u+ _. }" S1 d7 v3 W+ V   
5 Y5 |2 s; Y" k4 o" F, s7 G5 f) L   
   
   
, s$ J2 \; F" y3 ^0 k+ l, s   
   
! `0 o1 n8 f- H5 [        gen pmaxplddy240=`tempmaxphddy240' in 1
       
3 g8 c  Z# I" ~4 |        *-------------------------------------
   
5 t4 Z) M* H! `# @   gen price5d=v5 in 1/20
) O9 Z9 L' B7 N6 Y: A/ \      
3 \3 g; d" l! e: n   gen price5yy=_n in 1/20
" w  {; ~0 [9 A# z: E; B! o   
5 u2 u: N- y% L+ L& P! x$ J3 x  W# }   gen price5y=price5yy*(dy240max-dy240min)/20
, \$ r4 H6 f2 g1 L: F   
   *price5y &price5d
" e, B/ P9 v2 P3 }4 P   
& w' S; B9 u# g) r) x   gen priceny=price5y in 1
. d" t  P- h1 T9 ~. |' z   
5 i# @% ]& a6 x& j) _, z: j) F   gen pricen=v5 in 1

     
4 R* h1 i4 k: ^5 S: [9 Q   *priceny &pricen
   
   
   *-------------------------graph------------------
% [9 ~& I/ p2 u  
' v  N1 d7 w$ K0 j  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)
% E6 N1 @& a- V3 _. w7 B
6 u' \/ E! x( N* ---------可选项目----------------------
+ S+ t) M% A% F2 |*(spike yxis amean,lc(blue))(scatter yxis2 amean, lc(blue) msymbol(+) mlabel(amean))
- e8 E: s" }! P- H* c}
end

maodeng

awer 发表于 2013-2-5 14:40
+ c# p  @3 x8 R" v6 |; v' e/ v抱歉，最近太忙了

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

8 F- D% X6 T, T# z4 B哇，谢谢，我仔细看看，这些字母数字代码晕了
: `/ R) i+ S; `) v

喜欢喝冰茶

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

喜欢喝冰茶

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

0 X# K9 x7 }! c- P+ ^9 `" T  o
这东西是开源的？知道有个统计软件叫stata，不过要收钱。

汉水云潇

准备试一试
' x% [, N6 T8 M2 g5 Q. u% r& y6 z先找Stata软件

汉水云潇

晕呀
论坛复制程序代码
+ o4 z/ U* C0 u1 }, A) ?3 X+ a乱码夹杂在其中
很烦
对了两次
3 q: P- T: C+ x以为没问题了
2 K( x$ g0 k; Z1 e" p
$ d( a* w; ^: y; z5 k结果
- S6 N2 p. M& P% Z& I执行后
/ p& O& n. G' x还是出现
r 199
, g- j; S% j5 H/ p的错误

汉水云潇

试验成功，下面是三一重工的图
9 y% q+ e% U3 [$ d
* ^9 k! g2 B3 ]

peterbev

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

中秋不圆

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

bknite

这个帖子值得学习。

jeanny

感谢！！！

Hubert2000

明月回春 发表于 2011-8-13 20:17
4 B( ?% C+ A- g! _程序写的不错.
+ m+ d4 _" Q) n) n' f3 \: v不过有一点要注意.通达信导出的数据,只对价格复权,并没有对成交量复权.所以数据需要手工处理 ...

9 S& E3 k4 b2 G& o; s老哥，也是山寨这篇文章找不到了，不知为何？