设为首页收藏本站

爱吱声

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

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

    [复制链接]

该用户从未签到

61#
发表于 2012-11-7 11:46:47 | 只看该作者
本帖最后由 墨轩书童 于 2012-11-7 14:37 编辑 - D8 h& z, y  S- c4 B" L; U  g. c. @

* ?% D& o1 ]  E0 o" r感谢楼主!!
  • TA的每日心情
    开心
    2016-1-12 14:27
  • 签到天数: 1 天

    [LV.1]炼气

    62#
    发表于 2013-1-20 17:03:12 | 只看该作者
    本帖最后由 海阔 于 2013-1-20 17:08 编辑 + y# N# _' {; q+ y( i, U

    + g: Z+ k6 i6 b& v按照版主的步骤操作到了最后出现% \6 |  d6 u$ a$ G/ C( ~& e$ {

    5 R( p$ T' u# l+ ]' X* g. hello 600031
    8 v9 v; F, `. a# V+ XavolDo not found
    4 G5 Z% V) e; i3 u8 @: ?% \1 \---------------------------------------换成其他代买,也是一样;总是出现这个问题,怎么办呢?奇怪

    该用户从未签到

    63#
    发表于 2013-1-25 13:15:07 | 只看该作者
    大师,那两条曲线代表什么啊) C: v6 M) X  o& @( K

    该用户从未签到

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

    该用户从未签到

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

    点评

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

    该用户从未签到

    66#
    发表于 2013-1-28 13:13:51 | 只看该作者
    . hello 300236. ?* l: ~2 J/ P. g( Q
    statement out of context% c0 F( u, |  ^* a, |: P. x' ~0 q* p
    r(119);, m* j; l- k1 |* S8 i0 @

      B0 E" l6 Y8 g& I& N+ D大师 就是这个

    该用户从未签到

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

    该用户从未签到

    68#
     楼主| 发表于 2013-2-5 14:40:59 | 只看该作者
    maodeng 发表于 2013-1-28 13:13
    3 y. d+ u) w! n& J; w$ M2 o# c. hello 300236
    : C7 b$ O5 O6 }' u4 |2 n/ estatement out of context2 P  N& c/ B. z# s7 g
    r(119);

    ! d# G" n4 ~# V$ d抱歉,最近太忙了& W# Z+ x# o% J! e. ~9 O: ~$ l8 O
    5 F# x+ K! B  y. j2 _4 U; z
    300236 我用下面的程序试过,没问题
    4 i9 ~7 `: X* s4 G1 u6 d) F* M. R% ]  O7 h+ x
    你把 D:\stock\test\source 改成你导出的数据存放的目录,再试一试( t. u: a3 k! W, `" |* I
    ! T+ _2 @, s! T; G
    capture program drop hello
    . O  M7 X5 D* |6 }  L- J/ k8 g- S$ P  ^' O+ i6 p
    program hello
    * a) ]- B5 Y* d, d6 L& L* G
    8 {! k1 u' d9 B% o  q& {: C quietly {6 R5 |/ g; J9 P0 b$ n
    ) j. K! w2 r* o1 U; K! ^
    clear6 k. U$ {5 U  w/ E( u

    3 T, t9 K3 C9 T- P; G; ~cd D:\stock\test\source, d# L7 D: I9 P! j; O6 m
    0 l! {: p& I8 s# V
    insheet using `1'.txt- s$ @( A2 D$ O, }' W
    , V  A5 d! l: u  u9 e- |
       keep if v5>0
    . L9 w( g  N3 v% u$ @9 r) N   
    ! e1 r' ]2 E1 z/ n6 Q* s! k* S   gen vt=_n
      f8 y! A6 G2 K$ }- b   
    " U- y, h3 y" h0 g% }7 `   tsset vt  {$ [& M# \: i* J# t( p! t
    , b0 V* q) \  t1 n8 C
       gen zf=D1.v5" I4 {& S2 r; F7 \

    + C. K$ @6 \8 @, A& J   gsort -v1
    * [. s8 R. a2 @# C) U7 b
    5 C, a2 Y+ t5 J8 ~5 O+ D6 z# B   gen avol=v7/v55 V7 G+ Y+ I# j* B' {
       # F. v" a$ J+ b7 c! e9 I$ [3 s
       keep if avol>0
    ) u. {# ~0 z3 ^   
    7 \7 z# c3 M7 R% X5 j% o( p; y   drop vt6 l) h4 g) a0 F9 j
      : G( S. k9 k) u
       gen vt=_n3 }& y2 y) s2 Q
       . H3 x) g% U( s* |( U) b
       tsset vt
    2 O' a( A' ]8 ?: d; F. S4 V   
    & q/ c8 J: j8 F. \* k& T    sum v5 [aw =avol]3 Y3 r' }) }( E) h6 U! ^  O
    ( k( \/ e6 H8 j) _! f7 T9 t
       gen amean=r(mean) in 1
    8 y6 U1 T$ g8 h5 v! {+ V   4 W+ o# @/ A, H  k
       * 总加权平均值 amean9 _4 l' D* h2 S% |1 K+ x3 H
       
    0 P# P5 }* i9 f) z/ J   
    % j$ ~0 W0 m  J   keep in 1/240
    7 U- D. m! _+ f' W$ m6 ^
    ; H$ z- ?+ E' V  }7 y0 }   sum v5 [aw =avol]
    $ m7 C, T+ E) ^( y
    ) O0 r' R# a/ ?6 ~3 I! P& Y# b- t9 Q   gen mean240=r(mean) in 1
    " O1 _9 W! Y- G. i8 f   
    / V$ L) k) @  T   *240加权平均值 mean240, o& f! i- K) w& `" b' _+ s: J
       
    5 P4 `$ E3 X0 k  b   set obs 500
    " B& f' ]! n" Z   9 T& m! T! Y8 v$ C2 ?
       kdensity v5 [aw =avol],gaussian n(500) generate(dx240 dy240) nograph
    . Y+ ~8 T& ]7 F   8 f5 i# f8 `1 z$ K" g! n
       *240 kdensity(dx240 dy240)( Y1 g- c6 c2 Y1 v
       
    2 n; w. ?2 C! Q$ c   sum dy240
    6 ^7 {) W% Y; J# z2 h1 ^   ) a/ B3 r" |5 N9 j! B3 n
       gen yxis=r(max) in 12 ?# A9 \2 G# `0 D2 A5 h
       
    % I* A. C/ n, A& q   gen yxis2=yxis/2. U$ z* G0 p3 L6 \5 \$ B# E1 I+ z
       4 y4 I# t2 w) _
       *yxis* \, q/ `$ O' B
       
    & n. h3 ?9 L% v% ]1 K   gen dy240max=r(max); J! i9 i( V5 a4 i; G
       3 ?5 q9 A* U4 f0 ]; @! `
       gen dy240minp=r(min). Z4 o, t; Y  D) X

    6 q* ?8 E5 ~) \    gen vtt=_n
    5 i" v9 ]9 ]& ]+ V
    ; t: B2 P/ v( J' w) K    tsset vtt$ k" P$ s  h$ e5 x3 [& }( L! N
    ' p3 i, n# G1 z& c% p' J( V9 w
       gen ddy240=D1.dy240/D1.dx240; ^% y5 s  F# u* ~& C# u
       
    + I5 i& T) f5 ^0 r8 Q   sum ddy240
    ) D1 n+ |3 a/ W   
    $ H9 q3 b8 N& e   gen ddy240max=r(max)% a' c2 O5 I. O0 U4 E$ z' U, O: d( e6 c
       
    ' [" B: x" G  M: G5 ^   gen ddy240min=r(min)5 g0 J8 I6 B7 ?1 Q. B" Q2 }
         8 Q3 g# l; |- s0 |3 i& @' ^
       gen ddy240g=ddy240*(dy240max-dy240min)/(ddy240max-ddy240min)
    * n- |& s; Q, ~) Z! N8 m! w   
    8 r. w# x; F% `3 b, B1 E# l  J   sum ddy240g$ f( i0 {+ {; s) r/ W: v
       7 T- L+ `( ~' _2 i
       gen yxis3=r(min) in 1' C0 _! A# I) z: K
       
    + a0 R  o+ g' p5 k7 X   gen yxis4=0 in 1$ k. \7 p3 W  K  U; L1 i  m, g1 |2 q
      # g- r( D) P8 h; z9 [3 d* T- y
      *ddy240g
    9 F. `8 E2 Y5 s: S   & C' \% {8 m0 W
       
    ) S2 k9 ]5 K# [) {% t+ y! I  preserve
    5 i9 i4 Y; P; n3 K+ K' y1 X( _  
    1 q: U0 M5 ~% }$ ]: C. j  sum v5 [aw =avol]
    " H* R, G( S% z/ r3 d) k  
    $ |0 @+ f2 `" L6 R( X) q3 I  keep if dx240<r(mean)
    ! N( Z+ |0 V: V' D  
    : T. j  L. `% I6 s% u" ]: V4 u  sum ddy2406 E; }) ]: h7 t5 s2 m7 H
      
    # P2 b+ z& q  m$ p; z6 w& d6 M  if r(min)<0 {+ O' @0 q0 z- f6 q7 G1 n( I% h9 N
            4 O2 g: P/ o: {9 `
                    keep if ddy240<=0$ V2 H, o4 F; E/ s/ D6 ^9 m" x/ Q8 _
                    % M$ D) J0 `0 ^4 ?  a4 z7 N
                    count
      m/ u4 l! u: ]                3 D# t/ e8 h. x. u+ |% N
                    if r(N)>2 {
    4 X- A6 k3 x4 o, c* n                 
    + P: C6 @8 ?! Y3 E; C                 pkexamine dx240 dy240
    3 \: a" A2 K- d* K  K% p: H               
      _& B: S: D$ P! P6 q$ F                local tempminp=r(tomc)0 e+ w5 J9 E: J
                    restore 9 E7 ?5 J: r, v2 e
                     }. c  f% a# p" i( c+ k4 m
                     
    ! O+ [- Z4 \( C+ g" t* R2 E        else {
    & _$ w& Q2 g0 E$ {4 t  P; F8 \$ @                & H6 ?: Y. c, s  i
                     di "cannot find minp"
    2 C$ i( h, u" a2 u6 M1 |                 3 m! H$ e! u) M! L5 T+ \8 w
                    restore
    " G  e: v' V% }  r) ^                 sum dx240
    . Y" i2 u! v1 c+ @2 z8 T                 local tempminp=r(mean)2 c( U9 M/ H* |0 ?7 _; y. o
                     
    + O: ?  n9 y0 I2 L( D$ ^                 }
    % ]9 P1 r; o& f% i" b; I+ A. I                }
    * T4 {, q$ F# O' y   else {
    % a. t( ]- S; w) g& M                 5 T# h' _% T/ v0 ?
                     keep if ddy240<=r(min)
    & F6 d) Q0 T2 p2 w9 _                 ) t# E9 w9 u6 g6 h+ I1 `! ^$ x
                     sum dx240
    : ^' C' ?/ e4 M1 {6 U3 q3 y                 
    8 ?9 A& Q# @$ K+ z* @                 local tempminp=r(min)4 Y/ A( p- A# [
                     restore 2 m3 ]; h0 ?( q" _- |5 |
                     }6 e0 d6 J5 J( f
    ( @; P7 P. {2 T& Q6 \/ s: @
    ' J: \; X. }/ D/ p  R
    4 z6 F2 I8 U# s/ f! }! p$ [
       gen minp=`tempminp' in 1
    + E0 l& q$ a) n3 @+ H/ W   *240 minp
    8 O; L. w& r3 j. J$ R. w* e; C, r; h7 v

    , ~: |0 B! V6 |5 A* f     V+ H$ T6 b0 e
       preserve
    : B4 ~( a, n0 ]9 l   % k5 L: H* B% D# h4 F- ~7 V
       sum v5 [aw =avol]
    & h1 Y7 p( T) a" A7 S   
    / q9 G8 Z+ X9 a   keep if dx240>r(mean)
    * t$ t9 V* Y7 c   2 A+ Y* A. G) H" u2 k
       sum ddy240
    " ]7 J  A5 A% x1 Q+ q; }2 U2 k   $ b5 b# S* u' i# t1 ?
       if r(max)>0 {
    0 y3 C. l. w, I: \                9 i1 q' Y, i3 R! {. W/ D) f
                    keep if ddy240>=0' j  O! I/ B1 d5 ^$ @
                    ! m' q( N* ~' L# e8 w; C' I
                    count
    , v; @& M9 g( K; T                3 p* Z; o! G+ }8 Y
                    if r(N)>2 {, z# y0 a6 o7 I( F; k8 e8 W
                     
    2 {( p& u; [3 J- }  E$ \; s7 F                 pkexamine dx240 dy240
    ( |6 n- d1 t8 j                + k! |6 T- b! J, l
                    local tempmaxp=r(tomc) 4 g- O( E2 L' T' D- T/ R
                    restore 0 ]( a0 F0 P" b- n
                     }
    7 Z: w8 \* e9 G5 [7 @, t                 
    ) J- b3 @' P; F' t3 ~- U9 f                 , k( A& z' ?' c" z0 A: E
            else {* c$ c1 L5 d4 H; d! r
                    restore 3 {6 I5 I7 S* `4 h! A0 Q% a
                     di "cannot find maxp"
    * Z4 s. |- j8 b               
    5 B. w" D' @. }  Z3 l  H               
    * v$ `& v/ ^' j                 sum dx2403 C6 X" W) D0 U" x' P. w/ K
                     local tempmaxp=r(mean)
    9 c4 K' N6 L. p8 q# x; V+ e                 " B- H: }6 a5 L: z. u
                     }
    , [" i! P: T; `; R                 ' M: P0 |3 v& L
                     }/ w/ f% N! w; i
                     6 I  ]. O$ x9 X6 k  a8 V
        else {) w3 N, q" \1 [* P% s0 x
            6 @7 q" E+ W9 S7 Z8 n
               keep if ddy240>=r(max)
    - _  Q+ f1 Q% x- r5 a           ( o5 D9 a5 p8 G, ~* |
               sum dx240
    . F: B; U8 @# M; H' F/ g          
    ! @/ i2 K  @/ f; r5 |- V           local tempmaxp=r(max) / I. F! i: O& s. Y6 u2 D: `1 h
               restore . p7 m, |- u9 R6 M
               }- S. x8 W6 }0 U8 V- k
            , d3 U! t7 G# U& K  T! Q6 r

    ) ~4 w- U, u1 R( ~; S& B7 G6 \  
    ' Z' V' k9 ?% D; J
    9 D6 U4 n5 n( \( j   gen maxp=`tempmaxp' in 1+ C) T: m! _( _9 i2 S" V/ I
       - @/ u/ i/ Z% f1 B8 U
       *240maxp
    4 V1 g2 T3 b3 N" P8 [5 \+ y1 S   , H3 V2 i1 y; y9 n- q! G$ M, E9 B
       
    . }+ Q3 b9 @) }( ]# Q   *-------------ddy240pminp左侧高峰值-----
      r  a$ o+ H$ S0 G1 }+ [- r3 K   : P% j1 f* c6 m: Y
        preserve
    " s: i# v4 H5 I) |" A9 @& b   
    " M/ X( z2 j4 \: Z+ P; Q: |   keep if dx240<`tempminp'
    0 U. p7 u& P; {- b2 ~, f   
    " q; L+ D* K: y3 n   keep if ddy240>0% b, u6 I& r; Q- h
       
    / J% c1 @: Z% `3 k: G( q4 o  G   count3 }5 R+ s, q/ j) f
                   
    2 g8 Z: {7 A5 u( ~. x                if r(N)>2 {, S2 \5 B6 |, D+ k
                     
    ! y+ E8 n( q* A3 ]& c7 O                pkexamine dx240 ddy240
    ( B, j' e  u3 G& r# G# V; h! m                . ?+ S0 }& i: x
                    local tempminphddy240=r(tomc) 6 h& F7 X9 C# N8 U1 U  O* {3 `& G
                    restore
    1 ^; L: w" X; {                 }
    # F. f4 N1 V( S: y8 M                 % a6 I( G9 l5 l* W0 n: K
            else {/ d( u3 ]: o/ r5 e6 [
                    restore & E, ^; W+ c0 |, j
                     di "cannot find pminphddy240"
      q3 O2 H- S& W2 H                 2 j, @- D4 I, Y8 `5 q
                    , D/ n8 ?4 t2 ?( t! A0 q7 e
                     sum dx240
    * W2 U' S0 u- q0 f: `% A5 p                 local tempminphddy240=r(mean) 7 I, R* I: y# Z9 U. K& O
                   
    : e; h, b2 L4 c, a! c                 ) n0 D/ B' [- r/ E0 m0 ]: P2 w
                     }
    8 n3 f- C8 \* W1 Y+ Q( D' i4 T                 
    6 B9 @9 L! g9 c- f/ D   1 v, r/ p& j2 c! l
        & G6 ]% J4 p- f) f
       
    9 L/ U! F6 K4 ~9 A" x9 Q1 ]   
    ; a6 E8 p6 r, J   # v6 U. J. U, t6 O! ?, j" z
            gen pminphddy240=`tempminphddy240' in 15 t: ?# ~+ k4 A* m5 V1 G
           
    / D) d* v8 K  W0 ?% z/ n8 t       
    8 f/ _( [7 }/ s$ L0 l        *-------------ddy240pminp右侧低峰值 -----------------
    . N/ w4 L# q3 r/ Q( ]% t, z    preserve
    & M8 ?3 |  E8 u  Q% N
    ( h' Z# r2 w# P9 K$ {4 V& c. Z   keep if dx240>`tempminp'
    4 q, S6 t3 r8 h   , P) t# B9 p$ P3 o5 Y) B
       keep if dx240<`tempmaxp'# e2 v: H: x1 a* D
       + n( p/ a2 I$ \/ c8 m/ l
       keep if ddy240<0
    ' j* u6 |' i' E6 @+ }5 B% u, o& c   
    # u1 c2 h1 t. J* }1 i9 I9 Y+ ?7 H0 v   gen temddy240=-ddy240" E0 v7 y9 P1 S. j2 n
       
    # I& U  v/ G! D4 u/ P9 k5 q& n    count7 w. b& z- Y6 J  F; m
                    ' N& W9 S3 b  a' B* l8 Q
                    if r(N)>2 {
    . z9 u) g. i7 ]  k# w, D) C- H5 ]                 
    # c" b/ p% A; S& O6 ~; o                 pkexamine dx240 temddy240" [& w) K3 T1 M+ h1 G2 @, i
       - U: j# J2 v' J' N3 z5 }
      local tempminplddy240=r(tomc)
    4 j7 A9 E2 P' Q                restore
    ) l0 s( n2 v, w. p                 }
    9 x6 M$ ^2 }9 B  I                 ' ~& R6 F+ Y3 r, \
            else {
    8 G9 k, L! n5 s' V0 B4 d6 \                2 f4 g% j& r0 o. \
                     di "cannot find pminplddy240"
    & A( [; }* N; s  `8 s                restore& p8 `  n7 J6 ~, Q
                     
    - k  }  Q  M- [% P2 h                  sum dx240
    $ f' |! p1 u6 M& r- R( q                 local tempminplddy240=r(mean)2 @# i& [& M% E$ C  O& R
                     
    8 T2 z: o1 A$ X/ o4 U" H$ ]$ w                 
    6 B! E3 P) B# J0 a+ K                 }
    % k& G. d5 H' R& J8 ^   
    7 K: k) s% f. ^2 F; H    / q1 D7 Z# R4 v: @; s
       # Q( K6 C, X8 ~5 p- s" l$ |7 M
         }% O9 P" l- b5 X' O% V0 T* f
       
    . u0 F  b3 X1 ?9 N& ~- r/ [% Y        gen pminplddy240=`tempminplddy240' in 1$ V* g9 A9 I5 h
            / T1 M" _% Z( d% w+ t& g' g/ E
            *-------------ddy240pmaxp左侧高峰值 -------------------
    1 F* d2 q9 Y0 _$ c: h        0 R3 c7 X! m9 J( S
             preserve
    . d* u. ^4 `$ I   7 n' R% m7 O: T+ y- o0 i* i; ?
       keep if dx240>`tempminp'4 S; u$ S: K: Y7 [' C
       / j, U! F4 M5 _. N. ~
       keep if dx240<`tempmaxp'
    6 V/ s/ x1 w8 Q& O( l3 C+ O   8 B( \4 o( T* r) k+ F2 N7 F2 U
       keep if ddy240>0! t& ?& ^9 M1 b$ R$ X- r
       count3 u7 {8 v$ X+ K
                   
    $ q, B8 s0 J8 N  a: }5 W                if r(N)>2 {$ R4 E' \% o4 G9 d3 I
                     
    4 K* p/ \! K/ F' }. g+ n& t9 G1 z                 pkexamine dx240 ddy240
    9 w$ y* c; ~" X4 |+ |   
    8 e8 q- V6 O" z- x- \" i   local tempmaxphddy240=r(tomc)
    2 k% D4 h, v, U5 y   * F. {/ E1 \0 j: W8 P6 b
       restore
    8 |/ s7 C0 a' N( O4 w5 K                 }
    1 b" Z- H! Q# J                 6 Q1 d& [6 G/ j8 c; R+ `
            else {; c, J! a# R0 [0 D# F
                   
    + M* P$ \! ^0 y2 V7 _( n" g                 di "cannot find pmaxphddy240"
    ) A' q* ^7 j, B+ k( v                 restore
    & X2 A9 b0 d& @. M6 P% r               
    : T' V( j% v# l9 J# [* m2 A0 }, a6 V& n                 sum dx240
    3 J" r- [8 s2 u/ W# @                 local tempmaxphddy240=r(mean)
    5 x$ W4 E1 A+ s+ V3 b8 c! B4 a) P2 x                 * E" a: \# _1 M3 [5 B# m) c( ?
                     
    . J+ c8 F* i+ s4 L8 u                 }
    1 {7 p5 r/ u2 I   2 f$ J- c8 I# j7 O
       " z: p" B# N- ?
       
      U+ h2 V  }7 {: C8 H% p. c   + t( F2 V9 c: m* W/ ]' p
       ! N4 x! p0 a. J# v/ r
       
    , T8 k* [7 I0 M/ w+ B0 `3 t        gen pmaxphddy240=`tempmaxphddy240' in 1
    ( I  k+ n% b$ h       
    " ]7 g: t. M) }! B        2 N: |" O* Z, S  w4 k
            *-------------ddy240pmaxp右侧低峰值 -------------------
    , W* m. X- E; P& i# @9 b        # y7 v" V9 f1 ^( V+ f) b0 }8 H
             preserve
    & |0 ~! W  s( I0 I$ c) i+ @! W9 Z6 X  
    ) \2 P8 Z3 B/ Z  V; V5 Y+ B: ~   keep if dx240>`tempmaxp'3 U. k0 n% w, q3 L
       ' k6 H- Q% @6 c) z
       gen temddy240=-ddy2405 ]  \' G8 `' @2 a( x/ U
       ( b5 }! X5 g5 y' X* P0 {% J6 D4 T
       keep if temddy240>0
    0 f6 A3 F6 W. V/ H7 b" _   
    - L1 [- H7 d! |5 L1 @+ h9 S   count
    7 K: F* Q  y: v               
    8 d' B, E7 ~' R. r  ^. b/ J                if r(N)>2 {
    / X2 i4 I$ c/ Q( X6 q! }  l' K                 
    # m, C& Z' d, q2 v& h9 U$ l                pkexamine dx240 temddy240
    ' O7 \0 [( x$ ^; u, \   
    ! g1 a1 D, |' t, {   local tempmaxplddy240=r(tomc)
    ) R" G4 W7 j+ G; r4 W   restore1 F6 Y  x- r  @
                     }
    - g" q/ T- q+ ]4 H7 u$ J                 
    7 Q5 ?. \4 d6 V* v! W        else {
    8 w  m( M8 w9 O9 a- {                restore
    / H2 V% R' _: A) h                 di "cannot find pmaxlhddy240"
      _2 I( G7 O7 ~1 z- l                 
    + J' K& M  I! Z# E5 y8 b# |               
    3 f5 I  H3 X+ t+ P  N7 i. h+ P. R$ N                           sum dx240
    - @4 F4 q/ z- M0 p5 M% H; @7 Z                 local tempmaxplddy240=r(mean) in 1* e8 ?% k- E8 j# \2 W
                    - O/ n* g0 V+ R# E
                     ' [3 d( R; q5 W; B  \
                     }
    $ M3 G, w1 ]  r: E   : t3 d: g+ v5 s1 H1 b! l, @9 k
       % P+ z8 U2 o0 u; W# q0 n& @9 H
       
    8 g& l5 {- W4 |4 ^+ a2 ~, z' W   / C* K; g$ }3 ~
       
    . m5 j$ c  P0 x* W) t5 i   
    + `( _9 u% l- _        gen pmaxplddy240=`tempmaxphddy240' in 14 H6 n2 W# x$ [5 X" n. W
           
    7 G2 e7 u% G1 x1 @/ S        *-------------------------------------: j* J8 C& `; f$ q& W
       8 v4 x3 ]8 O0 T8 T
       gen price5d=v5 in 1/209 G4 y8 K( O! @' H* p. W) f
          
    , A9 [. B& j* U   gen price5yy=_n in 1/20: }: g5 ?, C3 e: Z
       ( _' L. D4 a2 O; |& _6 H/ W5 G; n
       gen price5y=price5yy*(dy240max-dy240min)/20# C! o8 I8 W6 h& y1 N/ @9 T
       
    3 \4 r) D; B7 M, q; \   *price5y &price5d/ L5 g3 X4 D5 M  Q1 T3 E
       
    ( g( `2 |% m7 y" ^- |   gen priceny=price5y in 1$ I7 P" u+ B* i6 F7 @
       ; v$ J; O1 [2 Y
       gen pricen=v5 in 1( _/ z# z: e. `- u0 c
    ) E6 z8 z* `* a
         3 I! C7 H) O: T2 p# Y& l, ~. h9 t
       *priceny &pricen
    ! y% [6 t- T5 P/ ]: w   
    + Q0 }' L$ U) x8 N8 v1 o# \   
    $ c: j" K  q) y6 U4 o   *-------------------------graph------------------
    # x  u. w; }( e/ i, s0 s  ' o- d! F3 G* Z4 @: H" q, R# b
      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)
    $ V7 {7 c+ `1 E& J' i: \9 J/ @6 g6 v
    * ---------可选项目----------------------
    - ^- U- [1 Q2 a) _*(spike yxis amean,lc(blue))(scatter yxis2 amean, lc(blue) msymbol(+) mlabel(amean))" ~( ~4 A3 I# H3 r) J
    }- ]5 x- v. p, `( u# U
    end
    - U* p3 p$ ^2 N! ~& g  o. b

    该用户从未签到

    69#
    发表于 2013-2-7 00:03:44 | 只看该作者
    awer 发表于 2013-2-5 14:40
    1 y% e+ M8 B; V! @抱歉,最近太忙了
    0 {% s' ]2 _' z1 C! z- j0 j7 L$ m' u
    300236 我用下面的程序试过,没问题
    $ D% i; a# n8 m8 r4 @
    哇,谢谢,我仔细看看,这些字母数字代码晕了, D& q6 I6 n* u8 d* z. q

    该用户从未签到

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

    点评

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

    该用户从未签到

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

    / [& u$ _+ @! D* h; ~& q( L+ O( u* N3 u0 x7 U
    这东西是开源的?知道有个统计软件叫stata,不过要收钱。

    点评

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

    [LV.3]辟谷

    72#
    发表于 2013-4-17 17:14:44 | 只看该作者
    准备试一试
    % z% n8 ^: n( B+ ?/ B- v$ c先找Stata软件
  • TA的每日心情
    郁闷
    2018-10-21 12:19
  • 签到天数: 7 天

    [LV.3]辟谷

    73#
    发表于 2013-4-17 22:21:02 | 只看该作者
    本帖最后由 汉水云潇 于 2013-4-17 22:27 编辑
    ) O1 K" `; p$ W* t3 ]6 V
    7 w) B" q+ k+ f# u! m8 q晕呀7 d4 n+ E; v  s# f
    论坛复制程序代码
    # B' b0 Z$ m1 }* ~. z' A1 S乱码夹杂在其中
    " H, r" x3 h$ H! A: a. |/ |很烦
    5 p& B/ @2 M9 _. m9 k4 M对了两次
    ' i) Q" i, m" p以为没问题了. w4 ]( z0 q6 t* B
    + v8 J# M3 s5 t- f
    结果9 y5 B, t2 W. I$ A3 [$ W
    执行后
    + V$ a3 t7 y2 I- m还是出现
    + p& p0 k6 i/ [# E2 b+ a* b' m* _r 199
    * k5 I& ]( q! r- H6 N的错误3 P/ t6 M  B$ c7 K( O( r6 c4 `
  • TA的每日心情
    郁闷
    2018-10-21 12:19
  • 签到天数: 7 天

    [LV.3]辟谷

    74#
    发表于 2013-4-18 07:50:56 | 只看该作者
    试验成功,下面是三一重工的图$ e/ _0 R. T) m8 @$ Q
    5 W7 ^9 ~" c6 |0 w4 O  r* f

    点评

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

    [LV.1]炼气

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

    该用户从未签到

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

    [LV.6]出窍

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

    [LV.1]炼气

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

    使用道具 举报

    该用户从未签到

    79#
    发表于 2016-3-23 10:23:05 | 只看该作者
    明月回春 发表于 2011-8-13 20:17
    2 t" V/ n: t/ k3 j9 G/ r! d程序写的不错.* k0 `! r' H0 Z4 `9 g8 x2 e1 L
    不过有一点要注意.通达信导出的数据,只对价格复权,并没有对成交量复权.所以数据需要手工处理 ...
    1 {, d% \/ m5 ^- F$ J
    老哥,也是山寨这篇文章找不到了,不知为何?

    点评

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

    [LV.9]渡劫

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

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

    GMT+8, 2026-1-12 08:06 , Processed in 0.049799 second(s), 27 queries , Gzip On.

    Powered by Discuz! X3.2

    © 2001-2013 Comsenz Inc.

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