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

maodeng

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

maodeng

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

maodeng

. hello 300236
! f7 [9 ^- D7 h3 }& ?* J* hstatement out of context
& D- I2 Z4 `7 R" R7 ~4 cr(119);

4 O- Q7 ]/ l' ]8 F4 x1 K) O  M+ @大师 就是这个

sergio

这贴实在！

awer

maodeng 发表于 2013-1-28 13:13
. hello 300236
statement out of context
) j# U$ _& [! x3 i% vr(119);

抱歉，最近太忙了
3 H" N4 G& V  J; B5 i7 [
" J! y" }$ I; z300236 我用下面的程序试过，没问题

你把 D:\stock\test\source 改成你导出的数据存放的目录，再试一试
4 X. E+ V1 n: e1 N% N5 f
capture program drop hello

program hello

' c2 p3 V' X9 d6 v. U quietly {

clear
: s/ K) t9 w1 {) k% R
cd D:\stock\test\source

2 U9 L7 O  v/ {$ l# Finsheet using `1'.txt

   keep if v5>0
- Z" c, s; x# A# N! Z# Y% J0 H; z   
   gen vt=_n
  H# S5 {  d. U( U$ N- z   
   tsset vt

8 \: V. |9 ^* J0 @6 E4 S' }& t5 p   gen zf=D1.v5

0 a3 p, S6 k0 {3 ^) Y! _  Y   gsort -v1

   gen avol=v7/v5
   
7 O$ k, x6 }4 S# t; i   keep if avol>0
1 D& \9 D6 Q5 U   
3 h2 b( O6 A+ b( {# E  R   drop vt
9 ~0 [0 l* y" D- W  
7 b, I# s- Y0 [7 z  }   gen vt=_n
   
, g4 D& ^/ D/ V5 t" a8 e   tsset vt
/ O. t3 q4 s/ k" v9 T6 i) `   
2 A; E& X& L( V    sum v5 [aw =avol]
1 ^' H3 h+ P' s6 \2 j
1 X; F4 v( {8 ]% Z   gen amean=r(mean) in 1
; W& [' Y! |. s+ Y/ a5 J   
   * 总加权平均值 amean
. R" T- L* P: z; E   
   
   keep in 1/240
( V$ C. z* ]7 p) ]& [
" Y5 Q3 a% K2 K, i" e$ T8 E   sum v5 [aw =avol]
; W, }1 r6 W/ u8 k. E
   gen mean240=r(mean) in 1
   
; A2 ^9 r- o& G- c   *240加权平均值 mean240
% v4 w" c- s/ \* G- V9 x   
   set obs 500
   
   kdensity v5 [aw =avol],gaussian n(500) generate(dx240 dy240) nograph
" B- g) B7 b% Z   
7 w* I  P# `3 \( Q" E   *240 kdensity(dx240 dy240)
   
5 ^3 z) n0 i$ l5 b* M2 z+ n8 @   sum dy240
   
   gen yxis=r(max) in 1
   
) T$ S" y5 V8 {& d% E$ t   gen yxis2=yxis/2
( N9 D# x- T; S# E3 v& h# l$ o   
   *yxis
4 @( p; t& F4 y' E8 J   
' C0 i! J8 t1 Z! O/ j) n   gen dy240max=r(max)
   
: t- g( s) E; g   gen dy240minp=r(min)

1 q) H& w$ n( N: p# F! m, b    gen vtt=_n

    tsset vtt
1 T( z/ r- k+ P9 Y/ c
, y! Q. X  q. @  i   gen ddy240=D1.dy240/D1.dx240
   
   sum ddy240
. L" K* Y& U2 @  m+ [& J) C   
   gen ddy240max=r(max)
/ y& c% j) S+ Y( u- A; ]; p   
   gen ddy240min=r(min)
' P( O2 V1 C% A* U9 c     
0 G9 x/ o: M$ s( t. \) L3 B" O5 Y   gen ddy240g=ddy240*(dy240max-dy240min)/(ddy240max-ddy240min)
- q! ?$ _( ]) f- c( r% M* F   
& O  R' l7 J2 c* D2 H" C   sum ddy240g
% |5 r% R3 @% g2 s- n2 w7 _+ _   
% S: }& |+ P8 M7 W- V) E   gen yxis3=r(min) in 1
   
; b0 k2 x8 L' o% M3 t' |. }* K   gen yxis4=0 in 1
5 m; B3 g- N) R7 l4 Y  
( q) i: B: H& a/ f. j- d  *ddy240g
   
   
  preserve
6 o4 l: x$ W6 c- \3 T9 S3 y  
; }8 f  ~' P* y  g! O  sum v5 [aw =avol]
  
  keep if dx240<r(mean)
  
  sum ddy240
