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

maodeng

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

maodeng

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

maodeng

. hello 300236
statement out of context
r(119);
9 U) b6 e' r/ V  c; }- ^
大师 就是这个

sergio

这贴实在！

awer

maodeng 发表于 2013-1-28 13:13
. u# l8 x( H4 C8 W5 F; I4 W6 ~$ B. hello 300236
7 d# ?  [8 {& |. ~% Estatement out of context
r(119);

% X, g9 [$ ~# @+ l* k% ~# x2 m抱歉，最近太忙了
0 I" l, i" x+ `. V- l
300236 我用下面的程序试过，没问题
* i- @( N2 S- Q% Q6 B; g; ?
你把 D:\stock\test\source 改成你导出的数据存放的目录，再试一试
1 [, J1 S! @: |; Y
capture program drop hello

program hello
' I- F3 v3 m, P3 Q/ J* F
quietly {
: V, _5 |6 d& ~* Z
clear

8 _3 E: u" \( k, w, Pcd D:\stock\test\source

insheet using `1'.txt
  t$ J' e% c0 G( M9 f! i2 l2 u/ m3 |
; u) a* w0 l5 @& H$ G   keep if v5>0
/ a% q$ k4 R, s0 _/ U4 D4 R   
9 T' {  S* f* N. f0 i, {7 x, p$ W" t   gen vt=_n
9 p) q( O2 v4 e) h* {& }   
   tsset vt
' V& K5 G( d5 Z) n  E
7 l% }+ o. O( F5 Y+ A( r; Z   gen zf=D1.v5

) A  S6 n* h- `2 x- L/ n: K* m   gsort -v1
/ V/ H6 }8 ]. ^# {1 ]& S0 F+ W
1 R8 D3 x' P& S8 n   gen avol=v7/v5
   
: I: k5 z/ G/ u. @3 O4 U   keep if avol>0
   
9 ]- J: ?* ?3 C) G# O! S) d   drop vt
/ E. c% w0 k( T0 G6 g5 N. @  
) [+ i) n5 s, z! E- A8 T% I   gen vt=_n
% V, O6 Z1 c% T; B# O   
8 }9 d8 v9 n8 `$ _+ r   tsset vt
- M( }0 Q1 A5 w   
    sum v5 [aw =avol]

/ w2 B5 F4 J7 Q- R2 {& l   gen amean=r(mean) in 1
   
   * 总加权平均值 amean
0 r9 A+ q" M0 ^7 u8 _8 M$ p& z- C   
   
   keep in 1/240
  i4 n! ]( {- Q) _  s
5 B  l& x0 n" b: Q) x   sum v5 [aw =avol]
7 \8 K6 s6 u- U, |$ @2 `- a9 r* f
   gen mean240=r(mean) in 1
   
   *240加权平均值 mean240
4 g, E- i* g; ~9 \- Z( s   
   set obs 500
* v$ B% \* W* W6 V   
   kdensity v5 [aw =avol],gaussian n(500) generate(dx240 dy240) nograph
   
8 R# C# H5 G. ^' e" n4 M7 h   *240 kdensity(dx240 dy240)
2 u+ E* w0 t+ |. T& O/ S( s   
0 F- L  B! i$ D   sum dy240
   
3 h5 H/ K- [% H) ]0 W   gen yxis=r(max) in 1
   
   gen yxis2=yxis/2
   
   *yxis
7 l5 w$ m$ ~: L: }) ?. O   
; \; K$ e) p. A$ y   gen dy240max=r(max)
   
   gen dy240minp=r(min)

9 J, `: m9 s' z& Z# u/ C    gen vtt=_n
, r! o- W: q& p9 S
    tsset vtt
  b9 B" D7 M/ h+ d, ~
/ ?2 w& j0 X2 ]* h   gen ddy240=D1.dy240/D1.dx240
   
8 a/ J0 _" _! S$ W; l   sum ddy240
- Z. I! R) ~2 O1 G8 Y7 a   
   gen ddy240max=r(max)
   
# n# C' A% Y' o7 l1 i* a   gen ddy240min=r(min)
% B" ?4 z7 w; b! ~- k: ~     
$ u6 E: V+ w$ G  L! c   gen ddy240g=ddy240*(dy240max-dy240min)/(ddy240max-ddy240min)
) Z- X" j8 u$ f( f   
   sum ddy240g
   
   gen yxis3=r(min) in 1
