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

maodeng

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

maodeng

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

maodeng

. hello 300236
statement out of context
r(119);

' T1 S# D. a0 Y' Y& Z$ J大师 就是这个

sergio

这贴实在！

awer

maodeng 发表于 2013-1-28 13:13
# H0 i) ^1 ]1 g* g. hello 300236
; z. `6 }8 @8 a! O8 Y2 A; ystatement out of context
r(119);

+ A0 y4 `, T2 \( ^# M抱歉，最近太忙了

300236 我用下面的程序试过，没问题
0 _6 ?! w$ \/ b0 T; |
你把 D:\stock\test\source 改成你导出的数据存放的目录，再试一试
+ I% t+ D5 O* S4 j( t! B$ `) X
" M5 a: k& W  `capture program drop hello
* L4 u. |. j5 q8 w8 a- J3 Y
program hello
: f/ ?3 w( v/ v# e
quietly {

clear

cd D:\stock\test\source
$ C9 n8 q' B  f" A- W$ i' ?
insheet using `1'.txt
5 G  _1 w/ [- L) U( ?. a
   keep if v5>0
   
3 R+ o# _( ?/ e   gen vt=_n
   
   tsset vt
& }6 e/ z- G- Q" y; Z
1 \8 O- q: C0 [; g$ K   gen zf=D1.v5
: R- F  ^: F) c7 g& S/ x1 \6 |
   gsort -v1

   gen avol=v7/v5
- g2 W) Z" g7 w$ n7 m  ~+ }   
   keep if avol>0
   
   drop vt
  
4 h2 W- J/ I, A) A   gen vt=_n
   
   tsset vt
9 f' E# D! N! M8 M; ^   
    sum v5 [aw =avol]
) X% k, P# a" u; J) T2 N2 g
   gen amean=r(mean) in 1
   
   * 总加权平均值 amean
   
* B( n* \& \2 V! l$ t7 F9 ~( ]& C   
$ Z! s. ?6 h5 ^7 Y* ]7 j, W   keep in 1/240
, l. d4 j5 _' z+ S$ L$ g
   sum v5 [aw =avol]

# U. H- m$ H' ?; Y( U   gen mean240=r(mean) in 1
   
8 ~6 Y+ y$ b/ A. t   *240加权平均值 mean240
# L! l+ ?6 s# i, b; k' n5 P+ o   
   set obs 500
   
  z6 E3 N4 v. F! i   kdensity v5 [aw =avol],gaussian n(500) generate(dx240 dy240) nograph
   
2 F, L) c, V8 j% V; ?   *240 kdensity(dx240 dy240)
   
; l- k) p0 w+ E9 Y' S3 l& n1 V   sum dy240
   
- t; d, m. L9 Z: C2 ]   gen yxis=r(max) in 1
; W  d4 u0 h0 l   
   gen yxis2=yxis/2
   
   *yxis
   
   gen dy240max=r(max)
- x! c/ h. J4 P. y4 k& i   
& G% }& m$ R# e' B; T; b   gen dy240minp=r(min)

    gen vtt=_n
( g9 |+ I1 j- K8 L% D0 e) z2 }
    tsset vtt
) p+ a2 c3 I6 m- j( o* f/ L
+ N6 s- p+ {6 O4 X* \   gen ddy240=D1.dy240/D1.dx240
4 K. h! }2 u9 o0 o' Q6 |   
" [+ c3 D* A4 _) k; _0 o   sum ddy240
   
   gen ddy240max=r(max)
8 o9 u% y! s% ^8 ^: i; g   
   gen ddy240min=r(min)
! n, D7 u/ v* R7 ~     
1 ?' }8 Z( M3 L8 t, U4 o   gen ddy240g=ddy240*(dy240max-dy240min)/(ddy240max-ddy240min)
   
* a0 s! u3 N* X  ~6 o' Q6 I+ d   sum ddy240g
8 [2 x" h9 j9 q+ z) [# D! |   
3 x3 f8 H" E/ }! F0 Y   gen yxis3=r(min) in 1
   
   gen yxis4=0 in 1
  
  *ddy240g
4 L$ `. P7 ~0 n0 _- `4 B   
   
( b) Q+ m5 k/ o4 v' s% K- N( v  preserve
  
  sum v5 [aw =avol]
  
  keep if dx240<r(mean)
  
$ |& c1 \" \/ H5 J, X9 j& j  h  sum ddy240
0 _6 @4 {1 [/ W2 H& \  
1 D% @' H/ y. D2 L5 d/ \  if r(min)<0 {
  t7 G; Z" Q* g' C3 n        
                keep if ddy240<=0
               
' Y$ U) s- c  b9 C                count
               
. U" K  I% F& |1 P; u: d( \/ m                if r(N)>2 {
                 
& k% Z1 ]! `  ?9 k5 b                 pkexamine dx240 dy240
6 R" u- d( @! H, o* }               
- I/ K$ f( a3 T- v                local tempminp=r(tomc)
                restore
# G' c! ~6 P* T4 F0 R3 X6 v                 }
                 
        else {
               
                 di "cannot find minp"
                 
                restore
. K* A; v/ e$ a& t5 |( U$ Y  J                 sum dx240
                 local tempminp=r(mean)
                 
                 }
                }
4 `. @! U$ G+ e7 }# N3 }! \   else {
% P( `! M/ v+ L3 F% M2 }+ e                 
                 keep if ddy240<=r(min)
2 x# B1 A  k0 u                 
' @/ @" l3 L, |) x1 l                 sum dx240
+ w* g' n0 Z$ T# W: ?5 |$ J; z                 
, T4 q0 G8 K) i# x3 S/ t                 local tempminp=r(min)
                 restore
) W* A! F' C" O& D                 }

+ P3 e! t* R2 p/ S6 \$ Y4 }
8 g4 {. a6 K. G8 U
   gen minp=`tempminp' in 1
- A! H0 Z/ N# ~9 v% S; v/ |   *240 minp


   
" f' V+ `, l7 e; |   preserve
3 O. r9 p/ `" j# ]   
   sum v5 [aw =avol]
   
9 a3 p) C' q/ N. l% w: d# N) k/ m   keep if dx240>r(mean)
   
   sum ddy240
   
2 t! {9 @9 d$ S. r2 I! s$ A7 ]   if r(max)>0 {
               
1 T( s/ Y) e* z  F0 N2 c( B                keep if ddy240>=0
               
                count
               
( \  r6 }! q- r$ F9 B4 X                if r(N)>2 {
                 
                 pkexamine dx240 dy240
" A) M* m, x# [( J. w9 N0 y2 o8 H               
                local tempmaxp=r(tomc)
2 p5 o. E4 S  M) g& }2 ~( M                restore
                 }
                 
                 
        else {
                restore
                 di "cannot find maxp"
6 x% Z1 J; p; ?: O9 p$ m3 h               
. {# z' @5 M4 o  u' i2 J- `5 M               
5 d" Y9 s2 }, o, E' L                 sum dx240
# p+ u) h" U5 G5 H: M$ Z7 J  R                 local tempmaxp=r(mean)
# W7 Z( C5 L- i! a8 z                 
+ j6 R; J: T/ v3 E                 }
3 J: s  `. I3 m7 [9 L                 
                 }
                 
: o' z7 y. ]# Y  @' v: \3 {2 D    else {
       
           keep if ddy240>=r(max)
          
           sum dx240
( h, \+ ?1 M5 n( e          
' b; A, ^0 M6 k8 m: [2 Z1 }8 ]           local tempmaxp=r(max)
: A9 f2 q" }+ g. d  Y           restore
           }
& @' V, x  ^2 D$ z( u4 s% |' Z       

2 Q8 V2 d7 G1 N7 p2 |  

: Q0 ?( c- _* B+ f  d   gen maxp=`tempmaxp' in 1
   
7 u& g; P- ]% H; A$ y; [" C. O   *240maxp
+ K* R; N4 n& V/ i6 ]9 n+ V& C( p   
! |. M- G1 Y! X0 s6 y0 T   
! r' m- ]+ _; {# e  |7 x# _   *-------------ddy240pminp左侧高峰值-----
& j- C1 D- B7 R8 O   
& J+ ?0 P4 i3 _& n  _    preserve
$ b$ ~+ B2 z9 H   
3 O1 [. r+ F9 p( E   keep if dx240<`tempminp'
   
. I( h# J# C  o' j  _/ j) O3 _   keep if ddy240>0
   
   count
# i) l) V3 |3 h               
                if r(N)>2 {
4 P5 V8 ~5 c0 T$ F( ^) k; t                 
                pkexamine dx240 ddy240
8 }1 s; B8 y* J8 {               
                local tempminphddy240=r(tomc)
                restore
( X( D: ?3 z1 ?8 V' k                 }
                 
        else {
1 m- @5 m. Z8 x, U0 K; Y2 v) Y3 m                restore
& m( i+ C& ^0 l4 D                 di "cannot find pminphddy240"
                 
               
2 |- P, e- I  q/ ^, I5 F                 sum dx240
% D5 ~/ r7 O3 l2 H, r# ~; f                 local tempminphddy240=r(mean)
- X  b: \( i( p  x3 {; i& W2 v               
' m* r6 X) L1 _! {1 {) c* N- E                 
: ]9 N) R: C. P# g+ n; i4 W                 }
3 W7 ~, f) E& h1 W+ N! `                 
   
7 I' m# Z, U5 v! E% C0 O6 m) \/ k   
   
   
   
2 \( n; h! K% L        gen pminphddy240=`tempminphddy240' in 1
; l) ^* X3 {! s       
" W( ^9 h" r7 k# B. d       
        *-------------ddy240pminp右侧低峰值 -----------------
) F9 q' J0 V; ~4 f    preserve

1 _4 M% L- E) O3 j* R) r* D   keep if dx240>`tempminp'
6 s6 l$ \# @- s0 a' j4 q# U  j   
7 ~" v: K# b8 J   keep if dx240<`tempmaxp'
   
5 A; z, B  d5 ?) d  t2 L; K   keep if ddy240<0
   
   gen temddy240=-ddy240
3 N3 P& e: c+ T1 t) k& U* K4 ^7 c2 q   
0 B! }$ {1 `3 s! w$ z% A; s) o    count
               
                if r(N)>2 {
                 
                 pkexamine dx240 temddy240
   
  local tempminplddy240=r(tomc)
                restore
) C3 a  N$ d; m6 d7 m5 X, _                 }
# N7 E9 z2 e" N5 S4 W                 
        else {
- `7 ]% t: Z$ I- a  W4 o5 B               
                 di "cannot find pminplddy240"
8 }6 M5 k% K' L- c                restore
. k( Q: a! N9 l4 U( o                 
                  sum dx240
7 T! U% y) Z. }% D: B                 local tempminplddy240=r(mean)
1 f7 H3 o8 O- C7 r  U                 
                 
