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

maodeng

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

maodeng

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

maodeng

. hello 300236
6 B8 P% Z, M  |' kstatement out of context
6 u3 A0 W! E% U0 o+ J+ wr(119);
/ f5 M: T  ^  K1 {
大师 就是这个

sergio

这贴实在！

awer

maodeng 发表于 2013-1-28 13:13
9 h7 n3 i$ d& z8 R4 p0 l. hello 300236
1 A( l7 f. R1 u  K' astatement out of context
4 C0 A  d0 ?3 [! q+ B: i0 mr(119);

1 h4 M' S; x( i5 V) S1 v' b1 M6 c5 Z7 A抱歉，最近太忙了

9 ?  z; t2 X' d1 Y+ D300236 我用下面的程序试过，没问题

9 s4 k4 {1 y: b7 L2 i/ ]& b5 v5 E你把 D:\stock\test\source 改成你导出的数据存放的目录，再试一试
3 M: E- ~+ W) p" T% o) z$ i* B
3 M/ Q8 D% p" w8 }8 u! J3 ]+ \capture program drop hello

program hello
- Z6 h# W2 D- ]) r5 }; q
quietly {
0 R8 p* ^6 N& O  v' e/ ^
* \# ~/ g9 A( g- D7 ?+ Tclear

cd D:\stock\test\source

% D1 L/ }; e' X& \( {0 [/ d" {insheet using `1'.txt
  @& X3 @) n) K% C+ E. u2 \/ c
# E9 v" B0 ^& A; A! @   keep if v5>0
   
   gen vt=_n
   
5 i3 S4 c8 n0 L9 S   tsset vt

. x% E! r+ t' w/ o   gen zf=D1.v5

   gsort -v1
1 X6 A- q% x4 l3 ]0 m/ ^3 D5 s
   gen avol=v7/v5
   
   keep if avol>0
2 ?) s' L: x: W! G   
   drop vt
& [+ ]5 |/ c( [  
   gen vt=_n
   
* w% q/ B, ?. A% @; {' e   tsset vt
   
$ t  B1 k0 R" p' P! `# b    sum v5 [aw =avol]

, b: L1 a2 b+ }1 r' b& v/ h) f$ M3 U   gen amean=r(mean) in 1
1 a1 Z$ Z! R. R; q   
; ]6 |* V+ u+ b: e2 o+ G0 A   * 总加权平均值 amean
6 q/ _7 r1 c. |$ s3 a6 k   
   
   keep in 1/240

   sum v5 [aw =avol]
8 Y0 {+ a1 v$ ~8 h$ f2 @/ |# @& O0 X  M
   gen mean240=r(mean) in 1
) ^0 s6 S1 m9 O% P) R   
( A' K, h' ?# P7 O# ]7 F% z   *240加权平均值 mean240
1 s* N# ]) s+ ~. [% _8 G0 ]   
   set obs 500
9 a2 u2 Z. K( j) m6 |% w3 T3 V2 v   
1 T) }2 H1 z9 J8 \   kdensity v5 [aw =avol],gaussian n(500) generate(dx240 dy240) nograph
   
. e- p" P+ A- d( O: j) x) U   *240 kdensity(dx240 dy240)
   
) y0 T( ~. C* B7 W" Q' u& T   sum dy240
7 j! L: I0 e# Z' G  t0 S   
   gen yxis=r(max) in 1
* @; `6 o" ^4 q/ f* ]! `' Q   
   gen yxis2=yxis/2
   
   *yxis
   
2 I6 Y+ \* ^# E) K, C& o0 C6 E   gen dy240max=r(max)
& m9 @# y7 K: I   
5 Q4 H- l6 b: j( m0 n   gen dy240minp=r(min)

- D4 r% j: w: }4 j( v$ _* U    gen vtt=_n
4 S+ f  T) X/ b  r
    tsset vtt
  |9 Y7 }+ x1 z; t
' c2 l" k9 v! e5 s* k- w7 e% J   gen ddy240=D1.dy240/D1.dx240
# Y+ s7 w+ h, J' F   
   sum ddy240
   
