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

maodeng

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

maodeng

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

maodeng

. hello 300236
8 O$ W$ ~8 m- s1 L- ]( Zstatement out of context
' F+ h$ \2 b' ?3 S0 l7 wr(119);

大师 就是这个

sergio

这贴实在！

awer

maodeng 发表于 2013-1-28 13:13
. hello 300236
8 y. d0 h- V& G' e% i9 u. r# k: ?statement out of context
r(119);

抱歉，最近太忙了
  k) a% ?# C2 W: D) z# A
# d" `9 g; h! h9 H* o300236 我用下面的程序试过，没问题
6 f) }$ Y, u. L+ Q0 a0 C! l
: W* P9 c1 [/ Q/ f- L3 j1 b" K9 k: n你把 D:\stock\test\source 改成你导出的数据存放的目录，再试一试
0 m+ R, Y& J6 O6 p" J
' F0 D( \; a) I: O# p% s$ N% Ccapture program drop hello

program hello

quietly {
3 r+ A0 f1 Q% M- k
  u$ d$ D/ X8 ~9 O+ P2 _! qclear

: \, E8 E* t# ^5 x8 ccd D:\stock\test\source
) o( _3 t/ X' ]
insheet using `1'.txt
2 f: b" L+ d) |' ^: c" z" I
: Y4 o: s+ T( |7 Q5 ~7 k, p   keep if v5>0
   
   gen vt=_n
   
( o8 Y  t! m. `: }' O   tsset vt

   gen zf=D1.v5

; g4 ~0 D1 v, X; k; g   gsort -v1

   gen avol=v7/v5
' m' ^9 ~; U/ ]. {- }/ |5 w   
0 f* t7 s$ L- k9 ^5 v, v! a   keep if avol>0
% J% N) [: r, o' M8 G. ?   
   drop vt
  
   gen vt=_n
- i7 |; v. Z( ?- }" P   
! x6 p* X/ U, ?# c5 h" S* S0 A   tsset vt
* @( I9 c* b  |7 [% x. W$ F   
    sum v5 [aw =avol]
$ [* ^5 B9 A. Z
0 V" w; d" L6 y2 I, M   gen amean=r(mean) in 1
7 A0 ~, h; Y$ W) ?6 `% K9 F   
# R& N8 K$ T. P+ G& c   * 总加权平均值 amean
   
& U# r( t9 p  f/ M' C   
5 r' T" B& B2 q+ ~% k   keep in 1/240

   sum v5 [aw =avol]
8 c. ^0 ]" W* f% ?( @, D( `
   gen mean240=r(mean) in 1
6 G! N2 x! f5 f% U6 {; |% ]2 J   
% W6 ^% K7 }; F. F   *240加权平均值 mean240
   
! m* |8 T, r5 _5 E   set obs 500
   
   kdensity v5 [aw =avol],gaussian n(500) generate(dx240 dy240) nograph
+ ^6 U: W- [# G5 ^6 b   
8 U. @; m; ?& T* u. j. ]0 M   *240 kdensity(dx240 dy240)
   
   sum dy240
# ^9 ]2 y! R* ]* x& T9 V1 |9 V) g   
   gen yxis=r(max) in 1
+ Y8 Q4 ^6 [5 f4 r* D! E4 P   
3 u# O; s& E! ^5 T+ w, ]* [% }   gen yxis2=yxis/2
5 k& s" [# T2 `2 J/ z# _. l! Z   
   *yxis
   
; o  J4 G7 _) a! h   gen dy240max=r(max)
   
" d0 J  O, T# g, o  @) D8 D   gen dy240minp=r(min)

) h5 ]" y- v1 @2 d    gen vtt=_n

/ j7 r4 W) M9 f! m7 w: j    tsset vtt

   gen ddy240=D1.dy240/D1.dx240
   
$ D% ]% F1 n& L( i4 G   sum ddy240
2 F0 Q6 f9 \6 i' H, y- w- g   
   gen ddy240max=r(max)
   
2 ]2 E9 h$ I8 z% C6 a  y; L; B* B   gen ddy240min=r(min)
9 q$ u- J; @5 X     
   gen ddy240g=ddy240*(dy240max-dy240min)/(ddy240max-ddy240min)
