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

maodeng

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

maodeng

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

maodeng

. hello 300236
* {; C' q7 J5 `' |6 sstatement out of context
r(119);

大师 就是这个

sergio

这贴实在！

awer

maodeng 发表于 2013-1-28 13:13
. hello 300236
6 e' N9 {! Q  M' P8 ?9 Hstatement out of context
6 B* r2 k8 `! U8 qr(119);

( y* H/ r$ c' u; V8 j抱歉，最近太忙了
) h$ o2 F% R) |; Q
/ \& S9 O7 j$ x9 j8 Y( W300236 我用下面的程序试过，没问题

8 b* C+ w# P& Z, `( G你把 D:\stock\test\source 改成你导出的数据存放的目录，再试一试
& A  Y; X/ O, n( ~9 a8 U& a$ Z8 a6 w
capture program drop hello

1 Y5 u6 C* d8 }. Wprogram hello

) F  _! ~" f; b8 r$ s& A' t$ I2 ~ quietly {

clear

cd D:\stock\test\source
0 J. z% c. A" m9 S7 @
: O9 z* L& h7 ?insheet using `1'.txt
- {0 I5 D/ w) y$ F3 J
3 k0 b( y8 L# k& I   keep if v5>0
0 n$ _% v6 q1 N8 c   
   gen vt=_n
   
   tsset vt

   gen zf=D1.v5
0 K0 s& n' ]2 P9 X) z% _  j" s
   gsort -v1
; A. \/ ^+ J. A9 x1 o
   gen avol=v7/v5
   
   keep if avol>0
   
   drop vt
6 [" Z9 u, l4 d5 r5 o3 \! r4 L  
   gen vt=_n
  f* B1 ~: L' ~5 I7 u1 |. s) X   
   tsset vt
   
    sum v5 [aw =avol]
: l& s5 r6 x1 `* X, P2 t
+ o7 w( {* O" y( O   gen amean=r(mean) in 1
$ o& a* Q  z5 Q, T  g, o; u6 e   
& j: }$ n, l4 F% f9 V   * 总加权平均值 amean
   
8 {9 Z+ G: k) a+ v* q1 `) W   
/ q) c9 D- |/ n- v  }+ f   keep in 1/240

   sum v5 [aw =avol]

: r( `: W& n; M0 [  a  P$ \   gen mean240=r(mean) in 1
( r/ y3 [+ K+ k5 R1 \6 Q/ ~5 {$ D9 q   
   *240加权平均值 mean240
/ h7 U+ ?8 q- W" ^0 q/ ]   
% q% o. Y  {$ X- X   set obs 500
3 P3 I" {$ h+ i& B- a7 S   
   kdensity v5 [aw =avol],gaussian n(500) generate(dx240 dy240) nograph
/ ^2 S, [) v, U: E# i; r   
   *240 kdensity(dx240 dy240)
   
   sum dy240
. D# {5 F' X- k0 b# r8 R   
   gen yxis=r(max) in 1
   
   gen yxis2=yxis/2
   
' P  X4 t$ v! [! h/ x   *yxis
# B3 R- Z. p0 h( ]$ P' o+ f: @   
   gen dy240max=r(max)
   
   gen dy240minp=r(min)
! w7 b, N8 ~7 M/ S3 U3 D6 A
    gen vtt=_n

6 ^- K3 t" v! `: {& n& K+ j    tsset vtt

$ a& T( e) S3 V( P/ K2 t/ K- J   gen ddy240=D1.dy240/D1.dx240
0 P7 s+ K5 K; p# C7 J' |   
   sum ddy240
/ ^4 ]; W& d1 p0 P* D( m$ r   
- I. t* t& e( R5 t( l$ @" z   gen ddy240max=r(max)
, E* d0 `7 r0 V- m7 j   
) b0 z, @' s8 ~  k- _6 C* r+ o( Q" X   gen ddy240min=r(min)
- ]; z; s& \5 T4 s0 z+ O) O     
3 n" W% m5 r- ~% W8 ?" C   gen ddy240g=ddy240*(dy240max-dy240min)/(ddy240max-ddy240min)
   
0 K; D0 @' U  ~; \4 O1 ?4 Y   sum ddy240g
   
+ m. j2 Q: ~5 S5 x& w3 ]8 P   gen yxis3=r(min) in 1
* F' G; r. f  H% O. B   
   gen yxis4=0 in 1
0 c8 J8 ~, E( B- G  ^  
  *ddy240g
   
) x3 \+ o; h% T- [& K   
+ p6 C5 I7 I. A# h' V  preserve
  
  sum v5 [aw =avol]
  
  keep if dx240<r(mean)
) l0 ]1 G2 w9 ^0 N6 k  
5 C# T# X* K1 w6 ~  sum ddy240
5 h2 k* c6 B1 }! @" L  
  if r(min)<0 {
        
) z# J  o$ M7 F% U% Z                keep if ddy240<=0
               
7 E1 j% T4 V1 ^( y9 r9 u1 |                count
& L2 h" L. G2 c; V2 R               
6 V) }$ p/ @+ T                if r(N)>2 {
                 
2 A: J2 n0 m! S3 f: V6 i: D                 pkexamine dx240 dy240
               
( B- q3 v- s6 M# _( v9 b                local tempminp=r(tomc)
: r8 U7 N: I; F/ x/ {                restore
                 }
                 
        else {
: a. \" u/ o% H* ]6 W. C6 B               
9 E( K  @- B' T' z                 di "cannot find minp"
) F7 V5 R) @4 t: \6 O                 
                restore
                 sum dx240
" E* Q5 g# l6 ~3 [                 local tempminp=r(mean)
                 
                 }
                }
   else {
                 
                 keep if ddy240<=r(min)
2 ^/ }/ @3 J$ ^                 
                 sum dx240
                 
* E' g6 W0 |7 d) v                 local tempminp=r(min)
* P3 _, ~  d" P                 restore
" l4 M; X1 x+ T( E                 }

0 p/ N" [4 S+ |4 Q

   gen minp=`tempminp' in 1
   *240 minp
/ F# d. i5 W3 Y

! [% P- n3 y9 r* o) T8 T6 K: R( G   
" k9 z- X2 J. W' r5 ^* Y   preserve
; c% C6 {& G) \% P   
   sum v5 [aw =avol]
; U" b8 e( R7 B' H) V  p( ]- q3 o   
   keep if dx240>r(mean)
6 r$ i6 u) M& u4 {1 D  A4 F$ ~2 U4 s   
   sum ddy240
   
% u4 ^4 Q; {) W5 Y. v   if r(max)>0 {
( I3 a" \" a. {               
                keep if ddy240>=0
2 d8 W5 c1 V: G3 C1 g               
8 ^' K) ~) G; U# t( C* K                count
, i2 f' c7 H& ^: M               
( P" b* x6 x1 K- e! f+ \                if r(N)>2 {
& r$ A, s9 ?' C2 {: m4 @! Q                 
                 pkexamine dx240 dy240
+ n6 f9 n  _4 X* [               
                local tempmaxp=r(tomc)
                restore
                 }
                 
" z+ d- s3 S, m! Z! z, A  _                 
6 v9 m; x1 E, M. O1 @" f' N* Y& B7 @2 R        else {
- L  v3 G* y* ^                restore
$ D7 p3 l2 E3 G                 di "cannot find maxp"
               
+ b! A6 W+ i  s+ m/ U% B0 C  ?8 k. o               
                 sum dx240
, A/ S, p& z* }% g9 T) ]6 j                 local tempmaxp=r(mean)
                 
                 }
                 
! ]( I7 V  B9 q0 c                 }
                 
/ G  x- L; G3 b: ?" `% a    else {
       
5 m' J, k9 b# M           keep if ddy240>=r(max)
, {) B+ i1 Q% F& P          
           sum dx240
. R( s1 W8 m) {  m: C& b% E4 N          
! @6 M6 ]8 ]- R5 Q8 c           local tempmaxp=r(max)
% r9 P1 Q: v6 q7 o" P5 F           restore
! K8 _4 r, x/ U           }
6 J4 }5 B  w9 l7 t! {       

  

   gen maxp=`tempmaxp' in 1
   
* _2 m" A; }5 q+ ?   *240maxp
   
   
( r% l0 A2 o1 T. s  d   *-------------ddy240pminp左侧高峰值-----
   
    preserve
   
, S) y& n* f$ ]3 u5 e  M   keep if dx240<`tempminp'
   
   keep if ddy240>0
( M( e8 j* x7 C! y$ o+ D3 f   
1 G+ h; B) i& X) n   count
               
8 [% u) z+ p$ d5 Y                if r(N)>2 {
                 
                pkexamine dx240 ddy240
( |$ r& r" \. }5 k2 w' ]4 v/ X* V               
                local tempminphddy240=r(tomc)
                restore
( W+ t$ J9 b3 d! z! ?                 }
! c+ J6 s! n5 j/ ~$ u4 @! s* U. Y                 
  `3 L, ]6 v9 r7 d1 ~9 J9 o' J  O        else {
                restore
                 di "cannot find pminphddy240"
& Y/ d1 {, ]0 j0 T                 
               
                 sum dx240
( g. C* A& w! d2 I8 Y  n- _                 local tempminphddy240=r(mean)
               
0 @$ D* _9 J8 g7 U/ N6 }                 
+ O" _& S7 S5 m* R* A5 T. X                 }
3 V4 j& @8 G* |  Q! s" T3 r) q3 N                 
" X0 @& k8 p6 G8 Z' X/ n. i   
. P9 H+ l5 h* _5 S# w0 R   
! r4 s( k6 Q( D+ t6 g/ E   
   
1 @% k; B/ v: ?& \6 e   
* H/ f; j4 k# l+ L        gen pminphddy240=`tempminphddy240' in 1
/ b, k: O; d5 `" {       
7 Q4 y2 O) z. g* N( \; C( B& N       
        *-------------ddy240pminp右侧低峰值 -----------------
6 U8 I# ~2 t5 o. y- K1 W    preserve
; p) s3 V: M9 \8 l
   keep if dx240>`tempminp'
   
   keep if dx240<`tempmaxp'
) j* c' p7 q1 m3 n4 G( \% `! Z) N   
$ f+ ^( |- L/ T+ L   keep if ddy240<0
   
   gen temddy240=-ddy240
   
    count
! l  Y  u( W; P! \" b% x1 P, X               
+ q* F, X, S$ [4 C5 Y% p" Y) g9 p                if r(N)>2 {
2 c0 s7 m% G+ i" L$ Z- B                 
                 pkexamine dx240 temddy240
; ^& x% {- Q; M5 y- S   
' h# i8 r# R+ [3 w  local tempminplddy240=r(tomc)
                restore
                 }
- s: z, {- s! H3 z+ U# K: s                 
6 A! H1 @( u/ P  I7 H4 A        else {
               
5 ?; O& Z( |- C4 r6 ~. ^                 di "cannot find pminplddy240"
                restore
. `, Z! w4 r. A                 
& J- V6 ~9 B; }1 \; q7 q                  sum dx240
: k/ J2 C: _+ M3 Z  J' w/ E                 local tempminplddy240=r(mean)
                 
) T7 c  z, ]4 ]4 H                 
+ a, u; O: S  L                 }
* Q$ e9 ^- r# o. }6 W/ e" _' M   
   
   
   
$ z7 Y2 Q: g2 X! k% p   
        gen pminplddy240=`tempminplddy240' in 1
       
        *-------------ddy240pmaxp左侧高峰值 -------------------
# ^% j( m6 n1 g0 W& K* V       
8 }: a4 Y, I9 e* x         preserve
" q, t5 ^1 w( A" f; ~2 m2 r: Z   
   keep if dx240>`tempminp'
( S- H/ v8 k( a6 U0 z   
/ J: e+ M/ y8 D. u! D. ~! `. E   keep if dx240<`tempmaxp'
) z9 H# W. K7 i5 i6 Y   
   keep if ddy240>0
! |" B* [$ u  `" r/ F& [   count
               
/ N" c( a( A+ D5 N                if r(N)>2 {
! J8 P6 {% y) `; _% `# P1 k                 
                 pkexamine dx240 ddy240
   
/ z* i, \  k6 i   local tempmaxphddy240=r(tomc)
   
   restore
. z) E- {" ~& d. t4 p8 y                 }
' P& Z: [( g' l" Q  d, X0 e                 
        else {
% ?! \8 j2 o% h6 x" _; j               
1 ~. F) {9 j2 G. t' b* i                 di "cannot find pmaxphddy240"
                 restore
               
( B& _& y. i; A' _' G6 X                 sum dx240
                 local tempmaxphddy240=r(mean)
                 
( ]- Q9 ~8 W! @7 ^0 \/ V( M                 
                 }
+ v( N. ?9 g6 D0 w- Q   
   
   
   
1 g  c. o  a7 m+ y0 O   
   
' E# P/ s3 U1 }: b% \0 A1 ?        gen pmaxphddy240=`tempmaxphddy240' in 1
8 {$ w* M+ X4 i/ Y7 J       
; f' i) x) ~3 b       
        *-------------ddy240pmaxp右侧低峰值 -------------------
       
         preserve
  
   keep if dx240>`tempmaxp'
1 G; G5 I7 n, g0 f+ G. }   
   gen temddy240=-ddy240
   
   keep if temddy240>0
4 J* J6 i1 u& d& s4 V, ?   
1 L( l2 b2 x  r' a) Q% q   count
               
                if r(N)>2 {
: W8 L/ I- y" }' _( f2 e0 C                 
                pkexamine dx240 temddy240
/ ^8 D/ o7 s5 b5 N0 x3 F   
, \4 t/ I, [! y$ n( c. S   local tempmaxplddy240=r(tomc)
   restore
; M, e- ^' z- I                 }
                 
        else {
- e& R. `/ d9 J8 n7 i) T$ {                restore
* f+ f. E/ c) q: v                 di "cannot find pmaxlhddy240"
                 
! L& B; A& N  d2 T3 N               
                           sum dx240
                 local tempmaxplddy240=r(mean) in 1
               
5 r! w3 w/ u- t! `6 N( }                 
2 f& i& q5 ~0 r- U3 y! x' n% J                 }
   
   
* @" L3 M0 a! P" U8 d: g   
5 Y- _3 }4 a6 H" J   
   
, R$ e% Z- \( e  N, o* U9 Y   
        gen pmaxplddy240=`tempmaxphddy240' in 1
       
. ^( t& t" B* w1 I4 y        *-------------------------------------
6 \5 o. N, }& P+ Y) s1 ]   
   gen price5d=v5 in 1/20
8 u, M4 G4 P7 O) q6 c      
+ O7 u8 v4 u( I4 ~   gen price5yy=_n in 1/20
: n6 a& ^) [* ^) |   
   gen price5y=price5yy*(dy240max-dy240min)/20
& |  x; y0 I3 u( @" Y   
6 P* w) {1 s# @/ g8 L7 }   *price5y &price5d
1 ?7 S: t8 R% T0 W5 T" _, M8 [8 [1 `! F( B   
5 Z2 S0 f1 i/ R8 T8 n2 c2 ]   gen priceny=price5y in 1
   
   gen pricen=v5 in 1
' h1 A# {+ @2 y9 T1 R& B9 _
     
2 l5 ?4 T5 r9 U- Y$ x   *priceny &pricen
/ P# j/ b+ b) S0 D, Z: b" h   
   
9 @$ R$ \; e' h! J. V/ ]) a   *-------------------------graph------------------
7 l8 ?% b1 ?/ }" h; B  
% s& L+ j+ `- I; ~6 g" \  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)

; n) F+ @) i# H: B% g+ P& y* ---------可选项目----------------------
0 r4 ~) H' a& S*(spike yxis amean,lc(blue))(scatter yxis2 amean, lc(blue) msymbol(+) mlabel(amean))
/ y. T, f0 Z6 z" _2 ^; y8 w/ ?}
, @1 k$ N! `& X" zend
1 S6 P+ h$ X/ O4 ~6 N

maodeng

awer 发表于 2013-2-5 14:40
. b0 T1 g5 {7 ^) t" L  i3 _' p抱歉，最近太忙了

* g: }; k  e, ~, ?+ ]0 G300236 我用下面的程序试过，没问题

9 |. U" h' P0 F+ J: b. Z哇，谢谢，我仔细看看，这些字母数字代码晕了
7 j  s3 j+ R( ~6 o

喜欢喝冰茶

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

喜欢喝冰茶

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

- v' t6 s1 l; b8 C! b
这东西是开源的？知道有个统计软件叫stata，不过要收钱。

汉水云潇

准备试一试
2 C8 d* T% U# n) C4 M9 P先找Stata软件

汉水云潇

8 d, u  K6 j2 d2 W- N7 b
) \$ r4 d8 s3 ]晕呀
论坛复制程序代码
乱码夹杂在其中
很烦
对了两次
3 B/ \* m, V* G# K4 ~! {7 ?/ Z以为没问题了
, W: b4 a6 `8 e' F3 ?4 S& }, |( x
7 H$ @5 W5 r0 T- `结果
- L, Y( J3 g+ V) j- @9 E" ~执行后
还是出现
r 199
的错误
/ v- V7 a/ V& Q8 D8 F2 {

汉水云潇

试验成功，下面是三一重工的图
/ ~: \& [% x0 m/ k- P

peterbev

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

中秋不圆

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

bknite

这个帖子值得学习。

jeanny

感谢！！！

Hubert2000

明月回春 发表于 2011-8-13 20:17
0 {1 Y2 l6 L6 S1 o程序写的不错.
不过有一点要注意.通达信导出的数据,只对价格复权,并没有对成交量复权.所以数据需要手工处理 ...

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