/ i+ e- Y( G. |* W% c   gen ddy240max=r(max)
* y4 J3 c/ T3 {   
: e0 ]! t$ g7 a# ]   gen ddy240min=r(min)
* q) `7 I6 J# P* T0 C     
   gen ddy240g=ddy240*(dy240max-dy240min)/(ddy240max-ddy240min)
) T2 s* D& ~3 V9 a8 M   
: l9 o/ m( v8 h   sum ddy240g
   
   gen yxis3=r(min) in 1
3 z- C) v1 U1 X9 _# y2 ?5 d   
* j1 j7 n7 b( t' A  q   gen yxis4=0 in 1
  
  *ddy240g
   
   
  preserve
  
  sum v5 [aw =avol]
& z; c" j8 K9 i1 j$ U* e; I  
% f( Y5 \- J1 o  keep if dx240<r(mean)
  
  sum ddy240
. k5 o* @: W: Y$ ^  
2 {4 m& W4 L) f8 j3 b  if r(min)<0 {
        
6 i, Q- x# t5 i+ P  T7 s                keep if ddy240<=0
" x) X( _9 _* k) k- i6 ?( O0 I               
                count
: O( \# u3 j( X: C               
7 L$ ]3 b0 ^( A5 Y                if r(N)>2 {
8 F' e8 l3 C2 |/ J: e& y7 [6 j- f                 
3 ?! V9 X; K2 ?. k                 pkexamine dx240 dy240
: I: b& y2 S6 T) J  e) r; h- i               
; i$ T" y5 w3 l' `                local tempminp=r(tomc)
                restore
                 }
' m; }' |5 U2 q, x                 
6 \8 F2 {$ ~7 f1 s" u$ E# G( T8 D* }        else {
5 L& \6 _% Q( `0 t               
                 di "cannot find minp"
1 o" w. g6 h" \! R                 
: `% r' Z5 D4 U- z9 c' _                restore
                 sum dx240
% e6 d7 p* P) z/ I" y4 G0 C& n                 local tempminp=r(mean)
                 
) V7 Y! G/ `  X- C* E* L+ e                 }
                }
   else {
                 
                 keep if ddy240<=r(min)
                 
                 sum dx240
                 
' I* i/ r& E* h4 V                 local tempminp=r(min)
( e; ^' q1 |: k# m& g2 G; a                 restore
                 }


; [5 u1 l/ h0 y2 U1 M  S. c
/ j4 K' R4 Q$ v/ ^   gen minp=`tempminp' in 1
" N# m+ k: }* D5 ~   *240 minp
+ k- I9 R/ O4 s, ?
- j7 _$ ~! d8 C0 W, V0 K# R0 c
7 ]/ q& j3 Q2 x   
; J; ]% c. n- P: m/ w; R   preserve
   
   sum v5 [aw =avol]
   
   keep if dx240>r(mean)
3 l# \* Z" G1 T. R   
3 I9 _' ^: I! e, m  k+ v   sum ddy240
   
$ V8 y7 Y% g7 U5 m; l9 i, E/ W9 R   if r(max)>0 {
/ N! ^+ }/ F* `* B6 }               
7 e* ]' g, M0 \3 c4 {                keep if ddy240>=0
               
                count
" c+ V4 ~! L" s. }9 u4 d3 O" B               
                if r(N)>2 {
                 
                 pkexamine dx240 dy240
               
) \5 j2 A( [+ @, `  a                local tempmaxp=r(tomc)
                restore
                 }
& J! f  `% I9 k9 O- d. {& `% R                 
                 
        else {
' H* m( |! Z; I+ |, q4 l                restore
% |, `2 H( m/ r) l- }                 di "cannot find maxp"
               
               
4 t' F& }0 {1 T& s) G2 U# X                 sum dx240
) Q9 O5 m$ B  M% Y  P/ c2 e  Q                 local tempmaxp=r(mean)
6 ?% _1 h3 m" K) T7 n                 
  l) B$ ]3 f. J  m/ C                 }
                 
: @; L6 ]4 f5 s3 f9 _                 }
& p7 `1 D8 }2 W3 i                 
    else {
       
           keep if ddy240>=r(max)
          
, A- d/ @  R# c3 Y           sum dx240
# `0 x9 h+ P( H6 t4 q' \5 `6 ^          
           local tempmaxp=r(max)
& o  y' G) h9 G% m3 R* u+ {           restore
           }
