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

maodeng

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

maodeng

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

maodeng

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

  R& s% z% U9 R* x% d) k9 Z* M4 V/ m2 f大师 就是这个

sergio

这贴实在！

awer

maodeng 发表于 2013-1-28 13:13
# p; ?) |) H0 {1 X9 x- j: B( D( R, J. hello 300236
statement out of context
" O' ~, o' M2 f! b7 g- f+ Kr(119);

抱歉，最近太忙了
" g- D1 A* _1 w1 E9 ~2 d4 A- u* r
6 m4 B+ J& }; m( D" ]300236 我用下面的程序试过，没问题

你把 D:\stock\test\source 改成你导出的数据存放的目录，再试一试
/ `/ A# J- a0 d2 E
; ]% ^7 U% T, v) n1 ~: Gcapture program drop hello
5 Z- A+ b) Y) Q' Z$ d7 K
program hello

quietly {
2 R, w# A4 Y/ O
% J+ E$ ?/ n1 X2 B! n9 Q- f" ]8 |clear

cd D:\stock\test\source

insheet using `1'.txt

   keep if v5>0
   
   gen vt=_n
% o% \$ Q, c3 @9 l4 W   
   tsset vt
0 h, ~6 S$ \, c! `& n1 m: a4 E# X
! Q9 Q( E& A* i   gen zf=D1.v5
: H0 n: J. m9 U& ~8 J! }
   gsort -v1

   gen avol=v7/v5
   
   keep if avol>0
* G' h  Z$ x! u: N5 Z" y/ T   
   drop vt
  
   gen vt=_n
   
   tsset vt
  j; ?4 B8 }- x0 A7 y+ g: H+ ]   
    sum v5 [aw =avol]
& G3 a* @$ a$ S8 ~, H* r
' U+ D, G( Y* `! Z   gen amean=r(mean) in 1
$ D: w; e% ~+ o# o7 t5 O  q( j1 [& P   
6 F2 O5 G8 _8 C, Z   * 总加权平均值 amean
1 K) ^- D3 q' }! a, x" h, ]4 Z   
   
4 P1 U2 [8 A% ~5 X1 H, x   keep in 1/240
1 E' `* }' `% H
   sum v5 [aw =avol]

   gen mean240=r(mean) in 1
5 `# W9 O: r8 [5 y   
   *240加权平均值 mean240
   
# h, `: I4 i9 L$ L/ p8 l   set obs 500
   
   kdensity v5 [aw =avol],gaussian n(500) generate(dx240 dy240) nograph
   
3 \) I  L. f! `# _) ]& d   *240 kdensity(dx240 dy240)
1 p1 B( f' }# e; h. K3 s0 c/ M+ E   
   sum dy240
   
3 Y% s# M8 p* Z  G6 Y   gen yxis=r(max) in 1
   
   gen yxis2=yxis/2
5 I  g% h# z- A  [( ?9 ~   
   *yxis
   
# y0 s! Z7 R- f# {7 n+ g$ ^8 J0 ~  e   gen dy240max=r(max)
5 R6 N$ c3 `2 w! U. m   
   gen dy240minp=r(min)

$ H5 R( t6 I  r4 C    gen vtt=_n
$ b* F# `  D3 \$ y
    tsset vtt

7 ?$ i" `+ O1 {- L$ I/ ]3 A   gen ddy240=D1.dy240/D1.dx240
! U5 c* A. R/ k   
   sum ddy240
8 X. I% |/ K' w% v9 y; f& b   
3 E4 Z) M# U, `2 _# Z   gen ddy240max=r(max)
   
   gen ddy240min=r(min)
     
   gen ddy240g=ddy240*(dy240max-dy240min)/(ddy240max-ddy240min)
   
   sum ddy240g
   
   gen yxis3=r(min) in 1
   
$ m3 ~0 H# v0 r0 f' X8 U   gen yxis4=0 in 1
  
