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

maodeng

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

maodeng

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

maodeng

. hello 300236
* V8 d" J$ h6 Qstatement out of context
. L+ F, N3 j0 e# _r(119);
6 n+ r! R0 J- G9 f: h# X. S
# Q2 G9 w  V3 G- T  j大师 就是这个

sergio

这贴实在！

awer

maodeng 发表于 2013-1-28 13:13
. hello 300236
8 o8 Z7 @6 `2 q+ G7 _; |' Rstatement out of context
9 e" X2 B. Q& n! ^% l, Jr(119);

抱歉，最近太忙了
$ k6 M# r1 V6 P( {' s/ a1 M
' y2 V6 k* m# v# Q300236 我用下面的程序试过，没问题

你把 D:\stock\test\source 改成你导出的数据存放的目录，再试一试
/ |+ k& E4 H0 a! a8 A
capture program drop hello

program hello

quietly {

clear
2 a: s' ]5 N( ~3 |1 P
cd D:\stock\test\source

insheet using `1'.txt
, D& p6 b/ R" p3 ?' Y9 i
6 ~5 H  y. a) l7 j& W) S& Q   keep if v5>0
/ N* n/ T# T8 {' `& |# n   
- Y& n, ]/ ~- e   gen vt=_n
   
+ E* T) t% F4 {/ C( n  f   tsset vt
% N1 U! O& M: T6 D7 d: G! M
   gen zf=D1.v5
5 y+ Z: q$ H: C0 S$ N/ ^. M  W
: y( H& F1 m" u0 w   gsort -v1

   gen avol=v7/v5
   
' e$ l) H% c  Y   keep if avol>0
   
   drop vt
  
2 [, }! Q3 O7 i5 q/ s5 @   gen vt=_n
   
* q! b) v" |+ n2 ?5 Q   tsset vt
   
    sum v5 [aw =avol]
, b& G) W, n" w2 h" |( s1 S  H( w
/ Y1 J* D- V. {! @6 |$ }7 M' ?   gen amean=r(mean) in 1
   
- u% |7 V8 F1 V& ?4 r   * 总加权平均值 amean
* G7 C( [, [+ Q   
   
   keep in 1/240
' f/ [7 P+ N7 k
9 Z2 \7 f; J  i   sum v5 [aw =avol]

& G2 e; ~. V4 n# T+ \) z   gen mean240=r(mean) in 1
   
   *240加权平均值 mean240
   
& y' z' Q$ d0 W* h   set obs 500
   
! C3 R% B( H2 ]/ F" }( d& V' V/ G   kdensity v5 [aw =avol],gaussian n(500) generate(dx240 dy240) nograph
1 ^# ^& j/ ?) d3 [6 u. {   
   *240 kdensity(dx240 dy240)
! B# K0 g0 C, g+ J6 t6 B   
   sum dy240
   
* a2 J# U' M2 [: d: Z   gen yxis=r(max) in 1
   
   gen yxis2=yxis/2
1 _% Z1 K# n" _4 s: h* H! n   
   *yxis
. F  z7 V/ K" B' H- `0 X' D( ]   
   gen dy240max=r(max)
   
( ?. ?5 A  S  u7 V, i& S   gen dy240minp=r(min)

    gen vtt=_n

    tsset vtt
3 I: I* x. y+ t. B0 E
5 r0 V1 P0 `. y   gen ddy240=D1.dy240/D1.dx240
( Z' M' k$ d6 K; K' p   
   sum ddy240
   
, ?# I, ^& M( h4 b+ e: Q   gen ddy240max=r(max)
   
   gen ddy240min=r(min)
     
   gen ddy240g=ddy240*(dy240max-dy240min)/(ddy240max-ddy240min)
( w8 F* y; X. n- q4 \$ s: ?' ?7 \( t7 ^   
: Q9 d3 H% x$ @2 O" A   sum ddy240g
   
   gen yxis3=r(min) in 1
   
+ d/ I9 j" c, X   gen yxis4=0 in 1
1 Z- O/ B8 S7 n5 I9 m  
5 e% {& ?8 u" A; q2 W4 y2 Q5 }  *ddy240g
   
   
  preserve
. V" z5 H. S( n9 [( Z  
7 B0 V9 U7 a2 @4 t) G6 A) q  sum v5 [aw =avol]
6 _9 j! |* A4 C6 q( A, Y9 {" d' A8 [9 G  
  keep if dx240<r(mean)
7 M6 v# n1 a: Q  ^& x  
  sum ddy240
  
) S+ b' n4 |& ~* V  if r(min)<0 {
9 F- Q" J4 u/ {6 ?* H7 X        
                keep if ddy240<=0
               
                count
& [2 [- Z- l! c0 v. P+ v/ v               
( G) c, Q8 g- I& O: ?1 |: S; |                if r(N)>2 {
+ C4 f5 a" `+ M4 y. w2 a2 U                 
. f( G) i7 J0 b1 |+ l0 r6 l                 pkexamine dx240 dy240
+ M# X+ k' `6 H% e( w               
. h0 S7 }& g/ C% D3 d7 ^5 ^                local tempminp=r(tomc)
                restore
                 }
* ^2 _1 g8 Q4 N4 l                 
8 A# r7 m9 A7 f  }4 Y  V& E% ~        else {
               
3 R7 o# s1 N8 C! Y; i0 l) t: c. m                 di "cannot find minp"
' }; J2 P, R+ F1 T- N# p- t" ^                 
- k- P0 N) M) o4 Z3 s2 D  s9 e' }9 x                restore
                 sum dx240
                 local tempminp=r(mean)
                 
+ s0 k7 E: M: t( U& p                 }
; M0 Y; F! M4 y4 t! [! ]                }
   else {
; J' S) o1 L3 R  C. j( j6 z/ {+ g                 
$ `. N' D& }$ A* j+ X* G! X# @9 _                 keep if ddy240<=r(min)
  B: K0 J9 u7 R" @                 
                 sum dx240
                 
                 local tempminp=r(min)
                 restore
  y3 r4 A, O1 Y                 }
5 b+ l' r( m- {


3 F) B/ l0 c, e$ P3 ?4 _   gen minp=`tempminp' in 1
   *240 minp


1 J6 e- \5 E/ x7 d7 F, y/ x   
3 f3 g2 n" D  u, T3 O  m   preserve
6 u' S2 J' O( y& A   
   sum v5 [aw =avol]
3 L2 _& v% S# s+ J' L9 E   
) B# S' O! ~) Z' m3 }   keep if dx240>r(mean)
; O" U/ G3 ~. m# F  v% N   
   sum ddy240
   
   if r(max)>0 {
. w" H  c, W' T% |               
2 }) U5 q+ F% T; N  U' f- ^$ Q+ y  R5 z                keep if ddy240>=0
               
                count
               
/ ?6 T" T3 t/ ~                if r(N)>2 {
& J  f' S9 j- L( ^                 
                 pkexamine dx240 dy240
               
                local tempmaxp=r(tomc)
                restore
; u$ F: z5 D; k, G                 }
                 
                 
  r; B* `, t' A0 T        else {
; o: c9 r$ `$ k) [! p1 I                restore
: x- o% v% R9 V. J- h3 W; d, t1 M& \                 di "cannot find maxp"
9 h( i6 q! P! b8 c9 J4 f5 D               
9 h) `! P& n' u: @% j. Y# k/ r               
) R$ B# b# B0 b- R/ C                 sum dx240
                 local tempmaxp=r(mean)
. H* G2 r. y3 i                 
                 }
  e! `% u+ c2 ]6 u* Z7 C                 