8 Y. i4 l. d- x( m8 ^# S' C/ {       
7 K1 M' \' y; v6 i+ A
4 ^! @$ d+ j3 z5 `$ I1 [/ ^  
0 M  C$ q4 y; i9 I
4 ]3 f/ ]4 ]( b8 u% U' [6 t   gen maxp=`tempmaxp' in 1
5 {) Q0 O1 H" O, j% Y   
7 p3 s6 \, K  _9 a   *240maxp
: v) O. Z+ \/ Z$ C6 \" `4 g3 v   
- K5 E5 m5 F( ^& r3 K   
- [0 J1 {6 k$ G, T8 k   *-------------ddy240pminp左侧高峰值-----
% U" A5 P8 Z* d1 c5 C: [   
+ V. @4 h- b$ ]    preserve
2 n+ A7 }9 l) l; `, ~$ `   
   keep if dx240<`tempminp'
( b# k) V/ R7 V# a0 e5 q9 F   
1 f' s" N8 L6 w3 m1 ~4 Z2 k" ]   keep if ddy240>0
7 C) x( M( t4 W+ ~6 L   
) q% S0 [7 K1 P6 ~- @8 ^   count
) I# Z" X$ o4 V5 r" T               
                if r(N)>2 {
/ _( W1 L# H$ I# k2 ^+ {                 
/ j6 a2 U6 }7 n$ U' R                pkexamine dx240 ddy240
               
$ H" B& T1 m2 v, C                local tempminphddy240=r(tomc)
                restore
& M0 P- ?8 a. e( D                 }
                 
        else {
                restore
                 di "cannot find pminphddy240"
8 _6 I5 x$ w" x8 L6 O1 X                 
5 ?3 I) F: m2 [% L  q* g               
6 t$ A" k9 ]( s                 sum dx240
) q6 m6 z# K( K' |1 K; i                 local tempminphddy240=r(mean)
+ t  a! O3 i& p7 @9 R2 a& b               
$ `: {. n3 S+ s* ?( \/ X3 R  m4 s                 
3 T/ b, V. b8 P                 }
+ ?3 G' G3 ^- o- P9 i+ t                 
   
   
" V) r' i) [1 ^% b; g* S5 {  {* J4 b   
& g  K0 H7 U3 ]. L   
   
        gen pminphddy240=`tempminphddy240' in 1
3 k. g( f: Y1 D       
       
        *-------------ddy240pminp右侧低峰值 -----------------
; W# H% O- f9 ^: O    preserve
( N4 p; j" Y0 S
- W2 J2 X/ j5 m  s0 D+ }   keep if dx240>`tempminp'
  h: E/ C# z4 ^3 M" e4 U   
   keep if dx240<`tempmaxp'
7 S  n5 }1 ^- m1 P2 }5 |   
   keep if ddy240<0
! q# h) g3 Q9 l: v; `4 @   
   gen temddy240=-ddy240
: f9 s' ~2 ~; J7 O! @   
) l! N' Y: p3 s: D6 [1 d' @& \% J    count
# O8 c7 f! j7 Q3 k& ]: H               
                if r(N)>2 {
9 l# p" S# x, R  ]- j# A8 z, R0 w                 
                 pkexamine dx240 temddy240
   
4 x5 w: [5 w& \, C/ ?. a  local tempminplddy240=r(tomc)
                restore
. A. g, ?  N4 m                 }
                 
        else {
               
                 di "cannot find pminplddy240"
                restore
: E! S' Y) e9 H' G                 
* A  T7 u5 c# f) e! m1 l" ~  x                  sum dx240
                 local tempminplddy240=r(mean)
                 
1 {: d6 D# b3 I# X. U( Q" |$ x                 
                 }
" V( d- R5 q" f5 J   
   
   
   
- b& ~0 Q' X' l" S, s7 i   
6 U: }4 X  |  f3 Y, ^# M        gen pminplddy240=`tempminplddy240' in 1
) S# S7 e( D3 B  p* D2 C9 `       
) l' _+ ~. F4 t1 _        *-------------ddy240pmaxp左侧高峰值 -------------------
5 C- G- L7 o; k3 P       
         preserve
8 b* t! c2 K& R- e% O" \  P   
7 d6 A- H+ b. l0 P. X& q( ~3 H   keep if dx240>`tempminp'
# d# S0 q: ]0 x0 K/ T. k1 J: }/ {( o   
   keep if dx240<`tempmaxp'
   
& v( X  F- }. j( b; |3 @   keep if ddy240>0
2 M1 O4 X8 X3 m% n6 Z   count
6 P4 o  I  S& \7 }5 W! @7 H               
; ]. S5 h  R" B                if r(N)>2 {
                 
                 pkexamine dx240 ddy240
   
   local tempmaxphddy240=r(tomc)
5 Z% w' Z1 y3 M8 {& I+ q   
6 n3 f7 R; f1 H: H+ d# t; F; y   restore
& R1 L1 U% j2 l* {' n3 N. V                 }
                 
+ D3 Z: i) r* }% l+ D9 w2 g        else {
               
                 di "cannot find pmaxphddy240"
7 b* p8 Q9 R* |5 e* b                 restore
, e; T  Q+ d) f$ ~: i               
9 i' v* m2 U$ ]' y. z" Q                 sum dx240
2 h3 c0 S- N  d* x2 h& Q; I+ R5 }                 local tempmaxphddy240=r(mean)
                 
                 
, \/ v6 m5 w8 }                 }
- V) Q8 F, w4 a9 s   
   
   
   