8 Y: R: ]4 W( _9 ~: C  *ddy240g
, l3 o7 d6 h9 j   
   
2 N* b! _& b6 Z+ Y% d% X3 d  w  preserve
' w5 r/ W2 c* p+ [7 R6 w  
  sum v5 [aw =avol]
  
7 E; }. K1 Q! y5 \; T$ c7 [' Q  keep if dx240<r(mean)
% W% w5 I) {- b+ N: g  
  sum ddy240
  
  if r(min)<0 {
        
* n: n2 w& `8 C$ S$ P3 `                keep if ddy240<=0
               
                count
               
' r+ L7 d6 I, G1 @( A. g3 g2 @                if r(N)>2 {
- Y6 k0 E6 v- w7 x                 
                 pkexamine dx240 dy240
               
                local tempminp=r(tomc)
( `0 k$ Z- _( G8 H& s: K8 t                restore
                 }
                 
        else {
               
, Z3 B3 P. [5 {& X+ }# Z- v& Q+ R                 di "cannot find minp"
                 
* Y/ ]( z5 w6 ?: [; A: q6 M                restore
                 sum dx240
                 local tempminp=r(mean)
                 
! `3 l+ U7 U& N1 K8 @% w0 ^                 }
                }
% J# U. q# Z, I+ n   else {
; d8 C, x3 q; |) V2 t                 
: Z! c6 r& I8 n; Y  k                 keep if ddy240<=r(min)
% E. H( C3 M( E; K                 
* L) x0 k) [- w7 m                 sum dx240
3 L+ w2 j$ P8 Z3 |                 
7 G* G. F9 H2 W+ I9 X& M- T                 local tempminp=r(min)
# j: X$ J6 C  @( {! \; a                 restore
                 }
7 h; G2 ^6 G2 E) c$ ~' v6 g

; N  A% F7 b" u: m# j8 c
8 J8 E( M0 t0 l* c  _6 D$ M) |   gen minp=`tempminp' in 1
. w4 y$ ?( k* }4 E% c& [( }   *240 minp
9 a9 D8 ~& k: I* Q

   
9 i" y! Z+ @3 T5 \- d$ H   preserve
( E- s7 d3 K/ W  h( P  ^7 I- \   
* R8 q0 _% ~% M9 W, p& D9 }) D   sum v5 [aw =avol]
   
   keep if dx240>r(mean)
; l& F" g3 }  Q9 P   
   sum ddy240
$ M* O% M% H: ?3 {, l8 Y   
( t3 A8 Q  W$ \/ @% [( V   if r(max)>0 {
               
                keep if ddy240>=0
; c: F' k7 g. A5 L( H* g" R1 n& ?               
( \+ W# k3 j5 w0 e7 E' j" F                count
               
9 A& U% f5 f; \4 m" N# k                if r(N)>2 {
                 
                 pkexamine dx240 dy240
               
+ O0 F0 {4 v6 X6 z5 P5 l                local tempmaxp=r(tomc)
                restore
# p- i; @6 H6 C* {                 }
                 
2 N7 H/ ?0 ^9 N                 
        else {
( V4 G; K  J7 @3 ~! ?                restore
# I8 \5 Q# ]3 _8 {$ }8 Y9 {                 di "cannot find maxp"
( D) J+ }- A" u5 {$ c$ L               
# T" S- K' o, Y1 X& _# }               
                 sum dx240
                 local tempmaxp=r(mean)
# R1 @: o/ F7 e8 {( J' r* r                 
0 n# f, K1 A1 `" P" ^3 O* q                 }
                 
5 f1 E+ C4 i0 p, k                 }
                 
3 k# S- _+ |" j: \8 @    else {
0 \. r7 p+ V: R' B" T$ z       
           keep if ddy240>=r(max)
