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

maodeng

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

maodeng

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

maodeng

. hello 300236
statement out of context
r(119);
& X8 }! x; ~7 j; H, t
: J# g/ _9 q( f+ {4 S大师 就是这个

sergio

这贴实在！

awer

maodeng 发表于 2013-1-28 13:13
. hello 300236
7 K4 ]7 _6 f% r7 \; fstatement out of context
2 y* O1 x. N, n6 t7 a# F* Zr(119);

4 |/ s' p  v0 s抱歉，最近太忙了

6 p7 \7 J: M9 ~+ \$ P: V( _300236 我用下面的程序试过，没问题
) A% E% c5 j8 L! @- D3 ~+ }8 ], A+ R
% X. V7 ^5 Z7 l4 G5 J; p; s你把 D:\stock\test\source 改成你导出的数据存放的目录，再试一试

9 N+ D/ U/ C9 w4 H& p7 x# Jcapture program drop hello

program hello

quietly {

clear
6 `- I; O% q, E: p* U: C# M& L
cd D:\stock\test\source

insheet using `1'.txt
- }; t6 q- |+ z: z
$ q- w  F, L) ^: W# K/ }   keep if v5>0
4 o1 P- E0 ?( W   
   gen vt=_n
   
2 ^) f* P6 r6 v1 k4 y   tsset vt

' A+ r2 _4 E  t" z, Q9 O2 T# a+ n   gen zf=D1.v5
9 Q% U8 L" x# Y% O* d# G
9 e7 N1 h4 }3 |7 X) |, L( P( W: U   gsort -v1

  L* E1 m% u" Z- P   gen avol=v7/v5
; G/ p2 K: w: t7 s0 l2 U   
   keep if avol>0
0 d" w- j. P, v% m% G) J. m   
# a; ^8 x; e/ q0 L2 O   drop vt
  
7 f$ ]9 ^+ }  b) N7 {. ?   gen vt=_n
   
   tsset vt
   
$ V/ [) |* o# X0 @" M1 Z    sum v5 [aw =avol]
2 F) U2 S+ g/ v* L* A
   gen amean=r(mean) in 1
   
   * 总加权平均值 amean
   
   
