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

maodeng

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

maodeng

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

maodeng

. hello 300236
statement out of context
r(119);

! H5 j9 E' @9 |5 F: y9 D大师 就是这个

sergio

这贴实在！

awer

maodeng 发表于 2013-1-28 13:13
. hello 300236
statement out of context
# B4 e$ C  m; i' }. |r(119);

抱歉，最近太忙了

300236 我用下面的程序试过，没问题

你把 D:\stock\test\source 改成你导出的数据存放的目录，再试一试

0 }( ?  t% H9 V2 e+ Scapture program drop hello
8 i, i0 x5 e/ I
program hello

quietly {
( B! n" d, p' d0 D2 \9 V4 Y
! Q' P9 `  r' m1 s/ P: {+ {: Eclear
! ~; [3 S- C* q7 J! s5 a
9 N6 s+ \" _5 s5 R* T6 H5 Jcd D:\stock\test\source

0 w, q9 W' R3 G. x3 p+ Winsheet using `1'.txt
  Y. S$ d) }1 q! m, P( R
   keep if v5>0
, t3 `8 w! r+ g; D& I   
   gen vt=_n
   
   tsset vt
1 c" B/ [* m* f/ O
   gen zf=D1.v5

   gsort -v1
, \, ]7 X8 s/ R/ G) J7 h
   gen avol=v7/v5
( X; w% J3 X3 h   
   keep if avol>0
   
   drop vt
: t* C) A" T* G. n  
   gen vt=_n
   
   tsset vt
/ u' l' @9 R4 @: g   
    sum v5 [aw =avol]

   gen amean=r(mean) in 1
   
6 r7 I) g2 `* ~4 a+ X7 t7 v   * 总加权平均值 amean
" H3 W% [7 D6 S1 o# I7 c   
   
* p: K! V7 ?. T   keep in 1/240

( \3 ^- Y& @$ B- g   sum v5 [aw =avol]

   gen mean240=r(mean) in 1
. h# Q" F9 R: F. ]( f& ]" z9 U" N   
  i  O3 ?: ]- Q8 F* T: ?' k: X   *240加权平均值 mean240
   
9 @: L  Z4 P6 h- q* }   set obs 500
   
* u8 V5 I4 z( A! R; C8 b   kdensity v5 [aw =avol],gaussian n(500) generate(dx240 dy240) nograph
& N* R( `5 r; z) `6 j+ @  h. F   
   *240 kdensity(dx240 dy240)
   
   sum dy240
   
   gen yxis=r(max) in 1
   
   gen yxis2=yxis/2
  X. |% A/ {: H% G$ a8 }& k/ k. C- X   
   *yxis
/ Z* U, W5 n, |/ v2 \/ d1 {   
   gen dy240max=r(max)
6 C: n" c: J% q+ X+ q   
) @2 H( X$ S! g8 W   gen dy240minp=r(min)
/ u3 r5 A2 ^  W4 ]1 H% X
+ ^# D: u! y5 Q" Y0 q, J    gen vtt=_n

    tsset vtt

   gen ddy240=D1.dy240/D1.dx240
  K$ Z, z0 C6 F: z   
: U4 h; N5 Z' j  f6 B+ T   sum ddy240
   
/ }; u$ g: X6 J   gen ddy240max=r(max)
& A5 F% F: {9 @) m   
& K0 A) c( r  F6 A! s: l   gen ddy240min=r(min)
3 O; Q7 r" f) L# t: [/ B7 b7 ?     
4 K+ c! @# E: u1 s) ?: R   gen ddy240g=ddy240*(dy240max-dy240min)/(ddy240max-ddy240min)
( }5 G6 c3 Y/ @# V" G8 g' d   
/ ~2 q1 }: Q9 |! g   sum ddy240g
   
* a6 k2 p8 g1 W2 Y   gen yxis3=r(min) in 1
   
   gen yxis4=0 in 1
  
. A* `5 l+ L7 V9 J  *ddy240g
   
/ h- G1 ~$ Y1 P6 u( h" g3 J0 Z' _8 I2 l   
' F. J6 t7 h; M5 I, a  preserve
- u" V) \3 Y- y) R- }; \  
  sum v5 [aw =avol]
" h/ r7 B5 m& E8 h! V/ [( A  
7 [  P* I6 P: j* i. y7 [8 R  keep if dx240<r(mean)
4 B2 s; J5 I5 u- }  
4 {# ~8 |- }- O) @( H# J8 n$ ?* t  sum ddy240
  
  if r(min)<0 {
        
! @+ y  N4 f1 C2 m                keep if ddy240<=0
               
                count
               
                if r(N)>2 {
                 
' T) V5 \, f2 h- C                 pkexamine dx240 dy240
               
                local tempminp=r(tomc)
                restore
" m' Y& T) k. j2 j) _                 }
                 
: b, i% E0 N+ Z! n! ^        else {
               
                 di "cannot find minp"
                 
                restore
                 sum dx240
% e& X6 ]8 x4 j0 }) }  y                 local tempminp=r(mean)
                 
- m! }! _8 S! `2 v- C7 a' i                 }
                }
, x  a0 {- B/ `* J, X6 d0 `   else {
( E" o$ w9 K' r- y8 L( ^; O( c                 
                 keep if ddy240<=r(min)
                 
                 sum dx240
  j4 l) j" Z% g+ y; `# i                 
                 local tempminp=r(min)
                 restore
" r1 K4 K* F$ l- l) @8 x                 }
% N" A( y; J) Z

+ R" H% g: }. o" p
6 _! G: p* R8 ^) J3 j% K   gen minp=`tempminp' in 1
' b+ Y" }% O8 a" M8 n3 I- @% q( s   *240 minp

0 \2 y5 x# @, G1 A. j
   
   preserve
% C. M4 J8 ?  T, x9 H7 F: ^9 g4 Y' D  j   
/ A# H# T; s( }+ x% Q- J. ?& P   sum v5 [aw =avol]
   
& C4 F* [+ n9 P1 ~: ?6 F% G   keep if dx240>r(mean)
   
) g! E- X9 a& g9 I# g9 Z9 U   sum ddy240
" [9 T8 O/ {0 K( q% Y0 J' V   
1 }2 t0 @: K- V" x6 l2 ^5 ^   if r(max)>0 {
. z+ K0 W6 O$ |               
& A- R. c( M" B; N' O$ Q3 ~                keep if ddy240>=0
) E6 L6 S: q* e/ P7 m% F               
                count
6 E# X( n* N6 L- Y4 C' e               
                if r(N)>2 {
7 ^! \5 l% c2 P4 X                 
0 Z9 z1 O2 S0 C$ |0 ^) a/ Y                 pkexamine dx240 dy240
               
                local tempmaxp=r(tomc)
0 }* \/ |2 g, ~' v  G5 P1 w                restore
6 u" F8 O" D+ ?: v3 g2 l6 L                 }
                 
                 
1 y% M' Z& D& S        else {
                restore
" x# g' U5 d% p/ O+ a# l                 di "cannot find maxp"
               
  ]$ I, z$ L$ i               
( r1 C! Y! A# F6 h; e' l                 sum dx240
* a9 f0 I3 ~5 ]1 }7 X" }7 K                 local tempmaxp=r(mean)
1 H: k( \: y( H" \" G" ?) D                 
                 }
; s+ K" ~3 c' z" X8 k: S                 
7 n. T- s7 z" V( {! x                 }
( ~3 l* t+ v3 G                 
    else {
, |) V6 K; s( e' X5 ]; B* J/ w       
           keep if ddy240>=r(max)
          
           sum dx240
          
. h3 Z, R% R; d) R% K           local tempmaxp=r(max)
           restore
' Z) u8 R" W3 Q. K2 r9 s           }
: _2 _. f. O2 a& u  J       

; q5 q% @# s' l- z- I2 J; E  

   gen maxp=`tempmaxp' in 1
6 g, H: {* @  y% b  q" i5 a   
# T, R1 d" r& _# ]! B& @1 m; ]8 Q   *240maxp
   
! _5 x# I% E$ C4 Q5 _# A   
   *-------------ddy240pminp左侧高峰值-----
   
    preserve
! y9 F. C. B: _5 _- ~7 V   
   keep if dx240<`tempminp'
   
   keep if ddy240>0
6 ]* |6 Q/ b  ]5 W; H; r# x  w   
   count
; ]/ q! f* z: M' v8 j; i               
7 z5 y& R0 g; S4 a& [% |0 U' |                if r(N)>2 {
4 V9 P3 ]$ i# l% o: |                 
" P* D0 |0 {5 ]- Z+ x                pkexamine dx240 ddy240
* R& `: ^3 a' X/ v4 l" J- R0 t               
" ^6 L$ N0 T, d% a* _5 U* V5 z& w: Y                local tempminphddy240=r(tomc)
2 W' }+ b8 l8 |1 R6 r                restore
                 }
- L7 @* b  [- ~                 
. d% E$ A8 U  ]4 O" [9 C- u        else {
1 }; G8 W# @, A, U$ d8 j( }                restore
                 di "cannot find pminphddy240"
, ?  G5 F8 b! z                 
               
                 sum dx240
; Q$ [1 R1 ^; Q+ T- w6 j5 q7 Y* O: E                 local tempminphddy240=r(mean)
4 d; D* h5 O: A               
; P. |: `! Z5 s  d                 
                 }
                 
! A, d, G/ c$ d$ S   
+ K6 V3 o$ x8 Q1 Q2 }8 X' U. Z, n# ^   
   
3 k$ j* h) x; [  g; c2 q, q! Q   
   
        gen pminphddy240=`tempminphddy240' in 1
       
( i( ^) j4 P6 I; H' ~  W$ _* g       
2 W! h6 O: f0 v. B        *-------------ddy240pminp右侧低峰值 -----------------
( ?7 t5 t' b1 Z' S/ ^    preserve

   keep if dx240>`tempminp'
+ x4 o& A3 i9 N; D. `   
   keep if dx240<`tempmaxp'
9 O( a; t  m. F# j. L* B/ S) ~   
   keep if ddy240<0
4 m6 v' p  x* W   
; _/ c7 M8 m( P. K0 L   gen temddy240=-ddy240
+ P2 Z& n6 @, Z- N   
, p; I9 C! p/ s    count
  g. d+ p- U) b) a. f. p  B               
                if r(N)>2 {
2 j4 z% P. Z" U9 ~. K                 
+ A2 f& r: U8 J9 i' m: ^/ S# o                 pkexamine dx240 temddy240
   
  local tempminplddy240=r(tomc)
( M1 W* k1 @1 k: H# P1 ~  V                restore
                 }
                 
        else {
               
                 di "cannot find pminplddy240"
                restore
: p$ Z8 \# y7 {3 V! m: v7 o8 |6 Q                 
! F" l% U$ F/ C, y9 J                  sum dx240
                 local tempminplddy240=r(mean)
- U8 s4 Q9 `, T- ~                 
                 
$ Z7 U6 r) `) K" s                 }
3 P& b  Y, ^6 D/ a8 l# p. G   
8 R& B1 t: z+ P   
7 h/ B; A9 W  a1 b   
   
4 v2 L* \; @- V: t   
" ]( B5 V- y( w, W( \        gen pminplddy240=`tempminplddy240' in 1
( p6 [7 J# W2 M7 f) a       
$ e1 e0 v  ]+ g0 t/ M- [$ l        *-------------ddy240pmaxp左侧高峰值 -------------------
/ ^* Z2 f; q. T+ x$ `8 _( c       
         preserve
1 ^, @" F7 Z- s/ R  X: Q   
# ~  y' w# `' u1 a- k! c. e+ `   keep if dx240>`tempminp'
7 B" [$ u9 Y( y2 U7 j9 p   
" m* H; z7 _; ]8 G   keep if dx240<`tempmaxp'
   
1 f- p1 A3 u9 k. C# v$ C) {   keep if ddy240>0
" Y, N1 j0 ^4 f' P- }8 c   count
: T: I. {+ e* }+ q6 f               
6 `) j# z, j4 ^. W6 j( u0 i! h& j                if r(N)>2 {
7 e& z2 }- P7 m7 f' I                 
                 pkexamine dx240 ddy240
8 w, t! I, X: J) e% J   
   local tempmaxphddy240=r(tomc)
) W6 j* I2 ]$ P/ }8 a   
4 n% X" B: g* H& W   restore
                 }
                 
8 X; ~4 E! H6 M8 F# H        else {
               
4 E" h& `7 J' I- B/ q9 g                 di "cannot find pmaxphddy240"
: f. Z3 i1 z+ b: V                 restore
               
                 sum dx240
2 w% [9 U# A7 h' Q                 local tempmaxphddy240=r(mean)
$ ?) v1 Q& v, A  D                 
- T5 `* n- A2 u$ I# D7 F3 _                 
                 }
! a% m/ a( _3 \+ W& X, v3 O; S: A* E   
   
* c  E# A8 L+ }1 A7 `& e   
; P0 \: G- k2 w# q) ~  R/ A9 w   
, H% o$ s) _+ T; D3 c   
4 e; @, }& |( \2 v! t4 K   
' d6 ?6 g( Z* ?        gen pmaxphddy240=`tempmaxphddy240' in 1
1 e, i. ^! v) u1 P. Z1 z       
' y' Z" N$ d7 R# q       
        *-------------ddy240pmaxp右侧低峰值 -------------------
       
         preserve
/ f+ M/ M; A( k* O9 c( p* t. c- H  
   keep if dx240>`tempmaxp'
7 U. U7 ~. z/ B- l2 O* p   
- w# K9 a$ Q8 i2 J% Z7 Y   gen temddy240=-ddy240
0 L" Q$ c- p& @& S   
   keep if temddy240>0
   
   count
               
* X* i8 S# `- e3 j6 |: p! ^                if r(N)>2 {
) Q# G1 i* l4 X                 
                pkexamine dx240 temddy240
   
   local tempmaxplddy240=r(tomc)
6 ~) C: j8 }$ i   restore
2 U; {0 H& s/ W& g8 Q" R                 }
                 
' A, G( p/ x+ f        else {
                restore
                 di "cannot find pmaxlhddy240"
$ g# l) j$ b6 I# a; _0 t                 
; x8 m% d# j2 i) d7 N$ ?; c               
; z# N' c; x. K6 [                           sum dx240
                 local tempmaxplddy240=r(mean) in 1
               
                 
                 }
; C- P" _  |, n   
% m& K8 T" }4 G' a0 I4 `: i   
) ]( B5 `$ I+ Y   
   
5 z3 u, ]: C8 ~# T6 k' u0 @5 d   
   
" w1 ~1 D# _7 k. g        gen pmaxplddy240=`tempmaxphddy240' in 1
5 y5 Q$ T+ d1 R; L. g9 O) X5 e       
        *-------------------------------------
   
4 y' H3 ~5 u9 L9 R; m) u9 B   gen price5d=v5 in 1/20
      
   gen price5yy=_n in 1/20
   
* f# O) e. @0 ^& N   gen price5y=price5yy*(dy240max-dy240min)/20
   
" w( d* \/ N8 J8 C7 ?' k" c   *price5y &price5d
- T9 }3 d$ b" U7 B* T, U/ C7 `   
   gen priceny=price5y in 1
   
0 l) }2 u6 b1 q: H   gen pricen=v5 in 1
! G  m1 y& h7 k
) \4 ]  e9 y9 a7 K     
   *priceny &pricen
   
   
   *-------------------------graph------------------
  
  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)