% m+ X4 t0 I9 |) c9 C% I: |1 x- b  
  if r(min)<0 {
7 T* o) L# I4 b+ i* W1 o8 C( g- v        
% Q# k$ g0 |$ V3 b8 \. o" q3 A                keep if ddy240<=0
+ V+ _! g, B) O9 L" I" ~               
                count
  m/ I/ w6 Y, @* s0 I: f               
                if r(N)>2 {
+ a) G) t7 \0 |) e; [3 R/ s6 F                 
                 pkexamine dx240 dy240
               
                local tempminp=r(tomc)
                restore
                 }
) L$ z) W  }7 o" P                 
        else {
               
                 di "cannot find minp"
: V8 O9 b( k% e* p8 y/ z' P% T6 i! [                 
                restore
                 sum dx240
* F& n. ]: U+ |3 b* u4 ?! i/ f                 local tempminp=r(mean)
                 
                 }
) [  N# @5 X  t  f: R8 u4 F# p5 k; h                }
# q+ a, K: j9 ^' d- T3 h7 P   else {
$ k7 A# B1 \" `6 `; ^                 
                 keep if ddy240<=r(min)
& I. C0 y- z/ M2 M9 W. N                 
                 sum dx240
% W: b) a1 W& N  B$ z3 W                 
9 G: m% n8 s1 Y- w6 s! i+ Z                 local tempminp=r(min)
3 v/ o( P5 V" O9 m0 Q0 ?                 restore
                 }
( i5 U3 I$ ?( V- h* }

; h, c4 {. X' g0 ]2 T4 w' @, Q
   gen minp=`tempminp' in 1
   *240 minp
+ e* g' A* l+ V

   
0 P0 j$ @4 Z7 i. `2 d: W   preserve
2 t( Q8 I& Y2 Y$ A4 w9 o$ h   
   sum v5 [aw =avol]
8 _# o4 a" Y+ E: n3 {1 l3 L6 Q   
; ?. n* N2 H6 l' b7 X   keep if dx240>r(mean)
   
   sum ddy240
   
- S' ~. _8 x. I7 x& t3 ]8 m   if r(max)>0 {
               
) y7 i$ W8 h6 X( N0 U                keep if ddy240>=0
  D1 Q* q& ]6 N               
                count
6 T6 T7 M) k% ^. s6 V               
; D: ~! \( a  e) h1 s: c                if r(N)>2 {
7 ^" n! ^' b% b: w9 W. X+ \$ [                 
                 pkexamine dx240 dy240
/ D# _% q5 |% y) @& G+ z5 T& T               
                local tempmaxp=r(tomc)
                restore
: I6 F  }4 j+ W  A' i4 `                 }
                 
: Q3 }4 q; j  [$ M# B, E* X                 
        else {
                restore
                 di "cannot find maxp"
               
               
$ h: B+ d/ K6 @# `, c                 sum dx240
                 local tempmaxp=r(mean)
. F5 p  q" L% j% O                 
                 }
                 
                 }
                 
$ y6 A( H9 x# R1 a    else {
       
           keep if ddy240>=r(max)
/ A8 b' r5 k: K! ~& i& o          
. f! h* ?- q  \! H/ q& a           sum dx240
          
           local tempmaxp=r(max)
+ A8 h& }( n/ {: s! m           restore
           }
, E# y! g! g) V* Y4 G# ]       

1 a# r/ r1 b; ^$ p2 n  

   gen maxp=`tempmaxp' in 1
   
   *240maxp
4 V/ o- n. ]4 X' P5 ?   
   
   *-------------ddy240pminp左侧高峰值-----
4 W+ F. d+ f8 X0 A   
    preserve
   
   keep if dx240<`tempminp'
   
   keep if ddy240>0
0 p, F  a: h: r$ ~2 E8 l$ ?6 \   
$ S7 [# \3 j$ p8 ?- p& H   count
               
# N7 l# ~' J' ^8 t+ D                if r(N)>2 {
                 
7 ^# r5 `8 L% T/ t                pkexamine dx240 ddy240
& J' m& Y' z% h; A2 }  m- j( |9 O               
8 _. U0 m, G# I$ V                local tempminphddy240=r(tomc)
- |9 B8 I! ]8 E9 e1 }" a$ I, ~% L1 a% G                restore
                 }
