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

maodeng

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

maodeng

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

maodeng

. hello 300236
statement out of context
r(119);
0 H5 y, O/ _" E1 u7 m3 L9 k5 b
大师 就是这个

sergio

这贴实在！

awer

maodeng 发表于 2013-1-28 13:13
! [' s) h1 ]( n- A8 ]6 f8 x* C. hello 300236
statement out of context
r(119);

抱歉，最近太忙了

300236 我用下面的程序试过，没问题
' |  j  I9 _0 K, x" c
你把 D:\stock\test\source 改成你导出的数据存放的目录，再试一试
% {7 L$ b0 W7 T% r6 i" \$ {5 V
+ G- }; X% K3 V3 }# [' Fcapture program drop hello
2 t* b: R$ b9 u3 b/ _
8 a8 y3 @2 r5 O% e' E% H) Bprogram hello

quietly {

clear

cd D:\stock\test\source

insheet using `1'.txt
/ w& }8 @( x& T6 w, I+ I
   keep if v5>0
   
   gen vt=_n
   
   tsset vt
; t' }, m5 F3 y
   gen zf=D1.v5
5 y' p/ ~, d* \: b& }* O& V
   gsort -v1

   gen avol=v7/v5
/ X) z. \) G# R7 T   
   keep if avol>0
- r9 \0 q0 ~0 P: {& Q9 ?4 J" |   
   drop vt
  
   gen vt=_n
$ _5 ^' D( z6 h- f# c   
+ s$ C3 m& e; f4 i* b: Y+ q   tsset vt
   
    sum v5 [aw =avol]

: m$ F; u" G- o& q$ F( ]   gen amean=r(mean) in 1
   
% b$ D2 l' V5 u& y) `" {0 E8 o   * 总加权平均值 amean
6 F/ p4 u7 r" |- {, R( V4 {% u   
4 d5 ~$ D$ ]9 h' {2 h   
: M' F) Y  i( B+ q   keep in 1/240
8 Z7 ?+ E3 k: h. e
   sum v5 [aw =avol]
8 r% i) L& L  {* B7 [, V; ]! r8 w, m
   gen mean240=r(mean) in 1
   
7 \7 `$ O$ f% r; P* c; {* L8 q, t' U   *240加权平均值 mean240
   
% W# q- F* z+ C' ], f   set obs 500
3 v7 k8 I8 q9 q' [# x( I' @1 W   
   kdensity v5 [aw =avol],gaussian n(500) generate(dx240 dy240) nograph
5 o& H4 [4 u7 N# w: h   
   *240 kdensity(dx240 dy240)
   
! [- @1 S. O! t: d+ }   sum dy240
   
' q( E5 Z6 n0 q7 |8 c: W- o   gen yxis=r(max) in 1
   
$ C8 b2 Q4 `  G$ E" A' B   gen yxis2=yxis/2
   
   *yxis
$ _( X* \7 I. q/ l: [   
   gen dy240max=r(max)
   
   gen dy240minp=r(min)
# K. F* I9 n4 v. n# K0 ~
/ q- ]. F% \9 Z! M  n/ W    gen vtt=_n

. \9 v8 Z" a. \5 a: F# w7 H+ c1 V; X; j    tsset vtt
' l1 w& u3 }! t
   gen ddy240=D1.dy240/D1.dx240
# l+ B" x, ]- m/ A- {   
' I2 G; c7 @8 {  n: \   sum ddy240
7 s% K1 s' g) d% o  c2 p1 t0 t, j, }   
* f/ {: _4 F) m   gen ddy240max=r(max)
   
   gen ddy240min=r(min)
) E" y# b: J/ S/ X/ ]     
   gen ddy240g=ddy240*(dy240max-dy240min)/(ddy240max-ddy240min)
   
   sum ddy240g
6 a/ n0 J' I# H) r) ^  e# G   
3 r' m+ a/ |" U! P& m5 Z0 V+ K   gen yxis3=r(min) in 1
" F9 r, V# K2 a   
8 R4 o* t2 d8 N" T. s   gen yxis4=0 in 1
  
  g( ]5 y9 f9 q  y$ K3 p  *ddy240g
3 k) t/ n0 W, p3 M; Q9 [   
   
  preserve
  
  sum v5 [aw =avol]
  
4 e6 U0 X2 \0 p. K' S' }+ b3 f  keep if dx240<r(mean)
  
