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

maodeng

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

maodeng

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

maodeng

. hello 300236
statement out of context
) o3 r: q, ]1 u# r" ]/ J! e2 Pr(119);
8 o9 v- x$ @. l4 O# K& G
大师 就是这个

sergio

这贴实在！

awer

maodeng 发表于 2013-1-28 13:13
. hello 300236
, u' F' l" \& u" hstatement out of context
" i1 a/ ], l3 K% S  z/ Cr(119);

0 u1 ?: e* Z  `1 ?抱歉，最近太忙了
/ O4 Z2 R; P' h- |- L" B  c
) h& O) a- `* K, C2 N) ]300236 我用下面的程序试过，没问题
9 L3 y: k+ [) a& _! y$ U4 i
1 U1 g3 c8 S. f( ~你把 D:\stock\test\source 改成你导出的数据存放的目录，再试一试

0 V5 t5 _9 `: Y9 a% q) T+ Icapture program drop hello
& a$ w" c& ^* J$ ~$ ^6 M
program hello
: r- M2 B) J9 W3 v
quietly {

- y5 M4 U8 ~' l: a/ A. kclear

1 L6 h6 r  N; \& p8 L# scd D:\stock\test\source
' n$ j8 M4 v, v% \: x3 Y/ N& }
insheet using `1'.txt
6 x; b' G( F& _4 ?! h% D
: ?+ ^' [' e3 ~. a. u   keep if v5>0
+ i( e% C& e" j! V   
   gen vt=_n
   
7 j+ I( U3 O9 K' I   tsset vt
- U5 J3 u6 y8 }" b
: z  X. g/ l$ l$ E% g   gen zf=D1.v5
# n. g; H& J# P
6 C2 I9 H0 X$ m& |% @) u   gsort -v1
  Y, M9 h2 H' |2 R6 l0 J
   gen avol=v7/v5
! T2 @& I4 C0 O  ~6 l   
   keep if avol>0
   
$ ~1 l+ ]8 W. L   drop vt
* P3 y9 b( d0 |  t  
$ \5 _0 @* D" H   gen vt=_n
   
   tsset vt
   
, G8 Q7 L5 O8 @3 }3 a    sum v5 [aw =avol]
" a% \1 l4 ~% P9 \& J) f- O) W
& M; [+ S: v! Q   gen amean=r(mean) in 1
: ?- h! D6 l0 K  e   
   * 总加权平均值 amean
   
8 }/ {. m9 r# ~3 u   
   keep in 1/240
% J' s+ W2 t% {8 N% |
   sum v5 [aw =avol]

   gen mean240=r(mean) in 1
   
   *240加权平均值 mean240
' L& {+ ^0 ^5 s/ l   
/ y: q+ n1 V" R   set obs 500
2 `! j5 n0 {% k! a& O, |   
   kdensity v5 [aw =avol],gaussian n(500) generate(dx240 dy240) nograph
   
( {: P5 }/ H' P- S# ~8 i   *240 kdensity(dx240 dy240)
: e; e  k( j9 O3 R; r5 P$ {1 ]$ I   
   sum dy240
# m+ G" l7 R4 Q; v9 r' o   
7 N6 K/ a2 F" b( {   gen yxis=r(max) in 1
* h, a* e  s! S2 d; W/ J- Y% F   
# K* @; E. T) G1 _5 A6 U   gen yxis2=yxis/2
$ Q8 A$ a8 y. S( O$ Q   
% }: T1 d, N$ ^   *yxis
" u! v& G- c8 p" N0 H   
   gen dy240max=r(max)
  W: o+ S+ z$ @( V   
# j4 V  T: @+ ]   gen dy240minp=r(min)
( ~7 D: \- Q* K
8 A! C% B4 s  _    gen vtt=_n

+ X8 V: T8 q/ C) j  s    tsset vtt

   gen ddy240=D1.dy240/D1.dx240
   
3 ?2 p' m: Z( w9 i1 h8 A. E   sum ddy240
  V$ M$ V% U- ~, K8 U   
   gen ddy240max=r(max)
0 @0 m6 Y4 G/ i9 }4 j  l& I% J' F  R! b   
   gen ddy240min=r(min)
     
! i# z) p! h9 r% Q$ q  r" X) ~   gen ddy240g=ddy240*(dy240max-dy240min)/(ddy240max-ddy240min)
9 h' |9 T* o' b1 B4 l8 `1 ^# q   
+ m: A# a9 V- N5 b$ Q. _( n2 z3 E   sum ddy240g
   
- a/ b  h) [& M& d( K( [% }- A   gen yxis3=r(min) in 1
% X" c4 h8 F0 U2 X   
2 H( _) l$ p/ M: y   gen yxis4=0 in 1
( Y3 @4 X. n( x2 }+ v) ]9 U3 B! h  
* d+ s: q6 M, x% O0 S: g  *ddy240g
   
1 e$ U1 ]1 P  {   
  preserve
  
, h1 J0 E' E, N2 ^# G0 l( L  sum v5 [aw =avol]
  
  keep if dx240<r(mean)
9 w: x$ N! r8 n1 u3 u  
  sum ddy240
$ }6 P8 a& W+ W$ Q  
5 f, x7 ?( \/ q" |! f$ D0 B  if r(min)<0 {
2 ~2 R! G/ v' |4 [( r        
- O8 O. N: y) d7 s) e7 C                keep if ddy240<=0
  u( @0 |2 Q$ g: w. V; p               
                count
% P" M( A) M! ]5 w               
                if r(N)>2 {
) L% i7 W+ E" `/ N' L                 
                 pkexamine dx240 dy240
               
                local tempminp=r(tomc)
/ ]' f( Z$ P1 s" t$ |                restore
                 }
