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

maodeng

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

maodeng

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

maodeng

. hello 300236
6 x+ C2 Y% B* i! ]& [* `! {3 Vstatement out of context
r(119);

大师 就是这个

sergio

这贴实在！

awer

maodeng 发表于 2013-1-28 13:13
8 f+ R3 [1 ]# n1 {. hello 300236
) e2 ^( k8 O. r+ Estatement out of context
) i/ S7 E8 k  H1 M" G6 ~r(119);

抱歉，最近太忙了
$ {/ \* r% t6 d7 `- l9 A6 _
0 Y$ W: L# x& e# U" j300236 我用下面的程序试过，没问题
* {) o& `1 F% z6 {
5 g  X2 J# L3 \! s你把 D:\stock\test\source 改成你导出的数据存放的目录，再试一试
. ], D6 |2 L% K/ z( K5 O% l
1 B. E4 E: h+ v$ }% p% Lcapture program drop hello

; [7 h( a) E& |2 mprogram hello
' t/ k4 m- w3 f* {# T& \, N, H
% L4 s8 x" o( d/ b5 Q; r$ D quietly {
5 w( q( {6 B% K: e" R
/ l& H0 Y; e7 a, X) `6 m) O/ sclear

cd D:\stock\test\source

insheet using `1'.txt

   keep if v5>0
   
   gen vt=_n
% K5 w3 ?! {/ O7 t+ V6 Z   
   tsset vt
6 k1 ^  D: m* A0 K6 `
, q4 x& z% l8 C   gen zf=D1.v5
" ?' @1 r2 L* U+ F, V, c
   gsort -v1
) _! l* S7 x0 s
. D, A& e7 w: y' E4 v- X* @   gen avol=v7/v5
, o; B5 x/ S: ~   
. q: j  Z/ O& q! L2 x% A   keep if avol>0
1 w) ~# a9 b; O  o* ^   
+ p) d% p4 Z5 _0 W2 ]4 W& U  o0 h   drop vt
2 m" R$ a9 |8 Z5 h6 s5 Z  
3 l; F5 t! D2 c5 J5 w   gen vt=_n
   
