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

maodeng

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

maodeng

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

maodeng

. hello 300236
1 k% D1 K* \& j' ^8 ?1 u! p9 S& wstatement out of context
# `; \7 r1 D* u" H) E; s( er(119);
7 k9 @3 T: W; T5 {
大师 就是这个

sergio

这贴实在！

awer

maodeng 发表于 2013-1-28 13:13
. hello 300236
0 g7 x7 c2 |% f2 E# o& fstatement out of context
/ L0 T+ s/ u, @2 V9 B5 Z3 |3 Nr(119);

抱歉，最近太忙了
) B) S% P7 T+ r- x# ]
3 B( ]3 p% z  m4 v% z" D300236 我用下面的程序试过，没问题
# ?* v  O' h/ d1 S8 m, k! z
  M- y0 c! ]( z: }3 g# M# c你把 D:\stock\test\source 改成你导出的数据存放的目录，再试一试

capture program drop hello

7 K: ]. X5 l$ f* Z" ?' tprogram hello

: u+ Y7 @8 o: [4 ~& I7 S quietly {

( G* F3 y" s# i/ X3 q% G+ S& |clear
( c/ d& p% M+ I! p/ I
" r- O! c) U; S. n' W7 D; l# scd D:\stock\test\source

" E: q+ F& L2 E+ v0 Winsheet using `1'.txt
! T' f4 _; @. z, d0 V( T, y
% q9 `% g6 |1 S) z( x   keep if v5>0
) ?2 h5 M7 u) e) `. s4 n   
0 ?0 t1 m: A7 s" P. a; }! ?   gen vt=_n
! E/ G# F3 s3 I( I  [   
   tsset vt

, i* }) f) m, b' w: S" A   gen zf=D1.v5
+ z- m3 h' r+ }4 l
5 H8 ~3 S6 l5 X. H0 b   gsort -v1

   gen avol=v7/v5
   
4 [# ]4 M) D  g3 {   keep if avol>0
   
* F# R& V6 {6 U; {- z   drop vt
5 q7 P/ |( J: x# x$ t+ A  
   gen vt=_n
   
( z, d) R0 t; B( s7 L, ^- N4 f   tsset vt
7 n9 t4 l) i+ t   
    sum v5 [aw =avol]

- ]" `5 ~/ Q5 b   gen amean=r(mean) in 1
   
, w: g" D4 L2 ]. y8 d; |   * 总加权平均值 amean
   
   
   keep in 1/240
: Q* {/ |0 q6 M) S$ n9 e7 l
: r9 ^8 \1 E6 L- V9 @   sum v5 [aw =avol]

   gen mean240=r(mean) in 1
! G5 s& p% j( Z  c* v   
: g2 t/ M# {0 [/ W+ y: i" I   *240加权平均值 mean240
   
7 [* i: ?. b+ ?, Z# _9 a   set obs 500
   
3 h4 X7 W; v# l% O" _3 k   kdensity v5 [aw =avol],gaussian n(500) generate(dx240 dy240) nograph
) q6 `; J8 ], m   
% L5 c! s6 S+ C8 e) p   *240 kdensity(dx240 dy240)
   
# A  |/ d% A3 K" a   sum dy240
0 ~& f- v+ E2 A( w0 x& [   
8 h% v( `! B- k7 B+ a   gen yxis=r(max) in 1
% T3 d' v" X, @3 m. A) \2 @   
   gen yxis2=yxis/2
! M' s$ P. ~. Q$ N; u8 f+ s! x   
! x) s8 b. {+ j0 A4 V5 R   *yxis
   
   gen dy240max=r(max)
   
& d& o: r' r9 F/ w( h+ j, J3 [   gen dy240minp=r(min)

    gen vtt=_n
6 I2 c* ]9 v- j2 @! `
1 ~, [4 X& j; t1 U6 u# O' V    tsset vtt

6 O! h/ M% w1 }* ]  C, v- J5 H% J# ]   gen ddy240=D1.dy240/D1.dx240
   
   sum ddy240
* z: M+ P# f" m' v+ o" r   
3 [2 ~- Y7 J8 @$ t2 G6 g   gen ddy240max=r(max)
   
   gen ddy240min=r(min)
     
2 I9 b0 B4 t# u/ Z+ p( \   gen ddy240g=ddy240*(dy240max-dy240min)/(ddy240max-ddy240min)
   
   sum ddy240g
   
   gen yxis3=r(min) in 1
$ B  u/ D" s( C# M   
   gen yxis4=0 in 1
  
  *ddy240g
. N9 j! S9 w7 F9 U   
   
  preserve
  
/ c5 ?4 l: |; y8 j  sum v5 [aw =avol]
  
  keep if dx240<r(mean)
. ^* A- i' p1 A4 d, T* x; n* q% V  
  X8 Y9 q" u; z$ g$ k2 v  sum ddy240
. U4 y; @, I' b$ K, e# E  
  if r(min)<0 {
( R* f) V/ E, F        
                keep if ddy240<=0
               
) V# e3 L% n2 g; }                count
2 r9 A9 w7 ^) {               
: q5 k* A; @$ [: v1 q                if r(N)>2 {
6 N2 w$ L/ T  r5 u6 P  V' p3 O                 
                 pkexamine dx240 dy240
# H9 R, M3 C3 x' X* }               
" `6 b9 E% p, r                local tempminp=r(tomc)
                restore
                 }
$ n2 w: h2 Q2 Y# e3 p5 k; }7 P                 
        else {
2 |. ~6 i5 `8 o1 k2 Q               
                 di "cannot find minp"
* i: f$ H8 I8 U/ F4 V: P& \# W1 T1 [* V                 
4 I' N, J+ n- E! c                restore
& }4 `# o  o0 Y  A; W0 M, v) o$ v! Z6 G                 sum dx240
5 f' y  g4 h1 q! k" n                 local tempminp=r(mean)
! F0 N1 M$ N% ^7 R7 X: H                 
( y' t. s! E8 m2 p) B! ~                 }
: x" `, W2 M7 v) N8 d5 K  W                }
5 H5 y0 L8 A. [0 \   else {
                 
' q) }4 R+ m) L9 |                 keep if ddy240<=r(min)
! `+ [0 i! f5 K. y) E. |                 
  N8 ]5 _8 P9 U                 sum dx240
                 
                 local tempminp=r(min)
                 restore
                 }
# @  W" Q- b0 ]% l# q9 B+ x
( X: m- ?1 ]% c9 V2 F+ q- x9 ^

/ H- ?" o$ q  H! z( A  e0 ]   gen minp=`tempminp' in 1
# Y) [$ |/ f9 M: M* s, [7 @   *240 minp

. t8 ^6 f+ E# m8 o% U
: ?# ]& L, Q* n3 M* R   
6 m! F$ i9 i/ _/ w9 ~   preserve
& N8 Y+ Z: f, [6 h+ c; g4 v   
   sum v5 [aw =avol]
! q, y; ?; S# r! A; K, S1 U   
' R$ P5 j- @5 o   keep if dx240>r(mean)
0 `1 a! s/ ~* m! F  K" Q0 L   
# a1 s; t$ b* J" j   sum ddy240
   
   if r(max)>0 {
5 o4 T  V7 ~* q6 Q: i: D% ^               
) d% r' l# Z& K( x                keep if ddy240>=0
               
. a2 K, j4 n+ G; ^: Q                count
               
                if r(N)>2 {
" {2 F+ d1 n* f: W' ~, \# A                 
                 pkexamine dx240 dy240
               
                local tempmaxp=r(tomc)
                restore
                 }
                 
( ^! y# o- H$ z& O' W4 {; r                 
) L- G# k+ q9 }8 |7 W# f% c        else {
                restore
' M- n" [) E6 z+ |* W                 di "cannot find maxp"
8 i$ R8 h! D5 U% t' X: N0 D7 L               
8 I) o- o9 ?2 c1 ^               
) D- u) I+ ^) R                 sum dx240
                 local tempmaxp=r(mean)
                 
                 }
$ m9 i5 e: X' h& u* G7 j! A0 N                 
                 }
) T3 P, O' _# E8 b. j                 
    else {
       
           keep if ddy240>=r(max)
! c! a0 [' T8 v6 Y" |$ ]* Z+ h- ?' t5 X          
           sum dx240
: P2 Z. r2 H9 P- u4 o' r          
0 _$ {% E. l6 d; m; m: Y( h! ^           local tempmaxp=r(max)
           restore
           }
3 a) }4 ?6 T* U: f: j: C       
1 j7 u) j. O2 _5 L# G1 k# I% F0 m
  

5 x) J4 [6 Z# n3 V* P; p   gen maxp=`tempmaxp' in 1
4 m. J6 y! t: B, ~1 X* b& _* n7 J8 g   
- Q* p: J0 X( E   *240maxp
2 a% T; N/ m2 M) X! z3 B0 I' [. Q   
   
7 H# i+ J0 D7 H/ t" w- r; w   *-------------ddy240pminp左侧高峰值-----
/ |; D5 [& {' a% Y$ U) B   
    preserve
! b3 k& |! ~9 ~- w5 m, u$ \7 {   
1 }: n5 R* |7 \4 T) ~1 J: @   keep if dx240<`tempminp'
   
   keep if ddy240>0
( P& h2 x( A$ S/ w. |   
7 e- Q. w. Z# L6 I   count
6 H% Y" A% M3 C9 r* L: `* r               
                if r(N)>2 {
; ]9 g# l8 Y3 O. P% B7 N; A7 q( t/ q                 
                pkexamine dx240 ddy240
5 M" C( Z+ X3 G% j0 l1 `* H               
  q5 O" n0 o7 [$ E                local tempminphddy240=r(tomc)
                restore
                 }
1 l# \- O8 `. G  A0 J2 z9 {                 
# D( X8 X4 e6 ~: _& d        else {
                restore
- I/ t1 R! ~# `1 g                 di "cannot find pminphddy240"
                 
$ |3 X- R5 ]$ x. O. x! V9 Z               
' f' B! [2 L* o4 \7 u                 sum dx240
8 }, a/ _7 y: F+ x: ]  |6 k                 local tempminphddy240=r(mean)
               
                 
/ r6 d& u7 |" K# n7 }                 }
* T7 ~) V% W3 l; ?                 
- P) N7 |" M) o2 Q& t3 b$ A& Y6 p   
5 l5 ]8 l2 A2 e+ x! l' C6 x   
   
7 ?9 \, r* e/ ^7 K9 m   
" K% y( \8 ?" z, @0 R' a/ W/ w   
. B# o* N' ~  j7 e/ v: Y        gen pminphddy240=`tempminphddy240' in 1
       
       
        *-------------ddy240pminp右侧低峰值 -----------------
    preserve

   keep if dx240>`tempminp'
; B; m: _0 I" t0 a   
+ E# V  ]! w* r) y$ e5 _$ [. ?   keep if dx240<`tempmaxp'
   
   keep if ddy240<0
   
- D" c4 J1 Q* v( ~3 ~; t   gen temddy240=-ddy240
" A8 {& }  G5 G3 \   
8 z% i( J+ @5 |( R. _8 g8 C    count
3 w3 B% |6 c+ Y1 @: L6 o3 w9 U               
$ d5 X) D0 ]/ y7 \2 o+ n                if r(N)>2 {
) c9 c0 C( M5 a5 ~                 
8 N. C$ Y2 O+ Z, X$ R. T. D+ t                 pkexamine dx240 temddy240
   
$ {8 r1 o8 d1 }, w* ]  local tempminplddy240=r(tomc)
8 Q( o9 e& U% ^4 n+ f# }                restore
                 }
2 c' j) D; U$ L( ^0 X                 
+ N! l5 z+ P" q        else {
+ m" S. x3 D* n5 S$ f7 ^               
. m- Y, p% a# J                 di "cannot find pminplddy240"
                restore
                 
, E; X' L1 W4 f# E                  sum dx240
9 B. Q: J* v, J1 |                 local tempminplddy240=r(mean)
                 
9 U( z+ y- o6 r$ U7 a3 A0 V) Q                 
                 }
