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

maodeng

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

maodeng

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

maodeng

. hello 300236
; C8 {) P$ m& }) j0 Ustatement out of context
r(119);
. q+ d$ b& f  a; a. [
& M$ |) ^1 D* {* `7 d大师 就是这个

sergio

这贴实在！

awer

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

抱歉，最近太忙了

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

; O- r1 X& I) Y- C5 M你把 D:\stock\test\source 改成你导出的数据存放的目录，再试一试
: r* l) x0 m5 ^; m
; Q$ b9 e7 X, b% r! t" f5 Vcapture program drop hello

0 ~1 E1 s. a% v) ]9 w. i" Xprogram hello

* f  U" z) F% L* U/ v3 c quietly {

clear

0 e% s* A: E& O: H5 m! u$ |cd D:\stock\test\source
% Z  K2 S$ [% |2 y( n( t! J
5 A, D' U+ D* ^2 o2 ~3 Y# S- C9 F. hinsheet using `1'.txt

   keep if v5>0
; D& [) |9 u# O4 Q1 v# `' A$ m: G   
   gen vt=_n
   
   tsset vt

   gen zf=D1.v5
6 M! ]0 z% Z, r; Q) n9 d/ B
   gsort -v1

& ?2 k6 e8 n4 A1 |5 |   gen avol=v7/v5
+ }7 E" p, W4 {" |& k  q4 ~   
" v  ]0 s4 k, M   keep if avol>0
   
! o2 K. K, ^, j8 w/ y$ o( P   drop vt
+ l  T; b  q7 {4 q7 _) J  {  
: X, j* ~9 ?3 C! h( q' w: v5 ~1 i   gen vt=_n
, [5 m8 M' }6 _3 x   
5 T+ @; W; s# K; I9 ^( a   tsset vt
( b" _' {' Q7 f3 l; [   
    sum v5 [aw =avol]

   gen amean=r(mean) in 1
   
, Z) t1 _0 O* o. @   * 总加权平均值 amean
   
   
   keep in 1/240

   sum v5 [aw =avol]

$ @, O2 U! z0 r% E7 U, L   gen mean240=r(mean) in 1
   
   *240加权平均值 mean240
: s' X5 Y' S! J6 l% ^) K0 K   
) \# ^  ]9 K3 M) x3 h1 K   set obs 500
# n# k, e4 b9 ]: R$ o6 j" @! B   
6 S& w3 l$ n2 Q   kdensity v5 [aw =avol],gaussian n(500) generate(dx240 dy240) nograph
, c* \1 \. R7 f0 ~- P) O: I   
' s/ \! _- j! Q5 b9 b: g# o; T   *240 kdensity(dx240 dy240)
   
4 `4 F1 b& P+ b$ }) ]   sum dy240
   
6 \2 [  }( [& Z5 g   gen yxis=r(max) in 1
0 O9 i* \2 P7 s5 V; H9 D4 c" q   
& H2 i/ K$ f% `/ X) v! x* O  P" w   gen yxis2=yxis/2
$ A+ j3 H4 r" b+ q   
: Y- Z+ I) y( J0 F. z. G   *yxis
   
   gen dy240max=r(max)
$ n: l1 K1 g: x: r# J2 s   
   gen dy240minp=r(min)

    gen vtt=_n

+ v5 |+ s0 q1 W1 T+ b( U. l    tsset vtt

/ {/ p. R. m0 ?$ \& _- I   gen ddy240=D1.dy240/D1.dx240
6 b: y) w1 b8 I9 [  Y   
   sum ddy240
   
5 b$ g9 [) A6 `4 ~' I3 @6 r1 d/ f   gen ddy240max=r(max)
0 I3 f/ u" |4 j: S; M. W- V3 |5 ~9 A   
0 O$ w/ q* s1 V; q   gen ddy240min=r(min)
     
; y% P4 v; R- |# h  V  A0 {/ y   gen ddy240g=ddy240*(dy240max-dy240min)/(ddy240max-ddy240min)
1 ?- V% ^( D9 `; f0 F9 D& q3 ~   
   sum ddy240g
) }$ T  s* j4 _! L* B   
   gen yxis3=r(min) in 1
   
   gen yxis4=0 in 1
% Z3 w( B4 P& H6 y9 v  
: [/ t9 [3 x( q0 w  *ddy240g
. ^; D2 O: W- \( O   
$ a( @% B/ _( `* [   
  preserve
  
  sum v5 [aw =avol]
  
  keep if dx240<r(mean)
  
; x) |% S6 y2 ^: i3 W  sum ddy240
  
" f) C5 f4 Y' x: @/ f) W  if r(min)<0 {
4 p* G4 R& K% ~9 A        
                keep if ddy240<=0
- t  c1 V- U. I  a  g               
                count
               
                if r(N)>2 {
& B; r* A- h' u( n+ J- w: Q                 
                 pkexamine dx240 dy240
% e: E; M5 P8 y, W6 H" H               
                local tempminp=r(tomc)
% k# V$ i2 x( q  v                restore
9 {. y% R7 u3 G) j* A) i3 H0 D                 }
% a% e  f7 g- }: A: f0 s                 
        else {
# j9 \6 E( B" c+ B6 C               
                 di "cannot find minp"
% d* x3 f' `2 P6 b                 
1 S: t% |9 H- @) c                restore
# r& k3 U* }& L: k1 R" F; q2 s                 sum dx240
                 local tempminp=r(mean)
% W% @8 o2 W) l$ }% @+ o                 
                 }
                }
. N+ c5 r) a1 Z1 Z# S   else {
7 f+ `& W3 ]2 b8 I                 
4 s; ?& u2 }  W3 w: }                 keep if ddy240<=r(min)
                 
                 sum dx240
                 
9 S! l% m% m. Z& }* H. U                 local tempminp=r(min)
$ @" z" x" U% n/ T: D/ S                 restore
                 }
5 c. N) y& }. F5 H' P7 u
0 p1 C: i% `) b+ N) P: p& C

   gen minp=`tempminp' in 1
   *240 minp
: C( Z- R4 b" A- e! l1 p
  k+ ?5 i8 j& s% D0 t  x
   
   preserve
   
) k1 e7 f1 s" |2 i; k4 Z0 S- _   sum v5 [aw =avol]
   
   keep if dx240>r(mean)
" i- X  O' ~& t, ?: h* C   
   sum ddy240
* K8 _! c; N- `2 u8 T! K, C   
9 k7 R. a, H$ H1 ], ^5 N   if r(max)>0 {
               
5 R! l# h$ S3 K: w. T                keep if ddy240>=0
               
                count
               
; n* A5 c. ~- q" E1 J                if r(N)>2 {
                 
; u; B' a8 G% h! @                 pkexamine dx240 dy240
$ g4 i0 g% t' o% p0 q               
                local tempmaxp=r(tomc)
                restore
                 }
                 
