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

maodeng

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

maodeng

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

maodeng

. hello 300236
statement out of context
# M" [+ u- g$ k. vr(119);

大师 就是这个

sergio

这贴实在！

awer

maodeng 发表于 2013-1-28 13:13
. hello 300236
7 t- A; H0 C2 m# istatement out of context
1 _0 Q7 L. Y2 j: j: L; W! Hr(119);

抱歉，最近太忙了

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

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

capture program drop hello
6 s5 |% E, ]3 L* ?
program hello

3 b1 r3 f1 @8 K. A quietly {
9 v% s0 ^3 l2 v: D
clear
- |: _& r1 h# Z0 J
$ `, |' z% c& ~5 z; ^5 U( u' lcd D:\stock\test\source
8 F) K3 u7 e) O2 }. ^
insheet using `1'.txt

- Y  ^1 a* X# |9 l) f   keep if v5>0
! o, E0 ^% J, y. A   
) E$ S8 Y6 Q0 s6 Y, j) M   gen vt=_n
   
" ~  v, \  d: e2 ^, J   tsset vt
! @' k/ A. T' I3 Z1 y) E
' A, g9 A* S: Y" [6 i) g3 t   gen zf=D1.v5

   gsort -v1
; U5 d% j  t! M
   gen avol=v7/v5
. k; l- W1 _4 E   
   keep if avol>0
/ k. ?9 ?" l/ h8 P) b1 m( b7 F   
   drop vt
  
   gen vt=_n
" d1 ]! K4 s$ R  W; ?4 ^   
   tsset vt
   
    sum v5 [aw =avol]
; ?; ]: T$ m' u/ p! G$ w
   gen amean=r(mean) in 1
   
: `; _" y9 A! h: E9 d- g   * 总加权平均值 amean
/ j3 h% X* u% k$ x   
   
   keep in 1/240
/ y0 E+ s" s" G! Z# D. W" S
2 a2 u, k+ l0 k% y# C( v   sum v5 [aw =avol]
& E0 L: k2 ~' ]( e5 q& A
* z4 w2 A3 ~/ Q* \; k   gen mean240=r(mean) in 1
+ l+ V# }: }6 b" [# q   
   *240加权平均值 mean240
* I) F. I& v7 k# w9 a   
8 e7 D$ Y' d; o. e5 J0 Y3 r* i$ J   set obs 500
) h) c/ X1 W5 }9 |, r# W   
   kdensity v5 [aw =avol],gaussian n(500) generate(dx240 dy240) nograph
   
$ R% u* I/ F. l   *240 kdensity(dx240 dy240)
   
   sum dy240
   
/ o+ \, ^4 C0 V' D7 Q   gen yxis=r(max) in 1
9 i5 x  J2 `6 n* [. L7 ?: ^   
& R% D9 |0 E) s   gen yxis2=yxis/2
! N7 z' F( I) ]( b   
   *yxis
& p- X! @  @  p5 @! _5 U$ \   
   gen dy240max=r(max)
   
   gen dy240minp=r(min)
' |$ F% P4 I+ O* o  Q' d1 r# P+ j1 U
    gen vtt=_n
3 z9 i+ C+ E, q1 J) [
    tsset vtt

% E6 B# M, M. @! L: m2 y% p   gen ddy240=D1.dy240/D1.dx240
   
   sum ddy240
+ B2 T: k& b  ?( S. i$ X1 F   
   gen ddy240max=r(max)
0 J% j* X8 }' m$ s$ {/ Y   
. C8 o  B" d4 ^8 C/ S, Q+ S9 D! {& y3 G   gen ddy240min=r(min)
3 e9 n/ h3 L2 |4 F' r. d     
1 `1 y% A3 Q0 V: O   gen ddy240g=ddy240*(dy240max-dy240min)/(ddy240max-ddy240min)
   
8 F- @" ?1 k! `   sum ddy240g
   
   gen yxis3=r(min) in 1
   
   gen yxis4=0 in 1
) H: C2 E) \5 k  
. s' e+ x. K8 O, M! H5 d  *ddy240g
3 ]7 [$ o$ }% N% b! @$ Z; ^   
   
2 |' D* y9 N# T2 s- |  preserve
  
  sum v5 [aw =avol]