5 ]" C3 C* s1 z& Q4 E6 V% _   tsset vt
0 u( C1 M; Z- i7 A# J  S3 y   
    sum v5 [aw =avol]

   gen amean=r(mean) in 1
( l; \: X& {- w4 `% S; ?0 ^% T   
1 ^" M1 ~( H5 Z; O/ H1 G   * 总加权平均值 amean
   
   
   keep in 1/240
9 ~8 ^; l8 a( j' Z
4 S3 w" ~& C0 ]7 ~& G   sum v5 [aw =avol]
% O  _* w  W3 W% B5 h$ V  N, B+ t
   gen mean240=r(mean) in 1
   
! A) p9 A% A  D; H2 F   *240加权平均值 mean240
   
: I+ i  u+ h! r3 a& T; L4 M, v   set obs 500
' t5 j1 h6 C2 q& w1 K2 [5 q8 W' I   
( w2 [2 |- b# ]6 b   kdensity v5 [aw =avol],gaussian n(500) generate(dx240 dy240) nograph
: X. A. i, J, t# L. t$ \   
/ q: M8 N8 C8 R2 T* u6 e% n   *240 kdensity(dx240 dy240)
2 I% E6 ]( p% d6 |  C/ @' F   
9 F6 w4 E" E: l   sum dy240
   
$ @2 x: l3 C6 P6 A7 J, O   gen yxis=r(max) in 1
0 }3 K1 F. @7 m% s+ T. Y   
4 ^7 F' b2 @! v3 ~/ L7 {. E" s   gen yxis2=yxis/2
  M( j: Q" t) p/ b" O   
   *yxis
  ^% {, B' W: w   
   gen dy240max=r(max)
   
   gen dy240minp=r(min)

    gen vtt=_n
: H$ O1 [# @0 |- d! L: W
    tsset vtt

9 H) L- a: Z0 c; a* f0 e) ^   gen ddy240=D1.dy240/D1.dx240
   
2 b2 o  T9 X2 l; b' j$ r- r9 \   sum ddy240
$ l: L1 S3 ~) z+ p   
   gen ddy240max=r(max)
   
( U) \9 Q1 W" n* J. z   gen ddy240min=r(min)
     
( O$ Z5 g. Q) p. K+ v   gen ddy240g=ddy240*(dy240max-dy240min)/(ddy240max-ddy240min)
   
   sum ddy240g
   
   gen yxis3=r(min) in 1
   
* o1 M* \/ ~$ s+ h5 B   gen yxis4=0 in 1
  
1 g+ w$ X: M% C( k' D' Z1 ~2 ~, o# c  *ddy240g
   
6 E& I, K! d$ Q( p6 w   
! x4 }3 o* q# U% D  preserve
  
  sum v5 [aw =avol]
+ `* j8 l0 v0 n% P  
  keep if dx240<r(mean)
# G& T* f6 @( B: E& P' R9 u8 ?/ v  
2 L0 e3 o. Z) u4 q  sum ddy240
$ ~+ G7 {2 ~1 a  x& D1 p% Q  
  if r(min)<0 {
        
                keep if ddy240<=0
               
" C, e( r* N9 n6 O. I' `                count
* @4 F! y& n. c3 [$ q               
0 r7 E7 b5 u/ Y2 ~% N                if r(N)>2 {
' V) {& b6 A7 s& [! X! r. s                 
# }4 A/ B9 ]2 w* d+ {  v6 O$ a                 pkexamine dx240 dy240
               
                local tempminp=r(tomc)
                restore
9 s$ L, P) }* {- ]- H                 }
                 
        else {
+ H( j# X5 y! n* \               
                 di "cannot find minp"
9 R2 T: y+ h, Y                 
: f: B+ N4 `3 d9 V                restore
                 sum dx240
                 local tempminp=r(mean)
3 I) @6 Y* o/ |- Z- c2 j# `, r                 
9 u8 d2 V# Q! E; L# r5 |                 }
: u9 @# s% q. T; Y6 r                }
# U: O3 w0 x, L7 C; h   else {
2 M" E' p) t; `5 L                 
                 keep if ddy240<=r(min)
                 
                 sum dx240
- c- d6 H4 d! U+ R6 W; o                 
                 local tempminp=r(min)
- U3 o' z7 L; F7 I* c                 restore
$ I1 a  F/ _$ S. o                 }
4 ^  n6 U! q6 T% ~- q5 r

. Q; I- }6 l, w  y7 n* C
   gen minp=`tempminp' in 1
! E- V2 h4 U, @8 c9 q- c8 z   *240 minp
1 z# ?! M+ Y) n/ q' x9 Q2 g

   
* a6 A/ M1 r9 ^& O8 ?( `; q   preserve
, R0 m0 h+ b- Z! _, ^   
* G& |) `$ \, o) l   sum v5 [aw =avol]
   
2 z6 B, G( D4 u( p6 j   keep if dx240>r(mean)
   
4 o- q$ i/ Z$ v2 ~2 F3 }   sum ddy240
   
   if r(max)>0 {
               
                keep if ddy240>=0
$ x$ z5 ?' I/ Z+ V0 X' q               
8 v& R- y" z8 Y3 {" Y. z6 R6 l                count
               
                if r(N)>2 {
                 
- Q8 P2 K, g+ [5 k                 pkexamine dx240 dy240
               
                local tempmaxp=r(tomc)
                restore
                 }
! j6 G; Z. Z$ |2 d                 
0 J# T. M6 t6 {% D# `                 
        else {
5 s: H* }- f1 D( Y5 l7 t                restore
                 di "cannot find maxp"
               
7 V) E( R: s5 L. P, s- \3 L" h               
! w$ a$ P+ q7 X. u                 sum dx240
- t7 q( j3 J' D# l, y4 M                 local tempmaxp=r(mean)
                 
% |- G; q( L$ Q+ x                 }
, q8 F) a, H1 d* L. r                 
                 }
, Z% m' j: E5 D4 e* S3 W                 
& ]5 p( k) B  p) a6 q' y( z* I    else {
       
           keep if ddy240>=r(max)
          
           sum dx240
# D7 j3 D* p0 s" V5 }          
. t5 @$ N4 K, w           local tempmaxp=r(max)
           restore
4 _$ U( n) _7 O           }
  p: ]) F/ O6 b       
; Z. l9 h* Y+ |: e
  
( t7 V5 Z) w9 B' b/ ]2 N
   gen maxp=`tempmaxp' in 1
   
   *240maxp
   
   
   *-------------ddy240pminp左侧高峰值-----
   
    preserve
4 d1 p1 r% ?, ~3 f, r! ]   
   keep if dx240<`tempminp'