+ K% \: y, T! s2 |" Q   
   
   
% H+ V  R* Q$ f; |   
   
        gen pminplddy240=`tempminplddy240' in 1
# U3 }( l) T* v3 s7 q       
        *-------------ddy240pmaxp左侧高峰值 -------------------
% {  s/ y# F* T) e: ^7 p       
2 X9 n2 S- _# |; m  b' p         preserve
   
   keep if dx240>`tempminp'
0 }+ q0 s- ]1 D  h; }4 u6 l7 [   
& M( f  v! B* z   keep if dx240<`tempmaxp'
( T2 R, S  \4 c( g( N+ F& n- e4 c7 ~   
   keep if ddy240>0
   count
               
                if r(N)>2 {
                 
; p) E1 Z; f% N; `- A# Q5 q+ M                 pkexamine dx240 ddy240
   
+ {7 o2 H* J8 Z& G   local tempmaxphddy240=r(tomc)
   
   restore
                 }
5 o* D/ y4 o7 J2 s7 ]( n2 U8 A3 p                 
        else {
' G) a  f- l+ i               
1 ^4 i& g3 M2 B; r1 d                 di "cannot find pmaxphddy240"
                 restore
  ^* M/ V! E' l, k7 B$ T               
# f2 p" F2 u# K. p# S( N' j                 sum dx240
, Z& _% L* Q5 K# ^! t! c5 N4 P9 m% v                 local tempmaxphddy240=r(mean)
' Q- z& q. i% g8 ^3 H  I# G6 F                 
* b. Z- h/ b; o                 
6 M& l  Y$ w# c* S                 }
   
   
, Y) [5 p$ b$ K) o) D, @! q   
  a6 K( A1 }  R7 V$ {+ B   
   
3 [/ m' O( Z$ c. \/ r: k( w5 i8 l   
        gen pmaxphddy240=`tempmaxphddy240' in 1
       
       
        *-------------ddy240pmaxp右侧低峰值 -------------------
       
3 q: f$ g$ \7 L5 o9 `) Z8 x# ?         preserve
# L& ^# X7 T! ]1 s0 P' i  
   keep if dx240>`tempmaxp'
   
& \8 U, A% y. D' Q5 a   gen temddy240=-ddy240
/ B5 R! x4 }+ ?6 N9 T# ?$ j   
   keep if temddy240>0
   
   count
. t1 h3 n3 U- W, I, J               
                if r(N)>2 {
" i' a( o+ X7 z0 |                 
2 q( G" Q. i! H: W                pkexamine dx240 temddy240
# h+ P1 ~) o/ g/ Q  E   
   local tempmaxplddy240=r(tomc)
   restore
                 }
                 
6 U9 E3 U  Y4 a1 J& m        else {
8 y( U  H4 b! D                restore
5 X  S& y+ ?$ _; u9 A                 di "cannot find pmaxlhddy240"
. l8 x& p3 v# B) E' Y1 X! e$ q                 
               
                           sum dx240
0 v) u  `, I( e4 |                 local tempmaxplddy240=r(mean) in 1
( z$ o2 n6 v# F5 j5 O6 h7 D' J               
6 f" E3 e6 w6 A' U$ p                 
3 M* t0 F! I) D: i) V+ i                 }
8 C" u7 c# u1 J1 d  |0 X1 ?8 W# G   
. Q% F$ _7 |. ^& \   
1 l1 h1 m$ U* r! ^# c   
4 t2 F8 {& \' T5 Y4 h1 S/ w! |   
   
9 b* z6 ^# A0 M& F4 C   
! \" C; V& ^1 ]( |. ^9 C        gen pmaxplddy240=`tempmaxphddy240' in 1
# m# G* z: F: T. l       
        *-------------------------------------
2 w  p5 k2 B5 T2 F   
9 n4 T  D: l+ `3 p7 G   gen price5d=v5 in 1/20
      
   gen price5yy=_n in 1/20
" Z# z7 t' ~3 E2 }- ]8 `   
   gen price5y=price5yy*(dy240max-dy240min)/20
   
" y" H, X! ~2 L8 g   *price5y &price5d
7 ]; Q- H* Q% X: g$ u' r- _   
   gen priceny=price5y in 1
5 G8 R9 J8 c* \, a8 I, `   
9 y; M- j8 R7 S3 E   gen pricen=v5 in 1

     
   *priceny &pricen
   
( \* h' P3 ?# e- S- V+ F! S   
   *-------------------------graph------------------
7 I! ^9 r9 \2 V, `2 \. u  
0 ], \& l- M* T" r$ ~8 [  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)