: B( E, Z( ^- N; h. }                 }
   
* `* E3 S, F3 l7 Q( m- y) K   
   
   
   
; n- R- c! e( q, {. ~$ J7 ]& B        gen pminplddy240=`tempminplddy240' in 1
) z6 L* r; S+ V; b- Y       
        *-------------ddy240pmaxp左侧高峰值 -------------------
       
         preserve
7 |  f9 G1 a  f: T   
: E$ x2 u5 L3 W# R$ E   keep if dx240>`tempminp'
6 {6 a$ j7 C& h8 A: g9 N$ ?6 l# |   
4 p# I  S# L+ x# @: ?   keep if dx240<`tempmaxp'
$ x& s! E6 n; \   
   keep if ddy240>0
   count
" l3 c9 b6 z- `1 k, K% K) A               
                if r(N)>2 {
, a+ b/ B- d$ _, v& G, _                 
                 pkexamine dx240 ddy240
   
   local tempmaxphddy240=r(tomc)
  C" A# L/ o# t5 s% p" Y   
   restore
$ q- M6 t3 [; j  g* H' p$ O. N                 }
                 
        else {
8 \  s" N7 w. C$ s: ?( `               
                 di "cannot find pmaxphddy240"
2 ^. J9 Q- l8 T                 restore
2 b2 c% e! \0 p3 L               
3 }! c# W! F7 z- o/ G                 sum dx240
                 local tempmaxphddy240=r(mean)
                 
                 
+ Y. r* F) z- Y9 y+ Y# u                 }
2 I( I+ z6 l7 x$ T2 S3 \' t   
' `' j0 `" N0 k2 ^) ~3 T   
   
+ Q6 ~7 J" \7 I$ ?/ M- h. N   
% j4 k: z$ e8 z  {9 f$ r. C   
3 W" w6 s) n" V' b# i7 P/ z  Z   
        gen pmaxphddy240=`tempmaxphddy240' in 1
       
       
        *-------------ddy240pmaxp右侧低峰值 -------------------
% z4 x0 a  I" a. T       
8 S. [! x. K5 y2 T/ `- u$ E9 e         preserve
; }7 o7 ?- S  s  
  i& [! f8 m9 ~6 t   keep if dx240>`tempmaxp'
, b# o. |( S2 F+ S% @5 Z   
   gen temddy240=-ddy240
) X" U; |; _5 j/ Y. J3 C  R; S$ D* U   
' y- x4 }) `+ U   keep if temddy240>0
3 d9 V% O: M; M+ Y   
   count
               
8 M5 e! `' o9 E5 T1 H                if r(N)>2 {
                 
                pkexamine dx240 temddy240
0 k3 v  N. i' \: x: T( c   
   local tempmaxplddy240=r(tomc)
& }/ J$ `& @1 m0 q, ]2 a" B/ _   restore
# Q7 F7 f9 Y- P5 {                 }
( |8 A' V0 @' {1 V3 i. M' z                 
  x. O6 s5 s' L2 q        else {
% a4 \' N$ u3 P6 T9 {                restore
  @2 X, A) l! A- h6 g                 di "cannot find pmaxlhddy240"
; v) O! N4 q. W! |  v$ r9 k                 
               
6 @; E2 |7 M: d* _2 u                           sum dx240
; p  B# v( c7 h% u+ M1 {# E  q                 local tempmaxplddy240=r(mean) in 1
3 Q+ ^3 _, |, }" y, t5 M4 `               
% l/ g, m/ {7 j5 |) i2 s8 g) s                 
                 }