8 ~: R5 X  e1 D! S- b4 I                 }
/ W/ d8 `. i8 Z: P$ i6 D  \$ H5 Y                 
! i+ h3 u9 h- B% _& I    else {
       
$ H1 I$ {3 Y( u" d, N2 a- {* H           keep if ddy240>=r(max)
          
           sum dx240
          
# u  g4 r: q+ l' i+ C' x* l           local tempmaxp=r(max)
& \- e5 {; e" u/ [: R7 `1 E           restore
           }
: d8 j# F$ R  V0 f" G       
( f, x; \! z0 H. K2 w) V2 I% u
  

   gen maxp=`tempmaxp' in 1
+ k+ _- E& C* L   
- d( q% ?5 H& ~+ ~( R' }3 ^( H: r   *240maxp
4 v6 u, B5 b) @   
/ G: f$ U+ A1 G. l& K5 t   
   *-------------ddy240pminp左侧高峰值-----
   
    preserve
0 d8 F( n- }( l& A1 \   
   keep if dx240<`tempminp'
   
   keep if ddy240>0
" L! j* ?$ C" g9 x; A) U) Y2 r   
9 F. Y( r2 ?' K   count
               
( H9 O; U" [1 a                if r(N)>2 {
8 p1 a% Q0 h- N+ _8 ^                 
                pkexamine dx240 ddy240
               
& N5 p) u: `2 _# R                local tempminphddy240=r(tomc)
9 u' q* r% E, i/ {# p9 p                restore
3 Y5 ^, p+ e  i: e- C$ M3 n                 }
! e! _5 B$ T6 i- q                 
0 Y2 @4 j3 b, l0 _        else {
7 ?( b  l0 p5 N                restore
                 di "cannot find pminphddy240"
                 
               
( \' c, N& ]9 Y$ j% ~+ y- }' T% i                 sum dx240
0 n# h+ g9 b- L% t3 m4 G' _                 local tempminphddy240=r(mean)
& G, j- T) C2 ]8 t+ j. L0 w               
                 
: ~" Z  _( f7 A: X" z                 }
1 r! i! ?* H, C                 
   
   
, x5 d' }$ T8 X* `8 \   
8 n& a  E4 a( D8 S! I6 `0 x   
   
' R3 H% ]/ k' y# Y6 t1 |* M2 U5 B        gen pminphddy240=`tempminphddy240' in 1
! c3 w& [5 ?; h$ `* a- T6 L, {       
       
3 j0 f) [$ Z' C6 j3 u        *-------------ddy240pminp右侧低峰值 -----------------
; ?5 Q0 a( T' h6 M    preserve

   keep if dx240>`tempminp'
8 B6 Q; u+ f) D$ ~   
- J6 U/ s" x, e# c4 O$ `   keep if dx240<`tempmaxp'
" A6 J  m; l( }2 z9 w! W& J   
   keep if ddy240<0
! G/ q$ y! Y9 j   
   gen temddy240=-ddy240
8 ^$ g% o6 S& G: c   
, M, I4 ]* U! G' G6 R    count
4 N& d% Y; Y) }1 P5 a               
                if r(N)>2 {
                 
                 pkexamine dx240 temddy240
   
  |$ l* A% x' U" l) h  local tempminplddy240=r(tomc)
                restore
                 }
                 