7 z! U1 ?4 B3 G8 j% H! {% j( W' B* ---------可选项目----------------------
*(spike yxis amean,lc(blue))(scatter yxis2 amean, lc(blue) msymbol(+) mlabel(amean))
# }- e3 U$ `3 s  y) W}
end

maodeng

awer 发表于 2013-2-5 14:40
$ y* h* w9 N# Y* o+ ]$ ]) {8 D抱歉，最近太忙了

. \& W4 v$ A0 P- M3 z8 z  X) t300236 我用下面的程序试过，没问题

哇，谢谢，我仔细看看，这些字母数字代码晕了
" y* ]5 l' @1 @1 e7 q& K

喜欢喝冰茶

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

喜欢喝冰茶

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

4 D* _: c0 w% e, K# F9 v$ |: f
  w" t8 l! D6 E* p' h这东西是开源的？知道有个统计软件叫stata，不过要收钱。

汉水云潇

准备试一试
先找Stata软件

汉水云潇

5 r. P9 I- \+ B8 E7 C) b$ d晕呀
论坛复制程序代码
6 a  C9 M' b9 T: T乱码夹杂在其中
很烦
对了两次
4 V/ ?2 V% s/ q' [% V3 F! q以为没问题了
- R8 A2 E- G0 a# P
) I8 c( a; e) h0 m8 ~结果
' L* \" M+ }: T& I9 f执行后
还是出现
r 199
* k+ k9 U3 l/ b; X的错误

汉水云潇

试验成功，下面是三一重工的图
4 @  Z+ b2 A  L

peterbev

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

中秋不圆

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

bknite

这个帖子值得学习。

jeanny

感谢！！！

Hubert2000

明月回春 发表于 2011-8-13 20:17
程序写的不错.
$ C. V7 Y! t2 H% B( c, A5 N4 z* ]不过有一点要注意.通达信导出的数据,只对价格复权,并没有对成交量复权.所以数据需要手工处理 ...

8 L: d6 L# f) A- q- d老哥，也是山寨这篇文章找不到了，不知为何？