设为首页收藏本站

爱吱声

 找回密码
 注册
搜索
楼主: awer
打印 上一主题 下一主题

[经济] 一个用stata 计算股票成本的小程序

    [复制链接]

该用户从未签到

61#
发表于 2013-1-25 13:15:49 | 只看该作者
如果我想看最近 一个月,三个月 六个月 的成本,怎么修改程序啊,谢谢大师了

该用户从未签到

62#
发表于 2013-1-25 13:18:55 | 只看该作者
我用同样的方法导出数据,发现有的能生成,有的就不行,不知道问题在哪儿,比如300236就生成不了,我的邮箱是19860225gaoyuan@163.com,麻烦大师了

点评

能告诉我stata报错显示的是什么么?  发表于 2013-1-27 21:37

该用户从未签到

63#
发表于 2013-1-28 13:13:51 | 只看该作者
. hello 300236, g8 f+ b" c! f6 ]
statement out of context
. ^5 n6 z) a5 M& b7 |$ {r(119);
8 F- e. R8 K. I! U& W
: p3 U( Z5 N+ y, ~大师 就是这个

该用户从未签到

64#
发表于 2013-1-28 21:11:43 | 只看该作者
这贴实在!

该用户从未签到

65#
 楼主| 发表于 2013-2-5 14:40:59 | 只看该作者
maodeng 发表于 2013-1-28 13:13 # B" i- v" m5 Q: ]* N4 M
. hello 300236% v, e& D/ b( Y% Z, I* {! g
statement out of context
0 G8 ?4 i; M$ Xr(119);

# z6 [3 q, f' K2 ~抱歉,最近太忙了
: p' ?4 G; w2 G; c1 x0 G1 `: ]/ M2 j- h9 s0 ]- ]' c4 P0 K
300236 我用下面的程序试过,没问题" |. B! B5 q: m
) ]7 r& G1 d0 w. L' p. g% G$ i
你把 D:\stock\test\source 改成你导出的数据存放的目录,再试一试5 h3 h  Y5 |& y( U

+ m8 {* V& [! ]) ?, tcapture program drop hello
) }- M$ P& J' z8 ^0 Q* \8 Z' m+ [0 K+ l1 N6 z# f
program hello2 ?: [2 |) Q5 ]; k# l3 R
5 ~- T+ ^9 b7 _/ R: B; }" H' D# g" j
quietly {
+ ~1 J6 N, T2 V0 j
/ O. Z7 V$ N3 `) F8 D% T# w: }0 Uclear
" W% e( z( y  ^$ X" c  p; k6 j0 ^1 p. e4 h# _' v5 c0 p9 x* a
cd D:\stock\test\source0 D1 o: E# I4 x' H7 P3 Y& S" S
' F! [; R" i' |) F" S- [
insheet using `1'.txt
6 M7 s& t* i) I/ T! d
; R4 ?; Y! x$ ?4 n& X' _  l   keep if v5>0
* q. j+ u. T' j1 b0 Q' h    " ]) ~( W; Q7 |; ?3 Q. r$ t" Q4 W- a
   gen vt=_n: |  X# n7 g+ }' ?
   
" b8 X  p# L( R$ L& p% V) G% p$ K0 H   tsset vt
% s. V$ j" _* d& l, V
# V0 J$ ?+ t) q( @1 d* u$ a' _   gen zf=D1.v5# I- a8 C+ e: J; g- y8 s

; c, Q) ~4 U  C- K) p   gsort -v1
  w9 z, T8 q, F% Y: s' }+ \, k2 v* z# c6 ?
   gen avol=v7/v5) C; L" O7 A1 V
   ) l, \% a( Y3 ]
   keep if avol>04 _1 [: @+ V, b& `3 U
   ) c) `% H* h1 t$ @7 T
   drop vt3 \4 l% w1 E0 r
  
$ T0 \' g) x4 x0 W   gen vt=_n
2 T& i/ D+ h/ o0 Y* i* {$ h   
6 h4 [! ]0 G1 q" n   tsset vt
3 g- d) a+ a! }6 K   
5 y" M, x1 j2 L    sum v5 [aw =avol]9 B( N( _5 Y  M/ `* J

  A' J1 T) E- v! j2 Q' S' z   gen amean=r(mean) in 1+ C0 Y+ l4 O$ }
   . p* X; e; E; C5 C3 H
   * 总加权平均值 amean, E3 D$ v" d' b2 ~% Y
   
  f# e, E) X. ]   7 _* M: n* `! }/ Z5 y* F
   keep in 1/240
