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

maodeng

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

maodeng

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

maodeng

. hello 300236
& r+ Y7 }2 v' I0 u# h4 X5 k7 \statement out of context
5 |# l$ a  z! ?2 i& I1 Nr(119);

% R# y/ }+ H( A* j4 H/ Z* m大师 就是这个

sergio

这贴实在！

awer

maodeng 发表于 2013-1-28 13:13
: I2 _( J. |2 N. hello 300236
statement out of context
: ]# E8 k, V- A9 K6 ~* Pr(119);

3 Z8 B6 G& D  [  v/ p9 t抱歉，最近太忙了

300236 我用下面的程序试过，没问题

! S. w7 F. r1 S你把 D:\stock\test\source 改成你导出的数据存放的目录，再试一试
- [$ ]) i% K9 m- m: {& w; a& g
$ [- [; a' F) `) Ecapture program drop hello
3 W7 y7 O% h7 d/ v  W
& @9 f5 k3 p0 ~& ^$ F4 s, K0 hprogram hello

- M/ |' ~8 x  k8 W" I% s quietly {
0 a. f* w6 G# `, y: o+ \
clear

cd D:\stock\test\source
* }: N' u0 V( y9 T; X+ q
insheet using `1'.txt

2 q  S* V0 r5 Y4 o# Q1 q; V" [   keep if v5>0
2 L. V4 p4 k6 H% R' k+ a) w7 v2 E   
. y( J2 H9 Q: `5 j- E   gen vt=_n
. g0 S* P- S% K* Y0 s: d   
( y' R0 @' V4 `( ~; g! I( [   tsset vt

; C$ O1 i. X1 p5 R2 N6 v2 y   gen zf=D1.v5
0 ~; a. E9 a4 e2 Y3 Y* r) Y* g; w
: |9 b# e6 p# m) t$ X   gsort -v1

   gen avol=v7/v5
0 E+ E/ D' \) o3 d6 {0 i( B# t   
   keep if avol>0
- k4 h" n8 ^+ D: w% F* ]1 v. T0 q; ]   
   drop vt
- H! m$ ]; I9 Q1 e  
   gen vt=_n
( o$ w0 e2 E  d+ |* u   
. b, F. f# k1 p3 N- o   tsset vt
   
3 t- B$ L8 R, k  o+ l& N    sum v5 [aw =avol]

- ?. i6 F* H% w. \   gen amean=r(mean) in 1
   
; h) P0 u  t% @; o& c   * 总加权平均值 amean
; t4 l# v8 z' @' W3 o   
   
$ d# ^3 G+ J4 ^   keep in 1/240
- S$ m6 v1 ^: i8 `3 B
   sum v5 [aw =avol]

   gen mean240=r(mean) in 1
   
   *240加权平均值 mean240
6 f- M. v. @/ B. m$ J1 X   
   set obs 500
2 n, T( g5 F6 `9 e" S" g. W   
  T- c# d/ x" e3 C; O   kdensity v5 [aw =avol],gaussian n(500) generate(dx240 dy240) nograph
   
   *240 kdensity(dx240 dy240)
   
   sum dy240
   
   gen yxis=r(max) in 1
8 m  H5 N' [$ a6 E* O   
   gen yxis2=yxis/2
   
0 u2 r- z9 [8 t5 ]1 d2 [+ }   *yxis
% p' x) x' W) L   
- |3 d7 J. e' x, O  t5 S7 r$ j0 x   gen dy240max=r(max)
7 `( Z, w% k; F' ]" x/ U! V   
   gen dy240minp=r(min)

* t$ y3 C* n- C, X    gen vtt=_n

, b! k+ l$ P4 s- R    tsset vtt
8 C+ Q$ g' m" B- H! \+ S7 _9 o
   gen ddy240=D1.dy240/D1.dx240
   
8 Z  @, a9 {% _! |; F   sum ddy240
3 c& p( ?$ I% T0 H+ k# A7 d   
5 F) R, @. a/ }0 L: [$ y" O   gen ddy240max=r(max)
   
   gen ddy240min=r(min)
* j6 w# r+ O3 R( L8 h4 z1 P" W     
% ^! K( b$ |. X$ p   gen ddy240g=ddy240*(dy240max-dy240min)/(ddy240max-ddy240min)
   
   sum ddy240g
   
   gen yxis3=r(min) in 1
   
   gen yxis4=0 in 1
  
  *ddy240g
   
' D& N8 h# G. H0 }+ K   
  P/ e  {4 n! Z2 d6 A% j- x9 m  preserve
+ z4 ?" G& P& \; l' s& b4 V+ g4 u  
  sum v5 [aw =avol]
; T5 l" H0 ]+ y# p3 Q8 w$ L* m  
. B0 f! }  {/ `! |7 u: J+ s3 ]- H1 P  keep if dx240<r(mean)
  
6 C4 }& @2 c% @5 G  sum ddy240
7 G$ k( C5 Z$ T4 ~6 V; m  
' n: Q% G+ E2 ~5 X' b! p1 E% y- b  if r(min)<0 {
        
                keep if ddy240<=0
               
1 e6 U( i& h) V# S' j! n                count
( J: V4 d0 G5 o' X/ `. O; B/ o               
                if r(N)>2 {
9 |6 M2 r' q7 ?1 O" r  U. Y                 
3 ?2 F4 o" n6 y' Z  }  v6 s                 pkexamine dx240 dy240
               
5 D* P+ R$ z( V% R" _( B, G                local tempminp=r(tomc)
2 K# T, B& r9 f3 @, G                restore
1 ]. X6 w5 D: X* A  r                 }
+ {1 L  [: N, U) |, J) m$ t                 
        else {
               
& e, V$ E/ `" m: n2 x) b6 A                 di "cannot find minp"
. G9 r- U1 t% |' Z4 ]1 H5 q                 
                restore
                 sum dx240
$ G" G& v( C4 i- O8 u: Q. ]$ A$ v                 local tempminp=r(mean)
' Q) \: ~/ z3 s$ f9 T  z/ i                 
                 }
                }
   else {
                 
                 keep if ddy240<=r(min)
( S2 F. X4 ~% w* B1 O* }                 
0 E9 K$ X3 O5 {8 Q/ R                 sum dx240
                 
                 local tempminp=r(min)
                 restore
" o. }- |8 Q# K# W  s                 }
" u4 N! j# a3 d) f+ m' a


& H5 u, i$ p- M; G- V# @   gen minp=`tempminp' in 1
   *240 minp
5 ~: v$ K  t# q
' i' `, a0 l/ a6 }7 m/ o2 W
6 V& k4 K, r7 A0 ~/ T1 @3 V   
1 {# z& _2 }* F' j   preserve
   
   sum v5 [aw =avol]
   
   keep if dx240>r(mean)
   
   sum ddy240
9 \( w  I. ?  W5 R( p# F   
   if r(max)>0 {
. m/ g1 v0 w; C- h) i& ?               
/ O2 a8 u; k- P, M                keep if ddy240>=0
               
4 A& d/ q6 Z+ o                count
, Z) ^- S" ~* x5 x, ^1 \5 p' |               
                if r(N)>2 {
% ?+ C+ |) T# Z' }: A                 
" A9 t- T" L7 g' Z$ @; r2 n                 pkexamine dx240 dy240
+ C$ J/ |2 i4 G) q: A7 G9 j; y' R               
$ V' v& b& A7 t  q                local tempmaxp=r(tomc)
                restore
                 }
                 
                 
        else {
                restore
                 di "cannot find maxp"
               
9 N5 f' V! n' z: p; S) u  g               
, t" h' E2 ]- C7 ~                 sum dx240
( [+ |: C9 u# f                 local tempmaxp=r(mean)
/ f. v1 J# m: J# U, W. S9 k                 
                 }
, w* `$ Q) e/ |) a$ Q9 [                 
0 y& b$ [' N4 V6 ?9 S# j1 U                 }
8 W* z) g* ~/ k  U8 \                 
7 N) {6 _6 C& U    else {
# V# d/ }4 _6 W( l( M# y       
           keep if ddy240>=r(max)
          
           sum dx240
5 @6 g, Z8 r6 z1 m8 f$ D          
           local tempmaxp=r(max)
           restore
- I! Q; I$ f) `' ~           }
5 I4 x0 z. C1 q% W5 _: e/ H; v1 O       

3 L3 m2 |. l1 `8 o; v2 _" _: Q* E3 w  
* @. P6 |0 b8 e9 u
% t9 q# N. v; {7 Y! j0 M: p% ^  s   gen maxp=`tempmaxp' in 1
$ F4 Q+ X4 k* Z   
+ E* H( h/ J/ `* {/ q- k   *240maxp
5 ~2 p# U% [; \5 k# I5 ^% m   
   
3 W2 |# d( x7 F: K" W5 _   *-------------ddy240pminp左侧高峰值-----
' ?9 X8 W/ }8 B3 y& \- j   
: o7 a4 H0 k2 P) f9 |2 A8 p    preserve
$ `5 ?/ O  j' X9 @/ f: L( Z   
6 C0 P+ O8 f6 b; H0 ]% X   keep if dx240<`tempminp'
0 Y/ I- q) g* w5 T& F  W   
. p* X0 v4 o7 u; }# D4 W   keep if ddy240>0
! b$ ^' X9 c( M8 K3 x) F   
   count
% \1 b5 H1 f" w0 R) \               
& |( ~/ t( ]: P9 l                if r(N)>2 {
& F  O9 k9 x: h0 b9 g7 U+ X- X                 
8 }' U8 Y8 V4 ^) e& F% @                pkexamine dx240 ddy240
- U: N! ~3 X8 g! B               
                local tempminphddy240=r(tomc)
                restore
* k" f6 W: @5 `                 }
                 
        else {
                restore
$ B, H6 r; A0 v" @" t                 di "cannot find pminphddy240"
                 
' Y1 L0 A& \5 A! O               
3 ]! f* Y/ c# W; N                 sum dx240
/ h' V8 r# _0 J5 n' V1 R                 local tempminphddy240=r(mean)
6 s) A, _& `% b5 N3 Y4 w               
                 
                 }
                 
   
   
8 ]6 }/ X- P, C4 P   
4 a# |8 [) x/ P5 @   
- j4 q9 k( Q: {2 E/ W   
: d) o. z# ]8 `. m6 e        gen pminphddy240=`tempminphddy240' in 1
% ~2 U% g6 w( V* N       
1 |4 r! ~- i8 _1 J2 H$ P# K       
        *-------------ddy240pminp右侧低峰值 -----------------
0 v6 A) f* \5 m2 L/ n$ n# T' i0 a' ~    preserve

9 R) D; P* V" R0 Z. j5 ?) `   keep if dx240>`tempminp'
   
   keep if dx240<`tempmaxp'
   
   keep if ddy240<0
- Q/ M# u! v$ T1 s   
' b: T7 o6 H4 i* [& i8 [   gen temddy240=-ddy240
   
% @5 Z- n8 \4 _2 M% S9 J* [; n% z5 W    count
               
! I1 D+ h1 o: k4 }8 N3 l                if r(N)>2 {
$ S" Y) }/ A6 M9 f, I* k                 
                 pkexamine dx240 temddy240
, Y0 Z# G' z. u7 E  P5 b   
  local tempminplddy240=r(tomc)
& p3 Q2 O, g: r7 b                restore
  i% M8 g/ ~7 |6 T9 \                 }
                 
        else {
, ]- h) l) j( l. R5 l/ Z" u' D# N2 G               
                 di "cannot find pminplddy240"
' r3 `/ u0 J9 Q                restore
                 
                  sum dx240
- \/ Y* a$ g! @  F9 j                 local tempminplddy240=r(mean)
                 
1 I, g5 x2 D  N' t& P                 
                 }
) r' a' N6 o4 t' K   
% a8 f- L$ o& m9 w% p9 ?" Z2 H% s   
   
2 _0 M- I2 B9 `6 L. e! U   
9 I7 f3 ^! U6 l2 C  u5 ]   
8 i4 B* n8 V4 z9 y' |) y        gen pminplddy240=`tempminplddy240' in 1
; N, I* U+ u# {2 _       
! V5 H6 e9 o. V7 ^# `        *-------------ddy240pmaxp左侧高峰值 -------------------
6 G, i2 c2 ^, w$ ~% C' S       
7 ~# z; p- M  D5 Q3 {0 i2 C0 {         preserve
( F( B4 i  [! l3 A4 v   
2 o' ^8 i: O5 I$ I5 S) a   keep if dx240>`tempminp'
   
& g  i* S, X: K8 T- H   keep if dx240<`tempmaxp'
   
1 x/ P9 r- C- r8 M7 V   keep if ddy240>0
+ z/ m% A9 G! l+ a   count
( F# P; V2 l* W               
: D3 {8 F/ M9 ]0 {& {( B; B                if r(N)>2 {
6 \# l2 X2 i& i) a. n                 
: ^$ y: s3 _' h, W, W. Z8 W                 pkexamine dx240 ddy240
   
4 d4 K0 C$ {+ J" a: a! U   local tempmaxphddy240=r(tomc)
   
   restore
                 }
- x# L/ }7 P& r' I" b                 
        else {
               
                 di "cannot find pmaxphddy240"
                 restore
$ Z$ U/ d2 e! @- w$ E& ~/ q2 P. U               
: P+ D. x0 |6 L0 _                 sum dx240
/ f2 E: X+ @3 ~1 r( K3 u                 local tempmaxphddy240=r(mean)
* b. T$ r: J( I- Y# ~; H2 |; G                 
                 
                 }
. B% ^0 z6 s( Y3 Q+ S   
. q* Y5 _4 u0 H& y. s   
, |0 G; u" \8 _2 l1 u   
+ _3 u* V5 C$ a( J3 m   
   
$ f* h* M. s+ y8 Y   
        gen pmaxphddy240=`tempmaxphddy240' in 1
       
       
        *-------------ddy240pmaxp右侧低峰值 -------------------
       
% S! u; A6 A$ l( o+ |         preserve
) q2 U3 N8 V1 w: j* h: H) ]9 S  
! c0 |( [" B& i* ]+ ^   keep if dx240>`tempmaxp'
' M2 `% g/ V8 C   
   gen temddy240=-ddy240
' c: A- R/ V- W4 w, n# v- E   
   keep if temddy240>0
. Z! g5 H7 `) d; m   
  B0 [5 l8 N" x3 ]! E: [& K1 @   count
               
                if r(N)>2 {
  `" ~+ m" D" `                 
                pkexamine dx240 temddy240
   
   local tempmaxplddy240=r(tomc)
& ~: p0 k6 L8 R5 a) i" H. S: o5 B   restore
                 }
                 
        else {
7 l3 i/ Q. h& i                restore
                 di "cannot find pmaxlhddy240"
. ]3 l% A( I, D4 C: N) t7 m6 q                 
               
% ]' `/ N7 E1 a* F+ U) y                           sum dx240
                 local tempmaxplddy240=r(mean) in 1
! w% m, |& Y3 w, f6 D. M+ L               
                 
4 P* K% k. Y0 m) j9 }                 }
3 r- R4 @+ q/ P/ y, ?1 ^   
   
   
   
