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

maodeng

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

maodeng

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

maodeng

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

大师 就是这个

sergio

这贴实在！

awer

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

抱歉，最近太忙了
" X8 A3 V+ W& t* i4 ~+ |
300236 我用下面的程序试过，没问题

你把 D:\stock\test\source 改成你导出的数据存放的目录，再试一试

* U, T- z0 R2 ]8 U! ?( Wcapture program drop hello

: N3 ?0 V9 f! eprogram hello
: D6 N7 c4 n# Z# j
quietly {

clear
, n' k3 D/ i4 |/ ^
cd D:\stock\test\source
. o0 l& X: ]+ h7 q2 R4 _
" e; K9 h/ [8 jinsheet using `1'.txt

- f) m- K1 c: D   keep if v5>0
   
' m; i! G& R0 H/ c   gen vt=_n
' |  @. Y' r9 A8 e' k   
   tsset vt
  n5 s3 {3 g, L$ y: `0 g
   gen zf=D1.v5
  {, \1 e" Y7 R% V' i# s
   gsort -v1
3 u1 u! R* \; w6 ?
/ N/ ]; X; H$ `7 }   gen avol=v7/v5
   
   keep if avol>0
0 v( L; i0 |, n' a% T   
   drop vt
* u. r+ [' }+ n- H; |: [: x  
, _. x: X! r& I' J/ Z, r% n, I   gen vt=_n
   
7 c2 R, c/ l% C   tsset vt
7 ?& `5 A) Z9 i$ Y/ \2 v   
    sum v5 [aw =avol]

   gen amean=r(mean) in 1
  `* Z- E' m6 s# n$ M5 p# Z   
   * 总加权平均值 amean
   
/ {- M" k* M3 X   
! G; q3 O3 l; u% Q: y/ w   keep in 1/240

5 v! V$ B0 i" i7 \   sum v5 [aw =avol]
' W" V4 C- z( A  \" p
   gen mean240=r(mean) in 1
0 B5 C% m4 U0 V* G: ^   
. d1 M% k3 B" J: ~   *240加权平均值 mean240
   
/ @3 }4 L( T* d% l; Y6 Q: b   set obs 500
) s+ z  H' g1 d' h& {   
) `" F" N% r( J+ k4 H   kdensity v5 [aw =avol],gaussian n(500) generate(dx240 dy240) nograph
   
   *240 kdensity(dx240 dy240)
   
   sum dy240
   
% r- U6 i) ?) N. f   gen yxis=r(max) in 1
* A5 |, t+ Y/ ?, I+ N   
  V" e) Q' a1 o( H: L8 x   gen yxis2=yxis/2
& i6 w( |3 U* h# S' h   
' C, z" H0 k/ }6 f   *yxis
% m7 e/ q0 y6 Q9 R* w   
   gen dy240max=r(max)
   
   gen dy240minp=r(min)

& i* O% F4 B* ^/ x0 ]& [+ J    gen vtt=_n

2 p4 N) F3 Q9 |9 G' V6 ~$ U3 x    tsset vtt

/ _- `" y% x3 T; Q   gen ddy240=D1.dy240/D1.dx240
7 B/ X# y4 P9 N/ M   
3 R3 `; s; N/ j& A( h   sum ddy240
9 @5 F: y, H0 j' P% k   
   gen ddy240max=r(max)
/ R" P# Q" o% n6 k* i5 X7 a   
   gen ddy240min=r(min)
     
   gen ddy240g=ddy240*(dy240max-dy240min)/(ddy240max-ddy240min)
) M! n* m; L( l0 H1 T6 S4 l   
' I3 s9 V3 M# x; M   sum ddy240g
   
   gen yxis3=r(min) in 1
7 M. @( ?2 [2 y" P   
   gen yxis4=0 in 1
  
( e6 H: l3 d6 M  *ddy240g
   
   
1 P! Y& P0 D: V* |3 s  preserve
  
3 u! `; ?& H% U2 R. J  sum v5 [aw =avol]
  
  keep if dx240<r(mean)
3 r. y% ]! P) ], `+ @% {0 H# t  
  sum ddy240
* z# D8 H1 h) a" Z  
. y" Z- w& K! |" v  if r(min)<0 {
; v+ Y* `" ]+ G( u+ V        
1 `% j) K; ^1 i                keep if ddy240<=0
               
                count
* W- h* }2 j9 _' D4 y3 ?) @               
, J4 L* P, c& V                if r(N)>2 {
. q$ O/ \: P6 O, F                 
  [) n% Z- [+ p5 k: v" L  `                 pkexamine dx240 dy240
