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

maodeng

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

maodeng

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

maodeng

. hello 300236
statement out of context
" U- p. f2 h( Dr(119);

大师 就是这个

sergio

这贴实在！

awer

maodeng 发表于 2013-1-28 13:13
. E- I9 p: v( Z' B' O. y) C& b7 a! A. hello 300236
statement out of context
- T! S5 h2 v0 n2 f' j. p7 G# xr(119);

抱歉，最近太忙了
6 Z0 r( {4 ~1 u+ S8 ~5 P
300236 我用下面的程序试过，没问题

你把 D:\stock\test\source 改成你导出的数据存放的目录，再试一试
9 k' y2 _' J* M6 o9 }
capture program drop hello
- E, G) V" B$ e
program hello
! q" V. E+ f2 _
9 e9 @+ I: @* X9 Y quietly {

clear

9 A; l+ B/ H5 C3 G6 M0 Rcd D:\stock\test\source

insheet using `1'.txt
0 N  j( i, p5 K' _
   keep if v5>0
   
   gen vt=_n
   
4 J2 M& n) E1 _% t$ A& ^" j   tsset vt

   gen zf=D1.v5

   gsort -v1

$ F" y5 ^$ U, N; y# ?   gen avol=v7/v5
+ L, ]: L1 M) w1 _. E   
   keep if avol>0
6 K; A" f) r/ s" e  X. r   
; P6 K9 G7 ]2 |$ V8 F   drop vt
1 Y9 h% y7 ~6 c+ h' O9 w  
   gen vt=_n
, O* x" O* N- d; a( R2 {. b   
: K2 m. p! i5 H( U4 s* r- N$ K   tsset vt
   
& {5 T' d, S2 w5 y5 W3 f    sum v5 [aw =avol]

   gen amean=r(mean) in 1
/ N8 |9 W+ N( }/ z% s5 W8 K   
   * 总加权平均值 amean
   
# v6 n3 x4 }- i- g7 j6 X   
# T- ?! Q% A$ n, K6 v- f   keep in 1/240

   sum v5 [aw =avol]

   gen mean240=r(mean) in 1
   
! @1 N# F" C. ~( v   *240加权平均值 mean240
$ z- N' d5 Q( j   
0 S) e. x1 T* K4 [  [& G. g   set obs 500
   
9 c* f$ g" c0 _! t0 m   kdensity v5 [aw =avol],gaussian n(500) generate(dx240 dy240) nograph
   
   *240 kdensity(dx240 dy240)
   
   sum dy240
   
   gen yxis=r(max) in 1
2 F% N+ J1 `& X: o1 F$ ~) A* N   
   gen yxis2=yxis/2
   
   *yxis
$ [$ U6 D1 R0 x  d1 M+ V- o# i   
   gen dy240max=r(max)
   
, Z: Z( h/ `- R! y7 l   gen dy240minp=r(min)

    gen vtt=_n
: r6 I( ]9 |- L* }$ _5 m
. ]. X  P6 A0 j" Z1 G* B4 }6 ~    tsset vtt

   gen ddy240=D1.dy240/D1.dx240
   
" \8 y: S" E' H- G7 P6 i   sum ddy240
4 i+ Z2 \6 q# t+ j' }   
9 Y0 l8 ^! m2 u% j% q& S6 H   gen ddy240max=r(max)
2 f7 |7 C. [* A- c* j( ~$ P   
$ B/ G7 b1 w* W   gen ddy240min=r(min)
! a8 E5 h5 g/ m; x1 e+ F     
   gen ddy240g=ddy240*(dy240max-dy240min)/(ddy240max-ddy240min)
   
$ t$ T* m4 |. a/ V' p- A   sum ddy240g
   
   gen yxis3=r(min) in 1
5 u3 X% G: @& l& j$ W) _   
6 r: U9 \! y: ]6 t2 t   gen yxis4=0 in 1
  
  *ddy240g
6 L/ }0 o) Z6 Z: R* o" v' t   
1 h2 c+ G, U) Z( q   
  preserve
, Z: x& V  K& \3 E3 B2 r; b2 l  
  sum v5 [aw =avol]
& r, A* Q5 Z, O4 t/ T  
) [9 i! w/ X* w4 b" v! W% z  keep if dx240<r(mean)
7 G* q. z, B; |) D( P4 S  
/ `- o0 V' e' U% B8 [3 @! V  sum ddy240
  
, J2 D7 `1 T, T# i1 ^  if r(min)<0 {
" ?  u% {7 M! G$ C4 m        
                keep if ddy240<=0
               
+ L+ k9 v! R3 z0 m  m" Y% \                count
               
                if r(N)>2 {
3 v1 l$ B( I5 m                 
                 pkexamine dx240 dy240
               
$ Z9 K! M, W. Z                local tempminp=r(tomc)
                restore
                 }
. e! }8 c3 o% V# [  J                 
        else {
               
                 di "cannot find minp"
                 
                restore
+ E- U) K0 o: P  c                 sum dx240
                 local tempminp=r(mean)
* O! X8 |8 B; [3 `2 m5 u, ?2 C- f                 
% ]6 y, v1 P5 A% M' Z/ B$ b1 t3 Z                 }
) _$ b: w+ _6 C0 v9 S$ H) I                }
( E9 q; ]# _& E! y, H  r+ t   else {
# G5 ~5 G( w0 `( H. A# T                 
                 keep if ddy240<=r(min)