8 k/ w/ @$ s: ^2 c5 H. @: [3 T3 Y9 |                 
6 v4 @5 ]3 }" }1 M7 p& Q  B' I' E        else {
) K0 D9 P* O4 a' Q               
, u% C( b+ e- K+ ^                 di "cannot find minp"
! @. P: r9 S* s2 \8 n                 
: y$ M/ ]7 m. z1 ^1 l9 o8 q0 L4 s                restore
                 sum dx240
( }  B# O# p7 }9 [0 [; e2 S0 \/ _                 local tempminp=r(mean)
                 
                 }
                }
: Z! _3 F7 f7 v# t0 U0 Y5 ]   else {
9 \: ?7 c; L8 a9 M0 G                 
4 F. {" q  P, T7 e' C                 keep if ddy240<=r(min)
$ D. F8 n2 P4 ~                 
* |$ ]' o1 p( d4 m3 O. S                 sum dx240
                 
- E0 [0 ^1 S- B( _; [: Q                 local tempminp=r(min)
7 E: v5 d0 J. ~, m+ _, }                 restore
/ P/ v( P- J- v                 }
0 V7 c/ T" z3 o: z
$ h/ y! ^% ~& v, N+ a& Q

) a' E/ w3 J7 V  C5 B7 N1 L( v   gen minp=`tempminp' in 1
   *240 minp
" ]2 ^  q# W7 F/ ^7 Q5 x
& H4 a% \9 }2 Y8 V- h
& }2 z0 c: @7 ~   
   preserve
. h+ N& ^) P6 P. ^# f: r   
   sum v5 [aw =avol]
. T! h# v8 {7 b5 _   
   keep if dx240>r(mean)
   
. c+ P0 i! b. m2 w2 K; o2 Y" Z   sum ddy240
! `# E. b+ ~" B3 X9 \; @6 Z   
   if r(max)>0 {
               
                keep if ddy240>=0
               
                count
               
" G$ u: U- q7 Q                if r(N)>2 {
                 
+ o& g& ]- n9 }6 q! T  |                 pkexamine dx240 dy240
3 Y* k* e& F% \4 _  Q               
                local tempmaxp=r(tomc)
                restore
3 [, {7 }; V! q2 V9 E& e1 ]3 P                 }
                 
                 
% F8 c$ H# @; z2 ?5 F        else {
                restore
                 di "cannot find maxp"
               
               
                 sum dx240
3 l% l% d2 C$ v- Y" {9 j. q) N                 local tempmaxp=r(mean)
; x& F( A" `- C7 K; T                 
# a0 k& p$ B1 Z* f+ m                 }
                 
: e* f/ w) V, b% A* F                 }
                 
    else {
! b) y  X6 z5 \% z% H' f6 V- K       
           keep if ddy240>=r(max)
          
, S9 T: Y/ i2 w6 d0 `( l  W           sum dx240
          
           local tempmaxp=r(max)
" K* B2 x% H' r( P           restore
0 E) r7 R: m7 @' X2 o0 N( X           }
2 w. B" _8 u& W2 K/ ~2 \8 e/ j       
/ [- h" ?$ ]2 a/ c, ~
: g% t5 _( r5 A8 B1 ?  
* N9 I" p2 @; u
* y; p& n& Y. L) q   gen maxp=`tempmaxp' in 1
   
2 q# F# R) D* x, ?1 Y! M+ @5 w   *240maxp
   
. V6 m& f' P. j. T' }6 e+ y  y   
2 ^* H! Y3 `* L   *-------------ddy240pminp左侧高峰值-----
& M' A6 l( X# b4 q7 j   
5 [; E; a' g8 D5 t    preserve
   
   keep if dx240<`tempminp'