! d7 }6 ]" [6 U* t: p  F" B; R, s* M$ h# N" o( @! g' [* N
   sum v5 [aw =avol]' W. T8 i4 t. j! [6 W- b# G9 ?' @) S

& C6 B2 u1 g# d- }! S- M   gen mean240=r(mean) in 1
" h+ h% G' \9 J- d+ Y' Y1 V# x   
$ D* e  J& G  q1 d0 ~$ Z0 A   *240加权平均值 mean240
8 g1 z( b7 v4 {2 o   
/ q" {# V& v" S3 ?8 @" f" h# ^- Q   set obs 500
) P% _4 K4 X7 E* C   
0 \, M( E5 g9 U, c   kdensity v5 [aw =avol],gaussian n(500) generate(dx240 dy240) nograph/ z( P" d0 c7 ]* K; k: Q
   ) D  L# g$ U  G9 U
   *240 kdensity(dx240 dy240)' B3 T$ I; G. {4 g+ m: W% J' d' M+ r- Z
   ( `, g  ?' A3 _2 ~
   sum dy240
$ Y) j  y' }& p1 U5 {5 Y   
/ t2 I7 G8 T/ K" F2 m# O2 Y   gen yxis=r(max) in 1
8 P. g; g0 T9 R' g: Q& c   0 b) I' M1 y/ c# `3 S0 R
   gen yxis2=yxis/2
, U6 ], H  m, E& O: T! X. _   ' p" Q# P* _0 G1 \$ i
   *yxis4 F( X! X1 ~: E7 A) U
   
# P4 b$ i5 V; D   gen dy240max=r(max)8 V" ^9 m) B1 o8 q+ r( r; y
   % s- v' {2 F- ^0 m3 C1 ]
   gen dy240minp=r(min)1 D, y. r  Z( S9 b  Y1 a0 e$ f

; Y. }% p, Q" t% n3 c    gen vtt=_n
) N+ k4 ]- i. J1 K
. s4 y) ?; G$ K5 j    tsset vtt' D) w6 G0 U; d+ K( |
+ N2 E! K6 H! ]; P
   gen ddy240=D1.dy240/D1.dx240