' M5 ~. F) a! M* j: |  J4 _                 
, `& u, Y: p: }! P' A# q/ w/ x7 P- V        else {
  V8 w6 p; t" f                restore
0 V2 j3 j5 i! U8 L                 di "cannot find maxp"
               
               
                 sum dx240
" i3 Y, m$ g/ F: m( ~) ?- c$ v4 a                 local tempmaxp=r(mean)
                 
& S: @* r+ \+ u: j                 }
                 
                 }
                 
$ Q& W* j8 K+ Z* X% B2 o    else {
1 y+ N; E  }3 a* j8 l       
0 s7 U$ g0 ?! a           keep if ddy240>=r(max)
3 j: ~9 O* N, {- I- m          
           sum dx240
          
* @: _2 ?' f1 Z. z+ y! T: S) `           local tempmaxp=r(max)
9 M& z" D! ]: W; C3 h. d7 c; c% c           restore
           }
0 a# L4 t8 J% B+ y       
6 [; t' @' d( g. n/ P
  

" X% N9 t3 P; K   gen maxp=`tempmaxp' in 1
$ Y' S! u# y0 U2 a2 {   
9 y. x9 ]' s" Z" {/ b   *240maxp
' u3 x5 O/ A& g, {/ d; ?# `   
7 m% L4 r# Y9 U* Y$ V8 @   
: `/ G7 i! Z! x4 S" l  L   *-------------ddy240pminp左侧高峰值-----
   
0 [" e# E2 O$ P; x6 g  d. B7 n+ Q    preserve
   
   keep if dx240<`tempminp'
   
6 k9 [$ P& \! i- a& v! e   keep if ddy240>0
   
   count
" t( h, A9 h  r- g9 U5 A  c               
                if r(N)>2 {
5 J7 f7 P- c/ b+ H: Q' |; ~3 B! B                 
: L0 K$ h! c: F+ I: e; b                pkexamine dx240 ddy240
               
* F; p& V: N+ u+ c! H1 E9 O                local tempminphddy240=r(tomc)
  a- X- v0 l7 v1 P* Z0 u5 c, R                restore
4 P3 E6 O$ D5 C                 }
7 L1 i5 j' q4 V& Z$ O$ N; Z                 
' I* s+ A+ E9 T* K3 ^' h        else {
$ l: T1 u0 d- y( S  Z                restore
                 di "cannot find pminphddy240"
) S# U8 i  H- }: E" a. Z                 
               
2 @# Z% k; B" u5 q* i( n- ~. m                 sum dx240
2 a+ c3 w0 `3 K9 }6 e                 local tempminphddy240=r(mean)
7 O5 T! g( z2 \9 Y               
% l  _: n: ~, a0 ^; c$ Y                 
                 }
                 
- l4 @3 Q4 A% {: A9 W   
   
   
/ j/ [* U7 A4 u4 }7 J' ?   
5 s+ W  V, Z% y   
        gen pminphddy240=`tempminphddy240' in 1
       
/ _/ e' ]- |6 ]2 m9 f+ F( Y       
7 [0 o6 D7 `9 L" C6 M* W' M; u        *-------------ddy240pminp右侧低峰值 -----------------
    preserve

2 V0 ]9 {$ X$ l3 c& y" N. C   keep if dx240>`tempminp'
   
   keep if dx240<`tempmaxp'
   
   keep if ddy240<0
   
   gen temddy240=-ddy240
+ ~5 o& i& w0 F" x+ l4 \$ E   
    count
& {3 [# K- D  M2 Y" O               
                if r(N)>2 {
  \" X4 X* F  P" x$ ~4 `                 
                 pkexamine dx240 temddy240
   
' e7 x1 D$ J7 f7 z5 D  local tempminplddy240=r(tomc)
, E7 X% u* i0 s                restore
* p9 t- @/ \  O                 }
                 
        else {
( J, T9 s) A; e% q3 e               
: V- J" _$ {7 W7 X% o4 `                 di "cannot find pminplddy240"
$ F% Q. g( [+ h, |: P5 X                restore
+ F# T4 ~3 [' h8 R) c( a% {8 z5 U                 
# j) n" i  I$ n, M7 e4 p5 j                  sum dx240
6 q, @3 E, T* j' S; o3 Y/ |                 local tempminplddy240=r(mean)
8 b0 Q$ @  s& @3 `& b# t5 N# Z6 C                 
9 p1 L1 B+ G" O                 
& O0 Z, \8 y/ C1 `                 }
   
   
; _' `7 S9 ?) {   
" ~9 z" I% R5 _   
   
        gen pminplddy240=`tempminplddy240' in 1
* z, X" d6 _" m' U& o       
" M0 d' Y% p5 [& p7 j' }6 p        *-------------ddy240pmaxp左侧高峰值 -------------------
       
         preserve
( g5 V- q0 y# q1 m; r8 d& [$ W   
3 n  h/ G6 K+ A   keep if dx240>`tempminp'
   
   keep if dx240<`tempmaxp'
   
$ w& @8 T' d+ `9 p# i   keep if ddy240>0
   count
               
8 u- H  s% e2 @9 t                if r(N)>2 {
                 
                 pkexamine dx240 ddy240
   
1 V3 U' v+ F- x* u5 M3 k4 C) m   local tempmaxphddy240=r(tomc)
   
7 s5 I  m( ]+ O% P9 l0 ^: b5 b   restore
8 X# J  p7 J0 _* E                 }
8 Z+ V% Z" C( _. P3 v% Q                 
( j4 c, h0 S: f6 e( v% J        else {
" T3 ^4 \  u- ?4 l               
                 di "cannot find pmaxphddy240"
! Q" J4 n8 M7 ~  L6 C1 `1 L9 Y                 restore
0 \( _8 z1 ]% i9 _# E               
5 q6 X- Y4 T8 s2 J" M                 sum dx240
                 local tempmaxphddy240=r(mean)
, e# X- D& \$ y% U6 u                 
                 
                 }
   
   
   
! t, x9 G, {  H. w- ]   
   
   
        gen pmaxphddy240=`tempmaxphddy240' in 1
: b8 k& W, a2 k" m6 ?! K) }       
4 U' K; g$ I* Y) X       
' M: u2 F$ W$ z( V2 V: L        *-------------ddy240pmaxp右侧低峰值 -------------------
4 ]% I7 b* y6 I& i" C       
         preserve
  
% ?  i+ D, O; w# [& n   keep if dx240>`tempmaxp'
6 G% a9 H/ \. S   
6 ^& X; b$ a+ n! U   gen temddy240=-ddy240
" A/ w. d/ y7 w$ x   
2 W( A+ o- |0 S. \   keep if temddy240>0
   
8 H: ?  X% U, ^, E- B   count
: h& r2 L  S& _) u               
0 F" f; z. w- e; C. }. e# |                if r(N)>2 {
                 
                pkexamine dx240 temddy240
$ L2 l9 @) a7 p* _   
   local tempmaxplddy240=r(tomc)
" I9 e) l* C+ W0 @: q: z" Z   restore
7 S# N2 G6 h4 J( |! S  f) X3 D6 o; ~                 }
                 
        else {
& u6 A2 G) B6 h: b( }; r5 X) Y                restore
! J$ ?" _# ^, P" T/ E                 di "cannot find pmaxlhddy240"
                 
+ N% k' Z1 f3 R8 k" S! Z1 S! E  E               
3 N/ {0 M% ?2 K- T- o                           sum dx240
                 local tempmaxplddy240=r(mean) in 1
+ s$ O8 C; G; _               
" g9 F5 Q# ?) I8 l2 X                 
                 }
# c2 O9 z9 C+ g; p4 ~$ \   
   
   
2 O7 h& M2 p7 e% P6 y   
   
   
2 d$ I) R. A3 V1 a        gen pmaxplddy240=`tempmaxphddy240' in 1
       
& ^/ t$ }1 G9 e2 _2 `1 Q        *-------------------------------------
" e- m5 y- X- N$ \. r- F0 B% Q  q   
8 B4 \$ l+ n$ A. q3 t. }0 n5 E   gen price5d=v5 in 1/20
      
   gen price5yy=_n in 1/20
) H0 Y8 V8 F, y- r6 M: P   
   gen price5y=price5yy*(dy240max-dy240min)/20
- |  F, C2 p6 D7 U6 J   
6 e5 o2 Z: r! S) n# I" I( y   *price5y &price5d
   
   gen priceny=price5y in 1
   
) R7 m4 f! A7 n' f- `# q   gen pricen=v5 in 1

7 `4 H2 t  T8 B) p     
   *priceny &pricen
- p' [$ m+ v8 m2 o9 ?   
8 B+ d/ p) j. Y# q0 b8 u8 f2 x   
; |: L) }9 b/ ?) R7 `. Y- ]: }$ B# U) X   *-------------------------graph------------------
  i; B* f) M) d1 n  