8 q4 m( d2 {! u* z' T          
           sum dx240
0 q1 d4 R: D; L' p          
           local tempmaxp=r(max)
0 o0 g; ~4 z) L1 R- E           restore
* J8 J: _% J( z2 J: O           }
       

: P" d: ?! S8 E  

# M3 h- v: L. X- ^, y& K' Q2 ^   gen maxp=`tempmaxp' in 1
9 E2 y$ i3 `7 z: q3 W   
& i; y$ A) D/ \: r   *240maxp
   
   
   *-------------ddy240pminp左侧高峰值-----
   
5 C* L5 d* L1 y7 ?! Z, L! J5 D    preserve
   
9 Y% n, |) g, U+ Z   keep if dx240<`tempminp'
   
   keep if ddy240>0
   
   count
               
( ^& p# U: t1 A$ a                if r(N)>2 {
                 
                pkexamine dx240 ddy240
' y. o; q( i6 W( E               
( m7 \9 Y/ l  V                local tempminphddy240=r(tomc)
                restore
, ?% c, n0 j1 S                 }
5 x  L, A: y' e2 l                 
" u, \, X4 u+ q# D        else {
( W4 o- U/ m5 j' [& F6 [" C, c; E                restore
( B4 P# Y. Y- L                 di "cannot find pminphddy240"
0 C+ X7 L# d2 [% y& s& c5 E                 
( l7 O# x- N, U( k3 D0 ?0 g               
                 sum dx240
0 r- H& n' ^9 a7 m# K) a& f4 g                 local tempminphddy240=r(mean)
               
                 
                 }
  O% b0 z/ D: ?                 
# m. v: N4 c9 h- V$ k7 r' |& f0 i   
   
% P8 h' r, m, T6 v   
   
& e4 ]" N: G- f+ O   
2 ?* ?3 n6 h8 m! H6 d# }        gen pminphddy240=`tempminphddy240' in 1
       
5 d  t. f) a. _! g) d2 O2 a5 b       
        *-------------ddy240pminp右侧低峰值 -----------------
/ i6 i3 E' T" z4 {. d+ y    preserve

% T) {: |, \% W2 r) M& B- x# h   keep if dx240>`tempminp'
+ _% R. Q* Y. z# c1 F1 p   
   keep if dx240<`tempmaxp'
   
( B# |  B& \- y7 w6 W9 I  v/ N   keep if ddy240<0
! C4 X$ v1 K- |5 U( v   
   gen temddy240=-ddy240
( {+ m) e1 F. ]% t! W; [   
    count
" a1 l. \- h! o& V& P" l               
                if r(N)>2 {
                 
                 pkexamine dx240 temddy240
' g  |* v$ Q% I7 ?' K3 k   
) ~2 S2 ?& A% @' k* E1 w% T  local tempminplddy240=r(tomc)
9 P; G- U' e! N* K1 n) u                restore
+ ~( B, F! {+ {7 X  F& z% C% c                 }
) P4 G* W+ Q' d, |" n                 
) {  u9 D5 d  t: E8 O        else {
* ^; c, `$ D0 U, v& }, j               
) V' B  Q  M4 c' c                 di "cannot find pminplddy240"
                restore
                 
) L; w& h3 }6 D+ Y& g& c                  sum dx240
  i9 B5 a# m( C( F2 Y& P                 local tempminplddy240=r(mean)
0 S- b6 [& t6 H6 a* Q3 V$ [                 
                 
                 }
   
- P5 a+ h* f  c% x, E9 y" g7 q8 u% ^   
  Q5 W0 c' B7 s! ]+ b   
   
$ V7 o2 |# I. g" V+ _; y   
# F# d  j  a: g1 x' D" _' i1 S5 O' P        gen pminplddy240=`tempminplddy240' in 1
9 s. I: j- V5 q) r0 X, d- u       
        *-------------ddy240pmaxp左侧高峰值 -------------------
4 a* I5 _& ^. x! m       
         preserve
   
   keep if dx240>`tempminp'
" @/ j! E. C8 C3 w7 W) ~   
   keep if dx240<`tempmaxp'