2 f. e+ I9 l/ J5 u, m; q1 c   
   sum ddy240g
  U8 |  |7 M' L6 i  ^* e7 K5 C  [   
( c- V1 `8 Z, Q5 M8 B9 b3 P   gen yxis3=r(min) in 1
% A" K3 B, g3 v( `" T   
   gen yxis4=0 in 1
  
  *ddy240g
5 q: h  E- Q1 C8 ~0 p, ~* }% N   
   
7 f( y, v+ g2 m# c! Z5 B  preserve
2 c" x, H. l3 R# v: U  
  sum v5 [aw =avol]
( `0 R: ~, |" P  s4 `" r  
; h* E% w( ~2 @/ H9 o% [; i  keep if dx240<r(mean)
) t. ~8 k0 O9 D8 Q" ]: Z# Y  
  sum ddy240
  
2 H2 N, {+ U3 C3 ?* b  if r(min)<0 {
% g5 }. Q% W+ l* w7 M# n) R, c* w7 |        
                keep if ddy240<=0
& |* f- `. M# F' W( Z# ^, b( p- y               
) T; d: s$ Y; L- w/ s                count
% A; K  S. J" A$ E               
                if r(N)>2 {
                 
                 pkexamine dx240 dy240
* y8 J+ u1 N3 ~               
$ _- I( B' L7 U: @) h. v1 F                local tempminp=r(tomc)
                restore
                 }
                 
        else {
: S$ _- _  W7 d. P, Z% U( q7 o               
, P2 M( ]$ P' b- A- v8 r                 di "cannot find minp"
6 X; F1 ^; w! K6 @                 
" v6 C4 f& r' v3 Z8 X                restore
                 sum dx240
$ T& ~8 S/ m% ?( n% m                 local tempminp=r(mean)
                 
' I: v- D( l1 {7 }                 }
                }
   else {
' a+ B# M0 M. I" `) O( f, f% P                 
                 keep if ddy240<=r(min)
                 
  h) l5 @( Q& H5 z! Q9 J1 ~                 sum dx240
1 a" E$ @- ^5 E. L6 F: L2 B, P                 
6 Q# @5 u# o5 m                 local tempminp=r(min)
% m! q. x9 I% v) w' v2 ?                 restore
                 }
! G5 Y+ h' D  F( ^9 z" ]
+ D, w9 r' x. ?' ?* y$ q

   gen minp=`tempminp' in 1
  s# I! G4 L* M9 Z' X3 D* A% Y& H   *240 minp
- o: u- M4 \* w4 w! ]) l

   
   preserve
   
1 @" R' K; z( ^3 g   sum v5 [aw =avol]
* C, P: {) I$ S  z. O) z" n   
   keep if dx240>r(mean)
   
   sum ddy240
, {- [! X& }$ q6 l   
   if r(max)>0 {
               
                keep if ddy240>=0
" `- [0 k4 p, b* z  o               
                count
               
                if r(N)>2 {
                 
                 pkexamine dx240 dy240
               
                local tempmaxp=r(tomc)
$ h; X  h9 p( O! T8 s                restore
5 _, `0 `6 N# `0 U* o; D                 }
7 R# h9 ^4 d: q, f: n$ V- f. b) L( x                 
                 
' w! a  ^1 x, l# D4 I        else {
                restore
) V# |. D8 s2 T' e7 N                 di "cannot find maxp"
               
0 t1 G$ m: u! L6 u2 }               
4 [: E. q9 C8 Q( T% m6 Y) x* I                 sum dx240
. t, D2 W9 B2 o9 l7 l  D                 local tempmaxp=r(mean)
                 
- I! z! c: _( p$ C                 }
/ Z. _8 `/ `/ ~* ?                 
                 }
                 
& X" O" ]* H) X7 w2 t    else {
       
           keep if ddy240>=r(max)
          
8 r" X: Y7 f- n; u+ s$ A           sum dx240
; Z; K+ B% [4 d% y: f          
/ F( b- p# E7 d1 i: r7 r& v           local tempmaxp=r(max)
           restore
           }
       
2 R$ D6 h) s; Q7 V0 K+ S- P$ L$ Y
( [$ F1 b/ s: p3 n3 W  

: G; S, ~+ K3 L" o   gen maxp=`tempmaxp' in 1
. B  W) k0 \0 Y   
$ b% u# q* C0 W6 j   *240maxp
5 [9 n' Q+ o) T; U- ~; W   
   
! [& }/ h) L' C1 h7 x% L   *-------------ddy240pminp左侧高峰值-----
  a- F8 X: _; O% g, ^   
    preserve
   
9 H) c% ]5 C) v: x   keep if dx240<`tempminp'
   
   keep if ddy240>0
" J% g  Z5 E- F- P- x7 p   
/ Z) T% s( [: r: Z, {, c, @# O   count
' Z% q4 n- h: t5 }               
                if r(N)>2 {
. p! f8 Q; b/ H. R& d6 I                 
                pkexamine dx240 ddy240
  u3 u$ K: u, e+ y, e# J3 ~  s               
                local tempminphddy240=r(tomc)
                restore
                 }
& g; E- E% e2 t; N                 
- g( v6 t6 h, H/ |. v) y+ ^        else {
                restore
                 di "cannot find pminphddy240"
                 
               
# \3 s! m6 o) [                 sum dx240
7 m" X  q$ ?# _, F8 [5 n/ L                 local tempminphddy240=r(mean)
# ^# {  A2 y3 K1 X" `4 }; ~& `               
) X. }  U6 q! I0 Q. b  e7 [: k1 U                 
                 }
7 F- n( X6 o5 ^( \! D0 ]' H" b; o                 
2 b0 h2 E$ \6 V7 Y   
' J! M6 J) h, n6 l   
   
+ I  n/ L' }% f8 L% o2 c   
- @: G9 r' J) L8 s5 i, A   
& t5 U8 g+ _1 S. A. N7 f3 J/ Y/ d        gen pminphddy240=`tempminphddy240' in 1
$ T) q; _+ R4 w+ v- ^8 X5 x       
       
        *-------------ddy240pminp右侧低峰值 -----------------