4 B5 Z" k) ?* W/ }   
2 A% ^* o/ x4 Z) ?3 U, ]3 L   keep if ddy240>0
   
  Z6 P, {; k5 j   count
" k, x" y9 x, v. `& g- B               
* B7 c, M# z, t0 M0 t                if r(N)>2 {
                 
7 _3 v+ C. {% b- ^/ y6 e. w% {# v  Q                pkexamine dx240 ddy240
               
                local tempminphddy240=r(tomc)
                restore
                 }
                 
        else {
                restore
' _3 D: }) `4 [                 di "cannot find pminphddy240"
                 
               
                 sum dx240
+ q; {1 s7 K# w; V9 k- j- L                 local tempminphddy240=r(mean)
               
                 
                 }
                 
- G% F7 j! D5 D# A: _   
* Q7 c, Z% h& q" T   
   
   
, g% G8 H/ G5 B& l   
! Z1 v' ~% E( M) \        gen pminphddy240=`tempminphddy240' in 1
       
' J+ v. e) n% V- Z6 U' y; F       
        *-------------ddy240pminp右侧低峰值 -----------------
    preserve
  n- M8 }) h- q: ]$ u. A
& k# ^, O. W* P% U2 R. q+ |   keep if dx240>`tempminp'
, j* H: L+ l/ p+ h9 J" Y   
" G8 {% {+ U6 I* }! r   keep if dx240<`tempmaxp'
   
5 ^# Z7 }9 R8 E7 ^8 y9 B: @   keep if ddy240<0
7 u  T, I# ?4 O$ i* N9 F   
2 U5 o3 g2 z* c; m   gen temddy240=-ddy240
   
    count
               
                if r(N)>2 {
9 r/ ~: A4 F2 n                 
                 pkexamine dx240 temddy240
/ E6 H' d% q: K' D2 H( @   
  local tempminplddy240=r(tomc)
# |7 ~+ c7 V! [5 U                restore
! J' E- y- |$ A" Z9 r7 N. a5 q                 }
. q- Q6 v) x* g# T$ p                 
5 _- l0 Z6 Y1 g0 g7 f        else {
$ A) s  n1 o9 M0 y- }7 |               
                 di "cannot find pminplddy240"
/ [% M/ M/ b, N) n  r) e                restore
                 
                  sum dx240
                 local tempminplddy240=r(mean)
1 _2 g0 C. R+ E( A' L; r6 R0 ]                 
                 
% q# g6 M1 c5 p5 @                 }
   
   
/ p, D7 R. k% Q2 Z% n7 I1 ^   
   
   
- c% o$ z. `& A  i5 L        gen pminplddy240=`tempminplddy240' in 1
; V/ k: t* h2 v8 K9 f       
% W1 o7 U8 ?  [; m4 u! W. w% b        *-------------ddy240pmaxp左侧高峰值 -------------------
       
9 E6 g4 H$ `" {6 x3 B         preserve
! _3 N# c1 m6 M- L. B# z   
. S8 ]- x% ?; U( B7 y   keep if dx240>`tempminp'
   
   keep if dx240<`tempmaxp'
3 ?5 F/ M& ?" N9 c% @6 s   
% U6 o1 @, Z  D* D  H2 o! ]7 r   keep if ddy240>0
   count
               
7 ?9 K0 L' i' L  S& y6 E" X                if r(N)>2 {
                 
                 pkexamine dx240 ddy240
   
5 n( u/ j/ n# F   local tempmaxphddy240=r(tomc)
2 S7 e- t8 n, h$ g$ G, n3 C- ]6 c5 h   
/ y7 c/ v+ `1 r   restore
. C5 ]3 U+ c1 o; g  k# d                 }
- x' U0 U2 y8 d8 A. `                 
& Q* P0 I, r5 s% @        else {
, x9 b  G# _" M1 b1 T               
( @# v; H# H: D- W1 X. t, G& ?                 di "cannot find pmaxphddy240"
                 restore
9 p' N' l3 f/ R' A: z               
                 sum dx240
. z4 H/ F/ P$ X                 local tempmaxphddy240=r(mean)
2 H$ ?/ Q1 X1 }' h( c                 
                 
2 L* k) F( |5 O7 R1 V" n                 }
2 h6 b( u/ u- o& p! A1 @- E   
# [- X' F: p. J5 ]: j9 R% u   
- L  g3 N4 ]3 A9 A7 b! _) n6 T   
4 K* B' P) C6 A   
   
% X+ a. @  ~% [8 _) {   
        gen pmaxphddy240=`tempmaxphddy240' in 1
% W, f+ `2 m2 A' j+ Y       
       
        *-------------ddy240pmaxp右侧低峰值 -------------------
1 ?8 o8 E& ]- r, G       
3 v: N' ^" i# T$ t         preserve
2 r1 K0 Y+ N9 a  R1 t" t  
( u9 z6 z2 Z# p. Y5 O! Z6 q  g" L9 I2 i   keep if dx240>`tempmaxp'
   
   gen temddy240=-ddy240
   
   keep if temddy240>0
   
; k5 V. I4 i; [, M) ?* z, U   count
  f9 e; i0 y# g9 J: Y' i               
2 I. G: i! J6 Y) j- l                if r(N)>2 {
                 
9 M8 _' w7 [! ?: A/ r' k) k                pkexamine dx240 temddy240
   
   local tempmaxplddy240=r(tomc)
- ^: @3 c0 B; T3 r0 I$ o( q   restore
                 }
                 
( r4 n7 t! c. s* F7 B' H+ N        else {
                restore
. |: D! g) V5 E' \( r                 di "cannot find pmaxlhddy240"
                 
               
                           sum dx240
7 z1 V9 ^% q" ~* g. u  b                 local tempmaxplddy240=r(mean) in 1
/ e" B4 ^1 x8 N               
/ b9 X; H7 Q) D+ [' x                 
                 }