% X  u, C4 @& C- O5 D0 t               
                local tempminp=r(tomc)
4 S( |& N+ D0 G: M1 ]  ^9 k                restore
# d+ Z8 N: X' H0 |                 }
' x& Y" ^( _/ d; K6 W                 
  A" L8 R2 [2 n. d        else {
  ~, m$ q6 O  O               
                 di "cannot find minp"
- h% V! r. a/ x/ j                 
' I" e8 J+ m& p& M$ m, O$ K                restore
                 sum dx240
                 local tempminp=r(mean)
                 
0 G, g/ k! U! h% a+ ?$ I+ {. y                 }
                }
6 n/ l8 R+ N3 b3 M  W- i# O+ C   else {
                 
, j% D3 r1 A! j                 keep if ddy240<=r(min)
                 
                 sum dx240
$ A1 F( S8 Z$ G7 k: e1 }                 
                 local tempminp=r(min)
/ i  Y" h  a" _7 S6 S; p                 restore
                 }
& h, p! |8 {6 |. G, d. r


   gen minp=`tempminp' in 1
. v$ M+ X  C" ], `1 H   *240 minp

- G" O1 F/ \; p/ R9 \
+ A2 A. ^5 Z& l7 b  M; N4 i0 ]   
( A. _  E% h) C8 A4 N4 Q6 _   preserve
1 w7 P- h9 d2 |7 a   
' T: N, L# `8 `( E* q, Q7 D   sum v5 [aw =avol]
   
   keep if dx240>r(mean)
4 g8 V0 v7 M9 z. ^   
   sum ddy240
" N3 N5 C1 E. @0 Y; y" \1 o6 |   
8 f, i7 J% h5 o" P. [3 {   if r(max)>0 {
8 t) T4 ?* y) W, J               
$ z- a  o8 {/ ]9 r# B                keep if ddy240>=0
' f2 \& V( K2 A: H8 d5 h               
7 {( i! K0 P% \  Q                count
' @$ `3 n5 N& E7 F7 o. R. s2 ~               
1 U, [- S: }% t! A7 X3 t+ M  Y& F                if r(N)>2 {
) ?9 r0 T$ d$ ]                 
2 [; R1 y7 f9 i                 pkexamine dx240 dy240
               
7 b! ]) I4 U0 L4 k+ D* E6 x                local tempmaxp=r(tomc)
) G6 N+ v# S" h                restore
                 }
                 
2 i% H9 _( Y9 P' [2 j                 
( P- Q) ?1 O, O: t, U* g        else {
: ]0 h# t, F& D3 M& ~, `7 Q7 D% ?                restore
  D- E0 Q  y$ X% Q# I2 b                 di "cannot find maxp"
* k; r. O$ P2 R+ p6 d- l               
               
8 E4 j4 Q- Z9 f$ A                 sum dx240
0 M( P- r- Y* k0 n                 local tempmaxp=r(mean)
                 
3 S9 y1 o, m9 v6 g8 w5 ~# @) G                 }
3 m% [: h1 i" q% ]                 
5 S4 Q8 Y  G9 h                 }
3 z; D1 W8 m( `: O$ ~( V) \                 
' W7 I+ r0 s+ ^: O- R    else {
       
6 q. O1 h) C; L" [: ]           keep if ddy240>=r(max)
+ |, T) y+ m6 d0 q6 ^, X2 x9 T+ w7 v. M          
           sum dx240
6 A0 D/ Z/ V) k7 h9 x5 Z. y          
( _( f9 b; c$ d4 t6 B           local tempmaxp=r(max)
: e- r+ Q6 |: K2 {4 v# ?           restore
1 W$ u: K3 T3 H: W# }           }
8 F1 M# b0 O) J+ b* ^1 @9 {, t       
* h7 c. i7 Z* Y. {0 r
5 B# n( v# T) Z: w/ a0 C  
; [0 B/ z9 W3 }: |& ~- U0 n, R
% P+ B7 E! Y% s/ h4 r   gen maxp=`tempmaxp' in 1
5 }4 S1 O* k  V2 B' ]; m* \) l   
  T# O+ y; d! l9 }; ^8 f6 s) a/ Y7 _   *240maxp
. N4 T. f3 i; F! W   
   
   *-------------ddy240pminp左侧高峰值-----
# }# A( b# M; f+ Y, J5 W   
" f# Z$ W; ^9 p% P* B6 A    preserve
   
   keep if dx240<`tempminp'
+ r0 o2 O# W* {; V3 p   
1 m" B* c# b* q; s0 t   keep if ddy240>0
   
( k: s) P1 ?1 t) m+ s! h5 i   count
               
                if r(N)>2 {
; Z& j4 R4 ?( c) N                 
2 K5 J+ Q$ O& I+ I0 ]3 }8 W                pkexamine dx240 ddy240
, Y1 d$ A! H: V: ], n2 O               
7 J9 `" J. y* K# b. i+ B                local tempminphddy240=r(tomc)
4 j( q" y& @' C+ b2 e- a8 ?                restore
( |- D8 ^6 J  n: s" S+ t" Z$ E                 }
$ Y/ j$ a$ O9 J  t2 t$ ?  d                 
: h2 \# v( _, U1 J        else {
                restore
6 u. y! r' ?# H! k+ Y                 di "cannot find pminphddy240"
+ C6 q4 g" L; J# R                 
               
& M" h  {. x0 y9 A) X  ~. t" E1 ^                 sum dx240
2 l$ }3 r. n: Y+ d                 local tempminphddy240=r(mean)
               
                 
                 }
                 
   
, ]) ?3 O6 n8 l" G# }   
   
   
   
        gen pminphddy240=`tempminphddy240' in 1
+ \* U/ F6 H& a! ?% v2 p       
       
+ M1 j; V, h) A; k" w7 q3 `+ R0 c        *-------------ddy240pminp右侧低峰值 -----------------
( v) K8 M4 I+ |8 X6 g8 N    preserve
# j( D- w* `* U* n2 ]1 E
   keep if dx240>`tempminp'
% }, ?" v9 M% X5 U   
/ z0 }% n# N1 i/ b% |   keep if dx240<`tempmaxp'
   