: z9 [' A# s% e* j: g  f5 b
) b) Y- w2 K% ~8 p& n2 s8 p+ [* ---------可选项目----------------------
*(spike yxis amean,lc(blue))(scatter yxis2 amean, lc(blue) msymbol(+) mlabel(amean))
! F5 a( o. s5 n" g6 C; x$ A}
% O& e3 R" \) K6 |2 e- t9 {$ d0 kend

maodeng

awer 发表于 2013-2-5 14:40
0 O' ]5 w. l& L抱歉，最近太忙了
! [& Z2 X; F# R- _2 h& D8 @
+ j( s1 R) Y6 {' N5 C4 ~300236 我用下面的程序试过，没问题

$ v  z9 }8 r5 b! ?- U5 g0 f哇，谢谢，我仔细看看，这些字母数字代码晕了

喜欢喝冰茶

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

喜欢喝冰茶

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

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

汉水云潇

准备试一试
先找Stata软件

汉水云潇

晕呀
论坛复制程序代码
; G' S8 W9 y% r$ ?乱码夹杂在其中
很烦
对了两次
1 G; ]5 D8 Q' S$ H* A# J" _以为没问题了

结果
执行后
还是出现
r 199
的错误

汉水云潇

试验成功，下面是三一重工的图
2 c  [1 J/ \! v' S- \3 t4 g) O

peterbev

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

中秋不圆

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

bknite

这个帖子值得学习。

jeanny

感谢！！！

Hubert2000

明月回春 发表于 2011-8-13 20:17
- m& O# ^. F" E% {5 ]) @. I9 P' g程序写的不错.
不过有一点要注意.通达信导出的数据,只对价格复权,并没有对成交量复权.所以数据需要手工处理 ...

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