5 B' ]$ T' u5 i2 d$ z5 ^: F   
   gen yxis4=0 in 1
/ W, {6 }7 A9 p6 Z* w; \' V  
9 {) {# K8 m: W6 n0 E  *ddy240g
) K# ~' B2 ]3 H2 }   
7 w! X" S% o8 F* h3 N8 V$ I   
8 {9 o" h+ e3 y: V- |" b$ Y  preserve
8 y8 Y( e  i4 i; r1 |4 f7 v5 h  
  sum v5 [aw =avol]
  
  keep if dx240<r(mean)
' f1 f2 f0 [8 C  
" ~8 M1 Z4 l+ @. b( b  sum ddy240
  
  if r(min)<0 {
1 a" R/ \2 c6 @( m8 m, B1 D        
- Y# f" V- I/ f                keep if ddy240<=0
4 V9 \# |; E6 u$ N! d' k               
                count
               
% S3 |$ E" b% f3 u                if r(N)>2 {
9 i* ]- Z/ t$ W, c$ b                 
                 pkexamine dx240 dy240
               
                local tempminp=r(tomc)
& U2 T2 r* T3 |! ?                restore
                 }
8 C+ Y6 I, h& H# B# Q2 c                 
        else {
- e) d0 @/ \' M) |$ V               
                 di "cannot find minp"
/ u0 p2 r3 g% B$ q                 
                restore
* b, G( F( V# g8 N0 H                 sum dx240
% f6 n) k. Q4 [  T; S5 \( W0 \                 local tempminp=r(mean)
                 
                 }
& X5 I, X+ z8 A$ `3 [; ]" d                }
6 [& c3 b  l6 w   else {
                 
                 keep if ddy240<=r(min)
                 
- a' U0 m. o& U5 o% z) a5 o3 m) @                 sum dx240
" F, L: P2 N& {7 }3 S8 K                 
                 local tempminp=r(min)
                 restore
                 }

4 n& K0 V6 T$ |8 ~9 T+ \
, c" _1 B' q& Y9 `  ^( F6 i% v
   gen minp=`tempminp' in 1
   *240 minp
  B$ v! c, h* k

1 w* _$ `8 t9 m8 K   
   preserve
   
   sum v5 [aw =avol]
. y8 B3 ~: I8 O$ ?" S) ?   
7 |- s, S+ z9 k1 J   keep if dx240>r(mean)
0 w3 B8 ~4 S5 J& F' Q# w: L   
; @* Q# X! R3 T$ d& g8 @   sum ddy240
   
   if r(max)>0 {
               
                keep if ddy240>=0
( M4 v: J* K& E( J& x, q- X; V: t               
                count
               
' }3 j2 l# V) ~" P4 i                if r(N)>2 {
0 h- Z+ `4 p! }2 |2 d% e+ z" V                 
                 pkexamine dx240 dy240
               
, ?7 B4 Z% j% n                local tempmaxp=r(tomc)
                restore
                 }
5 w8 G0 `; S4 H2 R: W                 
                 
, l2 `5 u; a: L7 t0 i        else {
  n# K: [7 v$ ?* y* f                restore
. A9 V" d, p- c" J. ^. ]6 D+ I                 di "cannot find maxp"
8 v7 {2 i5 ^7 h1 \! |* `) f               
               
' H( L+ d6 j0 e3 [  V                 sum dx240
                 local tempmaxp=r(mean)
                 
8 x/ g% w, p  Z/ f                 }
1 b2 m6 R8 J5 m$ a- j9 @2 z                 
                 }
) `  \0 y" a6 B0 f                 
    else {
       
           keep if ddy240>=r(max)
7 q9 C9 Z. F: h. V3 x7 `0 a0 B& k          
           sum dx240
          
4 Z2 Q* F$ i0 r5 D! D( S9 w           local tempmaxp=r(max)
           restore
           }
       

% ^& q$ b' e. H: \5 j  
- ^4 b( H1 v$ B! L% i
9 E+ G) s3 b7 d: q! i% N1 O/ r   gen maxp=`tempmaxp' in 1
   
   *240maxp
   
   
' L0 D" z3 z& e6 {9 s   *-------------ddy240pminp左侧高峰值-----
2 J8 t# ^( E& h/ s: n   
    preserve
   
   keep if dx240<`tempminp'
, [3 m7 G- C7 Z   
   keep if ddy240>0
   
. b9 t6 n6 u  W8 ~) ^2 p   count
               
& U6 i- Z* i  c6 J                if r(N)>2 {
' W/ F0 [4 R9 p3 v$ ~/ v                 
# \3 R6 K" G6 D) a3 Y1 Z7 J                pkexamine dx240 ddy240
9 n: c; @2 j9 w/ y; b5 \               
                local tempminphddy240=r(tomc)
" f+ L  i+ X# L4 b8 x                restore
; i6 H  b7 d( R+ j% `, q' T) H/ \                 }
* C/ F/ I. u" _; T8 u2 f( y                 
        else {
7 I" c% Q! x; G* j                restore
$ H9 f# ~. v% [& P2 ?                 di "cannot find pminphddy240"
                 
               
) H9 y+ `4 w; `- a+ |  Z5 {                 sum dx240
                 local tempminphddy240=r(mean)
: c! ]7 t1 e4 r. }               
                 
                 }
                 
4 x6 h! f! t. }% P* e* {6 R8 i   
   
   
) h/ ]+ W; k. W1 T% d7 k   
   
( D" i6 X7 \) J( P        gen pminphddy240=`tempminphddy240' in 1
       
: [% i7 N5 a9 X" j* I       
6 T' r8 Q0 B* m6 J% P$ @* J        *-------------ddy240pminp右侧低峰值 -----------------
% F+ ]$ {) ^# O  y- W$ u    preserve
0 I+ ?) K5 \+ t& O
   keep if dx240>`tempminp'
   
   keep if dx240<`tempmaxp'
   
- y5 b, d. H, q3 {. ]   keep if ddy240<0
   
   gen temddy240=-ddy240
/ ~" f( ]4 D* y; E& W! h- N! P   
% k& E5 l  `9 T/ O    count
0 f7 y3 a' [* ]& M9 U               
7 C. D4 e* T' c6 `- ?                if r(N)>2 {
* x& r6 x& z4 O                 
                 pkexamine dx240 temddy240
   
  local tempminplddy240=r(tomc)
& U3 Q# w5 i2 N                restore
                 }
                 
+ |8 ^, N" U+ ~! V        else {
               
# a. n" C) c4 H6 P' Z6 w- P' c                 di "cannot find pminplddy240"
                restore
                 
7 G! z5 L2 c2 c4 x) V- T                  sum dx240
# @) d; [! |$ Q9 c                 local tempminplddy240=r(mean)
" ^. ]9 k- D1 u5 |1 K                 
3 n$ D! M" K" |+ Z                 
                 }
  P2 w: M/ J. o% r4 f9 t   
9 V1 W: ]) m8 D: R- |# l   
   
