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

maodeng

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

maodeng

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

maodeng

. hello 300236
  s0 }. p! \. [; I& H: [9 ^statement out of context
4 L: O6 E0 T, \# }( Xr(119);
* f3 g1 B# K9 M
$ v  w( n; D( I" }$ I/ _大师 就是这个

sergio

这贴实在！

awer

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

8 a: i( B) R& F, Y* U5 p+ B. O抱歉，最近太忙了

300236 我用下面的程序试过，没问题
$ F+ f/ Y/ ^# @/ X% L5 F
3 m+ p, y6 `% S3 u6 d你把 D:\stock\test\source 改成你导出的数据存放的目录，再试一试
" x0 i& D5 d+ ^/ Q- j
capture program drop hello

program hello
7 G8 f" O6 \; h( d2 Z
quietly {
" r. p6 r9 X- n/ w+ t
3 ?! W/ d4 g: h0 X; |; D: m# Yclear

cd D:\stock\test\source
9 y0 P( w; t" h( \( D/ v
insheet using `1'.txt

   keep if v5>0
: D8 H: d# C2 p   
$ N, o8 R! n$ f! `" }3 D. B   gen vt=_n
   
   tsset vt

2 T" W2 L+ B+ G9 Y0 }   gen zf=D1.v5
+ U3 R7 {- X: V+ Z( L
   gsort -v1

; k' V, U# [5 D   gen avol=v7/v5
   
. K% o9 t: h- @# t   keep if avol>0
- Q7 ~9 n7 T( i* X* P   
   drop vt
: y! K: X/ f5 W* R5 M; w' D" I  
   gen vt=_n
' ~$ A9 X" @. ?3 b2 D- {   
   tsset vt
   
    sum v5 [aw =avol]
  V+ x5 D$ O8 A+ r& T8 j$ K5 H6 m
   gen amean=r(mean) in 1
   
   * 总加权平均值 amean
   
7 F- n1 U3 o4 h& M, v   
   keep in 1/240
- i  M; E' {  G  K# `: _. x! a
* m1 T/ T1 p1 W. S0 M: U   sum v5 [aw =avol]

; a( G# v$ ]: E" F! j; ?   gen mean240=r(mean) in 1
   
4 j* L7 f5 b+ T; I$ O: ~   *240加权平均值 mean240
# s6 O0 J5 d  t2 N   
6 [/ Y8 l6 E' v) Q! n: [% Z: |   set obs 500
; k8 e+ k+ r! i7 T. P   
6 |: _! H2 x. P, x6 e   kdensity v5 [aw =avol],gaussian n(500) generate(dx240 dy240) nograph
   
/ u% y+ g4 A' A- v, a   *240 kdensity(dx240 dy240)
, b0 E) S: U# D8 R# O$ ]   
   sum dy240
   
   gen yxis=r(max) in 1
   
: R7 n  C1 X+ u: s( W0 t. l   gen yxis2=yxis/2
3 l$ T) M! X' Z( [  Q" T9 H   
+ r  w% Q( a- v- z/ u9 }2 o   *yxis
   
   gen dy240max=r(max)
7 L' L9 T- z# i. I- B   
' m8 }, K" O4 P# T   gen dy240minp=r(min)

( E. w, F( R3 B' v: Y    gen vtt=_n

    tsset vtt

   gen ddy240=D1.dy240/D1.dx240
   
# D9 h7 ^/ @/ H   sum ddy240
7 |  e3 u9 B7 B! C   
   gen ddy240max=r(max)
# L! t9 H. h2 ^. X   
   gen ddy240min=r(min)
     
   gen ddy240g=ddy240*(dy240max-dy240min)/(ddy240max-ddy240min)
   
; S1 a9 Y' m+ q0 _3 [   sum ddy240g
   
, v' `- p5 i* ~  ?+ b   gen yxis3=r(min) in 1
6 z9 \0 [1 p; o% I2 ]   
   gen yxis4=0 in 1
  
+ W( K( z* _* r  *ddy240g
/ X4 [% A# d7 o2 I) O   
   
. r* @5 ^, H, T; i% r# ]7 p5 J  preserve
  
  sum v5 [aw =avol]
" t- U! F) a, p  
# c4 U7 y# l) {( h/ @% s  keep if dx240<r(mean)
  
* x, F  S' z. ]9 `! x+ b  sum ddy240
! x; i3 g8 @0 P( p/ S4 q( `  
  if r(min)<0 {
        
                keep if ddy240<=0
  W$ d& p5 B) ]2 L               
                count
               
6 b7 Z/ d: s2 w                if r(N)>2 {
                 
3 w7 N0 p8 ~* z" n4 [% A                 pkexamine dx240 dy240
               
                local tempminp=r(tomc)
# L% a- h4 N6 W% X3 G  Y' k' \                restore
* Z, E: z2 p: @9 \                 }
6 p2 ~2 d) I6 K. b9 H                 
        else {
               
9 }# A" m* S. C, s$ J7 [                 di "cannot find minp"
                 
                restore
) ~- ~+ L, R6 S6 K, {8 T                 sum dx240
                 local tempminp=r(mean)
8 C9 ^0 L# s# Z& G+ X5 t                 
% S% [( B) O" A1 ^                 }
  o- h" e8 J  M- o/ i  o. i4 q                }
+ I0 v" Y6 i0 K! Q% B2 N/ S   else {
, s0 \: j! _/ z1 ]* j0 N1 i* e                 
9 Z' v: q8 Q: n  P6 A+ H( f                 keep if ddy240<=r(min)
9 j7 I/ z! q4 h& F! _+ r                 
* j+ H% _8 w( o' c+ \                 sum dx240
7 u* D' I3 P" x3 r: r- l/ Q                 
7 D% _' p8 `& I! D& {1 ?                 local tempminp=r(min)
: |9 P8 x. W7 U/ l% y/ K- |( |                 restore
7 N. ^3 A. y+ `: p/ a/ B- G                 }
$ H0 \4 r' {: r, }
) \' }) T0 }# ?8 ~) ~3 t

: ]3 Z! S& r: r$ I" c  `. l, T! Z   gen minp=`tempminp' in 1
   *240 minp
; ]* z2 U: r9 w  j) q9 w  V
# i- q) W, R1 `1 e
1 `: d% U- W; X   
   preserve
   
   sum v5 [aw =avol]
   
   keep if dx240>r(mean)
7 d) N) ?4 ~! [7 W   
7 A6 E) d4 O9 @0 |" b' O5 O; d   sum ddy240
9 q8 A0 H5 F5 P( N& {   
; m# E' t& v6 w7 n3 {# o6 w   if r(max)>0 {
               
                keep if ddy240>=0
3 z1 [% v6 o1 W/ Y* q5 R0 o; h, W4 Z               
7 O1 A- B, N4 B4 A- x7 f1 m# d- ^4 P                count
               
                if r(N)>2 {
                 
                 pkexamine dx240 dy240
- ]& q) P: P. u& Y# G/ ]7 a: Z               
) E5 z, b# B/ B5 m. D+ C                local tempmaxp=r(tomc)
                restore
! g6 q* _8 ~7 [                 }
# W2 J# [0 a6 X. X! C1 v                 
                 
6 {2 N: i+ a; c- S  N4 P( F        else {
/ o5 j( L  a  O                restore
3 O: t8 m+ G' ~/ ?7 M                 di "cannot find maxp"
               
               
6 Z' f8 |/ w  a/ u                 sum dx240
+ u5 O6 L6 E3 ?2 m6 z2 a( e  Y/ [                 local tempmaxp=r(mean)
                 
, D& X. s2 J: x& g                 }
4 @! N9 F# d& ?: u/ n9 ~                 
                 }
2 u3 W# j2 z4 O2 q                 
    else {
       
           keep if ddy240>=r(max)
          
           sum dx240
6 ?: r: m' |- s          
           local tempmaxp=r(max)
           restore
           }
       

  

" s1 G' P4 B5 ?$ G: e  d+ A1 [  S   gen maxp=`tempmaxp' in 1
4 o8 [% K0 p2 i. ?5 d: v   
6 K$ u% J2 P- e9 N   *240maxp
' s& A& `$ H7 ?  `   
   
( D) T9 B, F8 C/ F) d4 I9 A   *-------------ddy240pminp左侧高峰值-----
   
    preserve
  s# w7 C+ K* K) P  |9 L0 F) e   
& {% k3 _3 J5 l5 V   keep if dx240<`tempminp'
/ }8 V4 H4 U4 T, q! K  B* u   
   keep if ddy240>0
   
, }: l6 f2 }* @; X1 p) j1 S" j   count
               
                if r(N)>2 {
$ h8 ~, ~% L8 f6 C                 
6 [" k" y* F1 [% ]6 S& C; {" g                pkexamine dx240 ddy240
2 S' ], [% ^& s; x               
                local tempminphddy240=r(tomc)
# J9 A/ `+ F- L                restore
6 L& S- N8 j% h3 L                 }
                 
8 _! D5 S* m& T/ W& V( w        else {
                restore
7 w0 ~  N: j$ j3 ^- {, I+ H                 di "cannot find pminphddy240"
                 
% `4 o2 h0 P6 l6 z; X               
                 sum dx240
5 r) [5 }6 z2 P. t9 i1 K                 local tempminphddy240=r(mean)
               
                 
7 p7 P) T4 q" e* _0 ^( J, X5 ]3 E+ v                 }
                 
1 Z; v+ D. [& m5 U6 O/ q* X   
+ e' n' S# I8 ^4 J   
   
   
" r# s+ ^) f4 W; {9 [1 b5 W9 l  J   
        gen pminphddy240=`tempminphddy240' in 1
       
       
  F  y. D6 E2 r        *-------------ddy240pminp右侧低峰值 -----------------
- R  Z, F$ z3 q9 \, L& y    preserve

: B, ~- R1 _3 r1 n- N5 M   keep if dx240>`tempminp'
   
; n8 M' a/ S5 e/ }   keep if dx240<`tempmaxp'
1 D: n8 N4 O+ z* m   
   keep if ddy240<0
" W. b6 ]* ?% O  f4 x4 c% E   
   gen temddy240=-ddy240
" ^+ [# O4 r% e   
    count
               
                if r(N)>2 {
3 }4 p" J6 |/ f2 h2 N& S                 
# r) H+ E+ Z# K/ v! X1 H: m                 pkexamine dx240 temddy240
   
, n5 q: N4 g! x% S7 J  local tempminplddy240=r(tomc)
                restore
) K' @6 ~6 k) N5 {0 Q- c1 P! B                 }
                 
- k4 L, u5 E. Y  Y  v        else {
- d0 B2 W. ]& h/ f               
( p% E$ X* M! z( w                 di "cannot find pminplddy240"
                restore
                 
! t! a4 X8 a! L' {3 Z; r                  sum dx240
8 {; Z1 C# s3 i& t                 local tempminplddy240=r(mean)
/ P8 P# b+ t: B% O2 q; \1 F& K                 
( z8 F$ k/ _; X( h2 b) Q. \                 
& t% s7 i" V( w- D                 }
4 s& L8 }2 Z: @# y6 G   
   
# v) k0 H6 Y8 G   
   
   
9 L9 k% [+ _- w4 Q1 v2 R        gen pminplddy240=`tempminplddy240' in 1
       
        *-------------ddy240pmaxp左侧高峰值 -------------------
( Z( f( p" C$ c2 m9 l       
/ S, s; E" o* Q% q* R! r9 }# S) Z         preserve
   
   keep if dx240>`tempminp'
   
! ]/ b3 T" s: k7 N. k! E! W* n   keep if dx240<`tempmaxp'
   
3 z; E( l* B# ?6 i# v9 J   keep if ddy240>0
   count
               
# o- I, {7 f- y* b, ^& L- o                if r(N)>2 {
                 
" Z/ o# }9 O; X6 R% R0 G                 pkexamine dx240 ddy240
, A8 h) O! e$ x1 n/ k6 p0 m   
% |6 o: ^3 F5 |( \   local tempmaxphddy240=r(tomc)
7 O+ D# l9 H- i- A* \5 ~) O! U0 j% v   
! i3 y* q; F7 U   restore
                 }
                 
- v$ _0 E+ y  z! i        else {
' Z( B+ W5 L/ E               
                 di "cannot find pmaxphddy240"
                 restore
               
: A7 E- b7 \1 _, ], C; @' W  s                 sum dx240
                 local tempmaxphddy240=r(mean)
% O- m2 E7 m6 m% t                 
) g5 o9 B9 N6 ~. V! k                 
: R3 c. m: T* t  w( e7 [: X) l6 Z  `( I                 }
. B  a! m9 C+ P3 V6 V) T4 S   
   
   
7 I) i) p0 g; Y5 q0 _! t   
# [  |/ z$ U. v" \   
   
$ _* J3 T, }* q        gen pmaxphddy240=`tempmaxphddy240' in 1
       
       
        *-------------ddy240pmaxp右侧低峰值 -------------------
6 ]/ A4 E$ b/ K       
$ Z, k$ Q  V  b5 t         preserve
  
- C+ N( }2 U0 H9 r( }( H   keep if dx240>`tempmaxp'
7 q  Z( K! F6 L   
1 Q/ a8 f0 [. S3 G% v   gen temddy240=-ddy240
   
+ r4 J  m# V$ j  C   keep if temddy240>0
   
   count
) E( e8 w# L; t               
                if r(N)>2 {
                 
6 l" N, B* x! Q0 F4 ]) ]6 V$ {                pkexamine dx240 temddy240
   
; W* p: |$ l$ O1 C0 s/ d: _, }0 ^   local tempmaxplddy240=r(tomc)
' H8 O" ]# k4 Y2 f  G# M   restore
9 F) s2 E8 \& I, Q                 }
                 
) }8 _$ I% w' p# w        else {
                restore
                 di "cannot find pmaxlhddy240"
                 
1 X8 O( Y' M0 g8 q               
                           sum dx240
  M6 e0 f( d* S1 l8 s                 local tempmaxplddy240=r(mean) in 1
- W! W8 v, ]1 b8 l% D+ P               
                 
3 c' X# g, j1 F% H' b# O                 }
   
   
   
   
   
, i- m( R" A1 V   
        gen pmaxplddy240=`tempmaxphddy240' in 1
5 q1 o; N& m) L( z1 D- {       
7 Y4 Z0 d4 ^% V/ T        *-------------------------------------
9 @  [; ?6 d6 ]6 t2 T5 ~   
   gen price5d=v5 in 1/20
: [% a3 C7 @* s      
   gen price5yy=_n in 1/20
   
   gen price5y=price5yy*(dy240max-dy240min)/20
   
   *price5y &price5d
- P9 o  p) x/ |$ D/ t( F+ n/ r& j   
   gen priceny=price5y in 1
   
   gen pricen=v5 in 1

: V' U+ [/ A( s: w. J1 S     
   *priceny &pricen
$ m/ h' I* x+ v   
   
   *-------------------------graph------------------
  
  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)

* ---------可选项目----------------------
; ~4 p1 b) I- s: O9 t4 J*(spike yxis amean,lc(blue))(scatter yxis2 amean, lc(blue) msymbol(+) mlabel(amean))
8 N- B$ y, T! D2 T# S9 \}
end

maodeng

awer 发表于 2013-2-5 14:40
抱歉，最近太忙了
2 o' J, L7 m9 V9 }5 ?( Z
300236 我用下面的程序试过，没问题

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

喜欢喝冰茶

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

喜欢喝冰茶

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

1 F5 \1 ^7 i, [" {1 J% A这东西是开源的？知道有个统计软件叫stata，不过要收钱。

汉水云潇

准备试一试
先找Stata软件

汉水云潇

1 @" k& B' b7 |" n# I
; H4 D$ L6 L8 Y  ^0 M晕呀
论坛复制程序代码
9 \9 D+ c2 I0 @乱码夹杂在其中
9 n  J1 S) o; N很烦
对了两次
以为没问题了

结果
执行后
$ e5 ]0 U0 W- X, Y; E$ z! d3 q还是出现
8 f) [$ z* g6 o# Jr 199
的错误

汉水云潇

试验成功，下面是三一重工的图
/ G' j2 t# E% K+ Q: ^7 |7 z) \: U

peterbev

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

中秋不圆

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

bknite

这个帖子值得学习。

jeanny

感谢！！！

Hubert2000

明月回春 发表于 2011-8-13 20:17
程序写的不错.
# t" y6 K8 ?& m% t9 c" w不过有一点要注意.通达信导出的数据,只对价格复权,并没有对成交量复权.所以数据需要手工处理 ...

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