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

maodeng

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

maodeng

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

maodeng

. hello 300236
4 v2 N9 U& a8 B& o  U3 Ustatement out of context
r(119);

& B3 _) j) T5 X大师 就是这个

sergio

这贴实在！

awer

maodeng 发表于 2013-1-28 13:13
1 o) k% ]6 e7 Q5 k, Y. hello 300236
& Q" b0 H, n# V9 \statement out of context
; A5 ~7 g2 {# G4 g# ]2 n/ {r(119);

抱歉，最近太忙了
5 e$ R/ b  @+ @& F( J# ^
$ t8 t! _# V! R. ~300236 我用下面的程序试过，没问题

你把 D:\stock\test\source 改成你导出的数据存放的目录，再试一试

1 l) S& F- g4 Z+ E& Ccapture program drop hello
3 b8 `3 O4 ]1 r7 S8 R  v1 q
program hello
; v+ `$ b. V1 I, z
1 e( W5 M* Z7 k- z! Z0 w! u quietly {

" G4 I% G2 t* Y3 v  J* d1 eclear
0 }% @  W% |$ @. N3 @# A- S( n- f
cd D:\stock\test\source

7 {& z  {4 g  ]- W4 f1 ainsheet using `1'.txt
8 z& V0 k. c* L+ o% w2 |7 S
   keep if v5>0
# l0 H1 l7 R% d) p% w6 x   
: h& K4 n: y  D6 \' `; {: q   gen vt=_n
   
   tsset vt

, B- W. w5 T; y3 a/ R   gen zf=D1.v5
$ n( u6 K8 U( ~* |
/ G; P7 b* ~1 _. s$ f8 \; M   gsort -v1

7 ?* |* S' k! F, n% q( o   gen avol=v7/v5
/ H$ t5 |* O1 X9 _! ~   
  ^+ c9 b: [9 ]) w; L$ F: W; a   keep if avol>0
   
   drop vt
8 Z+ ?5 o2 G9 Y& v3 ?  
   gen vt=_n
   
   tsset vt
5 h: h: x6 S, s: o   
7 w3 j- k: P" H, ?0 ?$ Y    sum v5 [aw =avol]
! K% d0 F6 ?' H& Z& E5 D# g/ E
9 B8 p: M) z. V1 k   gen amean=r(mean) in 1
' l$ }" {5 r+ O; M   
   * 总加权平均值 amean
   
   
2 p# q/ {( [! V3 R   keep in 1/240
. _4 H. M& j; K, r/ k4 f* Y4 M
   sum v5 [aw =avol]
+ J% S9 ?# E$ U: e
9 a+ L) R& {4 W   gen mean240=r(mean) in 1
   
; Q2 q& u+ z$ ~4 {# j+ u7 A   *240加权平均值 mean240
, G, Y9 W) [9 S0 z9 s* ^   
& k# u: b( n# `6 z   set obs 500
; z. T3 ?# h! u% i5 W   
* B/ Q. J* e; Y% t4 P6 u3 {7 N   kdensity v5 [aw =avol],gaussian n(500) generate(dx240 dy240) nograph
   
   *240 kdensity(dx240 dy240)
   
& b2 n3 k0 V4 C: O   sum dy240
   
5 [& P+ n/ l0 a1 w# N( c# ?1 _  g   gen yxis=r(max) in 1
   
   gen yxis2=yxis/2
   
8 S4 G- m. ]8 s& j! g7 T   *yxis
   
7 l: f* [6 z7 p1 \9 S' N& E% R   gen dy240max=r(max)
   
* i4 Q1 y- |# d: e( k( Z$ k# n& t0 \* l0 l, a   gen dy240minp=r(min)

% W7 q$ F! `" M) h$ U    gen vtt=_n
6 P8 z9 o$ V" O
& ]9 Y6 L" D" e( t, c5 _; x1 N1 p    tsset vtt

( z( F7 F( ^' [9 q, e3 g' e- w7 j   gen ddy240=D1.dy240/D1.dx240
   
7 K+ b' m$ a/ U1 ~+ ]: P   sum ddy240
   
# a$ B: V: ~: f* u$ ~8 a8 }' _   gen ddy240max=r(max)
& c3 P/ A. {* T$ U8 z   
' ^% a! G! _, V. k/ r   gen ddy240min=r(min)
* i4 q8 u1 i( [7 i     
   gen ddy240g=ddy240*(dy240max-dy240min)/(ddy240max-ddy240min)
   
5 u4 v3 [+ K  c- M0 W9 X   sum ddy240g
   
   gen yxis3=r(min) in 1
$ n7 M0 V6 ~; |5 }" A+ d   
   gen yxis4=0 in 1
  
  *ddy240g
   
   
  preserve
) P. u# h$ I# e  
  sum v5 [aw =avol]
5 ]9 e( e4 M$ v  
  keep if dx240<r(mean)
$ b. l* _0 ~/ n) k/ f  
  sum ddy240
  
6 ~3 E2 w# X9 U) y" z* e2 p) F  if r(min)<0 {
        
                keep if ddy240<=0
2 u5 j7 L6 e  \               
& Z0 G9 d6 A* |( V. T6 |. P$ W; ]                count
1 U3 a! s( i! o8 j               
! P) R/ f: _  L1 _* z                if r(N)>2 {
                 
9 [+ O; Y$ b( e2 ~- X: K0 _                 pkexamine dx240 dy240
               
                local tempminp=r(tomc)
                restore
                 }
0 h2 V: i7 b( F3 r, U                 
        else {
               
                 di "cannot find minp"
. b* [: f  H" Y                 
                restore
                 sum dx240
" O  q& D" B9 c                 local tempminp=r(mean)
                 
                 }
4 p4 N, w' f- j                }
   else {
                 
0 Q3 q! a$ H0 ?5 \                 keep if ddy240<=r(min)
# ^) N8 r# [% H1 D                 
1 M2 w! C: ~7 \2 G$ w                 sum dx240
                 
; j, @) r8 Q) F1 F& Q                 local tempminp=r(min)
/ A9 T8 N4 V9 \( \) v# U. c$ {                 restore
                 }

  G' u- g4 E, |. H3 k: o
: ?4 I0 q- k, A. N
7 l9 w) v2 Z9 }* }  A   gen minp=`tempminp' in 1
# @; S% w& t1 j( |" c  g% K- O' k6 u   *240 minp

9 C/ }/ c8 p7 v* Y
, [3 p4 I/ f6 F) P, O   
) w5 E9 s, L$ J9 Q/ R   preserve
   
7 W5 F: {) g( y$ M9 n   sum v5 [aw =avol]
; U  ]- n) d% Y" Y, J5 v' J5 c- d& ]   
# e; E4 G" S. `   keep if dx240>r(mean)
   
   sum ddy240
! \" b' O; v3 W4 L# v8 l   
4 A$ a6 N; j: i   if r(max)>0 {
               
                keep if ddy240>=0
               
                count
               
                if r(N)>2 {
                 
                 pkexamine dx240 dy240
               
                local tempmaxp=r(tomc)
                restore
                 }
4 j) [( _1 R# q9 k# Y( n1 R                 
2 E/ n- [! j' L( R& R8 E5 Y                 
0 h* A& w" O6 p  u        else {
- Q9 w: \* b' g1 b% ^) [8 ^: ^                restore
                 di "cannot find maxp"
) f" K# S7 p3 s! G$ B1 @& W               
               
                 sum dx240
3 W* r' \& J& N3 F' q0 C! n5 t1 F                 local tempmaxp=r(mean)
$ l3 _- e& Y! j& h. L4 I$ v- O                 
$ t' T# _8 B/ p5 T& f                 }
) O3 d5 y9 U' ?8 Y7 y3 C                 
                 }
8 y8 I7 U% `* C& k8 y: Y  V2 j                 
$ o. Q( Y" a, ~6 `* ~: i    else {
: v6 n) ~; N. ^. P1 \. u/ u       
  D2 ?$ }7 H; W; O" K           keep if ddy240>=r(max)
7 O. x9 `& [, S" S; D& g  s          
( ^7 r) Y% u3 a9 \  _           sum dx240
          
+ W& }% b) M) x! t! I' C/ [. V           local tempmaxp=r(max)
           restore
           }
       
# Z/ f9 \$ e$ |1 N
  

, Q3 |( b$ }; g& c8 U" [   gen maxp=`tempmaxp' in 1
, n+ J& b" C. ]* E   
- [$ q7 t1 T3 m" S1 S5 p6 t, X, {   *240maxp
9 L9 T: e0 G5 n  ?   
   
  o' i/ ]$ {- y) x0 A   *-------------ddy240pminp左侧高峰值-----
   
% Y- w1 [6 R7 X! ^5 }8 R2 |; v/ n5 E    preserve
1 v0 k+ }  j/ R5 `( i/ I0 `   
* ~5 T+ y6 p5 J3 H. l   keep if dx240<`tempminp'
/ t* F2 B% y4 D) f& q. Q4 y   
   keep if ddy240>0
   
/ l7 I$ U0 I) N   count
: \# ]- B2 Z  f$ _; A6 f! a7 d               
                if r(N)>2 {
, @' x/ p* L6 |/ p- l* W                 
                pkexamine dx240 ddy240
               
                local tempminphddy240=r(tomc)
! [" S/ h3 G1 a7 M2 z* \                restore
+ u+ [2 [: u7 \                 }
                 
        else {
                restore
; S' r  m, [* @                 di "cannot find pminphddy240"
                 
               
                 sum dx240
( W9 S- D! s* ]5 X7 q- ^                 local tempminphddy240=r(mean)
: b6 b- T* `3 L8 |               
                 
                 }
                 
! w- z, {! ?' {   
   
! h7 V- ]$ n% W. s7 K$ M   
: i5 h+ K7 l7 y) y. z   
0 F8 c9 T7 k) ]   
        gen pminphddy240=`tempminphddy240' in 1
       
: }  w1 t+ N3 R( ?  n! e9 y       
( g% s2 }9 L: T  `" V* B: Z        *-------------ddy240pminp右侧低峰值 -----------------
    preserve