. D$ |" f0 P) A2 W7 n! A; {& z% j- W( e                 
                 sum dx240
                 
                 local tempminp=r(min)
6 f0 P9 }: W! A2 P0 n- ]5 m                 restore
                 }



& S1 T: D& B, k1 F: u   gen minp=`tempminp' in 1
   *240 minp

% M4 y& L7 f1 S  O; R5 F
   
   preserve
   
   sum v5 [aw =avol]
- G0 }$ k2 ~' P! F   
4 H. r& |& Q  a' _$ D0 h* K   keep if dx240>r(mean)
   
   sum ddy240
- }; y. E3 q9 B1 F   
   if r(max)>0 {
               
                keep if ddy240>=0
               
! |/ a/ F# {2 X9 D" ~4 y: d& u: c; [7 H                count
) @, g; I$ g# |* b. {2 k* G$ R               
/ D9 F# x! B! ~9 Z7 l$ z; ?                if r(N)>2 {
/ T* R5 |, G" i* ^! u                 
0 L8 `" Q, Q+ Y+ O                 pkexamine dx240 dy240
               
                local tempmaxp=r(tomc)
                restore
                 }
                 
                 
        else {
0 G8 P6 e/ N% w! _8 x- u5 G, F                restore
                 di "cannot find maxp"
               
               
                 sum dx240
                 local tempmaxp=r(mean)
+ W" t/ X1 }, D6 Q/ M9 _( t$ |/ a                 
) _8 N: Q* s0 R# r# B) p6 M9 R$ y* V                 }
                 