! g6 W; v5 y  U! K    preserve

   keep if dx240>`tempminp'
- E2 ?% u5 [6 ~8 q! R   
   keep if dx240<`tempmaxp'
4 S' q6 i- |) }9 G1 P2 n( \   
   keep if ddy240<0
; t. J' t' P; G4 u3 j2 F0 e   
   gen temddy240=-ddy240
   
1 @$ [; f1 n4 H- y# A6 K9 U    count
               
; o) b% p7 w5 [! H- z% d                if r(N)>2 {
                 
                 pkexamine dx240 temddy240
9 Q+ }8 K1 @- {0 O   
  local tempminplddy240=r(tomc)
                restore
" N" h7 \2 o( c+ X& c* `                 }
                 
        else {
               
/ @. i4 H. k9 m5 U' P                 di "cannot find pminplddy240"
  q% K" f" X) V) v* g                restore
                 
                  sum dx240
                 local tempminplddy240=r(mean)
+ R, X6 S2 |7 o9 w: [$ ?8 r                 
& [' Q6 K; B, C/ K                 
                 }
- i$ c3 A' K4 l7 _* c, F  b# J' I   
$ U; B; w8 C7 i# v6 k   
   
, ^, h4 S* p" {   
8 k- l: l5 G+ D$ k9 W; @; j9 `% g   
        gen pminplddy240=`tempminplddy240' in 1
       
        *-------------ddy240pmaxp左侧高峰值 -------------------
% G( c6 }; @8 E       
: g6 B5 d4 V8 V3 ^, H" b0 O         preserve
: u/ ?& L" ~7 {6 d& C   
4 {. D! C: P) z; j$ n0 p  e   keep if dx240>`tempminp'
' F' E3 m4 Z& n" l+ p   
   keep if dx240<`tempmaxp'
   
   keep if ddy240>0
   count
               
6 q; k/ H) l$ w0 l7 y$ d( p                if r(N)>2 {
- p  ^8 ~* r  t9 R                 
# x3 J& L0 ~1 s                 pkexamine dx240 ddy240
   
6 h2 k& n# K" z, N- L8 b   local tempmaxphddy240=r(tomc)
% L1 F, [& c( ^% M   
   restore
                 }
                 
* ^. ]& T: G3 x* {! s+ |, }4 I        else {
               
, ~( N* @4 R( }3 C1 |6 f/ O" f                 di "cannot find pmaxphddy240"
                 restore
               
                 sum dx240
3 A5 e2 M; c0 }1 ^                 local tempmaxphddy240=r(mean)
9 q  k4 \& \, j( f+ _                 
                 
                 }
/ X1 g  i6 o$ u+ b+ u7 y" Y$ ^; B   
( G& I+ M) X' G% Z' b   
   
   
   
   
        gen pmaxphddy240=`tempmaxphddy240' in 1
9 L3 @% t" G4 ]0 x9 s8 ^       
( a) f( i0 n3 L( R9 M$ f       
        *-------------ddy240pmaxp右侧低峰值 -------------------
       
! G, R6 x- H2 S* G         preserve
  
   keep if dx240>`tempmaxp'
& j0 [; N7 R/ A   
   gen temddy240=-ddy240
   
/ n/ \8 N- F! K   keep if temddy240>0
- Y( l! O4 h5 _! @) y   
   count
               
                if r(N)>2 {
8 v! h+ C' @0 E/ G4 q- X( k                 
                pkexamine dx240 temddy240
   
   local tempmaxplddy240=r(tomc)
   restore
                 }
                 
        else {
                restore
9 g8 z6 \: i) F; h/ h$ B                 di "cannot find pmaxlhddy240"
                 
               
% E( \  Y7 Y: z' l( }- n( p                           sum dx240
                 local tempmaxplddy240=r(mean) in 1
               
& O9 o5 o" I: A2 R9 R                 
                 }
   
5 ^& _* z. z" M+ b$ n   
   
9 J: N( U/ l2 Z  ^   
# ?, z6 }2 \3 i8 D* K   
! c* W; V8 T- V+ p0 K   
        gen pmaxplddy240=`tempmaxphddy240' in 1
       
        *-------------------------------------
   
. p9 L" @" ^5 }& E' U5 g* ?0 V; l8 u   gen price5d=v5 in 1/20
& D7 G$ f* b3 ]% r      
   gen price5yy=_n in 1/20
8 _. d! o0 e3 N4 G( o   
   gen price5y=price5yy*(dy240max-dy240min)/20
( E: I' w- X0 l% w   
9 o0 ]; W+ @$ S" u$ Z/ o   *price5y &price5d
   
2 v# `3 z/ N. b* L   gen priceny=price5y in 1
" z( X- D9 f* S3 \   
$ M+ K& s* D, X1 V* r4 ~   gen pricen=v5 in 1
. M; I, Q5 x" u2 |; l
! d/ ]3 A- f. p' P; a     
: f8 `0 P* y, }' {/ e2 S   *priceny &pricen
   
; H. p# Q$ [; t5 [   
   *-------------------------graph------------------
  
& x. W" Z9 }& c( P* P  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)

2 x7 Q' N2 y" s# L* ---------可选项目----------------------
: y& L4 ]5 M: k. a0 q*(spike yxis amean,lc(blue))(scatter yxis2 amean, lc(blue) msymbol(+) mlabel(amean))
}
end

maodeng

awer 发表于 2013-2-5 14:40
- Y9 p& F/ J; Y% I" o抱歉，最近太忙了

+ G$ ]$ G- H( e2 m9 w! E300236 我用下面的程序试过，没问题

哇，谢谢，我仔细看看，这些字母数字代码晕了
( g( {+ x' ?. ~9 m: ?

喜欢喝冰茶

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

喜欢喝冰茶

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

# E( U) v2 p9 d& P* e5 H, S这东西是开源的？知道有个统计软件叫stata，不过要收钱。

汉水云潇

准备试一试
先找Stata软件

汉水云潇

. s3 m- }0 G* t3 R晕呀
' U3 Y" v3 r% L0 a$ M3 s' l9 J5 k论坛复制程序代码
乱码夹杂在其中
: t' H7 l' x; k7 g$ x: N很烦
对了两次
以为没问题了

$ N; l! Q/ n# a0 b8 V% a! g8 e结果
执行后
还是出现
2 Z% i0 t: L7 y, `r 199
的错误

汉水云潇

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

* s" N6 u+ g0 v0 p9 X8 o4 W

peterbev

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

中秋不圆

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

bknite

这个帖子值得学习。

jeanny

感谢！！！

Hubert2000

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

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