: [0 N; p6 n1 K6 ]) q: [9 o1 r   
   
        gen pmaxphddy240=`tempmaxphddy240' in 1
       
       
5 [9 d# E1 s+ D' w6 ]+ Z) t        *-------------ddy240pmaxp右侧低峰值 -------------------
8 {0 L' f* ]" h$ S6 W" h, a       
         preserve
$ |$ c4 L2 V+ J' `: |# ]1 B  
) z9 M" X% H( k1 q( p  z8 ]) Z' U: [   keep if dx240>`tempmaxp'
   
   gen temddy240=-ddy240
7 S% J, f2 x" s   
   keep if temddy240>0
   
   count
               
                if r(N)>2 {
                 
8 @! p- ^% E9 N5 v                pkexamine dx240 temddy240
   
   local tempmaxplddy240=r(tomc)
   restore
% {% i: a9 J( Q                 }
9 h' \, V5 F- n# Q! C4 E8 q                 
! \( `$ c  y" H" d" R! d        else {
5 C9 y0 d+ g; R! x! s6 {                restore
: B& w4 F2 `; f                 di "cannot find pmaxlhddy240"
4 A! \' L6 Q$ Y& U1 b                 
               
                           sum dx240
                 local tempmaxplddy240=r(mean) in 1
               
3 o4 {5 H  d( u" \, e# L                 
) W; q/ _- _2 k' X' D5 {$ `                 }
   
' Z+ u2 }+ V& M: i: @   
4 w8 {0 c7 t9 S   
6 L$ S- D# B6 D1 Q8 V: P0 f- ~   
5 y( ^: N, n! m0 @+ o' W   
+ q/ p  D6 C* t% S   
        gen pmaxplddy240=`tempmaxphddy240' in 1
' n: a/ `4 b% S! L8 A4 r$ i- F" e       
3 Y5 C& `7 b! t! e$ e) \: t        *-------------------------------------
   
( M, [: n* N! M7 o( h/ r   gen price5d=v5 in 1/20
      
  A4 X+ r& I; p   gen price5yy=_n in 1/20
   
   gen price5y=price5yy*(dy240max-dy240min)/20
   
   *price5y &price5d
   
   gen priceny=price5y in 1
2 o- A8 K0 P. c4 o- Z   
6 T4 m+ c' q7 H  H   gen pricen=v5 in 1
6 y1 e) b( h3 \' [5 e
     
   *priceny &pricen
   
4 s! W- ?" L: M   
   *-------------------------graph------------------
  
0 U+ H0 h2 Q$ k6 u* |/ n/ I' e! N  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)
6 ~& D) G  `  o! M5 O3 L) _
* ---------可选项目----------------------
*(spike yxis amean,lc(blue))(scatter yxis2 amean, lc(blue) msymbol(+) mlabel(amean))
0 _1 d% _2 O" Y. a$ K, x6 E2 A& d}
" O/ {, _; `4 }3 T0 gend
# B# H1 j- z; I; s6 p. L' _; r0 y

maodeng

awer 发表于 2013-2-5 14:40
0 i, _# t' K3 F" C* q2 A抱歉，最近太忙了
2 X% O7 _; G4 ]8 m9 \
, |. M. t, s! L1 |# p300236 我用下面的程序试过，没问题

6 h* h! D6 q! p哇，谢谢，我仔细看看，这些字母数字代码晕了

喜欢喝冰茶

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

喜欢喝冰茶

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

( B$ A3 o9 j" l8 y
2 N/ [3 j7 h$ S$ i# l0 f这东西是开源的？知道有个统计软件叫stata，不过要收钱。

汉水云潇

准备试一试
8 \9 Q8 Y* s7 ~& R% S& |先找Stata软件

汉水云潇

( F& D+ X# V3 J$ L/ l晕呀
论坛复制程序代码
乱码夹杂在其中
! B- {+ j1 r  F7 E1 a7 N1 b: L很烦
对了两次
$ z4 \+ O5 P! ~  B3 f1 x* g以为没问题了
( A9 x4 s4 ~8 ^( P
! Y1 _/ r& z1 L* G) K结果
执行后
% E/ j* N  Q# p2 ^. r2 S% M# c( w还是出现
r 199
# F1 |+ [/ e, F4 b的错误
4 [+ g  J6 U, s) u$ c' e

汉水云潇

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

peterbev

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

中秋不圆

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

bknite

这个帖子值得学习。

jeanny

感谢！！！

Hubert2000

明月回春 发表于 2011-8-13 20:17
程序写的不错.
/ M8 _$ v1 ^* {& B4 p' j# d不过有一点要注意.通达信导出的数据,只对价格复权,并没有对成交量复权.所以数据需要手工处理 ...

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