/ y& M$ E: r) K0 J7 V  T                 
. r1 u) \5 D" e% c% p  W        else {
                restore
                 di "cannot find pminphddy240"
                 
& g* E# d) S0 p3 H& ]8 J+ q4 U               
                 sum dx240
7 O4 S( @4 A( Z9 c/ e                 local tempminphddy240=r(mean)
6 W& j1 j7 i) M  a' K               
8 [" h8 s: ]: s& |                 
                 }
# {8 z+ o$ i' u                 
   
   
1 H3 ^$ m. S$ Q7 U   
: k) B$ ?( W$ h& _% _6 E6 H  ^+ H   
' C7 L0 `& @' Y! W   
. x5 T  _1 t" G/ j9 s+ c        gen pminphddy240=`tempminphddy240' in 1
2 K- F* P. S; y7 j- c, X. C. E/ c       
       
        *-------------ddy240pminp右侧低峰值 -----------------
. A& ?+ T' k* V( a    preserve
. f$ J2 n) y3 X& K
1 a, }( ]# w) I   keep if dx240>`tempminp'
   
% J. ?8 _# e1 Q  d6 |0 f   keep if dx240<`tempmaxp'
   
   keep if ddy240<0
+ k/ q$ Y3 {; t' ?$ x   
) Q$ C% g/ i( l7 f; j# A   gen temddy240=-ddy240
* [& R: T0 ?! N) q: b' S8 c: v, L   
' F4 [8 @. R5 L: H) C    count
               
) L5 t; v4 e- Q6 Z0 A2 y                if r(N)>2 {
                 
                 pkexamine dx240 temddy240
3 o: w% l3 V9 H4 S( @7 `# |   
" Q8 r; k3 `+ w8 ^  local tempminplddy240=r(tomc)
                restore
: r- U3 h( g: U& V$ n# B                 }
( d; b! b/ x) w1 J. O3 n                 
        else {
               
* ]; W  _3 `( o2 A% T1 [                 di "cannot find pminplddy240"
* t( ?* q4 U( R. E7 ~- v' U                restore
                 
& K( V7 U( ?( C+ p7 Y                  sum dx240
+ h0 K% n9 |+ ^& E8 @                 local tempminplddy240=r(mean)
* X& }3 g; d1 |. j; k                 
7 B9 h0 p& u$ W- _, J: O0 l3 \$ r                 
                 }
   
' W- ~8 j1 u: v" B1 [# e8 v   
   
   
   
' Q% k5 J3 G8 ]; _/ h3 U* c        gen pminplddy240=`tempminplddy240' in 1
       
$ M- z" H/ D# @$ z% s& }0 z        *-------------ddy240pmaxp左侧高峰值 -------------------
; M2 A1 }/ p! e2 p7 [       
         preserve
6 ?& E/ m! z/ Y; `( a+ F- [" z! A, S5 }   
/ q. Y* M1 W( H& ~" G3 l5 K   keep if dx240>`tempminp'
   
   keep if dx240<`tempmaxp'
" ^. N! V; w- g* x) s   
* ~2 b8 d& L4 _# t   keep if ddy240>0
   count
: J/ w2 v) h; p3 k) S) X               
7 E+ k2 R) ]6 ^$ \9 _' r0 K1 h0 w' \7 Y                if r(N)>2 {
                 
                 pkexamine dx240 ddy240
   
   local tempmaxphddy240=r(tomc)
7 E+ F( w' P' t, j3 o+ b   
   restore
1 k, g  Y% A# Z                 }
) ]3 z0 n8 O( V( C/ p" z/ q) ^                 
        else {
2 A- _( z& D& ]               
                 di "cannot find pmaxphddy240"
/ ]; T0 e2 d' Z                 restore
  c) N, t* b  ~$ G               
                 sum dx240
! T- n7 ?/ H8 e) B* |/ l# `                 local tempmaxphddy240=r(mean)
) E, O0 m/ H' W# G4 C. @2 I                 
5 ~& r3 S3 s8 S% K: ^                 
                 }
   
   
+ E8 C+ ~: R9 ^4 \$ q6 n0 I   
   
/ V4 N- j4 @+ |$ k   
   
        gen pmaxphddy240=`tempmaxphddy240' in 1
       
/ l7 A, k$ f( q       
3 }: X8 S! o' G% J  c        *-------------ddy240pmaxp右侧低峰值 -------------------
       
         preserve
  
! i9 |: ?0 X# u* h2 d2 N. X   keep if dx240>`tempmaxp'
   
   gen temddy240=-ddy240
5 S( s6 b5 y, D" c4 L   
% B  `7 }2 x2 o& g7 D8 F   keep if temddy240>0
   
! i/ p* }' P5 R2 ?  }0 L  U7 v   count
               
. H( l: P8 s! E4 H/ G  q                if r(N)>2 {
: @5 j7 R1 s1 `# _0 B- U                 
                pkexamine dx240 temddy240
2 z. e, ]/ K- }   
   local tempmaxplddy240=r(tomc)