1 l+ n* n! c. y, q# q- ^6 N   ) u% t% W2 L3 q" E
   sum ddy240( j9 M; S. h" h& L, N5 x
   
- Q1 g4 k8 J4 X6 ~9 l& H   gen ddy240max=r(max)
' v1 M" Y9 f& W! V( N$ P3 m/ a) u   9 f" k. H! [) Z8 x# D( U2 D2 o' d
   gen ddy240min=r(min)
: L0 p( N3 K! B; h3 I. ], g     $ k. k3 }( g7 z6 Q. D1 E
   gen ddy240g=ddy240*(dy240max-dy240min)/(ddy240max-ddy240min)) N/ X( }0 e/ W4 Q% n3 K* r
   . H9 b0 }6 ~+ K) h/ g0 D# m& H
   sum ddy240g2 Q9 {$ [$ h& Y
   
% H. ?1 ^1 |" x/ q; s( q   gen yxis3=r(min) in 10 w  U7 U+ B0 u4 v
   
6 _$ I3 B: \' M   gen yxis4=0 in 1( `1 P: w) D) C4 X
  % q. l; `: Z3 Z: G  O
  *ddy240g * h/ l; l) b. W- G6 N7 P
   
7 E7 B$ Y' b. b( v; Z/ T' \; f   ) f0 v) V* i7 n* V0 m
  preserve * M0 o0 v- S) S# k; ^+ m% e* c( t0 _
  
8 h5 k- t+ G( ], z- ]! Q* O" `/ Q+ _  sum v5 [aw =avol]
/ V$ i" O/ [' E* T. h, L  9 Z/ K/ }& g% w( \; R& c
  keep if dx240<r(mean)5 T+ d7 }5 c6 U, z4 _
  % |& @, e8 X7 Z) ~: [
  sum ddy240
5 x. p8 p  O& D  3 p0 a3 b: f, n2 v( C& s. @
  if r(min)<0 {
1 c) V- i( d' {; P0 J        
8 }' S1 W: Y, u                keep if ddy240<=0. v- m3 v1 b2 y5 l4 p- u& Z' q( n
               
& z, |$ B* E  `4 c4 Z0 g1 z                count+ M/ L! u9 p* p
                ! z$ {; ^' L( m' _0 v; V
                if r(N)>2 {
% A, o$ y0 Y2 F! t; a, u9 M1 ]                 
; q+ |: n+ m6 Y, Y( F4 v; g: a                 pkexamine dx240 dy2406 q( B' Z8 R) ]8 Z4 {& m; o& w
               
: o! [' T1 ]5 F                local tempminp=r(tomc)
7 _- t6 z9 t: y  Y8 F+ Z                restore 9 @8 x7 i; E  A& @) g+ r8 g
                 }
; A$ T' w; j# F  P+ x9 a. ^                 ; K9 ?0 Q/ j5 N' a& B* T& C  b
        else {+ d) q) P) }5 q, ^3 i( s
               
' w% t, f& ]$ r. t6 G; B4 K! T5 K                 di "cannot find minp"2 L& d- r: {9 E
                 
% s. E1 r$ t3 Q+ [! S1 W                restore
9 g% _9 E0 b* Z6 a( p: x' ?" w                 sum dx240! R% H6 D$ g. h: p5 }7 e
                 local tempminp=r(mean). _' c: Y4 q1 _2 e; S& K
                 
' E9 a. [8 E2 L7 R" Z$ R, r                 }0 l: n; ?+ K. \: O' D- Q% i3 e7 l: D8 e
                }
3 I# ?  e, J- u1 {# `   else {9 C& @. q' a( t% X. `
                 
: a' G- c1 t1 }: V# l4 t                 keep if ddy240<=r(min)  C8 ]" p$ d9 I* G, p8 l
                 
" E; X3 g/ x- Z, s9 w                 sum dx240/ Z7 Q- }0 J8 D. [" h+ _
                 
- r/ {" ?4 ]% N, L( T! }* P7 l                 local tempminp=r(min)
* N, G0 w, p) h9 o1 f* d( k) m! Q                 restore / S( U) \% f% M; Z
                 }5 e7 d; B- b' j( d9 K
. M" b, k0 S0 T5 L6 u

; [! Y& u8 q! x! D" H1 C* {% j7 e" A- z* F  M  ]5 x2 G. a% y3 \7 y% _
   gen minp=`tempminp' in 15 Q' W" V; Y* i6 v% G$ A  G! I
   *240 minp" j5 z. |' G2 L

3 ^2 |4 R% F- j% r' b$ d; [' x1 _1 W7 }6 g; W, G. P
   
% D7 j5 i+ T4 ?2 l0 ?( c& E1 Q) v   preserve & G: M; Y# K9 x$ }# z" I8 g' I
   + K7 K- \, g1 l7 o+ @, s$ d
   sum v5 [aw =avol]6 ?; B. _. L8 d! T4 r1 D
   
9 B4 x1 y5 n: ~( g   keep if dx240>r(mean)" m0 D* N* Z* F1 E2 r
   . Z6 T& X1 L( @4 J% S; K% I/ o
   sum ddy240; Z; T0 p0 n8 Y+ V: c: I/ p' A% o
   ; n2 z+ ?- R& r; D5 t; E2 b
   if r(max)>0 {
+ H' _5 [! X5 K) z+ Z               
/ b" n( k# j, |                keep if ddy240>=0
" v" `  V$ a, W# x) E               
' N$ h# R9 ~/ _5 F                count; R& Y' Y7 |6 p. d! l% p9 D# Q6 `
               
5 K/ |  O1 Z/ _2 z/ ?4 O- R  T; h                if r(N)>2 {3 A* m% R3 P1 A& A; Z2 W
                 8 f% J6 ~- w7 Q' t
                 pkexamine dx240 dy240* l. j  p  ^5 H2 H
                3 P1 W, U: ^# h( z0 h/ Z
                local tempmaxp=r(tomc)   O% x# D0 {9 }) A7 L
                restore
$ c3 f- I4 F* {1 W                 }
% V/ ?! Z" v4 {# Y                 & W# ?/ k( L6 T6 |2 i) e3 V1 W
                 
3 {2 e0 E$ R9 X' G7 [  f6 [        else {% w+ ?3 c' O  ~3 g; P
                restore . q+ y7 M% A) ~% a, }9 a
                 di "cannot find maxp"
$ B! q; I. {& \4 w! ^' |               
6 _, m3 v- i$ _1 F) b                ; k( w- B3 a9 w: {$ ]
                 sum dx240- C/ \2 Z  U6 G( r% v3 S. \5 ~
                 local tempmaxp=r(mean)
# q; Q  K: v' f8 ]" {3 x                 
) [4 I: I9 r; S' p6 p                 }
# [: A% o! Q  ~. N1 S$ N; U                 4 @+ L% ^/ Z/ U  M% ]$ c
                 }
9 R$ A& h, X% U) m% ^9 B                 1 v& h; ?' q2 ]) }$ N
    else {7 n7 h# r- J, @# A6 C& t
       
2 f- e/ ^+ }2 |           keep if ddy240>=r(max)! r3 _+ v3 q+ @9 ]
          
) s! V( T% u3 x7 Z  u2 S           sum dx2400 z' ^1 X( m  }' O! a0 P
           ) u4 B& M( h+ _; ~* q! a0 F1 @. c
           local tempmaxp=r(max)
- W# e) Y: G+ ^, _: @% R           restore / M- u4 ]; S0 T' O. ?& `
           }. [, z! z/ |7 i$ ~
       
& c& d- M9 \1 F4 b4 r) ~. L: P* p$ l6 \3 c" r1 C
  
. `2 i) A/ J% w* D. O( b: D/ Y8 H- O" \
   gen maxp=`tempmaxp' in 1( g2 V  `( A5 i1 \+ u
   , \2 w" c0 B7 E. C) Z& C
   *240maxp
! M! f9 ~* F' Y% I* H   
3 }+ q+ X1 q, C6 N# m8 A: x4 f/ u  L   
0 g" K; Y3 \& p' y3 H! a" w   *-------------ddy240pminp左侧高峰值-----
7 a$ l" D) M/ O2 k* F: y   6 }' B  L) F0 u& y
    preserve ) u3 C2 R' h% Q  Y" d, X
   