, d9 j; q* v- Q  sum ddy240
  
  if r(min)<0 {
' A7 Z; |3 c& J4 o) _& |: J        
9 t/ J. }( R+ z% ?3 o$ A/ G4 B                keep if ddy240<=0
* V* d* \4 {/ M/ S0 a               
# k8 |* }1 v: w' x/ b" s8 {) w                count
               
$ f0 Q. r0 ~# m( f6 O9 e4 z                if r(N)>2 {
                 
                 pkexamine dx240 dy240
               
2 l: Y. ?( ]9 i1 ?: ~3 A                local tempminp=r(tomc)
                restore
! j( z) M6 U, q! H                 }
) G& D7 z& |1 Z5 U! e* r6 m& E                 
        else {
               
; x! ~1 E& D) D2 `, U6 O/ x4 h% q                 di "cannot find minp"
                 
, A4 z/ P! K4 M8 U% W6 m                restore
- a4 {) `9 t2 O" i/ s0 Q                 sum dx240
* t; v2 t" D4 o" g, z7 P+ e                 local tempminp=r(mean)
                 
/ g/ p4 p" ?. r4 q, k; M                 }
' G  W0 k+ D" a                }
   else {
" G% M( X, v$ D& ^                 
                 keep if ddy240<=r(min)
6 `  {0 h* j& x& @; c- D" S0 ?                 
/ A  _9 J& D: A$ B0 ^* X9 v8 `                 sum dx240
3 d, x3 R6 u6 W* M9 e                 
2 @' |* q; W: b5 C$ _9 j                 local tempminp=r(min)
                 restore
8 B' h9 C; T+ ?) |) E: [: I$ r' y                 }
  R/ E. J, l, J) |/ b+ z$ N9 n& {6 u
" L. T: O! s$ x. f, r) G
4 e- e' @; G* o& ]
   gen minp=`tempminp' in 1
   *240 minp
- J* s! _$ A/ ~* h% i
9 L  A; h! R) m* U, c! }
4 u! e3 t; }+ t   
   preserve
   
   sum v5 [aw =avol]
   
8 ^# O- A; s  _: y6 B   keep if dx240>r(mean)
1 ?7 E6 B: Q! Q4 E   
1 t. g) z4 v0 X   sum ddy240
9 t3 J$ C0 s# y   
% N1 ^& u$ I8 t   if r(max)>0 {
               
                keep if ddy240>=0
8 y! i+ P* e: K4 h, l               
  {2 f4 V' K" y* F. k2 v: i                count
               
) n2 s! U+ F0 S) _                if r(N)>2 {
                 
2 g7 w9 x" ]  P9 E/ U* A# w                 pkexamine dx240 dy240
( D; N$ _' e4 U! z* K2 v6 d: u               
                local tempmaxp=r(tomc)
                restore
                 }
0 S5 D& a' o/ M, v* U& G                 
& L3 f; i/ ]& ^                 
        else {
                restore
                 di "cannot find maxp"
               
               
                 sum dx240
                 local tempmaxp=r(mean)
                 
                 }
                 
                 }
: p. _2 s. k' ^/ C6 v* J                 
" J* q1 T& X  i3 q6 d, ?    else {
       
           keep if ddy240>=r(max)
          
' S, C7 L: v/ d8 ~# w7 o           sum dx240
          
1 G4 j8 _* l; J; c' n( ?: h6 l           local tempmaxp=r(max)
           restore
           }
       
1 W  C, w! o/ n2 G
5 A0 B5 V- @, p7 r7 O  
" _: W& p: j+ U+ _1 w% o; M
5 i( z! m0 t. x5 K' Q6 \1 `   gen maxp=`tempmaxp' in 1
   
& w( D& m) p7 S5 ^8 Z( q" t   *240maxp
) ]+ H% N7 ~# H) p0 d   
   
   *-------------ddy240pminp左侧高峰值-----
. A* n$ n$ E& d6 D0 C! t   
: J8 t# k! M5 U% w3 k    preserve
6 V0 ]2 ]: ^9 u- X) p5 {. o! ?   
   keep if dx240<`tempminp'
; U" n# w; F! u' r$ q3 \' g   
   keep if ddy240>0
' Z& o2 G0 U8 m8 V8 S   
) u8 H( a# D- J  {* ?   count
               
% S( w+ X- _; b* j                if r(N)>2 {
( E: m. C) g, ^2 ~5 L                 
( c/ o3 a9 G! X                pkexamine dx240 ddy240
0 o' w4 x2 z# e               
  J$ L6 O  ]( S                local tempminphddy240=r(tomc)
0 z- V( \  ^- p: O% Q- W                restore
                 }
& ?8 g$ M, G4 c+ k9 p( t                 
        else {
+ s5 y, L% L, o                restore
                 di "cannot find pminphddy240"
                 
               
                 sum dx240
8 y! n8 s0 v6 e" X8 q& o                 local tempminphddy240=r(mean)
               
                 
                 }
                 
   
: |9 C# g& ?5 [, p0 c   
' n, I% [. t, Z0 B4 Q* a/ ?   
4 C5 B) b3 D0 i   
. Q+ w! _7 a5 }4 l   
- G  B& w. E6 M        gen pminphddy240=`tempminphddy240' in 1
       
: K) I4 ^7 W' c4 c0 J       
- O7 G2 z8 y2 U1 a2 X2 T2 `+ m        *-------------ddy240pminp右侧低峰值 -----------------
    preserve
  L6 p6 E1 e! ^" u+ u
   keep if dx240>`tempminp'
! W; i# d) `8 r5 S) S! W; {7 W   
) W# r. S  {- D5 Q# B   keep if dx240<`tempmaxp'
   
& s! m- [: O# v0 p2 h+ S* l   keep if ddy240<0
   
   gen temddy240=-ddy240
   
    count
               
! N; g- N; A$ K3 l7 {8 V; V                if r(N)>2 {
                 
' Q1 ]. a0 G" Q: ?% S/ [                 pkexamine dx240 temddy240
" t, D% L& A4 U0 U   
  local tempminplddy240=r(tomc)
                restore
8 R0 O1 M2 `- r6 z) E, V2 B* R                 }
                 
        else {
               
. }& V7 p8 W9 T                 di "cannot find pminplddy240"
' H  U; J! ]# _' L( {                restore
& h* t, M7 O) o7 B6 w1 x0 H. l& |                 
                  sum dx240
3 x" w) j7 @) |9 [7 x1 ]2 L1 K                 local tempminplddy240=r(mean)
- z) @% h' Z1 J3 E% }* E# ^; A, Y                 
; W* P( n$ }/ {4 w' y                 
2 ^+ w, {% }! [- b                 }
% v$ ^' a) x7 u& {1 J7 ^   
  D6 c* P1 I$ v& E6 r   
0 C. T8 V5 w+ r   
$ q! T1 n* B& o$ {   
2 v' z$ d+ f  {3 U0 J2 U   
7 \5 h: b; r+ w3 D+ l        gen pminplddy240=`tempminplddy240' in 1
+ D" u+ h8 e4 ~; s+ X# O       
        *-------------ddy240pmaxp左侧高峰值 -------------------
       
         preserve
   
$ C% r/ ^0 _9 f/ M8 p$ k   keep if dx240>`tempminp'
% k0 n, u: {- b2 U   
/ y+ S. b( p# }7 l2 N8 R1 Q   keep if dx240<`tempmaxp'
   
   keep if ddy240>0
$ V5 \- Z( @; k0 ]9 w$ D/ @8 D   count
               
                if r(N)>2 {
                 
5 I7 D3 z1 x" x: @* X) J$ P$ H                 pkexamine dx240 ddy240
   
7 ]" c9 m* b  k. ~2 f7 k' J" ^: q   local tempmaxphddy240=r(tomc)
! _/ \. {7 T' L" s$ y   
- G/ _( S# N2 W: r; ?( y   restore
5 ]5 c" j) F% ?$ ^* Y% Q0 A* e                 }
                 
+ R" ]' Q7 M3 r8 u" D        else {
               
                 di "cannot find pmaxphddy240"
                 restore
               
                 sum dx240
, o: P/ _, @& J. k/ N                 local tempmaxphddy240=r(mean)
6 |9 q# {* @4 x4 A+ R4 w                 
                 
9 Q, ?: o1 B& x0 v                 }
   
   
/ D/ I3 I4 p- u3 a# C; t" v   
! E2 R& g5 w  Z   
   
% P0 X6 A: ~6 s7 B8 U7 |" `, K   
        gen pmaxphddy240=`tempmaxphddy240' in 1
) s) X( N# T+ @% L' g7 ]" ?: l       
  A, m; S6 s3 o& t6 Q1 Z/ B8 B: T       
1 G7 H$ N1 a( b        *-------------ddy240pmaxp右侧低峰值 -------------------
& C8 ?- t/ K- T& \$ W       
         preserve
, r* p- n1 U; s, l! u& n  
   keep if dx240>`tempmaxp'
8 L3 z+ q2 y% c! b: b   
   gen temddy240=-ddy240
   
. Y, G0 Z7 c! `, @9 P  P* d, y   keep if temddy240>0
4 ~& L4 N6 m7 f) R% n! b   
   count
% o: V5 j" x$ S. W               
& o  j: h3 q  g* F* n; \% o                if r(N)>2 {
                 
                pkexamine dx240 temddy240
   
5 t1 a* W; Q' X" G4 D   local tempmaxplddy240=r(tomc)
7 U- W% @  V0 t+ N8 u   restore
6 }3 s' C: C. ^7 {3 M4 m- {                 }
/ h; S) [. H: U! u                 
+ V  \' W( s* q' }        else {
                restore
                 di "cannot find pmaxlhddy240"
                 
               
                           sum dx240
5 N+ r1 p4 A6 s4 w0 r( P! T' i4 B                 local tempmaxplddy240=r(mean) in 1
               
                 
- x- p# }  E& \: X* j                 }
   
* T4 l# }  q* Z0 X0 |9 s   
3 a' s" z( f3 o; ?' L/ S  o& G   
   
   
   
        gen pmaxplddy240=`tempmaxphddy240' in 1
       
5 v, ^  W# [% x6 R9 {. x5 \% h8 {7 W        *-------------------------------------
   
   gen price5d=v5 in 1/20
: f: U2 ?8 I: B$ v  `: a      
   gen price5yy=_n in 1/20
   
   gen price5y=price5yy*(dy240max-dy240min)/20
; j% L1 J5 i5 A0 Q3 F0 ?   
   *price5y &price5d
, K8 y/ q0 G3 o! r  p  B) Z* c   
7 |: ~' n0 S* C, q   gen priceny=price5y in 1
! g, c, g+ b+ o0 S! q  ]   
   gen pricen=v5 in 1

; }0 b& f& c' M     
   *priceny &pricen
   
   
   *-------------------------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)
: e# p0 q7 w  v0 v+ K. P. Z
* ---------可选项目----------------------
*(spike yxis amean,lc(blue))(scatter yxis2 amean, lc(blue) msymbol(+) mlabel(amean))
5 ?: Q: Q" @7 T+ r$ n' p}
# Z0 a9 g# Y. m9 N) K3 o! Aend
* Y0 A" E2 U" v7 ?3 M) Y

maodeng

awer 发表于 2013-2-5 14:40
抱歉，最近太忙了

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

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

喜欢喝冰茶

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

喜欢喝冰茶

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

+ B0 u( n' S" `这东西是开源的？知道有个统计软件叫stata，不过要收钱。

汉水云潇

准备试一试
先找Stata软件

汉水云潇

晕呀
  A6 c# a8 w* c$ a论坛复制程序代码
乱码夹杂在其中
1 f! }: c2 P$ L很烦
! s& F' {; b6 Q% J! m对了两次
. R: v7 {( S, M% F4 S7 `4 Q以为没问题了
0 H0 z1 T# S# [% i
结果
8 ^* Z  I% q9 Q1 {执行后
还是出现
r 199
的错误

汉水云潇

试验成功，下面是三一重工的图
# H% l6 L9 h. `8 J  l) c: R

peterbev

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

中秋不圆

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

bknite

这个帖子值得学习。

jeanny

感谢！！！

Hubert2000

明月回春 发表于 2011-8-13 20:17
0 y: [9 Q" z  H8 q程序写的不错.
不过有一点要注意.通达信导出的数据,只对价格复权,并没有对成交量复权.所以数据需要手工处理 ...

6 ~9 n% n* J' P  u; v* ]: W! h  }4 B老哥，也是山寨这篇文章找不到了，不知为何？