$ N- p" r' Q' x/ N; a4 S/ ~
5 Y- f5 P) G9 w1 d; }9 M   keep if dx240>`tempminp'
   
   keep if dx240<`tempmaxp'
* x' ^3 ~; k' |- x   
   keep if ddy240<0
   
, r5 b& \5 s3 G. J, G: a3 a+ y, e+ q   gen temddy240=-ddy240
- b, q0 u" u- ?   
, o0 c' B5 ^, g% s4 s% N  q    count
. N. B  d) F4 y4 C& b: r               
                if r(N)>2 {
                 
8 Q2 q2 K" t( |" g$ m- [                 pkexamine dx240 temddy240
   
  local tempminplddy240=r(tomc)
' x% q0 r* M0 o                restore
  R1 l. ^4 m5 l4 F                 }
                 
: o0 h. T% X! ]2 V: N" {5 A1 C' C        else {
: A6 `( g2 K  P               
                 di "cannot find pminplddy240"
) K) a/ @4 e( i1 o* H( L/ s) i8 L% s                restore
                 
& `1 d, g/ K3 {+ }                  sum dx240
1 k0 B1 A$ |( d# K  C                 local tempminplddy240=r(mean)
                 
                 
( p$ ^! \" ?7 L                 }
) j' I, P& L( Q( J3 D6 D   
& z( C# Y7 q/ o* }( l   
   
1 Z* W4 F# U& _6 G   
1 i$ T8 K: [5 P+ ^   
: L+ @0 o7 N& y        gen pminplddy240=`tempminplddy240' in 1
( B9 z' p& P5 Y8 t- Y       
        *-------------ddy240pmaxp左侧高峰值 -------------------
+ T, n/ N" S3 @6 }- g% x       
+ B" A3 }1 l( T; r         preserve
1 E: n2 C" f; Z. U' j- |   
+ S4 U8 K' D$ |' s# B7 W   keep if dx240>`tempminp'
   
% v  T$ S+ b" k/ N, ]   keep if dx240<`tempmaxp'
( N  g4 d- s2 n0 b( H4 t, d1 z   
9 g4 l# ^. J7 O8 V( s   keep if ddy240>0
   count
; I+ ^2 }2 p" z9 K) q) P               
                if r(N)>2 {
  e- Q, A( Y0 {) i( i                 
                 pkexamine dx240 ddy240
/ x# j$ G" f3 S2 Y   
. q, }( K. I6 u! n/ v0 h, M% a   local tempmaxphddy240=r(tomc)
7 C9 v$ h* G2 d/ H) ?   
   restore
                 }
                 
' h+ w8 F9 M( z8 Y5 x        else {
               
                 di "cannot find pmaxphddy240"
                 restore
" K2 A4 |! B& s6 F               
                 sum dx240
# A/ |5 X2 _% p( X                 local tempmaxphddy240=r(mean)
+ ~, z; E8 b, \+ K8 {/ E                 
                 
7 y" _6 j/ `- {" T1 y                 }
   
/ o4 ~' F  i( {+ Z( r# |   
   
5 A+ b  {4 ?* \4 Q2 \, _   
   
   
        gen pmaxphddy240=`tempmaxphddy240' in 1
" \0 U0 N; {) m! D: z8 V       
. O  `8 d& p) }; r2 v       
        *-------------ddy240pmaxp右侧低峰值 -------------------
       
         preserve
; A9 [' K1 {& r- J( i9 N3 j  
* v9 u2 E" S% _% k) f2 l' M   keep if dx240>`tempmaxp'
# \# p, Q/ n: V' G   
   gen temddy240=-ddy240
   
$ Q- N- l* b" b& r( p   keep if temddy240>0
0 O, q/ Y( O1 D0 _. h   
* o# d: N( r- {- n  o+ E4 ]   count
               
* \  x$ `3 w* a$ v4 i4 j                if r(N)>2 {
$ t8 C4 `. \5 k8 U" M+ k                 
+ O" c' t. Q4 T5 J+ d+ X                pkexamine dx240 temddy240
9 E( c& ~6 E( k7 U& A   
   local tempmaxplddy240=r(tomc)
7 I  E  Z% c9 k! L; _   restore
                 }
) P& p: ~* j5 a" {2 S3 ~2 O6 E6 P. h                 
( `/ R4 k) f0 @, L# W+ k        else {
% P4 o7 b, s0 Q0 @( P$ r3 I                restore
7 H# G4 h" \; h) O7 [- p                 di "cannot find pmaxlhddy240"
                 
/ W0 a% U3 S, B) d+ `- [( n4 V! E) j               
                           sum dx240
                 local tempmaxplddy240=r(mean) in 1
0 g8 l0 L, ~5 B+ t               
* |9 I% q3 o! A! U                 
: D6 _! b( o) y                 }
   
6 L; R+ w" ^% {" x' V   
; @3 g/ P6 i3 V, I! F   
   
   
   
        gen pmaxplddy240=`tempmaxphddy240' in 1
" \5 l/ z) ]! |       
        *-------------------------------------
   
( C  l  W* V( N7 a2 Q+ X   gen price5d=v5 in 1/20
      
: A% ]) t/ d/ s$ D0 f# C1 z   gen price5yy=_n in 1/20
, X8 v8 k& {/ }   
   gen price5y=price5yy*(dy240max-dy240min)/20
% i& O. [7 g! B; h7 w   
& Z+ O5 K/ A0 ^: h   *price5y &price5d
/ I, b) N4 [7 u   
   gen priceny=price5y in 1
   
3 z) T! d' J& r" X- ]   gen pricen=v5 in 1
2 [4 I, _1 c  B' \
     
   *priceny &pricen
   
/ u. z7 s- k7 j& o8 B4 H* w   
   *-------------------------graph------------------
5 v9 S# R9 n' ^  a1 t! s. d) D7 i) o  
  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 L7 d9 e, y* ^( M*(spike yxis amean,lc(blue))(scatter yxis2 amean, lc(blue) msymbol(+) mlabel(amean))
; e1 d3 n! W& v; j* K4 @}
# a; e% _7 L! A2 Q4 r& x! Jend

maodeng

awer 发表于 2013-2-5 14:40
抱歉，最近太忙了

5 H( B% y9 e' U7 d300236 我用下面的程序试过，没问题

, Y8 o) P4 r0 Z2 l: J  T3 n哇，谢谢，我仔细看看，这些字母数字代码晕了

喜欢喝冰茶

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

喜欢喝冰茶

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

& A  p, u% O: z8 X5 Y5 E! L7 U# k这东西是开源的？知道有个统计软件叫stata，不过要收钱。

汉水云潇

准备试一试
+ h, W) t+ i8 N5 q7 F先找Stata软件

汉水云潇

/ Z/ |& g0 V8 r4 J4 U- w晕呀
论坛复制程序代码
乱码夹杂在其中
很烦
- {  c: g/ \8 t5 Z9 H( e对了两次
以为没问题了

结果
" t9 L6 n: P6 g执行后
还是出现
8 h6 i% a8 d# T5 |6 Or 199
8 l8 T" @3 n" L: M( C0 G的错误

汉水云潇

试验成功，下面是三一重工的图
& D& k0 S+ l) o6 l) x0 i5 d1 @4 f) i
/ K3 m" m& u/ K& c( L2 @

peterbev

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

中秋不圆

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

bknite

这个帖子值得学习。

jeanny

感谢！！！

Hubert2000

明月回春 发表于 2011-8-13 20:17
程序写的不错.
; C' p$ ]% R. H' K不过有一点要注意.通达信导出的数据,只对价格复权,并没有对成交量复权.所以数据需要手工处理 ...

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