2 ]6 M0 |+ ~, d1 a2 B$ G7 U   
- F* d4 j7 o) }* B- B   keep if ddy240>0
   
   count
               
                if r(N)>2 {
* y- H3 X1 S# _& I/ j                 
+ @' Q( \& h5 u) x, R                pkexamine dx240 ddy240
/ F, U. o8 \) r8 I7 _8 C1 F               
8 i, }. F+ A5 ^; q                local tempminphddy240=r(tomc)
' _, l1 r: v( O$ E                restore
9 j) k, J+ L- l) h; E3 a3 |7 g                 }
                 
        else {
" z7 d+ k! u7 j3 i4 f                restore
                 di "cannot find pminphddy240"
2 a% u5 G5 [" l3 c& j* x9 E# l( A1 v                 
9 H9 b8 [+ y. P0 x5 j! p               
                 sum dx240
                 local tempminphddy240=r(mean)
               
1 |) _9 [: W) M4 o$ ]' H                 
2 L% e: C3 u% k                 }
                 
1 z' K6 c( {6 p2 @   
   
" H! t1 f' y6 D) K% N. L2 Q   
   
1 `/ b$ g! l4 \4 `! g1 E) X   
        gen pminphddy240=`tempminphddy240' in 1
- n6 S( ~8 w% ]8 |+ z8 P       
       
        *-------------ddy240pminp右侧低峰值 -----------------
9 q! l6 |+ r- H( P% Z; f    preserve
9 z7 ~2 w3 f" t8 k( ^9 S
9 u0 c  B/ k7 f1 y/ U   keep if dx240>`tempminp'
   
   keep if dx240<`tempmaxp'
   
1 N4 \0 c; X0 B6 q5 R   keep if ddy240<0
% R0 e! I1 P- e2 ~' Y7 x8 O   
' N# Z% \$ k3 S% y' V   gen temddy240=-ddy240
   
    count
               
$ y7 G" E7 t# w& l- J8 S. {8 m                if r(N)>2 {
; i/ L5 M$ V, R3 P" t                 
# w, V' d' J! c                 pkexamine dx240 temddy240
9 `; X/ ~6 }. C$ B5 d  j) z   
8 u. I* S! i' \5 X  X4 q) F$ h( k  local tempminplddy240=r(tomc)
                restore
                 }
                 