& c2 P" K) a9 X% x        else {
: }; A) }7 Z5 E! ?+ y* z               
                 di "cannot find pminplddy240"
                restore
                 
                  sum dx240
                 local tempminplddy240=r(mean)
                 
) o# S6 N/ a. }3 g0 T                 
5 g% l9 Q+ ]" q: E                 }
   
: X2 y% d5 ^$ d& x* ]9 G+ r4 _# b   
   
5 q6 ?. m4 g/ B  r* }. x6 `   
( Q+ _& f% @9 k& p) }   
" a! w! J2 T4 r9 D2 g        gen pminplddy240=`tempminplddy240' in 1
% v$ t2 l6 B4 A  }       
        *-------------ddy240pmaxp左侧高峰值 -------------------
. b# T7 }0 o/ q- I+ Q       
/ u- @2 t% J  J# U7 O( u         preserve
   
  [5 N, Y4 V  M6 x   keep if dx240>`tempminp'
   
   keep if dx240<`tempmaxp'
2 s5 a2 e+ W4 q+ f   
   keep if ddy240>0
   count
               
                if r(N)>2 {
2 k/ [+ \9 z8 N- z" y* E! ]1 _                 
0 R9 E# y: U, l; ~3 T3 W                 pkexamine dx240 ddy240
- y7 z( ~  ^  F4 I   
   local tempmaxphddy240=r(tomc)
9 S( i, j4 r+ W& k* {$ W   
' }& v9 a3 z/ G- L4 Q" L2 W. S  F   restore
, @1 D( ~/ c! `1 g                 }
                 
3 [9 m; U( X2 m# B' I' T        else {
9 T' t1 F+ n' u0 U) j               
+ U$ \) u, M! _5 w                 di "cannot find pmaxphddy240"
                 restore
7 \3 j/ D2 x3 l" V) g               
7 q; ?9 r. ~, Z                 sum dx240
                 local tempmaxphddy240=r(mean)
                 
5 s( A# ~3 X: h4 F0 }6 T7 {4 x                 
2 X9 v2 `7 M. ]$ z7 ]9 r                 }
   
   
   
/ v$ O% J  I% H) |0 a4 e   
( u* @; L" D* C1 u   
   
        gen pmaxphddy240=`tempmaxphddy240' in 1
       
/ O# x( G+ d- Q) F% P- {5 n' z% X       
% c6 \# m9 p  W$ T  }' K9 Y& p        *-------------ddy240pmaxp右侧低峰值 -------------------
0 K6 n! N: S( @4 K+ r" D       
5 F' o, E( h2 P5 ^- u+ j9 K8 @         preserve
  
  e5 }8 N8 d- h/ f   keep if dx240>`tempmaxp'
* O) K) G3 a) p3 A% p  |* t   
   gen temddy240=-ddy240
. g* a: d  ]$ h. i- n9 s   
; c) J$ n' Z4 {6 a9 W   keep if temddy240>0
   
   count
               