; o( o7 f1 W6 `   keep in 1/240

   sum v5 [aw =avol]

   gen mean240=r(mean) in 1
! e* T/ p9 B% a" \. T7 t1 h   
   *240加权平均值 mean240
   
   set obs 500
   
! v+ n/ P+ l1 D4 r6 P+ t   kdensity v5 [aw =avol],gaussian n(500) generate(dx240 dy240) nograph
' i  h2 D7 R2 G" A   
   *240 kdensity(dx240 dy240)
8 c9 q$ z7 k" q0 b2 x2 b: k2 u   
2 P' B! @& {) e/ f+ s0 F   sum dy240
   
( f* r& n3 v- x   gen yxis=r(max) in 1
/ H4 V/ k1 ~6 {- c  w" D$ O! o# B; }   
, r* ], H& x/ ]   gen yxis2=yxis/2
   
# W0 X' O; q9 i% F7 G8 _8 Q   *yxis
1 _' K7 J; n, b- R   
. g) H4 T! j* {$ B2 c   gen dy240max=r(max)
   
( q5 [! J" {* ]' o: ?   gen dy240minp=r(min)
& E/ e2 V! z: @9 Q- w& z
5 s( e* e- k! W& X$ H    gen vtt=_n
5 ]( c# S+ Y, J
3 s0 L; t% g' Z    tsset vtt

9 |7 q3 [" V8 G  E& @   gen ddy240=D1.dy240/D1.dx240
7 C9 c  H( d7 K  o' V   
+ D4 z# ~. t2 N7 j1 B   sum ddy240
   
/ K# `0 C* O* D( ~   gen ddy240max=r(max)
   
   gen ddy240min=r(min)
     
' j5 h3 H0 z  `, i   gen ddy240g=ddy240*(dy240max-dy240min)/(ddy240max-ddy240min)
0 i- F7 b' U% B5 I   
$ ?6 e' o$ m: W* Q1 V% T( c: U   sum ddy240g
   
   gen yxis3=r(min) in 1
* @. w1 R' N4 V- @; T) w   
   gen yxis4=0 in 1
7 K1 w+ G! G  Q3 ?  
  *ddy240g
   
   
  d$ R' ]: X0 c" i3 Z  preserve
1 b% ?3 V5 n- ?' U  
  sum v5 [aw =avol]
4 n" t- V/ m& ~/ q5 T" A: k. t  
  keep if dx240<r(mean)
8 G3 W- q0 [% {: x  
8 s1 k' N! n$ d* Z1 E( g- b  sum ddy240
8 a  O3 [. b+ n0 t4 Q  
  if r(min)<0 {
        
% q/ c# d1 K/ a* I+ i                keep if ddy240<=0
+ ?5 f% e* u. h( A* x               
- b# T/ z, \4 n  B6 f' M* _                count
. X1 z" B- s6 P' w) `               
0 S, F% V/ |! @% I+ g                if r(N)>2 {
                 
                 pkexamine dx240 dy240
               
                local tempminp=r(tomc)
                restore
                 }
                 
        else {
" P  R: Q& J+ C$ G6 N6 k- R               
                 di "cannot find minp"
                 
8 `0 _& W0 h6 k" F! k                restore
                 sum dx240
                 local tempminp=r(mean)
- `: C' l( b- [5 O2 U, n0 M                 
                 }
* j8 V. H* T0 Q: d                }
7 L4 R8 o+ {2 W7 t6 c2 Y3 V+ M   else {
                 
5 ^! r5 a/ _' B) `6 r' n& y4 Z  R                 keep if ddy240<=r(min)
                 
                 sum dx240
: ?! ]1 r* j; n! D9 c7 O                 
" Z4 T* f" x+ t. u/ V( ]                 local tempminp=r(min)
                 restore
+ Y; v* x- ]% e; j% v" c  o( Y1 B                 }

: h1 ]: h( m2 T, X6 z2 A
: Y% T) l* A( z8 T
+ g* k5 i" W9 M* Q) q+ C0 X   gen minp=`tempminp' in 1
   *240 minp


   
, v8 ^% ~6 _' [2 @4 P) ]# W   preserve
   
   sum v5 [aw =avol]
   
   keep if dx240>r(mean)
   
" l2 w6 T, L  q- X8 _! Z% o' X   sum ddy240
   
/ j" D5 p) M9 Y* V  y4 \   if r(max)>0 {
6 U; V' C4 c. P( U+ B, F7 x               
                keep if ddy240>=0
               
                count
               
! X; g2 d! C0 B& N6 x                if r(N)>2 {
                 
$ b/ m1 U  W( z- x1 L- X                 pkexamine dx240 dy240
               
  n9 W; V: O& J+ c7 u& O$ Q                local tempmaxp=r(tomc)
& b' Q2 d) u  k! D7 r$ z! o) ^( O2 f                restore
9 d" v6 J7 y; ~) K3 s* H: N7 S                 }
* F5 j1 o+ F- F. Q' u# J                 
                 
+ j# Y8 |( b- U. M        else {
                restore
2 ~9 J! g/ F' S6 V6 a                 di "cannot find maxp"
' F$ T& a. n% x/ I4 x3 M               
( x* e5 q" S0 h* w7 `; t* r* I- I               
" w3 `9 a; ?. d; h                 sum dx240
/ G9 B* ?8 }. z" h' S0 I0 D                 local tempmaxp=r(mean)
* M- Y4 d+ Y' Q& i/ e6 h                 
                 }
                 
                 }
4 X# Z: Q6 _; W2 F- `4 t  I5 K8 \  _                 
6 d  H: m" k+ x# S: a# N4 ^) R, T    else {
( i- q9 x7 e2 Y; `; q+ M+ p       
           keep if ddy240>=r(max)
- h- s' U7 g( Z$ n          
           sum dx240
- [4 y7 q- U3 W. K# u          
, M2 A3 [/ E7 _$ w' g4 T* k           local tempmaxp=r(max)
           restore
           }
/ V7 d! f( s/ z. b1 j+ f, r7 Q2 i1 R       
2 q) {8 O2 h# `% a1 u' N# B
  
0 o5 |% A1 K9 Y9 S2 D
   gen maxp=`tempmaxp' in 1
! k0 L: _! F2 r/ L5 B   
   *240maxp
   
, ^% F! ?" V, a& U, i" Z   
   *-------------ddy240pminp左侧高峰值-----
: u6 G5 _0 N  {: x   
    preserve
   
: V7 t# `4 ?8 C* z* u   keep if dx240<`tempminp'
) @* u4 l) T: z# N' L7 o9 y4 S   
. S6 a, b6 u% Z% v   keep if ddy240>0
- S. D0 ?$ E8 H4 z2 C" i0 @; M   
  {7 C) t' y1 A0 r   count
               
" B% j# U/ B9 j5 {- K                if r(N)>2 {
( t; w5 B/ h6 K: P                 
+ T/ |3 `+ K) s* N+ p$ r( c                pkexamine dx240 ddy240
               
                local tempminphddy240=r(tomc)
! P' Y  e* M" }                restore
                 }
                 
7 ^8 b# p! e4 V3 W* k        else {
9 e0 m( S' I0 ]# _! o5 u8 s                restore
                 di "cannot find pminphddy240"
. g' Z( h& g+ H; e8 R* ]+ q                 
# Q$ x0 G2 i, q               
                 sum dx240
                 local tempminphddy240=r(mean)
8 i- @4 a4 @0 b3 ?0 b1 X               
* c9 @$ g5 u' k; \                 
5 ^( q' A8 i# K- `/ _: I! {3 _                 }
                 
   
   
  C+ x. j6 p/ ?$ N   
   
2 u  o% u' p" H1 K, P1 e* v   
9 l; a' @4 }# [        gen pminphddy240=`tempminphddy240' in 1
       
       
        *-------------ddy240pminp右侧低峰值 -----------------
- e$ G0 D: [5 \2 o    preserve
, A' E% y) ]* e, ^  P# U. A6 {! D
   keep if dx240>`tempminp'
7 j% x3 p! S9 W# L* ~# L$ g  w4 R   
2 |+ C, N% }; _   keep if dx240<`tempmaxp'
7 x. a5 o! a4 `7 Q' P+ K) C3 R   
5 S9 W# |4 ]! P' L: w7 z( w+ J   keep if ddy240<0
   
4 p# N, P+ k$ T   gen temddy240=-ddy240
2 [) Y( u* r& Q   
+ [6 [5 v4 b" m5 S8 t    count
9 L/ v$ a0 R9 n; y8 Q( R2 a               
                if r(N)>2 {
- H+ R! _' u9 d" m& [& C! ~4 y# _                 
                 pkexamine dx240 temddy240
9 @( o0 [6 w: o) w   
$ O( ^$ Q/ c' ?7 m9 e4 _% h3 D  local tempminplddy240=r(tomc)
+ K0 `5 Z" ^2 M" t$ D( x3 c- N+ t                restore
# h7 H  q* ^( W# @7 O                 }
/ W7 x' w% V- y2 Z                 
3 h& `! G! y5 Q+ l7 O, r/ {        else {
, f! x) W* r% v5 Z               
/ _2 {4 x  Y6 \. A4 l  Z8 I                 di "cannot find pminplddy240"
  H# t5 |# ^0 R( L                restore
                 
                  sum dx240
                 local tempminplddy240=r(mean)
                 
                 
' w+ N# s9 ?- {$ p                 }
" `5 _, h  ]4 c/ i   
   
   
: `' S# S/ ]' l: @7 X! L   
, y# Q+ l1 _' S; H' x) o# P- z7 t6 q   
- Z4 l/ E3 D. E+ o; G% R        gen pminplddy240=`tempminplddy240' in 1
& y6 \) b! U9 o8 z. u$ }       
        *-------------ddy240pmaxp左侧高峰值 -------------------
6 D3 a- ~, c+ N6 w; k8 P8 W6 C       
         preserve
   
7 _  A& R8 K! R$ Y" a" L' K& @   keep if dx240>`tempminp'
0 S+ T9 h8 O, Q! L4 }   
6 X, D$ J. Q6 k2 c9 i9 @, |7 U   keep if dx240<`tempmaxp'
   
   keep if ddy240>0
   count
0 T" x+ G, B; f& L0 h" B* v               
                if r(N)>2 {
! a+ I6 h% O1 K9 m. x. p8 c                 
                 pkexamine dx240 ddy240
   
5 H& H3 n% M' }! c   local tempmaxphddy240=r(tomc)
, b- n3 y, p2 S7 ]   
; W% p' D" P# [/ z   restore
                 }
                 
        else {
, Q2 x6 S# r/ @               
) J! O; n' g& P" @2 j# q                 di "cannot find pmaxphddy240"
) @% P8 \0 Q& ]* a( U; _& b( m                 restore
               
                 sum dx240
                 local tempmaxphddy240=r(mean)
6 M, b5 {% n" L: O0 ?                 
                 
" |8 N, Q( N2 [' ?9 t$ N                 }
   
* N$ I  y' p5 Z   
" E3 R4 E& b" Q0 }* E, ]   
   
   
   
        gen pmaxphddy240=`tempmaxphddy240' in 1
. y; i# `- w1 Q; \       
       
        *-------------ddy240pmaxp右侧低峰值 -------------------
       
         preserve
  
   keep if dx240>`tempmaxp'
7 \' j( v( W) W3 t  o0 ?   
2 |/ n+ d5 T1 r0 w% g" P   gen temddy240=-ddy240
   
   keep if temddy240>0
! S: ~+ O9 |1 K0 u0 E- c   
   count
               
! t3 u3 }( s( [' q                if r(N)>2 {
+ ~8 f. q3 t. z6 v$ K                 
0 x7 ^- Y+ g4 P% J& |% a                pkexamine dx240 temddy240
0 t$ U/ w- K/ a: d9 [   
: k! z% H& s# z) M( \6 v9 |8 R   local tempmaxplddy240=r(tomc)
   restore
                 }
                 
, |( i; L, }; K8 I* N- q* e& j        else {
: D2 Z6 Y  r( f5 g5 f0 G2 M                restore
                 di "cannot find pmaxlhddy240"
                 
; D& B  n9 K. b1 j2 B0 G               
                           sum dx240
2 o) D% ]9 ~- M% w& Y8 G                 local tempmaxplddy240=r(mean) in 1
               
9 k) Q6 Q0 p. A- @0 {                 
                 }
0 J- _6 u1 D& ~  T  @  Z8 ~0 v   
5 X4 I1 _. [$ }0 P9 D   
   
# w0 A" P) G# K/ b   
   
3 t1 [" L( z( P! C- n+ _# m   
0 E7 v7 E1 a: p" }        gen pmaxplddy240=`tempmaxphddy240' in 1
       
2 K6 `: {! ]$ l% `        *-------------------------------------
: p# L: `* _# r# M1 u* x7 X   
, T; v1 V1 G" I5 ]   gen price5d=v5 in 1/20
$ r$ }8 e- g3 A+ e( ~      
   gen price5yy=_n in 1/20
   
   gen price5y=price5yy*(dy240max-dy240min)/20
$ O) \% P& N* p1 L  k  S0 F, d% z   
   *price5y &price5d
( P5 O; C2 u2 H. ]  L; J* C* m8 {   
9 d# ~0 N- r' I3 m5 t5 j- y1 b   gen priceny=price5y in 1
   
   gen pricen=v5 in 1
- T4 _0 e7 ?2 T' X/ a& S5 g$ ^7 G+ Y
+ f1 w+ V( |1 L+ f& _     
- o4 B1 G; i! w) H# m& W   *priceny &pricen
9 J& ?( n! e2 c0 S   
$ ]) d9 J8 i* W   
' F% R- h! M7 W4 Y0 @   *-------------------------graph------------------
2 M. |$ [, V4 ]1 j# e" }+ J: ^  
  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)
/ G: c% n9 S! c' j: J/ ]: L
' G+ m: ?7 y% W; B+ m* ---------可选项目----------------------
6 h3 b& l( C4 x0 X1 b! W*(spike yxis amean,lc(blue))(scatter yxis2 amean, lc(blue) msymbol(+) mlabel(amean))
* W% g0 n. _' d- \+ ]' g}
end

maodeng

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

" E" m' I# q& ^8 F; Z) r300236 我用下面的程序试过，没问题

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

喜欢喝冰茶

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

喜欢喝冰茶

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

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

汉水云潇

准备试一试
先找Stata软件

汉水云潇

. G0 T; C9 V  ^- L2 D% |晕呀
论坛复制程序代码
% o+ F. `. T) x8 q, v% e% T+ E2 a乱码夹杂在其中
$ E# z  z/ l: {% ^6 d& a很烦
+ h8 Q9 s5 D0 o$ \4 g! ~( [对了两次
以为没问题了

8 p/ u7 I; X! M* ^. ?结果
5 I5 H& `" W  \% x5 Y& X执行后
还是出现
r 199
的错误
2 B1 @- l) f- N3 G: ]3 e0 x

汉水云潇

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

peterbev

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

中秋不圆

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

bknite

这个帖子值得学习。

jeanny

感谢！！！

Hubert2000

明月回春 发表于 2011-8-13 20:17
5 V* C/ E% D! \3 K  B" M程序写的不错.
# ^, Q0 A, T( y- ?& l不过有一点要注意.通达信导出的数据,只对价格复权,并没有对成交量复权.所以数据需要手工处理 ...

, Y- S* ]1 a+ S! m' q4 D9 W/ f老哥，也是山寨这篇文章找不到了，不知为何？