1 U! V7 v- z; Z2 u- _4 e" x   keep if ddy240<0
9 U4 g! o9 }) r7 a0 S. g   
   gen temddy240=-ddy240
   
; O/ P( p* j# N, O& `7 ]" e5 @    count
- r, E/ `2 O$ T2 B) E& q5 b               
                if r(N)>2 {
: a8 i. Q7 c  o' ^1 h8 E6 T3 V9 f                 
6 `+ W5 v0 o9 }/ J                 pkexamine dx240 temddy240
   
. b5 Y8 z7 w. V# m( H. Y  K* D( l  local tempminplddy240=r(tomc)
                restore
                 }
                 
3 p$ h$ N" Y* c7 c        else {
- n5 r) ]3 \% }% `" z" ~8 r+ m               
                 di "cannot find pminplddy240"
                restore
                 
                  sum dx240
3 I. z8 z! N9 B$ v$ M- F5 n% j. p0 ~! m                 local tempminplddy240=r(mean)
                 
# L% E3 P. ]& t& X9 f2 f                 
                 }
   
   
! c% ~, X/ D/ O4 H7 z" [# w   
+ o, M! P* N0 P+ J  n   
   
5 n$ _/ J- i( i/ W        gen pminplddy240=`tempminplddy240' in 1
$ e7 `) z0 \6 U. e) j* M, M       
        *-------------ddy240pmaxp左侧高峰值 -------------------
( B  D1 q. I- @8 m8 ?% ~. e9 M' P       
         preserve
6 f2 y# E, {; K: |1 b   
   keep if dx240>`tempminp'
* \/ o) @/ z; e/ Z, z   
   keep if dx240<`tempmaxp'
+ J) z7 j- Z2 s3 [4 O5 k   
   keep if ddy240>0
   count
6 B; g& F/ t. C% O' }               
- u- k9 m0 s+ |1 ^                if r(N)>2 {
                 
/ `3 y: f* {; d                 pkexamine dx240 ddy240
   
: R& B8 l" ~; C7 V8 t   local tempmaxphddy240=r(tomc)
   
' S  Y: Q. _! w" {. A) B   restore
                 }
                 
1 J! o9 u( Z" B" V5 L' d5 U$ J" }        else {
               
                 di "cannot find pmaxphddy240"
$ [6 @% Z  G' M$ m4 m/ b7 y                 restore
" U& \0 o! u* }  j* h0 m               
$ P/ ?; P7 l8 I  q( B* N% ]9 \! o                 sum dx240
                 local tempmaxphddy240=r(mean)
                 
: ]: q% r- k# d) ?9 m                 
; T) ]) ?  M. o  M1 q! \                 }
   
   
   
# r4 R) `/ ~8 l6 f; |   
1 x0 F6 d& x$ v2 a2 N, \! D   
   