" G8 V' M! Z5 w8 {   
1 Q  P2 `5 a- l; ]8 T   
   
   
   
1 a3 ~5 j' i$ u   
        gen pmaxplddy240=`tempmaxphddy240' in 1
       
        *-------------------------------------
1 o: y* }2 v3 b! C0 }2 k4 G   
. O+ ~0 i: P, S& i2 H2 o   gen price5d=v5 in 1/20
      
   gen price5yy=_n in 1/20
) C$ C& e9 ]3 r   
   gen price5y=price5yy*(dy240max-dy240min)/20
  g7 D* L3 u( |2 g) i. |   
   *price5y &price5d
4 N3 l' J! F5 m/ f   
   gen priceny=price5y in 1
   
* r! I$ _0 ~5 {1 O) k3 s   gen pricen=v5 in 1

1 }) `6 G; K. a9 q% w% m" P     
% z0 a" A; h  y& T$ @  F9 {   *priceny &pricen
& X# l. m6 x0 Z4 |& K" O   
   
8 Z: j+ m' _; g% {   *-------------------------graph------------------
2 u6 D* e6 x% G5 o) Q8 F  
  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)
& p  H; j5 `6 L- d
* ---------可选项目----------------------
*(spike yxis amean,lc(blue))(scatter yxis2 amean, lc(blue) msymbol(+) mlabel(amean))
' X7 c- Y6 o+ W, s1 U+ Q% C3 c}
end
0 D& i+ w+ T) F. k8 S

maodeng

awer 发表于 2013-2-5 14:40
9 S& \/ F$ ?" h0 b# ^* I/ h抱歉，最近太忙了

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

; H7 W0 p! t- q  M" I" z5 E哇，谢谢，我仔细看看，这些字母数字代码晕了

喜欢喝冰茶

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

喜欢喝冰茶

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

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

汉水云潇

准备试一试
先找Stata软件

汉水云潇

晕呀
2 b# d4 N" k# n  I; B论坛复制程序代码
乱码夹杂在其中
很烦
对了两次
0 R( X) K4 i6 v以为没问题了

# h& k- A( G5 ~. Q6 h结果
+ H2 m% s( e& ?. \4 J5 i% d( K9 z执行后
还是出现
r 199
的错误
8 L7 j" j/ ?1 E

汉水云潇

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

  |/ d% v, T# D/ q

peterbev

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

中秋不圆

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

bknite

这个帖子值得学习。

jeanny

感谢！！！

Hubert2000

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

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