$ H0 U% N9 N1 P/ w3 {; B   restore
# b- l+ f0 y" E$ u: r                 }
. h6 ~) ]& t2 ?9 d                 
2 t. T* K+ g; {- y9 m' ^        else {
: \8 E4 {" T/ _2 [' ]                restore
) ^5 {7 p2 P# g! g                 di "cannot find pmaxlhddy240"
                 
( s  _. h0 u# G& t               
                           sum dx240
                 local tempmaxplddy240=r(mean) in 1
               
                 
                 }
   
   
   
& W5 C2 \7 t7 @   
% H0 X# W9 ~/ `* A   
. O' o, M5 c3 i& q   
        gen pmaxplddy240=`tempmaxphddy240' in 1
       
5 a% K; Y2 h: T& r/ k        *-------------------------------------
   
   gen price5d=v5 in 1/20
/ f' e* E* ~& J7 W      
, y( v- D. R$ n9 [. S   gen price5yy=_n in 1/20
8 Y, {9 K# @. ~7 f7 V   
   gen price5y=price5yy*(dy240max-dy240min)/20
. y4 E2 o" b! u' T   
   *price5y &price5d
   
" Z: A3 {: i/ a   gen priceny=price5y in 1
   
) w$ _0 v$ h0 H2 N1 H/ L   gen pricen=v5 in 1

  c5 I, H" W$ K& O     
   *priceny &pricen
   
& `, [9 X& t5 f   
   *-------------------------graph------------------
7 v$ }" q* V" _! p  s  
  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)

* ---------可选项目----------------------
*(spike yxis amean,lc(blue))(scatter yxis2 amean, lc(blue) msymbol(+) mlabel(amean))
$ Z1 g! E; [: F}
: e: k( n- k  K4 o$ b' _end
- b; W9 D$ {- X

maodeng

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

# z* t& Z: y2 o0 t9 Y( F300236 我用下面的程序试过，没问题

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

喜欢喝冰茶

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

喜欢喝冰茶

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

( `9 f  Q+ ^& M4 [# V9 }
这东西是开源的？知道有个统计软件叫stata，不过要收钱。

汉水云潇

准备试一试
先找Stata软件

汉水云潇

, ?& F  ^" f( k1 S- `8 F晕呀
论坛复制程序代码
; E$ _* [8 Y; W* i( e! G乱码夹杂在其中
很烦
对了两次
以为没问题了
, I% P+ m$ v" V8 ?% f, s
结果
0 ]" z( ]6 Q: [/ g4 W/ z& K执行后
还是出现
r 199
% P: ^1 g7 I. v) K8 r, C的错误

汉水云潇

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

peterbev

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

中秋不圆

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

bknite

这个帖子值得学习。

jeanny

感谢！！！

Hubert2000

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

, V; g7 q- `* I  j, B老哥，也是山寨这篇文章找不到了，不知为何？