: P4 w* r% @  Z2 z5 L                if r(N)>2 {
                 
8 U5 B! H4 ?$ [                pkexamine dx240 temddy240
: d3 C" \) ^7 P: u8 _) {+ H   
* v- }# I1 ]7 P1 [/ }   local tempmaxplddy240=r(tomc)
   restore
, V$ F' T: S" t( v                 }
                 
# T% }: x) J( p9 S6 d- J        else {
                restore
                 di "cannot find pmaxlhddy240"
                 
* J" l" M9 S' n  G+ `' F* z9 P( a               
8 R0 ?& G' S1 o$ T& ^: p+ s0 v( J                           sum dx240
5 S4 h1 n0 K- [* {7 Z8 g& R                 local tempmaxplddy240=r(mean) in 1
               
                 
' C2 [7 j1 g- W0 J/ ?/ k7 c                 }
  x+ _, d3 N7 L   
/ P: Q- \/ F( y: f+ ^: S   
   
   
   
   
        gen pmaxplddy240=`tempmaxphddy240' in 1
% e% b+ S0 k( K( I: e       
( P' x# `+ B, D3 |4 u) L0 u) P        *-------------------------------------
1 g( j7 ]  ?" }" v$ g* o+ A   
1 c3 L: p; ]! y; s. f7 {   gen price5d=v5 in 1/20
& k( ]5 B/ y4 _  }9 H3 w      
& Q* K" I; y8 e1 y% V2 L; e: _* X   gen price5yy=_n in 1/20
) k* G* \7 N+ r   
$ V2 P! U% H* g   gen price5y=price5yy*(dy240max-dy240min)/20
- `0 M! g# E  i. K, l   
   *price5y &price5d
; `4 L' r- I3 x+ @& P5 o6 g   
1 h+ e! a8 ]6 l% q2 [' r   gen priceny=price5y in 1
   
( W, j. x: d3 K) ], @   gen pricen=v5 in 1

. V4 g& x8 B$ y! Z     
   *priceny &pricen
   
   
   *-------------------------graph------------------
  
8 o; l1 \0 I: G1 A) m  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 f- O' a4 X# p; l" X3 {! f/ w
* ---------可选项目----------------------
6 p  X# f: o6 X& C) I8 _*(spike yxis amean,lc(blue))(scatter yxis2 amean, lc(blue) msymbol(+) mlabel(amean))
}
) O+ X2 ~3 q7 k9 _+ l* K5 zend

maodeng

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

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

4 o* b8 b4 _& Y( ~哇，谢谢，我仔细看看，这些字母数字代码晕了
2 v" b% V$ l9 M% g4 O; ]" i

喜欢喝冰茶

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

喜欢喝冰茶

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

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

汉水云潇

准备试一试
先找Stata软件

汉水云潇

3 @+ k; d! C3 x; ]$ G
4 M9 K3 |3 i8 P晕呀
论坛复制程序代码
5 K8 M8 }8 m& w3 T$ \8 j乱码夹杂在其中
很烦
对了两次
以为没问题了

8 L+ n8 e; l  t9 E7 z+ K结果
执行后
还是出现
r 199
的错误
7 J/ _, @! I7 Q" u" J8 C

汉水云潇

试验成功，下面是三一重工的图
5 E( d- A5 D" \; C6 Y9 N
! n; }3 K3 i# T: o9 ~

peterbev

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

中秋不圆

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

bknite

这个帖子值得学习。

jeanny

感谢！！！

Hubert2000

明月回春 发表于 2011-8-13 20:17
程序写的不错.
% R9 q& q: F; b+ K8 u+ s1 H不过有一点要注意.通达信导出的数据,只对价格复权,并没有对成交量复权.所以数据需要手工处理 ...

' @+ b5 h( [8 E$ _老哥，也是山寨这篇文章找不到了，不知为何？