9 r, i9 f0 J% ]5 w2 o, T/ ^' H  
* n: }+ f- `( F  keep if dx240<r(mean)
. P9 l; o% i4 M1 H' V% S  
  sum ddy240
  
0 H4 u+ e- S0 e  G% H1 s' z  if r(min)<0 {
+ \: J/ v2 \& f6 L        
4 p" V  ?& ^2 w! l  a                keep if ddy240<=0
               
                count
6 i* y: }  n- B$ \$ x               
                if r(N)>2 {
: |" ~! G& l/ P                 
+ L5 b1 F; p. J# Y                 pkexamine dx240 dy240
               
, N# k9 \4 |( G9 Q! ^                local tempminp=r(tomc)
  b; y) @" x- n% q) r. K                restore
  r: n1 B: g3 O. X                 }
                 
8 x+ _, m+ D6 N1 F* m# Q# A        else {
# }; m. }8 {. n+ O% X               
/ p7 f) v) ^$ z                 di "cannot find minp"
  T, G- m. b$ a& P8 m                 
2 a1 a- `9 m/ \- E, m                restore
7 p5 s# ?: U* \; G/ D/ e                 sum dx240
                 local tempminp=r(mean)
                 
7 P5 R) h/ P# B3 Q                 }
                }
) G  a9 K# j) H( j6 W; K   else {
* B  s. J  V( M' ^, l8 @  I% T                 
2 _! M9 L4 C7 P" j- Q+ i) a1 H                 keep if ddy240<=r(min)
3 }% G# E+ a/ G6 s                 
0 F3 `: H. f5 s  T0 B3 U4 c# O                 sum dx240
                 
5 M7 I5 Y. {" R2 ~                 local tempminp=r(min)
                 restore
                 }



   gen minp=`tempminp' in 1
' N3 x! W3 J* _   *240 minp

9 `  z+ s! P0 f5 N1 b  e5 h
   
   preserve
   
/ N& I* W1 ]) i. V6 Q6 D6 t5 a   sum v5 [aw =avol]
   
   keep if dx240>r(mean)
   
$ ~( N/ B; k8 E4 t, S   sum ddy240
$ m  t; H- m; T3 j; R7 k! K   
' i7 `7 a5 S9 z4 a: }8 }   if r(max)>0 {
; ?1 E* V4 ]0 v5 F. [; a# h               
                keep if ddy240>=0
# a( A) L  q! f4 E  A               
6 N- U* _) J" B0 W4 z+ x                count
               
                if r(N)>2 {
                 
) L) V5 K$ g* L' v0 D& G( A: J                 pkexamine dx240 dy240
" S% I( A3 k7 I- J0 i  I               
! y) R8 a+ t: f                local tempmaxp=r(tomc)
                restore
                 }
6 I5 e: X0 a$ l4 o* X7 k7 G+ b% O                 
: O# M+ ^# D9 l8 C$ @- U                 
        else {
- N' [6 X  Z. [5 @' x* o                restore
2 z, @* t1 [! ~3 j                 di "cannot find maxp"
               
               
  v: c7 g1 M. l                 sum dx240
                 local tempmaxp=r(mean)
                 
                 }
                 
                 }
                 
    else {
7 f% u3 O6 C, b       
7 g3 q4 V% i$ W; K) x0 a( h           keep if ddy240>=r(max)
- I; g  ]1 L* p3 c( z: n' L, t          
           sum dx240
          
           local tempmaxp=r(max)
/ k! E( Y% A3 U$ |$ I" @0 P5 V           restore
! z% y5 ~! Y9 p           }
3 l4 e# Y9 C( H0 s# D# j       

  
- X  n+ z4 h) o+ j( G
   gen maxp=`tempmaxp' in 1
   
   *240maxp
   
) y" F9 C6 |2 U2 ~6 a& s   
   *-------------ddy240pminp左侧高峰值-----
, M, a6 D, e& ^  A- V* j   
    preserve
: e8 N2 N9 Y& K3 G6 I7 ^   
   keep if dx240<`tempminp'
) g4 ?  \0 k  M( t   
   keep if ddy240>0
   
   count
- ?% G" v+ |& q$ F  C               
) L- z1 W% t! c/ ^0 b4 [                if r(N)>2 {
                 
% H7 b% E+ ~. L& z! w* I0 n                pkexamine dx240 ddy240
               
' L/ y3 e$ @# {, v7 N1 m                local tempminphddy240=r(tomc)
/ K* Z2 g4 x. R/ n' a( l                restore
                 }
1 x$ }) }, y0 y1 _7 q, G7 H: c+ A                 
/ H0 D# M7 q9 i3 P# q  p        else {
                restore