1 ]6 A4 ^  _2 x, O5 J. G) G   
( M$ v+ @- e6 i1 _9 b   
        gen pmaxplddy240=`tempmaxphddy240' in 1
       
" A# k# A8 L- i        *-------------------------------------
; j( L+ h7 F& ~' q+ W! l- j6 e( P   
$ T8 k1 C$ K5 G+ V4 z   gen price5d=v5 in 1/20
      
   gen price5yy=_n in 1/20
9 t2 L" X$ x" ~5 C   
   gen price5y=price5yy*(dy240max-dy240min)/20
   
   *price5y &price5d
2 A4 i2 l( H/ g3 @   
   gen priceny=price5y in 1
   
   gen pricen=v5 in 1
! b- ~# T0 |: p+ l4 X; }
     
   *priceny &pricen
   
3 O6 c2 K- H+ b! x/ N4 N   
   *-------------------------graph------------------
* x9 q# q) z+ d. A5 p8 R; ?+ f0 @  
6 ]" f1 ^3 E8 g) q  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)
' d" T9 w( L' c3 ?
* ---------可选项目----------------------
*(spike yxis amean,lc(blue))(scatter yxis2 amean, lc(blue) msymbol(+) mlabel(amean))
2 N0 y0 \/ Z; d$ B, Z( z- ~}
& M" J* D( g9 C6 r8 O, s' g9 uend

maodeng

awer 发表于 2013-2-5 14:40
抱歉，最近太忙了
8 i% {! m0 |0 Y, t4 J
, V! m8 m* z( |3 @! ?% g$ h! x; T7 y300236 我用下面的程序试过，没问题

哇，谢谢，我仔细看看，这些字母数字代码晕了

喜欢喝冰茶

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

喜欢喝冰茶

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

这东西是开源的？知道有个统计软件叫stata，不过要收钱。

汉水云潇

准备试一试
) P* T! y0 ^; J. ]9 ~: T先找Stata软件

汉水云潇

+ a4 N6 M, s- N1 B/ I8 Z/ e
: {* x/ T9 y  `- S" z0 {晕呀
# v, G3 `! x7 _2 z- X' H: I论坛复制程序代码
乱码夹杂在其中
" j) C. a3 H1 ~4 t% x很烦
* W9 k0 i( x/ O5 o/ ^' s$ H: `对了两次
  N4 m  h* a4 Z以为没问题了

结果
& ^! Z9 ^7 _: R  b# v执行后
还是出现
$ _- U7 K6 e* s5 O$ n4 Rr 199
! L9 W+ a4 n* p" j8 `8 S的错误
( j  Y7 N1 q; T/ {$ g

汉水云潇

试验成功，下面是三一重工的图

$ _% E6 p# T/ [% c0 p

peterbev

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

中秋不圆

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

bknite

这个帖子值得学习。

jeanny

感谢！！！

Hubert2000

明月回春 发表于 2011-8-13 20:17
, \4 o: G  T. d5 ^! Y程序写的不错.
不过有一点要注意.通达信导出的数据,只对价格复权,并没有对成交量复权.所以数据需要手工处理 ...

* g8 D6 s. |) r* {8 O1 a老哥，也是山寨这篇文章找不到了，不知为何？