: [; R3 `2 `+ `        else {
               
9 E' @& J' u. `4 o9 ~* E                 di "cannot find pminplddy240"
0 d3 W6 z: N0 \/ H! c5 t- h                restore
5 I4 r; T/ n& T* b8 Q2 y                 
                  sum dx240
; H. \8 b* B* F5 R5 P* z                 local tempminplddy240=r(mean)
                 
                 
3 w% Z7 {( L8 Q+ W                 }
* H1 G9 S3 h6 Y" I8 P( @   
: O$ H+ t/ a+ M7 X, W; n% b. U9 N8 v   
   
# [/ T% {, u+ q2 F" P( K   
  l' C% H6 H5 J! j* @: H   
        gen pminplddy240=`tempminplddy240' in 1
       
# R0 ^! D* E; _4 P4 b* |& a        *-------------ddy240pmaxp左侧高峰值 -------------------
! B5 J- c& ]. N4 q7 {       
         preserve
   
   keep if dx240>`tempminp'
- J' V/ p2 G8 I" ^* k! {   
   keep if dx240<`tempmaxp'
   
3 Z. ~$ M: T- `9 ~7 @- G   keep if ddy240>0
" \2 p9 z! M8 k* m; Z   count
               
                if r(N)>2 {
& y  X6 W6 a. X/ x: @                 
& Y' `; e- I. v7 t6 v9 J                 pkexamine dx240 ddy240
   
$ p4 _! [2 I$ P( P3 Z. z( S. p   local tempmaxphddy240=r(tomc)
$ G+ B) K6 K2 X# J( E$ @) l   
   restore
                 }
: h+ m+ f9 v% M4 _; j- E& Q                 
! O; U" }- s7 L: N+ j+ y" q9 `) ]8 R        else {
               
                 di "cannot find pmaxphddy240"
                 restore
, W% ~1 D9 b( p9 n               
6 s4 O$ z, ?3 i+ s                 sum dx240
8 D# a( X( q. h; ?" t. N& X3 p2 o                 local tempmaxphddy240=r(mean)
( @, M* U0 V5 T4 d! F7 q                 
                 
                 }
   
   
   
   
/ b+ F: m' X; H. o" r" q   
  x) d: ?& L( q' V   
        gen pmaxphddy240=`tempmaxphddy240' in 1
       
' j: J2 [! z) v9 u- N1 `       
' e, r7 u+ N, _        *-------------ddy240pmaxp右侧低峰值 -------------------
2 E  d3 u4 \' }       
         preserve
  
: {$ j* l9 Q( L! {, Y) U+ J) u   keep if dx240>`tempmaxp'
4 Q' y4 E- ?; `+ X0 x+ ~( g, F   
, d6 N; q& r% W# m5 b   gen temddy240=-ddy240
9 I3 C( A. n% [   
   keep if temddy240>0
   
+ Y! T5 w$ v* E6 g   count
               
- C' C; ~9 f7 w& }2 O0 ~* `                if r(N)>2 {
9 {8 u5 G0 ]7 U- [                 
                pkexamine dx240 temddy240
9 ^% j) S& h9 ]$ \$ V3 z   
% ^) ?+ ~3 B2 \! L   local tempmaxplddy240=r(tomc)
4 @4 Z, x6 `1 r2 ?/ x1 s, D; A   restore
( j6 T% f+ a2 V% F8 \' y                 }
                 
7 e" e6 c+ E! O) l" g$ ?( S0 P        else {
                restore
                 di "cannot find pmaxlhddy240"
$ K0 U- r; X% B7 M; X                 
8 H1 F9 \% ]6 Q               
                           sum dx240
* g: Y+ x+ T3 V+ X/ O1 \                 local tempmaxplddy240=r(mean) in 1
               
% I; R, l3 [; ]+ D                 
                 }
2 @7 J% U3 P% {4 O0 }- z   
) [) N. P/ M* m8 ~' P   
   
9 B: i1 F0 p0 O   
   
2 V9 `, O3 ^) C# L   
' Z0 l, p5 D2 x% V        gen pmaxplddy240=`tempmaxphddy240' in 1
2 R& x5 t; x$ X# u* C7 k' W       
& A# a% C" ]% f        *-------------------------------------
   
   gen price5d=v5 in 1/20
' D' F6 z  j4 q  O      
; i2 P5 |: G9 M   gen price5yy=_n in 1/20
6 e. {- E& _9 @! q- P( t: |   
   gen price5y=price5yy*(dy240max-dy240min)/20
" u3 L9 D, G4 Z/ U+ H! {0 V8 Z   
   *price5y &price5d
' N& S8 N. w  L% l   
! ^/ g0 l+ a) J4 ]   gen priceny=price5y in 1
' P" I* [: t' L# V+ W3 ]  r   
/ ?% [& s1 J* ]" @! R   gen pricen=v5 in 1
* a4 x: B2 D. G; f
& z: r& }  k/ D7 c8 f( ~9 }     
9 b" A9 X6 P1 p7 s) o& H) Q, B   *priceny &pricen
   
   
   *-------------------------graph------------------
6 `+ F; q) x6 {' P, a) Z' a  
  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))
}
end

maodeng

awer 发表于 2013-2-5 14:40
% `/ j7 q+ Q, q抱歉，最近太忙了

7 k1 b( ^* K# R300236 我用下面的程序试过，没问题

哇，谢谢，我仔细看看，这些字母数字代码晕了
! c. E, I+ m5 q8 |

喜欢喝冰茶

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

喜欢喝冰茶

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

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

汉水云潇

准备试一试
* U% B  g; L; I% M% n* a/ z先找Stata软件

汉水云潇

晕呀
# B4 Q4 p: u$ p# s1 S! V论坛复制程序代码
+ _5 E! \1 N' S: T8 U+ s9 q乱码夹杂在其中
很烦
0 n, K4 d0 z8 d. O1 \/ e对了两次
以为没问题了
1 E7 u+ C# Q, \- l
7 [: f7 K5 v% j: p0 |结果
/ M6 j1 I) E0 ~$ q. O执行后
9 r0 V5 ?- n; _" X7 S- P还是出现
r 199
9 n) f: |( h1 J的错误
. [1 Y+ F* H4 @8 u+ g5 I* _

汉水云潇

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

$ ?% v! J3 n) ?2 q9 |' I0 |* ~

peterbev

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

中秋不圆

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

bknite

这个帖子值得学习。

jeanny

感谢！！！

Hubert2000

明月回春 发表于 2011-8-13 20:17
  k7 ~7 A& ]9 j% Z. f程序写的不错.
% C! B: y: I6 S) w, R2 l' U不过有一点要注意.通达信导出的数据,只对价格复权,并没有对成交量复权.所以数据需要手工处理 ...

, C# F8 G+ o* D3 t2 w老哥，也是山寨这篇文章找不到了，不知为何？