% [& f+ A* W" s! `1 f) S! o8 x1 e                 di "cannot find pminphddy240"
% p4 b) X% N$ S5 J0 @                 
# Q4 f, q2 g3 m               
                 sum dx240
                 local tempminphddy240=r(mean)
               
                 
                 }
9 O, l, B% y! _! Q                 
   
5 H4 Q4 o  C0 t; c1 }/ Y6 [( y; |   
   
# c/ b6 D# J2 P& J0 S   
   
        gen pminphddy240=`tempminphddy240' in 1
6 _( ]3 `4 T% L5 g$ t# o       
       
        *-------------ddy240pminp右侧低峰值 -----------------
& ~2 e0 s" O6 N) c( r& C) G- g    preserve
% L' D3 Q3 i( \+ @' P2 T: _) Q
: S4 s5 F- M" k; H* }& v- }   keep if dx240>`tempminp'
* l) g0 T- q' u5 \8 _& ]   
   keep if dx240<`tempmaxp'
   
8 `0 Q' i: k" }2 J   keep if ddy240<0
   
2 P0 w$ M! [2 d( L2 M8 l% t* I   gen temddy240=-ddy240
   
    count
               
, v; i( x$ Z  C* w. L. T                if r(N)>2 {
                 
+ W1 n/ b7 J# a% R/ {8 v: W; r. \# |                 pkexamine dx240 temddy240
   
  local tempminplddy240=r(tomc)
                restore
! h0 r: n* z0 p" O% @8 b0 s                 }
                 
# K& T2 o/ p' X- z! W' H2 n9 N        else {
               
+ m, ^/ [. @5 ?- [                 di "cannot find pminplddy240"
                restore
0 Q$ a7 S0 r/ Y  o4 i  N                 
                  sum dx240
& M6 R/ Y8 j& \% i0 t                 local tempminplddy240=r(mean)
6 F6 d1 g- U( z5 B, P1 e$ ~( z                 
                 
1 N9 l$ p: ^5 `                 }
$ N6 `8 z5 a9 {5 ]1 m( _   
   
   
   
4 Y# Y( @) y- X7 K8 F* l   
4 x0 H8 M4 g; t8 `        gen pminplddy240=`tempminplddy240' in 1
       
8 y. ^0 P- H0 P8 _  z        *-------------ddy240pmaxp左侧高峰值 -------------------
5 F- c- J2 H3 k9 J. E       
         preserve
5 ]) y6 G  n$ y# W0 E6 B# ]0 D3 y   
   keep if dx240>`tempminp'
; C' w) n8 j- v1 a4 x, c   
. u4 [) e: L- W8 U  X   keep if dx240<`tempmaxp'
7 K3 T& s  w$ B2 R) e0 {   
   keep if ddy240>0
   count
               
4 o3 @" @1 _% J9 E                if r(N)>2 {
( c# x' f/ T$ G7 l                 
                 pkexamine dx240 ddy240
   
   local tempmaxphddy240=r(tomc)
   
1 r  M% n6 e1 q9 `; Y0 q   restore
                 }
                 
        else {
+ S) O3 D  I- k  q5 P  H) W               
                 di "cannot find pmaxphddy240"
                 restore
+ t+ o& g+ {- P+ c               
                 sum dx240
# N1 l- O8 v- i                 local tempmaxphddy240=r(mean)
                 
1 I5 O4 S6 J* O, r- g( A1 H0 d0 C5 t                 
' \, T& `4 _( X2 E, o9 b                 }
   
   
   
5 e8 i9 q9 q5 y   
3 z; y/ X; `  g& w/ i  ~3 F  ^   
4 e6 n/ s- s" g( {2 _   
! y: _8 Z3 u" s) K7 O' m        gen pmaxphddy240=`tempmaxphddy240' in 1
+ _  v8 L4 C  K/ ^) T* b; q       
- R$ [( t+ ]: j0 D& T# x       
        *-------------ddy240pmaxp右侧低峰值 -------------------
8 `1 y5 y  T/ j' _' W  L. U+ p       
         preserve
  
   keep if dx240>`tempmaxp'
6 S" m# {" I5 ~8 r/ `( B9 }   
! A8 q' m) t4 I! k$ q4 V   gen temddy240=-ddy240
   
1 F% Z2 B6 o6 Q1 j/ u. G   keep if temddy240>0
1 f2 A" H* G, n" Y   
4 b% c; A* W) P. `, f, }/ n   count
               
                if r(N)>2 {
" M1 B- O) D' i% A                 
                pkexamine dx240 temddy240
   
   local tempmaxplddy240=r(tomc)
   restore
                 }
                 
; L0 p: C3 ?: ~        else {
& t5 a' f5 f, T. l1 [- L                restore
$ N3 ]5 W& Z  ^$ S, B7 @1 A                 di "cannot find pmaxlhddy240"
                 
; S3 [1 I; \( e. G/ z  r. k% {               
                           sum dx240
                 local tempmaxplddy240=r(mean) in 1
               
2 W$ N7 e* T8 V; l  @2 E$ r  n                 
' W5 P" V- l, R) U5 S" R) G( @0 t                 }
   
3 C- a+ X9 M4 X5 z  g; a" ]4 T   
   
   
   
   
/ f6 [5 L+ M* D: z2 h! X        gen pmaxplddy240=`tempmaxphddy240' in 1
       
7 r6 B% Q/ x: e, y. u        *-------------------------------------
) t) Y! ^2 Q2 f8 F' C8 @3 J7 ?" X   
+ |- z* i5 a- r6 z- a2 v* X# K: T   gen price5d=v5 in 1/20
" D, Z8 d/ W% o1 C/ }1 J$ o- i      
   gen price5yy=_n in 1/20
   
% C; z2 M5 _7 u/ Z" L# z3 A   gen price5y=price5yy*(dy240max-dy240min)/20
   
5 L! p/ a7 j" x0 F( k) \   *price5y &price5d
, [0 O2 L' n) Q, w3 N: g   
   gen priceny=price5y in 1
   
# o, W/ [; t9 j: W   gen pricen=v5 in 1

% ]3 E' [  y2 d" l     
2 c; P. W4 O1 n# q6 @   *priceny &pricen
   
   
   *-------------------------graph------------------
  
  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 s# U0 @0 v/ r6 U& L* ---------可选项目----------------------
, s& ~" ~3 U& Y7 h+ h- Z/ d! h*(spike yxis amean,lc(blue))(scatter yxis2 amean, lc(blue) msymbol(+) mlabel(amean))
- n# k: ~- ?. W7 b1 N) F}
end
1 [7 F. Z2 _& F) C+ D

maodeng

awer 发表于 2013-2-5 14:40
) O2 ~. W. u5 z6 s, u( G. ]/ }. k抱歉，最近太忙了

0 I5 d  q$ B- h% B8 L300236 我用下面的程序试过，没问题

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

喜欢喝冰茶

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

喜欢喝冰茶

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

+ h& k! d! s* o, h0 R* g+ t  Y& i这东西是开源的？知道有个统计软件叫stata，不过要收钱。

汉水云潇

准备试一试
先找Stata软件

汉水云潇

晕呀
0 i1 @. ?' V0 ^$ G& v$ N7 J7 `论坛复制程序代码
乱码夹杂在其中
很烦
5 ]- n  {, g, W( @- w  Q* A6 D对了两次
以为没问题了

结果
执行后
. v6 E& c: S8 o还是出现
9 t; H- j9 r1 k7 Ar 199
7 d9 W; r  _) o# d& e+ Y的错误
9 f/ A$ M: ^2 M

汉水云潇

试验成功，下面是三一重工的图
" H1 X4 z! n6 W7 b8 W

peterbev

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

中秋不圆

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

bknite

这个帖子值得学习。

jeanny

感谢！！！

Hubert2000

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

+ ]% i9 k+ L# y% a+ a5 M0 Y+ D老哥，也是山寨这篇文章找不到了，不知为何？