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

墨轩书童

) s, J# ^) f8 `感谢楼主！！

海阔

2 J( Z$ k3 P3 X' `
按照版主的步骤操作到了最后出现
) L7 h. c! {& A2 {
. hello 600031
; r+ [7 X  |0 U4 y$ z" F7 I7 ZavolDo not found
---------------------------------------换成其他代买，也是一样；总是出现这个问题，怎么办呢？奇怪

maodeng

大师，那两条曲线代表什么啊

maodeng

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

maodeng

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

maodeng

. hello 300236
statement out of context
3 C$ ?9 O* O' `' H- M" Ar(119);
& s9 \2 c0 R+ W! ?2 q
大师 就是这个

sergio

这贴实在！

awer

maodeng 发表于 2013-1-28 13:13
. hello 300236
( `; t3 a& o% X1 S1 a7 Rstatement out of context
' j- ?* k) F( {' H: X4 r$ br(119);

5 ^/ b% |+ c' W5 f! i2 J& p: c抱歉，最近太忙了

300236 我用下面的程序试过，没问题
4 i5 n6 K( ~" ]# I* m5 U, \
你把 D:\stock\test\source 改成你导出的数据存放的目录，再试一试

+ h$ I- v1 n* w( Ecapture program drop hello
  ?% A- G9 M1 M0 H% N. M
/ d6 s7 ^0 Q' y, {2 G9 L) J/ u, Q+ uprogram hello
" K( R8 K) U+ T$ Q7 h# q# e
6 B7 p" w/ S3 D( f7 A& }: n3 v6 B; @ quietly {

. o: a( @9 \8 K+ `9 gclear

cd D:\stock\test\source

insheet using `1'.txt
( u6 M# U0 |7 e% H
* y! F5 z3 f# U) v" P. }   keep if v5>0
   
   gen vt=_n
   
   tsset vt
. t; f: J" t; k" U0 `
   gen zf=D1.v5
* l+ [. ]% J! J
: P( E5 \9 E0 |" I   gsort -v1

; \+ P, z6 R( v; o: v: ?7 B4 t   gen avol=v7/v5
   
   keep if avol>0
   
   drop vt
+ V4 x9 D% l7 m: u, }6 |  
   gen vt=_n
   
( N# H+ j$ p& Q; ]/ u5 H3 S   tsset vt
4 T" G) H! L9 Y4 U' o   
    sum v5 [aw =avol]
7 c/ [! t6 e; X1 w. }( r
   gen amean=r(mean) in 1
   
   * 总加权平均值 amean
! f7 C: L2 \9 E5 i" e& ^   
   
  R7 I* i  A  q+ v' |   keep in 1/240
& Z" q8 |1 n% r. Q) W
* `- @6 y+ [3 q' r+ j) k( P   sum v5 [aw =avol]
0 n7 n7 T6 N1 g. p8 a3 q9 }
2 I' Z) C! K5 f   gen mean240=r(mean) in 1
) c8 t( @) N/ `/ Y; k  u' `% s! S4 r   
   *240加权平均值 mean240
   
   set obs 500
6 i$ e9 [, e: e! {$ x+ D! y9 J   
- u7 n: ?9 C/ ?8 H   kdensity v5 [aw =avol],gaussian n(500) generate(dx240 dy240) nograph
, v; D2 C6 t* X7 h   
   *240 kdensity(dx240 dy240)
   
1 P% L' a0 a: {7 a: h9 b7 A   sum dy240
   
   gen yxis=r(max) in 1
* E  z$ v* u# a   
( z/ S- v+ ~' h1 P7 e. P   gen yxis2=yxis/2
   
1 K+ S/ b2 n: L+ {0 B$ n; t; T* X   *yxis
   
: ~, E6 S+ f, U( z( o7 [' n   gen dy240max=r(max)
   
3 P8 n& X2 D* L   gen dy240minp=r(min)

4 w3 _+ r3 _( \! B. F/ v    gen vtt=_n
$ C* @6 k5 }7 v
; B/ Y4 _" P# }% [7 A+ @    tsset vtt

   gen ddy240=D1.dy240/D1.dx240
+ J* ^. b7 A9 [! M   
   sum ddy240
   
   gen ddy240max=r(max)
/ L. i2 l1 \# B& j  q   
   gen ddy240min=r(min)
2 {* V' c7 ?" q4 t     
   gen ddy240g=ddy240*(dy240max-dy240min)/(ddy240max-ddy240min)
, j* W( C% X  G+ w% n5 W, ^   
   sum ddy240g
   
   gen yxis3=r(min) in 1
   
   gen yxis4=0 in 1
. o& f6 d2 _5 S* J% S& Q  
  *ddy240g
   
   
  preserve
8 N# h" `) X  Q: L$ J  
  sum v5 [aw =avol]
8 _9 Q4 n6 y( u% A* o  
  keep if dx240<r(mean)
  z' [, F4 g! S4 ~  
  sum ddy240
1 k$ n* P2 U5 ?+ p4 I5 \6 E  
0 d3 R4 H4 ]& R  if r(min)<0 {
8 _3 @: u- e1 F; v% k0 K: G4 z        
- W5 f2 o3 h; Q. |" K) ]                keep if ddy240<=0
               
                count
! z1 O4 ?% X4 F# ^( e  a; J               
                if r(N)>2 {
                 
                 pkexamine dx240 dy240
$ ?4 U9 h- Y* @& x8 g               
                local tempminp=r(tomc)
# h2 I3 D( T# |8 U                restore
; M$ g. T- s5 [% `2 z                 }
                 
        else {
               
                 di "cannot find minp"
3 r0 V/ a% ?* M/ q                 
4 N7 G7 q$ J; A7 k* K7 |. A                restore
+ ]3 Q, f, X" ?+ ^                 sum dx240
. h- v1 c- z( B3 e! G4 P- f                 local tempminp=r(mean)
                 
0 e  [3 c) v! i: \) {& Y% U, V: w                 }
- i2 l" C' Q8 I* _                }
+ \0 p6 a# }5 G+ z   else {
# S) c% Z) A& i: K                 
! T  [: o) B+ u; k7 _2 Q                 keep if ddy240<=r(min)
                 
$ K5 `1 @- y8 U# T) O                 sum dx240
                 
                 local tempminp=r(min)
5 q$ g( X1 u8 ]! F8 T                 restore
                 }

+ |5 I" T" u$ N' p
8 W9 M# S* q' L. q/ N
   gen minp=`tempminp' in 1
   *240 minp

9 ~3 M9 x! V' P
   
   preserve
   
9 G  Q1 G+ M$ e% S0 E   sum v5 [aw =avol]
   
   keep if dx240>r(mean)
! |, _- M" a6 |   
   sum ddy240
" L7 g# c. _4 W1 u  B   
   if r(max)>0 {
               
                keep if ddy240>=0
               
                count
               
                if r(N)>2 {
* }' v7 h7 R4 P, J3 V% m; y$ [                 
                 pkexamine dx240 dy240
6 x9 q6 b' ^0 A" w4 A2 F+ i               
                local tempmaxp=r(tomc)
1 h/ g8 Q3 {7 C$ Z                restore
                 }
                 
8 {" ]% `/ d1 J& I1 w# P) f/ ~3 y                 
- M3 }7 r5 }: i1 n* }        else {
                restore
                 di "cannot find maxp"
; @+ @3 g5 m9 k. X" Y: I4 i               
               
% S& x4 ~  j% W! Z9 v6 j- E, w                 sum dx240
* g' p0 {! b. p# {/ |7 ]! w; R                 local tempmaxp=r(mean)
# Q! w5 ^2 i9 u& Y4 C                 
                 }
% A! J0 |) ^+ x5 Z8 L                 
                 }
0 T* }. ]) \5 J) h) Z                 
. e3 d' P$ a5 ^! |0 g4 ~    else {
- b+ S6 H9 U; \: @       
           keep if ddy240>=r(max)
          
           sum dx240
          
- d! o6 n( r- O8 H) F4 G           local tempmaxp=r(max)
           restore
           }
       

4 c/ @6 b2 L6 L, s: L( P$ |8 k  

& l) m+ S7 k$ s# j4 z3 b   gen maxp=`tempmaxp' in 1
( `: C  s, E6 T: G1 n   
   *240maxp
, c6 E4 h7 X! f$ @' l2 ^   
   
) I, g& s3 `0 L& t   *-------------ddy240pminp左侧高峰值-----
   
8 N, v9 _6 x" t/ ^* v) k    preserve
   
   keep if dx240<`tempminp'
   
3 P& A/ z# I2 I' l   keep if ddy240>0
9 L8 w0 u, ~1 \) {5 m; u$ T   
   count
               
6 T+ K- ^9 k: f6 h- w- @& x                if r(N)>2 {
                 
$ m* e; h- Q; e: X( \1 N8 t                pkexamine dx240 ddy240
9 L& `9 k' u. _0 X+ F4 g               
# n7 k, Y% u! L0 l# I. G                local tempminphddy240=r(tomc)
. ~0 e( ^) L* z+ h7 k                restore
7 Z; q" K1 Z  X" w                 }
1 M8 E3 E% T; P: A                 
        else {
; H+ v0 e8 Q' Z1 f+ B                restore
                 di "cannot find pminphddy240"
8 {; ]$ r0 D; k) E- d8 F                 
               
$ g! c* v8 _5 V) m+ _                 sum dx240
                 local tempminphddy240=r(mean)
' F, d3 b; r$ ~5 a               
                 
                 }
                 
   
   
   
   
' Z1 r3 [4 H, Q: M4 ~% ~7 k   
. R3 _2 |4 Y; p6 A+ |7 [# K        gen pminphddy240=`tempminphddy240' in 1
- i$ K$ \& l! O+ f       
       
        *-------------ddy240pminp右侧低峰值 -----------------
/ ?) B4 [2 o0 d    preserve

6 {8 e8 e" A8 C8 D6 V7 R1 J$ I   keep if dx240>`tempminp'
# K: Z* S' ~/ ?   
   keep if dx240<`tempmaxp'
2 ?1 J# S& P1 d# O+ }, ?1 ~   
; c. `3 O( S, Z0 ?: B   keep if ddy240<0
   
. v8 O( e8 x. [+ x   gen temddy240=-ddy240
   
    count
               
1 a  }7 @0 X; l) d+ i                if r(N)>2 {
* c. C& r( q1 ?# `3 t                 
                 pkexamine dx240 temddy240
4 V, f' b- z2 K9 s1 Y   
; w8 h# o# d1 o2 C. n& e1 K  local tempminplddy240=r(tomc)
% Y  `  j6 F: s* `, b# y                restore
                 }
6 P8 D/ D* G$ @/ \' B9 O" o                 
        else {
% O) m1 @& s/ Q! y( B+ C               
                 di "cannot find pminplddy240"
5 I4 Q' Y  q* c0 n: M( b                restore
' E. y$ Q6 p6 T                 
! D" I8 g0 s) F6 n                  sum dx240
                 local tempminplddy240=r(mean)
/ v& s) ~/ O  _6 S                 
/ I1 C3 o9 U  l( `                 
                 }
   
   
3 }  ^+ G4 m3 ?5 }   
7 k9 r+ N3 J# A9 w, P   
   
        gen pminplddy240=`tempminplddy240' in 1
3 Z" Z6 h* S9 Z3 r, W       
        *-------------ddy240pmaxp左侧高峰值 -------------------
       
         preserve
3 K4 s! Y& J" x6 d. ?5 W   
% o, u( A2 B" j  M& R   keep if dx240>`tempminp'
  M6 ?3 t! s2 k! S7 ^   
   keep if dx240<`tempmaxp'
  x$ ~6 d* d+ d* a: i# X   
3 J* g0 G! ^. E6 n  @/ W/ H5 L   keep if ddy240>0
   count
$ z: }9 Y7 _5 g0 P# H# }" D( A               
+ W0 T. }* e. K% y                if r(N)>2 {
5 u6 v( Y* P# r8 G$ P" Q                 
2 e7 t, _# @  x$ t! X9 ]2 }                 pkexamine dx240 ddy240
   
   local tempmaxphddy240=r(tomc)
8 N2 u7 I2 G) Y; l/ {   
   restore
, U/ E8 y' Q$ y4 d                 }
                 
% O5 v2 t& X$ `' }5 s7 t+ b3 M( v' B        else {
               
                 di "cannot find pmaxphddy240"
                 restore
( z2 w! G/ e/ j, W. ^' z# t               
                 sum dx240
                 local tempmaxphddy240=r(mean)
                 
                 
( g+ B/ d5 x' l* T                 }
   
   
   
' `! D" P( F' D4 O' O   
% h  ^  W0 J+ l$ C   
, @. d0 [" M, t   
        gen pmaxphddy240=`tempmaxphddy240' in 1
* s( s4 T$ V2 u  \       
+ a4 [& H" V# b% s" B       
        *-------------ddy240pmaxp右侧低峰值 -------------------
2 U, F: `! T5 l, j7 {1 ^       
/ X' ~1 T8 q! O+ H  v/ W         preserve
! e8 u- u4 I' k; ~- V  
* }0 I# m% }( G0 B) {4 `0 E, k   keep if dx240>`tempmaxp'
   
7 v) J5 H, d7 ]" J' q   gen temddy240=-ddy240
   
4 u' D  S/ @  W; a% g9 t! Y   keep if temddy240>0
   
   count
               
$ t6 Y2 o3 r6 M% w5 m0 k                if r(N)>2 {
                 
' C3 j  l# v. {9 N$ ]1 ^3 a                pkexamine dx240 temddy240
   
4 \- Q6 m' R' ?- o7 T8 O   local tempmaxplddy240=r(tomc)
   restore
" I/ U& m( e$ j. N6 M  Z2 [                 }
2 \1 k6 w+ {& V: Y# p/ O0 A                 
        else {
                restore
; s" K+ L* o4 `! Q                 di "cannot find pmaxlhddy240"
2 F" L! s$ T9 i6 F% A                 
               
                           sum dx240
% g; G- [' U  ~/ l" E                 local tempmaxplddy240=r(mean) in 1
               
7 d& L5 o7 k* d. l# t) Y                 
                 }
   
   
: d! H! @" `; ]. c   
   
- U& w: H  v4 `# z/ K1 B   
   
0 Y2 ]6 v* I8 @8 j  R        gen pmaxplddy240=`tempmaxphddy240' in 1
       
$ j% \' `* k; |% e5 y) B- O- _1 H        *-------------------------------------
   
   gen price5d=v5 in 1/20
! l( [9 {, o+ @  x/ V4 q      
; ]; g) [) j# K   gen price5yy=_n in 1/20
   
   gen price5y=price5yy*(dy240max-dy240min)/20
+ v5 l% }8 R/ u/ W   
0 `2 t# V& f. Q1 l, j   *price5y &price5d
$ H$ I& o9 U8 N2 i# I' s. A5 X" W$ |   
   gen priceny=price5y in 1
   
   gen pricen=v5 in 1
. A; m7 g0 l$ g  v
/ H2 k+ w$ F* ~+ F% V     
   *priceny &pricen
  Q  t( g2 f4 o) i8 v$ u& {: q   
   
   *-------------------------graph------------------
, Z; s8 }; ^7 f. w  
+ o, k  G4 O/ w) p" h% i  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)

* ---------可选项目----------------------
6 f9 i% [% v' k1 I2 g9 }*(spike yxis amean,lc(blue))(scatter yxis2 amean, lc(blue) msymbol(+) mlabel(amean))
}
' l- ?9 p& t+ @. N! aend
% ~3 ]; Q- r  K- ]& G3 i

maodeng

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

3 ]3 p8 u4 K  o/ _5 ^2 B2 }3 Y+ G0 h300236 我用下面的程序试过，没问题

6 s2 G* @+ {6 T$ t* i哇，谢谢，我仔细看看，这些字母数字代码晕了

喜欢喝冰茶

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

喜欢喝冰茶

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

$ l+ {0 ^  o1 \* [" L0 R
这东西是开源的？知道有个统计软件叫stata，不过要收钱。

汉水云潇

准备试一试
先找Stata软件

汉水云潇

6 |6 _; y8 c2 ~; q  B. C+ M8 O
晕呀
; V3 Q) O7 a( }( V! Y1 b& [$ X论坛复制程序代码
- f( c; M# F8 m乱码夹杂在其中
3 ~6 ~: _0 J& x9 d8 P; ~5 |9 B很烦
! Q% w# k. O1 B对了两次
以为没问题了

结果
6 _# E0 v: I4 V7 ^3 d8 G; l8 A执行后
/ Z* `% _4 t. {( D, |+ G还是出现
1 `0 K( l7 A+ s* L/ _* e; N, _r 199
! _& ^) p/ Z6 ^$ b( g; P) J- u的错误

汉水云潇

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

peterbev

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

中秋不圆

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

bknite

这个帖子值得学习。

jeanny

感谢！！！

Hubert2000

明月回春 发表于 2011-8-13 20:17
! F7 G9 o- R6 j4 e) @4 y程序写的不错.
  B: Q5 S7 L' L不过有一点要注意.通达信导出的数据,只对价格复权,并没有对成交量复权.所以数据需要手工处理 ...

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