/ j8 \( m# {: ?' q# M- }0 V   keep if dx240<`tempminp'0 q# C% g3 H& W$ B3 c
   5 i- u) S3 u5 E+ a# A1 ?
   keep if ddy240>01 Q' R5 f0 l" {% z6 W& c& E% {
   + ?5 m5 ~1 w$ U( _! `& I5 ~+ S
   count
. m. e0 X/ P% S! W0 p1 E. Q               
; N6 F; ~' Z0 [2 o" q" |; `" {8 j                if r(N)>2 {( G" x) ~+ E! m
                 
, \7 ^* P  W( w" C; q3 X5 A! C3 A                pkexamine dx240 ddy240
: d/ U; i- ~$ o! C/ A                7 K( [6 t1 O0 D: n8 y( S. _0 t
                local tempminphddy240=r(tomc)
7 q7 t( T7 l7 t5 Y% g5 [                restore
9 b; h- m1 M3 G: J8 {                 }0 `# ?/ i7 Z0 r' X1 o) Y
                 
# R) [  d- C: d: T' c        else {& Z4 J7 m5 U9 R( ]# p- h& ], l9 U
                restore ' j4 u; ?& k- p8 l- y5 k! j
                 di "cannot find pminphddy240"
: N( M; f, x8 m4 F3 H  e$ H' C2 n                 0 y% f. c* M* C1 I3 S
               
0 `  H% {( o6 y7 k0 f5 i9 ?  Z                 sum dx240. s7 j8 k9 t. ]9 d
                 local tempminphddy240=r(mean)
) ^3 [; i& u( l9 `' c0 x               
7 [. U' _2 L$ @( v. w                 * R1 P  V6 q! a  m: A8 o
                 }
4 O5 `- c  Z9 |! D8 Y" t; G                 9 Z4 h9 D# k" f- W2 d2 |/ p5 P9 M5 D: f
   9 k% ]$ J9 b4 o6 M* Y
    1 N( D6 B2 d/ Q; m2 ~
   
; T6 @5 S6 \( s! q   4 {. S. d5 h( I" z
   9 B( A; D5 H+ z4 Q
        gen pminphddy240=`tempminphddy240' in 1
' R' L  b- j" C  q& Z: i- N       
# b* o  \5 ~4 N; u3 I0 |2 S) n       
# w4 t4 g/ ]' ^% @& q5 u- o        *-------------ddy240pminp右侧低峰值 -----------------
! U- t% K: Z- ^- }  {* ?% D- y* g    preserve
" s  F6 u1 r/ n8 f0 \5 Y # R* |1 a0 T# T; c, u
   keep if dx240>`tempminp'
' }. f7 Y( E" z: g6 H, g   
4 _: y- T6 d! h+ l& d; K8 w' F   keep if dx240<`tempmaxp'
  E; X9 u5 Z: |8 A% D. s. i   ( `8 H3 ]* W$ a4 Y4 ?( a
   keep if ddy240<0
. Q$ Z4 ^% N& c5 R! e   - a' U/ w+ M5 b. R3 t" J
   gen temddy240=-ddy2405 u+ C8 D+ y$ I0 I* S; r
   
) m7 y  U. {- R  c5 V& x    count
# T6 Z# m9 y) }+ D9 E               
0 c" `6 v" C- e3 O) k& R2 @. u                if r(N)>2 {- e* |/ d! m0 L6 s% x* ~3 J
                 3 R) ^9 R, A; M
                 pkexamine dx240 temddy240& y* Z( k  Q' c5 q( {
   
, J% D( l8 Q6 L/ V4 ]3 X0 P$ V  local tempminplddy240=r(tomc)
" a% p3 P- A7 Q9 z9 I                restore" S: q8 E  {, Q  M1 j7 ^' W
                 }' J* O! \+ i* O" p; X
                 
+ a4 |; q) _: z' i- z$ y        else {
; x1 w* H( u3 {& N( t% a7 D/ E               
0 C# X' y% s) ], }7 s$ M4 p( ?( v' u  T                 di "cannot find pminplddy240"
) B# N7 T  V& {! T( {! d1 {- D6 r                restore
; q4 a# I4 d6 f5 P- ?- ~8 Q                 
/ v/ M# a& u- }. b4 K% M2 h                  sum dx2405 w1 G6 C: W1 b3 A# M7 D
                 local tempminplddy240=r(mean)
: W1 \8 }4 z) b1 |) N                 + L8 a5 N$ K1 U# W
                 1 {$ j. `  K# [2 L
                 }% E7 b! C) f. ?; H8 {% q+ h* _3 g
   / e4 k- c* Z' P3 X( l' P
   
. Q+ `9 B# l6 P9 }   , B- P! [( P( a6 x% g3 O
   
3 E2 E" |, M1 H( y$ e) U- ~/ ~5 ?   . L. P: [4 B1 a* {4 P
        gen pminplddy240=`tempminplddy240' in 1
7 T/ v0 r( E& u& x- i        5 P( h" u; u. q0 a4 G% }
        *-------------ddy240pmaxp左侧高峰值 -------------------4 r& M3 {/ z4 x5 O& k5 l
       
  d9 |( w3 G3 d% o! L         preserve 9 U& ?8 e, q+ Y" L9 Q3 I9 n& ~
   
7 b( _1 j0 E/ \& D' Q: t$ Z   keep if dx240>`tempminp'
" }" i5 R9 W3 J+ r' J8 J8 V   
& t. h7 n" J& B% z* w% _. M* |% [   keep if dx240<`tempmaxp'
& O& z* ^! [3 p& J2 A2 r- i3 u$ @   + ]- n& ~4 u4 j# {& d( {1 F7 W
   keep if ddy240>02 H; U6 K( c0 C% T
   count" J5 W; E- L+ a/ O) G
               
: V" N+ ?0 e+ k5 t+ ]: a8 F# R9 O4 ]                if r(N)>2 {4 `7 D9 c' p; }. S# ^
                 
, C% H+ A: ~: J) V1 W+ t                 pkexamine dx240 ddy240+ ~# r$ s# C8 @3 @1 a) D
   : P* y0 K/ C- h' [5 ~9 E
   local tempmaxphddy240=r(tomc) ( f' ~* h" c- j1 A! c
   
0 j  a1 L1 Q" U. M7 R8 X   restore
; l7 D" }% ~+ }9 U: P* d                 }
( V; n3 _- _' E4 z                 
- ^7 |- X5 n6 B) S        else {3 z3 |) H. o1 X# D
                * ]1 S$ Z# g- x! |/ W* i
                 di "cannot find pmaxphddy240"# H6 L1 z  k4 [. Y; T6 J
                 restore3 i4 M3 c2 A# m% p& w
               
. k2 `+ R1 c! J! J! x8 r: h; p$ v4 r                 sum dx240( C9 g, }" `2 R6 q, f/ k4 ^0 g% C
                 local tempmaxphddy240=r(mean)
" o( }2 @% @! p4 r& F                 
" n" e& _1 T) _1 s& `: W                 
( z& u' X3 n  b                 }3 p) i% Z0 Q, n
   / O) z; y2 w" u0 r9 O
   # F! ]0 F0 A, T) E' {
   
5 ]5 k- A7 s; i. Q, D# c   
* r" }/ D) T" o+ x! M; q3 b   ' k* [% k6 F, C6 {! w) T( U
   
! r& r2 X7 ^; ^3 R5 K  S        gen pmaxphddy240=`tempmaxphddy240' in 1$ T1 R) M: b5 `0 s- x2 U
        * u* g/ \0 T8 y9 V' O* O- F' P+ j' L2 s
       
& D( G8 y3 I$ H8 X4 Q/ Q; [; F        *-------------ddy240pmaxp右侧低峰值 -------------------, @( R  E, K: p
        . V4 f" l) ^, Q# ]! e* {7 P' P
         preserve
; B! l% q! Z+ B8 v  * D/ a# M' W+ T" n
   keep if dx240>`tempmaxp'% i! s2 E  i0 D
   
- x; G  x" e4 v' z   gen temddy240=-ddy240
8 k1 o3 u& N2 E   8 T1 b& U4 g% x) D
   keep if temddy240>0
9 u. B6 s& w9 |0 e7 m   : D1 ?) N% I. N% ?: r  O
   count
# [$ P' d2 l" L# ?& Y               
& U4 M: `9 @' y$ U% c; F                if r(N)>2 {
9 \5 |8 ?) N" I# s2 f" D$ Q                 
3 P& k$ Q+ p8 J* f. r                pkexamine dx240 temddy240# m: v0 E" v: Z$ X2 q
   + e! r2 @; H0 l  K$ L% v; i
   local tempmaxplddy240=r(tomc)
2 J2 {# h( X7 q+ s   restore$ S5 P6 s: y  E( A! N# o8 r9 o1 }
                 }
8 w  w' ~/ R1 e; J, A* P* \                 " B7 ~1 n; ^1 i. f$ u% O! g
        else {+ ^3 q. n( u- N' C9 r) M# q
                restore/ Q- A' X  g/ Y+ T. g* X1 o& y: t
                 di "cannot find pmaxlhddy240"
7 t2 L# J* ~) c( V% H! }- w                 : p+ E5 {* \# T6 O8 b
               
3 @. X# I2 V& y  G4 X                           sum dx240$ s% m2 B9 y$ ^5 U: O) Y! [
                 local tempmaxplddy240=r(mean) in 1: A& E9 T; R3 P$ Z
               
, |# g( U  G. _: S2 k                 
" I; V5 G" C# h  N4 P% I                 }8 D& m! d3 Z5 v0 w1 X# I  w9 u3 i/ d1 e
   ! l- Y; X4 R, R* R3 ^4 T
   
( K, Y7 Y& z- J3 J    - a, s5 {  x* |* Z  K- @, G
   
0 J% O5 L8 P2 D" h* [! B$ q   $ c' R( g/ e1 m0 o
   , `, p. O( U  a# ]3 \! N) ~
        gen pmaxplddy240=`tempmaxphddy240' in 13 t: K9 c, N" @
        8 a% B, J/ m, [1 A$ f0 Z
        *-------------------------------------
( R/ P$ o4 T% l# c# `! x0 B   
  x( `; B7 X) W8 l: I' e4 p# l   gen price5d=v5 in 1/20
- W. @% d1 @1 p7 u" T+ K      
* T" ?5 _9 n5 Z% m   gen price5yy=_n in 1/20
' z9 q2 p2 Z" V   
1 o/ P$ D& y' R" V   gen price5y=price5yy*(dy240max-dy240min)/20
& t: u: p) u* E   $ {) \5 w1 {, t9 [( P
   *price5y &price5d" R6 b) s& M7 A
   
  U; ^3 k" f9 n0 M   gen priceny=price5y in 1
* q; g1 g. E; ~   : H3 P$ ?. [6 i: ^2 b
   gen pricen=v5 in 18 S3 x6 H, E4 e+ r4 n: ~2 D" V
+ [, g( U' K2 v
     
5 A) {+ K# z' Q1 D   *priceny &pricen7 G6 {- `3 z" q) N3 ^( ~! n$ ?9 N$ w
   3 E  O* Z$ p+ ^) W1 B0 d6 E
   : P3 s2 N1 u4 f2 b0 v
   *-------------------------graph------------------
0 X; i1 o# }6 T( j4 T  t  
' C$ `+ Q+ }5 Y- H7 K  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)5 g5 _/ }( P  M/ U2 g$ h

# ~" K: I  K% m. `& A; C* _+ m& c* ---------可选项目----------------------" ^$ w, E% A+ b$ c
*(spike yxis amean,lc(blue))(scatter yxis2 amean, lc(blue) msymbol(+) mlabel(amean))& _+ w& _8 y4 S8 B# {" ]. t
}
6 t3 M% b& z; g3 g6 Vend& d5 x$ k* m! U" ^* o, D' _! N% {

该用户从未签到

66#
发表于 2013-2-7 00:03:44 | 只看该作者
awer 发表于 2013-2-5 14:40
( ?- @( T( v3 C抱歉,最近太忙了) M3 |. ]6 y4 m. Y$ d

3 {5 I- o( T  A5 Z2 v& h( I* B300236 我用下面的程序试过,没问题
7 y. D0 X& M' e; d3 d
哇,谢谢,我仔细看看,这些字母数字代码晕了
8 K- a- A2 u/ B8 _% \& @

该用户从未签到

67#
发表于 2013-2-7 04:05:03 | 只看该作者
stata是不是商用的,R里应该有类似的东西吧

点评

老虎推荐的一个处理数据的软件,很小,只有10m左右,很好用  发表于 2013-2-8 01:27

该用户从未签到

68#
发表于 2013-2-8 07:16:33 | 只看该作者
awer  老虎推荐的一个处理数据的软件,很小,只有10m左右,很好用  发表于 5 小时前

3 l  w  t2 f: i* _5 c" I+ N1 L# z8 Q. J  H
这东西是开源的?知道有个统计软件叫stata,不过要收钱。

点评

就是那个,找破解版  发表于 2013-2-8 09:15
  • TA的每日心情
    郁闷
    2018-10-21 12:19
  • 签到天数: 7 天

    [LV.3]辟谷

    69#
    发表于 2013-4-17 17:14:44 | 只看该作者
    准备试一试
    & I- Z. H5 E  S7 ?先找Stata软件
  • TA的每日心情
    郁闷
    2018-10-21 12:19
  • 签到天数: 7 天

    [LV.3]辟谷

    70#
    发表于 2013-4-17 22:21:02 | 只看该作者
    本帖最后由 汉水云潇 于 2013-4-17 22:27 编辑 - P1 `3 E# K; {/ d! S9 }9 {, d

    6 c# L* i7 ^  A5 f晕呀( [) R+ c+ p& S: o; W
    论坛复制程序代码
    8 @1 Q, I5 o* H, N+ E: W2 O乱码夹杂在其中6 g7 x8 j2 |- x. I1 E) e
    很烦
    4 v( D  ~$ A# W2 ^) @对了两次' G2 Z7 C5 _. z/ I. Q% O& p
    以为没问题了. B% M+ X- F. Y4 K3 U$ P7 o

    9 v0 I$ p1 f% T4 v+ e结果
    5 y, j5 ~8 [9 c1 [# p2 s& W执行后5 y) {! t5 @7 C9 C! _
    还是出现
    : q+ c- o: Y; Dr 199
    # f* b# h9 n/ G4 X6 `0 Q% `  S8 R的错误! s& L% D1 Q/ }
  • TA的每日心情
    郁闷
    2018-10-21 12:19
  • 签到天数: 7 天

    [LV.3]辟谷

    71#
    发表于 2013-4-18 07:50:56 | 只看该作者
    试验成功,下面是三一重工的图
    8 P# l3 E! ?( l2 i
    ' p& s# D/ A& Y# y: W: W- {

    点评

    我的怎么没有红点? 三一上市以来的成本才1.8元  发表于 2013-4-18 08:02
  • TA的每日心情
    开心
    2018-8-27 15:57
  • 签到天数: 2 天

    [LV.1]炼气

    72#
    发表于 2013-4-26 12:46:50 | 只看该作者
    先用用试一试,不知道准确不

    该用户从未签到

    73#
    发表于 2016-3-12 15:46:06 | 只看该作者
    噢。。。这个帖子好,比较直观。
  • TA的每日心情
    开心
    2018-10-15 15:53
  • 签到天数: 80 天

    [LV.6]出窍

    74#
    发表于 2016-3-21 23:51:14 | 只看该作者
    这个帖子值得学习。
  • TA的每日心情
    开心
    2016-1-11 12:47
  • 签到天数: 2 天

    [LV.1]炼气

    75#
    发表于 2016-3-22 10:41:13 | 只看该作者
    感谢!!!
    回复

    使用道具 举报

    该用户从未签到

    76#
    发表于 2016-3-23 10:23:05 | 只看该作者
    明月回春 发表于 2011-8-13 20:17
    , E" O. V4 W. b# R程序写的不错.2 G% M3 c3 L; P# l# n5 `& o
    不过有一点要注意.通达信导出的数据,只对价格复权,并没有对成交量复权.所以数据需要手工处理 ...

    7 x. V" W: D" }6 {+ @& D/ H1 w老哥,也是山寨这篇文章找不到了,不知为何?

    点评

    日久年深,早就忘了。见谅  发表于 2016-3-24 09:24
  • TA的每日心情
    慵懒
    2019-4-17 22:04
  • 签到天数: 633 天

    [LV.9]渡劫

    77#
    发表于 2016-8-3 13:28:29 | 只看该作者
    提示: 该帖被管理员或版主屏蔽

    手机版|小黑屋|Archiver|网站错误报告|爱吱声   

    GMT+8, 2025-9-6 17:13 , Processed in 0.048879 second(s), 20 queries , Gzip On.

    Powered by Discuz! X3.2

    © 2001-2013 Comsenz Inc.

    快速回复 返回顶部 返回列表