/ K9 O7 m4 T: ^1 z                 }
6 k- L) Y( U# ?7 `7 q" V/ N( J% u                 
7 A' D* s, u# w% B# S5 N* c    else {
       
  a2 J# W6 O/ c  g/ c           keep if ddy240>=r(max)
  `5 f/ v$ V2 {# a% `5 p! D          
* R* B' n3 q- e           sum dx240
- F& h. @! A' `, T          
% J3 s6 j( n* g4 `# J3 t           local tempmaxp=r(max)
6 f3 K6 B- Q$ a$ T           restore
           }
2 e' T" m. t0 e; Q7 }       
2 e  u' r( ~5 j6 {! _0 J
# D; a, T- N! o% N( T, x* o  
- Z2 G4 P8 R$ ]
- r8 S8 T# i; [; N* E" S" D   gen maxp=`tempmaxp' in 1
7 ^$ l3 x8 m. O/ [. U   
. |9 B6 \$ M0 U7 m   *240maxp
   
   
   *-------------ddy240pminp左侧高峰值-----
' S' Y9 ]1 U) I- |3 R! d; h   
% K1 d9 Q& l4 E  W' a8 C" \    preserve
# P6 u3 V+ R4 S   
+ h+ o* e+ D! F) W- ?   keep if dx240<`tempminp'
$ m4 e: f3 D9 l# @( o   
4 \; L+ m+ _; T/ j: A/ T   keep if ddy240>0
   
: L8 N# y7 ?$ j" I! G% B. ^   count
! ]/ X( Q/ S( H/ q+ e1 {               
7 z+ S' d% `/ I7 T  u6 t                if r(N)>2 {
% k  _& L; K$ {                 
1 c! ?7 o! r7 v( G$ }- o                pkexamine dx240 ddy240
               
                local tempminphddy240=r(tomc)
$ ?3 X3 i3 q7 M7 b                restore
                 }
8 S% H0 D8 o+ i* h3 J9 \                 
        else {
                restore
                 di "cannot find pminphddy240"
                 
               
! p4 l: a. e9 j7 B                 sum dx240
$ E* r9 @9 R  {) C$ Q: B                 local tempminphddy240=r(mean)
6 }1 O0 a% ?: W9 c2 @# [               
$ u  c: K7 [# F) ~/ N, x- q+ r/ k2 b                 
9 w8 p2 m0 E0 y                 }
; L' H! h4 k8 s                 
& w7 E6 w% V9 M: Z   
6 t9 D4 s; E* w) C- w, a   
* @4 D- W" n2 N+ j* B+ I* z   
* B" }* J) J- n* g1 @   
   
  d3 N. h; q2 B        gen pminphddy240=`tempminphddy240' in 1
       
       
9 ]! F: f2 ~" W# r  b        *-------------ddy240pminp右侧低峰值 -----------------
' v5 R0 ]* d. O1 R+ K7 r4 ^    preserve

- f% m+ y5 h7 _! h% B; R   keep if dx240>`tempminp'
$ S& I1 b5 V9 v- w5 R   
1 h. p; D  S# n- [% `# A' ~   keep if dx240<`tempmaxp'
$ g, r0 a8 t9 W7 q, x4 U( S8 r   
   keep if ddy240<0
   
6 j! |( ~: @3 n& i; u% q& G   gen temddy240=-ddy240
   
3 j  ]. r+ N( w5 a    count
1 I( U; @0 D6 R% C+ M               
7 W, W& ?, u) A3 B1 }% d( r                if r(N)>2 {
                 
0 Q% x  A5 B. S6 Z6 x- b                 pkexamine dx240 temddy240
. h9 ?+ o; a2 t1 V6 ~" n   
  local tempminplddy240=r(tomc)
                restore
                 }
. e- c4 `/ L1 i( d$ q  g; }+ C                 
        else {
6 D5 o# t+ ~( U! m% g, u               
7 P0 A: b$ [' H: c- U0 t                 di "cannot find pminplddy240"
/ M/ f+ k# l$ w                restore
                 
                  sum dx240
6 ~. l. P' l# _) C0 K8 _                 local tempminplddy240=r(mean)
# u2 q1 m: T! h* w2 `                 
                 
. o; D0 T4 z* B5 \  J* R                 }
2 n8 h4 s$ ^! ], A+ }8 @5 k8 W   
/ c# \! g% X" U5 a' e: A2 q   
' [% j7 P6 J5 t) ^% s   
4 |! ?) q  t5 @4 U, f* M- b, D   
8 B$ G0 M$ X- \/ s   
        gen pminplddy240=`tempminplddy240' in 1
1 J  H% \4 E3 k: ~/ g5 L+ F: u  P       
        *-------------ddy240pmaxp左侧高峰值 -------------------
( `$ o# w/ O% k+ j0 d; T3 x       
. u! U+ ]- A6 }  U- I4 _7 I         preserve
7 X- u5 A( l+ @) k' b& d   
/ D3 m0 J: G; d: z( K9 U4 _' K   keep if dx240>`tempminp'
/ P3 D5 d1 }* u( E. t   
   keep if dx240<`tempmaxp'
$ @+ g5 V; U* N; @: ]. W7 P   
0 u0 C& `/ s! h1 t" M   keep if ddy240>0
   count
' A: _+ x1 }1 V4 ]               
* L7 `/ n9 Q2 n7 F* X; p* `                if r(N)>2 {
                 
                 pkexamine dx240 ddy240
$ n/ H6 J- |7 \# k& l+ @( w5 ^   
; q  @% I, E2 ], N6 ~! [   local tempmaxphddy240=r(tomc)
# f8 p) c) U% L# q8 x+ F! |   
   restore
1 u* [7 i3 \* x9 N) F/ g                 }
3 R* F! V. p# O+ ^" ]                 
( e. h/ Q1 ^; y        else {
               
( t/ ^) L3 l" M* f; I                 di "cannot find pmaxphddy240"
: }, z9 S$ H7 F                 restore
% ~# p+ M% b2 S1 m4 s6 J4 X8 i               
                 sum dx240
& r3 z) l" `$ V: `0 O5 c; G                 local tempmaxphddy240=r(mean)
& g7 f$ @* T2 I0 q/ ]' u                 
" r/ i9 r6 G/ b4 s; o9 E% V  W                 
# S2 U, B. q% W& }; H: O                 }
3 t- Z* _3 ~" R: t   
   
   
   
   
3 i' Y. w* ~$ a! L8 m   
        gen pmaxphddy240=`tempmaxphddy240' in 1
4 ?7 q# |  H( p       
0 {: c! l5 C. X       
        *-------------ddy240pmaxp右侧低峰值 -------------------
- p7 i5 s' v1 n8 F+ o9 e       
         preserve
  
% [. M2 s# y+ K* s  U* X0 _   keep if dx240>`tempmaxp'
! d, @3 H9 U" x+ O1 F   
& c* a7 X3 n$ T% l. C' _/ z& }: M- [6 |   gen temddy240=-ddy240
: u% {/ S3 y# c' q2 c   
   keep if temddy240>0
   
% O+ _+ y; i8 ~' K- x   count
               
                if r(N)>2 {
* ~# [: d# S9 Z                 
                pkexamine dx240 temddy240
# [9 z0 z  e  t7 {0 B7 X   
   local tempmaxplddy240=r(tomc)
& \: {+ d! ]$ y( B: [+ D$ C  a   restore
                 }
2 `2 w& k8 C& C, I# \                 
# N' L/ S+ u' X8 M% N  c. l0 V        else {
                restore
3 ]% C' t2 p3 u                 di "cannot find pmaxlhddy240"
                 
1 G, f6 i: Y* X1 V/ l, M               
) m; @/ T. M" k0 E& a1 w- [                           sum dx240
6 r! Y; F' d1 ?. c                 local tempmaxplddy240=r(mean) in 1
               
" g% K8 _( ^  w: O4 [% ~! u+ X                 
                 }
   
* W8 [" x; L* ^# A! B" Z; B2 T   
   
& w2 ]6 L: K6 x+ l9 L   
   
   
        gen pmaxplddy240=`tempmaxphddy240' in 1
       
        *-------------------------------------
3 ]! G' R3 C1 y+ i   
   gen price5d=v5 in 1/20
1 F) {  H5 c: u( U  T, W      
   gen price5yy=_n in 1/20
   
" V( K! M, z1 C  |; D# C2 F   gen price5y=price5yy*(dy240max-dy240min)/20
( M' d4 K/ ^' i# D9 a( _   
   *price5y &price5d
- F/ H# J' D$ c$ A5 z5 q   
2 P, `- f# r* V   gen priceny=price5y in 1
   
   gen pricen=v5 in 1
) [6 t1 _4 w% J
! d  K8 T% O* A4 t0 [     
   *priceny &pricen
   
   
   *-------------------------graph------------------
4 K$ l1 W* g. }  
4 `" ]8 k* F2 z$ s% {" j4 w) l7 h/ ]; }  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)