# s  f* w" `6 @2 l3 q        gen pmaxphddy240=`tempmaxphddy240' in 1
' @- _8 U2 i7 a; n       
       
        *-------------ddy240pmaxp右侧低峰值 -------------------
       
$ }  y- T0 b- m. S$ D* x: B  H         preserve
  
   keep if dx240>`tempmaxp'
   
" o2 x0 ]8 p; V" s2 n   gen temddy240=-ddy240
   
0 P" B9 I( k  G6 n9 Q* m   keep if temddy240>0
/ ^% K3 w& `, X   
3 h3 K% [5 ^( a0 z   count
) S$ G0 ]7 O# \9 d( \               
9 a3 l/ b  R  S: \+ M                if r(N)>2 {
1 ?: }, r9 z8 ]                 
. r! X$ K3 |' @" I+ X% e                pkexamine dx240 temddy240
   
   local tempmaxplddy240=r(tomc)
) d: \0 V3 x9 u: z   restore
9 B* `0 S4 E; e                 }
                 
        else {
- O, {8 s, K# n/ Y! ?                restore
                 di "cannot find pmaxlhddy240"
; ^; @5 M+ p  Q% ?0 D; j                 
5 ~1 K6 K7 Z- R9 J               
3 c/ \  i% D7 W; Z% R- D. Z                           sum dx240
3 [3 d: U) ?2 b1 K; ]. i1 P                 local tempmaxplddy240=r(mean) in 1
" b7 J( s& x" H3 v3 X5 o) b8 n/ t4 m% l               
                 
                 }
) m+ ?/ l% V2 M7 C. q: d0 \+ [   
   
8 H' r: x: i1 _1 d" U   
   
/ y& M- M+ a  e   
/ I" ^) z1 }' l! u1 F" Y3 |   
        gen pmaxplddy240=`tempmaxphddy240' in 1
, a# ?9 z) ?: o& p2 I       
        *-------------------------------------
   
   gen price5d=v5 in 1/20
      
   gen price5yy=_n in 1/20
   
   gen price5y=price5yy*(dy240max-dy240min)/20
   
   *price5y &price5d
   
   gen priceny=price5y in 1
7 @- e, I. w7 q- U/ m   
   gen pricen=v5 in 1

6 r% W8 m8 Y% m+ L7 [/ p  x- J6 b     
0 f) e! W) u2 ]1 j" c. R& Q1 V. f) K   *priceny &pricen
   
   
* N$ A7 \' _7 u# f5 c. l   *-------------------------graph------------------
. L8 T1 v: v  T: H  f' B2 y  
( X- m# x; W' v  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)
5 W+ O' `5 j5 e; k
* ---------可选项目----------------------
5 S4 V) }0 R9 b6 y6 @*(spike yxis amean,lc(blue))(scatter yxis2 amean, lc(blue) msymbol(+) mlabel(amean))
}
. C9 B5 e9 @5 ^% P$ Wend

maodeng

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

+ L1 g' r& o8 d3 q2 c/ o0 G300236 我用下面的程序试过，没问题

+ X! A4 e1 w7 o4 M- a0 t: H, W3 _! {哇，谢谢，我仔细看看，这些字母数字代码晕了

喜欢喝冰茶

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

喜欢喝冰茶

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

这东西是开源的？知道有个统计软件叫stata，不过要收钱。

汉水云潇

准备试一试
先找Stata软件

汉水云潇

( {" v! ?4 C$ f, E
5 K" ?, o  j) Z$ Q( d; I# @* v* k晕呀
/ h& R, w0 t  J% o  V3 Y论坛复制程序代码
乱码夹杂在其中
很烦
7 d: J2 K8 u2 X: b/ x' V对了两次
以为没问题了
! H. e* g0 P, p
结果
执行后
还是出现
r 199
7 I. V; `" P6 E$ E+ u) f6 K- M& X7 D的错误
6 B' e  s+ c" c* i7 h; y5 l& q

汉水云潇

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

peterbev

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

中秋不圆

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

bknite

这个帖子值得学习。

jeanny

感谢！！！

Hubert2000

明月回春 发表于 2011-8-13 20:17
6 B  U5 Q* f8 k, g# h) f5 c( q  ~% z程序写的不错.
不过有一点要注意.通达信导出的数据,只对价格复权,并没有对成交量复权.所以数据需要手工处理 ...

% F8 p% P% h) B" @( t老哥，也是山寨这篇文章找不到了，不知为何？