6 H. T5 R( _# A3 m   
   
5 N! X, Y8 X( Y: g6 ?" N        gen pminplddy240=`tempminplddy240' in 1
. y$ \; J" k% M0 [# b* N2 x* h8 }       
& M( f" G7 S. f6 `$ a        *-------------ddy240pmaxp左侧高峰值 -------------------
       
         preserve
   
/ q- i& t; Y" u   keep if dx240>`tempminp'
   
   keep if dx240<`tempmaxp'
   
' d  i- N" }+ C- o! k   keep if ddy240>0
- M- m( d7 k7 l6 U& H' C% x4 Y   count
" W# u! o& S/ A( m: M               
6 W( i: a$ r' H/ C8 Z0 f                if r(N)>2 {
- J$ ~$ o& p6 w! _                 
( \1 K; X' ^4 k) m& B                 pkexamine dx240 ddy240
   
   local tempmaxphddy240=r(tomc)
   
8 b% n/ m: I5 s& r% T; }- q( O$ M   restore
# q$ v. z( Z1 \6 @% B* [( `! ^  G                 }
2 ?  ^, w: m- z                 
        else {
               
                 di "cannot find pmaxphddy240"
                 restore
               
                 sum dx240
                 local tempmaxphddy240=r(mean)
                 
0 h+ ?/ I# r( }" ~( E" S" Y                 
                 }
   
: q$ v" z! T7 R6 S   
0 f+ M% [# C1 I; G* S3 T: Q9 Q   
   
2 Y; ~# \) d; F$ H7 O. ]   
   
* [1 f$ z6 u+ N- L4 Q2 F6 H: v        gen pmaxphddy240=`tempmaxphddy240' in 1
5 ^* j; ?% K9 {  o5 I, V  G& X6 V       
1 V; W8 G* Z' T* }       
        *-------------ddy240pmaxp右侧低峰值 -------------------
0 }; y# J) I$ C6 v" v       
         preserve
, V# k5 s) J0 {( g/ e- h  
+ g8 t! |) K4 O7 c8 h6 N  u   keep if dx240>`tempmaxp'
   
   gen temddy240=-ddy240
   
   keep if temddy240>0
7 P9 m3 A7 y7 R; u% Q* s. B2 \   
   count
* P1 v% E, C- R6 x  v: ?               
                if r(N)>2 {
                 
                pkexamine dx240 temddy240
& {, ?9 ?% H! _   
   local tempmaxplddy240=r(tomc)
0 C' c1 Q' K( s% d( K5 Q: i   restore
                 }
( \2 p' g$ T9 a: W                 
        else {
* B# e9 y0 M! m6 f' J, J: e. ]3 I                restore
                 di "cannot find pmaxlhddy240"
) G( a* g" C' I3 a                 
9 L) O1 G, |' ]# `, R               
                           sum dx240
                 local tempmaxplddy240=r(mean) in 1
               
4 V# v5 n( ?7 w! v                 
! H8 }) w8 N  y6 s                 }
   
   
& U1 {- R+ D+ ~1 P" A+ K4 I   
1 U, y3 C3 t. e9 O6 {1 K   
   
   
        gen pmaxplddy240=`tempmaxphddy240' in 1
7 s, N7 V3 w  b: M8 ^. J9 @; N       
/ Z0 y1 o7 o1 |4 }( Y$ z; A! C        *-------------------------------------
   
   gen price5d=v5 in 1/20
- ~- k5 f$ n# @: p7 ~; H      
   gen price5yy=_n in 1/20
6 @" J  j4 c7 g1 f8 N/ U4 _5 L   
" W' u4 e  A1 j   gen price5y=price5yy*(dy240max-dy240min)/20
; J% u( N$ N" {+ K" N; L   
( {8 q* h+ N/ U- G   *price5y &price5d
   
   gen priceny=price5y in 1
   
. A0 ^7 }4 L0 R( F+ B   gen pricen=v5 in 1

     
   *priceny &pricen
" O$ Q: |. ?' J. r# u! [) Y( x   
9 d! P1 D* M$ g; F! t   
3 o' k% ~# r3 c   *-------------------------graph------------------
' s! v! X; V0 K3 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)
( B! G% M1 J5 j
* ---------可选项目----------------------
*(spike yxis amean,lc(blue))(scatter yxis2 amean, lc(blue) msymbol(+) mlabel(amean))
, l$ F) X4 ~; k}
end

maodeng

awer 发表于 2013-2-5 14:40
抱歉，最近太忙了
3 c  Q5 D+ Z- U& U& }% h5 G* c. Q# {4 w
300236 我用下面的程序试过，没问题

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

喜欢喝冰茶

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

喜欢喝冰茶

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

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

汉水云潇

准备试一试
先找Stata软件

汉水云潇

9 a6 Y7 F" a& n3 \- e& n
晕呀
0 n6 D7 v/ [5 X论坛复制程序代码
. y- [% Y+ R5 i+ e9 W8 o  q乱码夹杂在其中
6 `" M$ \* T  Y' D6 r" i很烦
对了两次
以为没问题了
& u) b- s" q7 \" w+ \
! N# J& g7 e# _7 w/ M结果
执行后
' d. y% K0 z! g( `; \% O( K还是出现
r 199
! |3 X7 k  [% Q* ]8 ?6 h8 N的错误

汉水云潇

试验成功，下面是三一重工的图
2 V: e& u/ C" N$ j
( a$ Z) |7 n# N

peterbev

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

中秋不圆

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

bknite

这个帖子值得学习。

jeanny

感谢！！！

Hubert2000

明月回春 发表于 2011-8-13 20:17
* z6 @2 z6 z+ q( {! e5 `# R& j程序写的不错.
7 t/ f, b5 R5 Y) T$ D6 w7 @+ S5 O, g不过有一点要注意.通达信导出的数据,只对价格复权,并没有对成交量复权.所以数据需要手工处理 ...

( l7 V+ R3 q$ U+ T2 K6 ]8 X6 {$ d老哥，也是山寨这篇文章找不到了，不知为何？