6 j3 w. Q1 j' W0 i4 C2 ^3 |/ Z8 {   
   
   
- x) h. Y( x2 D% \! d# o   
" T' ]! h' n" o8 E0 z  H   
# @& e! r" c8 ]! s7 g: y8 j   
7 v7 B/ F% ?) @3 W+ D% O; f        gen pmaxplddy240=`tempmaxphddy240' in 1
9 U2 `% M% W9 u8 d9 ]: h2 J! r% [8 @       
+ W$ _- }" w9 p6 _  k/ J' `6 ?2 |+ Y        *-------------------------------------
   
! U* y3 q4 @' m; {3 y   gen price5d=v5 in 1/20
      
7 L5 u" z  L; k/ Y   gen price5yy=_n in 1/20
   
   gen price5y=price5yy*(dy240max-dy240min)/20
   
   *price5y &price5d
   
   gen priceny=price5y in 1
" r5 P9 T2 X2 ?( R9 O: w: I- h+ O   
- g" |+ G. U. o0 f: Q   gen pricen=v5 in 1

     
) Z" b, k. T* x2 X   *priceny &pricen
   
   
& B9 H/ E5 Y4 L' I0 c* n3 ^; l4 e   *-------------------------graph------------------
) k' i5 X4 a8 B, K5 Y# y  
  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)

! D, R5 c* _( J1 Q* ---------可选项目----------------------
- g: l( I( _) T*(spike yxis amean,lc(blue))(scatter yxis2 amean, lc(blue) msymbol(+) mlabel(amean))
}
( B- j- O1 ^) Qend

maodeng

awer 发表于 2013-2-5 14:40
: H" ^, ?/ G2 E0 X5 f, A/ n抱歉，最近太忙了

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

哇，谢谢，我仔细看看，这些字母数字代码晕了
/ i3 N! m0 N, o

喜欢喝冰茶

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

喜欢喝冰茶

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

* w# h3 X- X, I4 [
# O! q1 m: U2 e) o* t' R" k这东西是开源的？知道有个统计软件叫stata，不过要收钱。

汉水云潇

准备试一试
先找Stata软件

汉水云潇

1 J( g! J, L% X) o1 H晕呀
' m! I9 I) v1 }& q: [论坛复制程序代码
乱码夹杂在其中
很烦
对了两次
以为没问题了
/ ]7 d2 Q* {! X! K
结果
执行后
还是出现
7 c( Y: \6 Q0 ]+ J! j0 Wr 199
" P9 w* K: p( m7 Q+ ^3 N/ |的错误
. n5 h& \: ?6 q# T' P

汉水云潇

试验成功，下面是三一重工的图
  c$ a/ {5 X: I& O
3 n1 g+ u# H% ^

peterbev

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

中秋不圆

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

bknite

这个帖子值得学习。

jeanny

感谢！！！

Hubert2000

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

老哥，也是山寨这篇文章找不到了，不知为何？