( d& r9 ~' k* ~2 d9 o$ z( u   
   keep if ddy240>0
+ [& [4 a* ]$ D: q) G2 r   count
               
- }- s$ |( S: O- O0 R                if r(N)>2 {
& P; \. D( R0 v2 [3 B: B                 
$ K% T; t7 C3 N7 U% _                 pkexamine dx240 ddy240
   
- Z1 Z/ Z2 y9 J3 t   local tempmaxphddy240=r(tomc)
   
" G/ f, }, q" t   restore
* w- h& E' D* G  R7 o                 }
                 
        else {
/ s4 n+ V) X* u( f4 ^& m               
0 q3 {$ H( f" y/ ^& c                 di "cannot find pmaxphddy240"
3 j2 d, }& i/ [0 [4 o4 N                 restore
               
# w# q1 \. p% v. k0 F; h                 sum dx240
; l7 @$ p1 O% W5 j% L                 local tempmaxphddy240=r(mean)
                 
                 
                 }
9 z# i. ^2 z  t4 {7 h/ r9 G; g   
   
   
   
   
6 H- P) |  K& w- U8 L, Q* Z& H   
7 {- ]# @$ ?% X# ^: ~        gen pmaxphddy240=`tempmaxphddy240' in 1
       
+ x( w* n& d: C5 ?/ ?; O' ^# q4 u       
        *-------------ddy240pmaxp右侧低峰值 -------------------
       
         preserve
; \+ W7 c& Y$ Z9 ?4 Q; T8 `9 w4 b8 t  
   keep if dx240>`tempmaxp'
   
4 b* @, q/ m* c) y* c# ^& Y; [   gen temddy240=-ddy240
. P+ x% u- f8 k0 k6 Q   
$ o% S, o* R! M0 ]  \4 ]4 E/ ?5 C6 O' p   keep if temddy240>0
: x  d. [8 v/ `3 A) ^2 o, V   
+ I4 y! i2 N" Y; ^   count
               
                if r(N)>2 {
$ j3 q( b+ R8 U0 w. Y+ d                 
* A& y7 \$ n2 z  T2 i                pkexamine dx240 temddy240
/ ]+ V8 D' M6 A$ N  J8 f# x+ l& t+ X   
  x# A4 Y3 s  O% z% L& w& \, x   local tempmaxplddy240=r(tomc)
9 q( a# t2 D9 A   restore
7 w5 m4 e% L8 F                 }
1 y# I6 k& Y( I' a; {9 v                 
        else {
; t$ M4 u. G  ]% Q                restore
                 di "cannot find pmaxlhddy240"
: F( d7 v( {' v; P5 L, a                 
- T- Q/ L' Z/ v& H% z2 r               
2 z, L1 q1 w* P/ p                           sum dx240
                 local tempmaxplddy240=r(mean) in 1
$ x" \- W4 l  P) {# _1 {0 k               
                 
9 ^* u5 T8 ?. y& Q' e                 }
: h* M/ j0 N' }8 U3 @* X   
& ~: g4 s; k7 U+ i! q   
   
   
   
   
        gen pmaxplddy240=`tempmaxphddy240' in 1
* z1 i( }$ A0 u. I9 o2 R       
) r: b* o& t5 {6 N3 ~        *-------------------------------------
   
   gen price5d=v5 in 1/20
$ W) {, ]) M) n4 `  G% U      
. x. x8 g) m# o# q, J6 E' x6 `2 x- z   gen price5yy=_n in 1/20
2 D0 D8 Y, O% j4 t3 n   
   gen price5y=price5yy*(dy240max-dy240min)/20
9 d) B" U9 g) D4 X: K+ a, r   
" u% L$ V6 S0 A$ _) P   *price5y &price5d
. x! `2 o3 Y7 w( C, F* u2 u+ t6 n   
   gen priceny=price5y in 1
   
. Z0 c2 d1 |( _/ v   gen pricen=v5 in 1

     
   *priceny &pricen
% Z% y0 M' f3 X5 H& p   
2 [- Y7 z9 Y& {' _- z6 v: h; I   
6 b# y& }1 r: ~5 Y2 h   *-------------------------graph------------------
  
: ^( O! X) B# T+ p9 c) 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)

# i5 ]* \/ B! C' M+ f( A* ---------可选项目----------------------
( T2 M+ A$ i! |* _*(spike yxis amean,lc(blue))(scatter yxis2 amean, lc(blue) msymbol(+) mlabel(amean))
3 e& T3 J# i, F' A# Q# z}
' X3 l* U  Z4 `& m( c) Hend
" T: a" I- Q/ f

maodeng

awer 发表于 2013-2-5 14:40
, _! q3 [+ ]4 R抱歉，最近太忙了
" y1 g+ c% [% ~! a: u9 g6 e
- I# x* g5 @6 d. d3 w300236 我用下面的程序试过，没问题

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

喜欢喝冰茶

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

喜欢喝冰茶

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

& d! Z3 b4 R/ {# I4 i9 U这东西是开源的？知道有个统计软件叫stata，不过要收钱。

汉水云潇

准备试一试
先找Stata软件

汉水云潇

1 {7 Z/ {; v  t% W- G晕呀
论坛复制程序代码
" p: c# x: L# `乱码夹杂在其中
很烦
对了两次
以为没问题了

结果
执行后
/ Z* g' y; |4 A# H还是出现
& g7 d8 u4 d1 e4 i% Z3 {r 199
, X0 h' I$ W7 R: r1 M的错误
+ G! j9 w. G9 z5 V+ u3 F: V  _

汉水云潇

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

peterbev

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

中秋不圆

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

bknite

这个帖子值得学习。

jeanny

感谢！！！

Hubert2000

明月回春 发表于 2011-8-13 20:17
程序写的不错.
. I4 q' N; o& s. [% _1 l/ Y不过有一点要注意.通达信导出的数据,只对价格复权,并没有对成交量复权.所以数据需要手工处理 ...

0 t! E- T6 y8 M1 ^老哥，也是山寨这篇文章找不到了，不知为何？