7 s8 h& S) l) Q, k1 M5 m: M; k- T  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)

* ---------可选项目----------------------
: _: P; I# @; R, \- }*(spike yxis amean,lc(blue))(scatter yxis2 amean, lc(blue) msymbol(+) mlabel(amean))
: x% @, W9 I5 B7 W5 W1 w! e}
end
* J2 w/ i( v; P5 \3 _

maodeng

awer 发表于 2013-2-5 14:40
抱歉，最近太忙了
/ R8 z; `0 P5 r+ M
300236 我用下面的程序试过，没问题

哇，谢谢，我仔细看看，这些字母数字代码晕了
+ @# y( T! B9 s

喜欢喝冰茶

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

喜欢喝冰茶

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

8 i0 B0 j6 b' Q, t/ ?
这东西是开源的？知道有个统计软件叫stata，不过要收钱。

汉水云潇

准备试一试
先找Stata软件

汉水云潇

晕呀
论坛复制程序代码
乱码夹杂在其中
# D2 z+ q, r/ A) w很烦
对了两次
8 {5 x0 }: B0 |7 P; O6 ~以为没问题了
) Q, M- k) }* ]: v, y, B8 P/ G
7 ]# Z" Z3 M& Y0 e# Q& }  B结果
执行后
还是出现
- b/ n0 M7 S7 d, m5 r4 u0 k* }r 199
的错误

汉水云潇

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

' X0 \! N' f3 R& j0 d! T

peterbev

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

中秋不圆

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

bknite

这个帖子值得学习。

jeanny

感谢！！！

Hubert2000

明月回春 发表于 2011-8-13 20:17
程序写的不错.
* X: t3 D1 F4 X6 q不过有一点要注意.通达信导出的数据,只对价格复权,并没有对成交量复权.所以数据需要手工处理 ...

& a; K0 Q& _9 c( W老哥，也是山寨这篇文章找不到了，不知为何？