4 h* r, l, V6 A* L* ---------可选项目----------------------
$ I% b6 o% I8 f*(spike yxis amean,lc(blue))(scatter yxis2 amean, lc(blue) msymbol(+) mlabel(amean))
}
  W8 e1 S3 v8 Eend

maodeng

awer 发表于 2013-2-5 14:40
* E0 E! W) D, b+ S. q' |抱歉，最近太忙了
& b+ {7 j* S" U) R! d- @
300236 我用下面的程序试过，没问题

9 ~/ T# S" X0 I' F" ?8 P哇，谢谢，我仔细看看，这些字母数字代码晕了

喜欢喝冰茶

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

喜欢喝冰茶

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

& M- N$ b2 g: W) T. ~# F  m
$ B3 W( z; i* c6 R5 @6 W这东西是开源的？知道有个统计软件叫stata，不过要收钱。

汉水云潇

准备试一试
先找Stata软件

汉水云潇

晕呀
, k8 X$ u& H9 G. k% D4 ~2 R+ Z论坛复制程序代码
$ }7 o- P8 g. G8 q+ S- q  I乱码夹杂在其中
很烦
! @$ `; n: f: o2 X' @" o% A) B对了两次
- ?5 u7 z7 ]7 Z4 s9 D以为没问题了

结果
7 Z% }% R1 n3 q  N5 O" a执行后
# k  O9 o* A; U. n, Z2 B还是出现
r 199
/ u% Z- Y0 t3 S4 f6 v" J& t/ }# H5 ~的错误

汉水云潇

试验成功，下面是三一重工的图
# d- ]4 H3 Q4 Y0 g2 D

peterbev

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

中秋不圆

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

bknite

这个帖子值得学习。

jeanny

感谢！！！

Hubert2000

明月回春 发表于 2011-8-13 20:17
2 F4 e3 f5 J1 f: y- I0 q1 I程序写的不错.
; ?0 l7 k6 V" g0 c; C不过有一点要注意.通达信导出的数据,只对价格复权,并没有对成交量复权.所以数据需要手工处理 ...

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