爱吱声

标题: 一个用stata 计算股票成本的小程序 [打印本页]

作者: awer 时间: 2011-8-13 18:03
标题: 一个用stata 计算股票成本的小程序
本帖最后由 awer 于 2011-8-13 20:49 编辑

我是从西西河一直追虎大追到爱坛来的，虎大曾用stata 分析过股票的长期成本和一年的短期成本，但是直接用kdensity命令无法得出筹码峰的具体数值，因此我写了一个小程序来计算股票的筹码峰值、长期平均成本和短期平均成本。下面就给出具体实现方法。

计算之前需要有：

1）招商证券系统（通达信其它客户端也行，功能都差不多）

2） Stata软件（我用的版本是Stata/SE 11.0）

然后

1．新建一个文件夹储存股票数据，比如:

e:\stock\source

2. 股票数据导出

a 进入招商证券, 系统-盘后数据下载，更新数据

b快捷键“34”进行数据导出, 在弹出窗口选择
“高级导出”,

在高级导出窗口选择 “日线导出”,

将导出目录选择为你新建的储存股票数据文件夹的目录（e:\stock\source）! 切记切记

文件名为 XXXXXX.TXT

勾选数据精确复权

添加品种中选择想要导出的品种

开始导出

3下载 stata, 将本贴附件中的hello. do文件copy至stata根目录下

用Stata中的 Do file editor 打开hello.do 文件，确保hello.do 文件中 e:\stock\source 和你保存股票数据的文件夹目录一致

现在就可以计算了，以露天煤业（002128）为例

在Stata中输入如下命令

do hello

hello 002128

(002128为你想查询的任意一只股票的代码)

几秒钟之后应该会出现下面一张图

hello.do (5.25 KB, 下载次数: 13796) 其中蓝线为股票上市交易以来的平均成本值（注意!此值没有去掉一级市场发行和国家的0成本）

红线为一年内交易的平均成本值

黑实线和虚线为一年平均成本值附近两个筹码峰的值

黑方框代表的事20内每天股票的收盘价格

红方框代表股票最新收盘价

山寨完毕，收工!!

作者: 明月回春 时间: 2011-8-13 20:17
本帖最后由明月回春于 2011-8-13 20:19 编辑

程序写的不错.
不过有一点要注意.通达信导出的数据,只对价格复权,并没有对成交量复权.所以数据需要手工处理.虎大买的是Thomas One Bank.
还有就是虎大除了用了kde,还有摘要功能，你可以看看这个也是山寨

作者: awer 时间: 2011-8-13 20:20
这点我也注意到了，所以我用成交额/复权后的价格作为成交量用于权重计算

作者: wenty 时间: 2011-10-19 14:32
多谢，挺好用的。
通达信似乎无成交额一栏。导出的txt文件要删掉header,否则stata当数据为string了。

作者: 旅途愉快 时间: 2011-10-19 15:40
本帖最后由旅途愉快于 2011-10-19 17:44 编辑

没太看懂。。。

露天煤业最近20天的收盘价在20块左右？

作者: 旅途愉快 时间: 2011-10-19 17:47
仔细看了看，发现原来自己理解错了~：）

谢谢awer。有问题再请教哈~~

作者: 旅途愉快 时间: 2011-10-19 18:09
还貌似有一个问题啊，为什么我处理2000年以前上市的股票，都会出现：

~1 k% J% f+ Q# x4 x
hello 000425
# V! y2 Y( P8 R+ Q, w! N
time variable takes on negative values0 a6 h: @ n" g! v1 o% ~. o
r(459);
/ [1 _$ o& C9 m4 j) c" D3 D
7 m: w5 o6 T L/ s+ \3 O
. hello 0007835 ~5 t7 N# p% f5 v3 R8 c6 J
time variable takes on negative values
: m* Y K! T; e+ N7 ]6 e. F
r(459);2 e2 q: ^" A9 n

复制代码

2000年之后的没问题，比如露天煤业。

作者: 冷眼看世界 时间: 2011-10-19 19:25
不错，我也来试用一下

作者: 沉默的前哨 时间: 2011-10-19 20:45
为什么老是提示：type mismatch

作者: 沉默的前哨 时间: 2011-10-19 20:49
能不能帮忙看看：

do hello

.
. capture program drop hello

.
. set mem 1024m
no; data in memory would be lost
r(4);

end of do-file

r(4);

. hello 600063
unrecognized command: hello
r(199);

作者: awer 时间: 2011-10-19 22:03

沉默的前哨发表于 2011-10-19 20:49 5 j" ~8 k/ _$ P$ s, `
能不能帮忙看看：
9 C9 P1 w3 {. K9 x
4 p0 s+ k( Q" F7 s2 k do hello

你算得哪个股票？

作者: 阿雷 时间: 2011-10-20 08:52

awer 发表于 2011-10-19 22:03 1 C9 u0 R: [8 [5 a0 e$ r
你算得哪个股票？

那现在露天煤业可是在平均成本以下啊

作者: 沉默的前哨 时间: 2011-10-20 10:30

awer 发表于 2011-10-19 22:03 . R6 Q" h% J! ]/ T' ]1 T+ \9 Y
你算得哪个股票？

我是用通达信引出的数据，比如600663，在e:\stock\source 里的文件是以SH600663.txt存在，无法执行，后改成600663.txt还是不行，请指正！

作者: 沉默的前哨 时间: 2011-10-20 10:32

awer 发表于 2011-10-19 22:03 / d) l0 V) t. F! r- g. h1 ]
你算得哪个股票？

do hello 后变成
......

end

.
end of do-file

作者: awer 时间: 2011-10-20 11:18

沉默的前哨发表于 2011-10-20 10:32
( O/ Q- `2 ~3 |do hello 后变成; T) |7 [: Y: X$ A/ S
......

作者: awer 时间: 2011-10-20 13:32

沉默的前哨发表于 2011-10-19 20:45
' k# I6 x. E4 q9 U% }为什么老是提示：type mismatch

注意日期格式要选成 yyyy/mm/dd，要不然时间变量格式不对

作者: wenty 时间: 2011-10-20 14:16
"此值没有去掉一级市场发行和国家的0成本" 是指已包括一级市场发行和国家的0成本吗?
谢谢

作者: 琴心三叠 时间: 2011-10-20 17:35
下方的棕色曲线是代表什么的

作者: maqiangcn 时间: 2011-10-20 19:13

明月回春发表于 2011-8-13 20:17 + `, h% r+ K9 N5 L8 g9 Z
程序写的不错.
4 J) V# }1 P/ V+ @ y& n8 c$ v% S不过有一点要注意.通达信导出的数据,只对价格复权,并没有对成交量复权.所以数据需要手工处理 ...

你好，这个链接好像不对，能确认一下吗？

作者: 琴心三叠 时间: 2011-10-20 20:03

琴心三叠发表于 2011-10-20 17:35
/ r1 ]" p1 |$ W* G- w' E0 b0 i4 K下方的棕色曲线是代表什么的

两条曲线只是周期的不同吗

作者: 琴心三叠 时间: 2011-10-20 21:37
谢谢答复还有一个问题我导出的文件为何只有当天的数据而不是历史全部数据呢

作者: 偶然 时间: 2011-10-20 22:19
建议楼主给个数据格式模板

作者: awer 时间: 2011-10-21 10:05

琴心三叠发表于 2011-10-20 21:37 $ K+ q# y4 K [ ]" X; b: W q" b
谢谢答复还有一个问题我导出的文件为何只有当天的数据而不是历史全部数据呢 ...

有很多问题是数据不全造成的
如果想要全部数据，
系统--盘后数据下载--日线数据--日线和实时行情数据

把开始时间设成中国股市开盘的那天，下载全部数据

下载完全部数据后，在屏幕上使用快捷键 34

34--高级导出--

然后选中日线导出；选定导出目录；文件名设为：XXXXXX.TXT;勾选精确复权；不要生成导出头部；分隔格式为tab（）；日期格式为YYYMMDD，然后选择添加品种

我一般都是把所有数据下载到指定目录下，然后进行hello运算

作者: 沉默的前哨 时间: 2011-10-21 12:40

awer 发表于 2011-10-20 11:18

谢了！我自己还是没搞定，我再看看。

作者: 沉默的前哨 时间: 2011-10-22 12:31
搞定了，谢谢！

作者: warbrai 时间: 2011-10-22 16:38

沉默的前哨发表于 2011-10-21 12:40
% c7 O* ` b3 [$ \谢了！我自己还是没搞定，我再看看。

别选择 "导出数据头部"

作者: 沉默的前哨 时间: 2011-10-22 23:16

warbrai 发表于 2011-10-22 16:38
8 b0 W- @# }% P; N0 N' T8 p6 G& |别选择 "导出数据头部"

是的，谢谢了！现在可以了！

作者: 心静天地宽 时间: 2011-10-23 20:45
老早就想找个这样的软件和方法了，无奈见识有限，一直没有办法，感谢你的无私。

作者: 琴心三叠 时间: 2011-10-23 22:04
抱歉我还是没明白棕色曲线的意思黑色线代表筹码分布密度那棕色线你说代表随价格的变化率这我不太明白请说详细些不好意思我比较愚钝

作者: 琴心三叠 时间: 2011-10-24 18:45
hello 600085
unable to change to E:\fenxi
这是什么意思？

作者: 观风望月 时间: 2011-10-24 20:32
已经下载下来用了，非常感谢热心的分享。
因为不会使用stata，所以有个小问题想请教一下，我按照你教导的方法已经成功测试了露天煤业，兴业银行等几个股票，都成功出图，但是在计算600188.601999的时候，程序提示“
. hello 601999
statement out of context
r(119);

我查了一下Help的内容，对于这个错误也没有特别具体的解释，不知道是什么原因造成的，因为数据的导出和格式和前面几个股票都是一样的，所以我搞不太清楚到底是怎么回事。

另外，虎大提到过的国家成本和一级市场成本加权可以通过在数据中添加一行交易数据来实现，只不过这笔交易的成交价格和成交量需要自己加权计算，对于不分红的股票来说加权计算还是比较容易实现的

作者: 本拉风 时间: 2011-10-25 15:22
横纵坐标分别表示什么，2条曲线又表示什么

作者: 钱潮滚滚 时间: 2011-10-25 16:08
发个图验证一下和大伙的结果有没有出入.

收盘价到昨天为止的三一重工:

作者: kevenjames 时间: 2011-10-26 14:05
就楼上这个图，几条线，还有黑点分别都是什么含义啊？

作者: wenty 时间: 2011-10-27 04:17
总结一下使用心得
1）删掉set memory那一行
2）改路径
3）从软件导出数据时，要先下载，后导出，不选导出头部。
4）用notepad看导出的数据每列对齐没，用tab修正。
希望有帮助。

作者: 沉默的前哨 时间: 2011-10-27 09:09

老祖爱蓝军发表于 2011-10-26 22:50 a. g0 `0 I+ a8 h9 V% j
我出现同样问题，请问怎么搞定的？

关键点：
1、文件名设为：XXXXXX.TXT;
2、勾选精确复权；
3、不要生成导出头部；
4、日期格式为YYYMMDD；
5、如果不行，建议可以关闭stata再打开试试，我认为和操作系统的版本可能也有关；

作者: 范进中举 时间: 2011-11-8 16:10
不错，按照楼主操作完，最后是选择do命定，选hello文件，然后hello 代码就行。
几个问题

最边上的蓝色的4.473554是什么意思，红点13.84054是最近20天平均价？灰点15.3842，黄点19.09930是什么意思？

作者: 范进中举 时间: 2011-11-8 16:12

钱潮滚滚发表于 2011-10-25 16:08 2 C5 B) f( Z% E; E8 C+ ?: v4 Q
发个图验证一下和大伙的结果有没有出入.
; r9 k' F4 F1 ~7 i3 A7 h5 @0 w; t- m) u: o8 |- X/ j
收盘价到昨天为止的三一重工:

蓝色线上的8.520041是啥？红点13.85195是啥？灰点15.40402是啥？黄点19.09793是啥？

作者: maodeng 时间: 2011-11-9 21:47

大牛，小生这厢有礼了，啥也不懂，照着葫芦画瓢，我的数据应该是更新到11月9号的，图形如下，不知道对不对

作者: maodeng 时间: 2011-11-10 19:58

wenty 发表于 2011-10-27 04:17
W3 J2 M* E* D4 ] i, `总结一下使用心得( W) k; E5 g3 a! t/ a
1）删掉set memory那一行
7 ~! C( w5 m7 ~7 I& \3 g0 _2）改路径

记事本打开后，怎么用tab

作者: 冰火 时间: 2011-11-10 22:44
认真学习

作者: higer 时间: 2011-12-9 22:34
我也曾尝试过使用stata软件对股票平均成本进行分析，但是作为一个非专业的炒股人员，而且日常还有工作的人而言，这种方法似乎显得太专业，而且操作难度极大，不一定实用。

使用这种专业的技术分析方法容易造成这么一种现象：大家都对导出的光滑漂亮的图像感觉良好，然而在股市中依然无法进退自如。

作者: jway315 时间: 2012-3-24 20:12
谢谢楼主，希望继续给同学们挖宝哈！

作者: awer 时间: 2012-3-28 18:06
改了一下，这次运算应该没啥问题了

capture program drop hello

program hello

quietly {

clear

cd e:\finance\stock\20110101\source

insheet using `1'.txt

keep if v5>0

gen vt=_n

tsset vt

gen zf=D1.v5

gsort -v1

gen avol=v7/v5

keep if avol>0

drop vt

gen vt=_n

tsset vt

sum v5 [aw =avol]

gen amean=r(mean) in 1

* 总加权平均值 amean

keep in 1/240

sum v5 [aw =avol]

gen mean240=r(mean) in 1

*240加权平均值 mean240

set obs 500

kdensity v5 [aw =avol],gaussian n(500) generate(dx240 dy240) nograph

*240 kdensity(dx240 dy240)

sum dy240

gen yxis=r(max) in 1

gen yxis2=yxis/2

*yxis

gen dy240max=r(max)

gen dy240minp=r(min)

gen vtt=_n

tsset vtt

gen ddy240=D1.dy240/D1.dx240

sum ddy240

gen ddy240max=r(max)

gen ddy240min=r(min)

gen ddy240g=ddy240*(dy240max-dy240min)/(ddy240max-ddy240min)

sum ddy240g

gen yxis3=r(min) in 1

gen yxis4=0 in 1

  *ddy240g

  preserve

  sum v5 [aw =avol]

  keep if dx240<r(mean)

  sum ddy240

  if r(min)<0 {

keep if ddy240<=0

count

if r(N)>2 {

   pkexamine dx240 dy240

local tempminp=r(tomc)
restore
   }

      else {

   di "cannot find minp"

restore
   sum dx240
   local tempminp=r(mean)

   }
}
else {

   keep if ddy240<=r(min)

   sum dx240

   local tempminp=r(min)
   restore
   }

gen minp=`tempminp' in 1
*240 minp

preserve

sum v5 [aw =avol]

keep if dx240>r(mean)

sum ddy240

if r(max)>0 {

keep if ddy240>=0

count

if r(N)>2 {

   pkexamine dx240 dy240

local tempmaxp=r(tomc)
restore
   }


      else {
restore
   di "cannot find maxp"

   sum dx240
   local tempmaxp=r(mean)

   }

   }

else {

   keep if ddy240>=r(max)

   sum dx240

   local tempmaxp=r(max)
   restore
   }



gen maxp=`tempmaxp' in 1

*240maxp

*-------------ddy240pminp左侧高峰值-----

preserve

keep if dx240<`tempminp'

keep if ddy240>0

count

if r(N)>2 {

pkexamine dx240 ddy240

local tempminphddy240=r(tomc)
restore
   }

      else {
restore
   di "cannot find pminphddy240"


   sum dx240
   local tempminphddy240=r(mean)


   }


gen pminphddy240=`tempminphddy240' in 1

*-------------ddy240pminp右侧低峰值 -----------------
preserve

keep if dx240>`tempminp'

keep if dx240<`tempmaxp'

keep if ddy240<0

gen temddy240=-ddy240

count

if r(N)>2 {

   pkexamine dx240 temddy240

  local tempminplddy240=r(tomc)
restore
   }

      else {

   di "cannot find pminplddy240"
restore

   sum dx240
   local tempminplddy240=r(mean)


   }

gen pminplddy240=`tempminplddy240' in 1

*-------------ddy240pmaxp左侧高峰值 -------------------

   preserve

keep if dx240>`tempminp'

keep if dx240<`tempmaxp'

keep if ddy240>0
count

if r(N)>2 {

   pkexamine dx240 ddy240

local tempmaxphddy240=r(tomc)

restore
   }

      else {

   di "cannot find pmaxphddy240"
   restore

   sum dx240
   local tempmaxphddy240=r(mean)


   }

gen pmaxphddy240=`tempmaxphddy240' in 1

*-------------ddy240pmaxp右侧低峰值 -------------------

   preserve

keep if dx240>`tempmaxp'

gen temddy240=-ddy240

keep if temddy240>0

count

if r(N)>2 {

pkexamine dx240 temddy240

local tempmaxplddy240=r(tomc)
restore
   }

      else {
restore
   di "cannot find pmaxlhddy240"


        sum dx240
   local tempmaxplddy240=r(mean) in 1


   }

gen pmaxplddy240=`tempmaxphddy240' in 1

*-------------------------------------

gen price5d=v5 in 1/20

gen price5yy=_n in 1/20

gen price5y=price5yy*(dy240max-dy240min)/20

*price5y &price5d

gen priceny=price5y in 1

gen pricen=v5 in 1


*priceny &pricen

*-------------------------graph------------------

  tw(line dy240 ddy240g dx240)(spike yxis amean,lc(blue))(spike yxis mean240, lc(red))(spike yxis maxp, lp(-) lc(black))(spike yxis minp, lc(black))(scatter yxis2 amean, lc(blue) msymbol(+) mlabel(amean)) (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)

* ---------可选项目----------------------

}
end

作者: 旅途愉快 时间: 2012-3-28 23:23

awer 发表于 2012-3-28 18:06
$ _& V! e: I- {% X& Q/ s改了一下，这次运算应该没啥问题了

这段代码是重新复制粘贴到hello.do里面的？

作者: awer 时间: 2012-3-29 11:40

旅途愉快发表于 2012-3-28 23:23
3 K4 z8 @$ V$ \这段代码是重新复制粘贴到hello.do里面的？

对的，就是hello里面的那个

作者: David09 时间: 2012-7-2 21:35

awer 发表于 2012-3-29 11:40
, h% `/ o* c$ n; E对的，就是hello里面的那个

. hello 002680
Obs. nos. out of range
r(198);

我运行后，出现这个，不晓得是什么错误了，请大拿们帮助

作者: qqysc 时间: 2012-7-3 12:12

David09 发表于 2012-7-2 21:35 6 ^7 C, {0 U6 R5 W1 D4 y* n. I- y
. hello 002680
5 n! i# o" a3 g8 \8 [7 K9 JObs. nos. out of range/ R. m6 D3 A9 [5 o! l0 W+ j
r(198);

估计是导出文本格式问题。

作者: qqysc 时间: 2012-7-3 12:14

awer 发表于 2012-3-29 11:40
7 z8 t8 G( P8 y+ K对的，就是hello里面的那个

复制粘贴有很多乱码，请楼主把上面的代码做成附件以供下载，谢谢！

作者: David09 时间: 2012-7-3 17:13
本帖最后由 David09 于 2012-7-3 17:28 编辑

qqysc 发表于 2012-7-3 12:14 5 s- N$ A' Y5 ^/ D7 l' W7 G& ?
复制粘贴有很多乱码，请楼主把上面的代码做成附件以供下载，谢谢！

我清除过了乱码的啊....

作者: qqysc 时间: 2012-7-3 17:56

David09 发表于 2012-7-3 17:13 : N% v* P% ]$ h) j3 {
我清除过了乱码的啊....

说得不是一回事。你说的问题可以试试在 002680的文本中，把tab键补齐，导出时不要带表头

作者: David09 时间: 2012-7-3 19:45

qqysc 发表于 2012-7-3 17:56 * i F$ q5 F8 N& Q
说得不是一回事。你说的问题可以试试在 002680的文本中，把tab键补齐，导出时不要带表头 ...

很奇怪啊，我看stata中间的数据表，都整整齐齐的导入了成交价成交量。但是一运行就蹦出来那个错误

作者: David09 时间: 2012-7-4 14:44
搞定了，原来是我选了一个次新股的数据，样本数目太少了，引起程序的错误。

感谢楼主给的好东东。

作者: 碧水芳洲 时间: 2012-7-13 13:40
怎么没人问曲线是什么意思呢？ density指标也不是很明白呢，谁来解释一下？

作者: benhuashu 时间: 2012-7-20 15:16
请问，这一个语句：
gen zf=D1.v5
是什意思？主要是=之后的D1.v5，这个"." 是什么意思？谢谢。

作者: 海阔 时间: 2012-11-5 18:37
非常感谢；

作者: 墨轩书童 时间: 2012-11-7 11:46
本帖最后由墨轩书童于 2012-11-7 14:37 编辑

感谢楼主！！

作者: 海阔 时间: 2013-1-20 17:03
本帖最后由海阔于 2013-1-20 17:08 编辑

按照版主的步骤操作到了最后出现

. hello 600031
avolDo not found
---------------------------------------换成其他代买，也是一样；总是出现这个问题，怎么办呢？奇怪

作者: maodeng 时间: 2013-1-25 13:15
大师，那两条曲线代表什么啊

作者: maodeng 时间: 2013-1-25 13:15
如果我想看最近一个月，三个月六个月的成本，怎么修改程序啊，谢谢大师了

作者: maodeng 时间: 2013-1-25 13:18
我用同样的方法导出数据，发现有的能生成，有的就不行，不知道问题在哪儿，比如300236就生成不了，我的邮箱是19860225gaoyuan@163.com，麻烦大师了

作者: maodeng 时间: 2013-1-28 13:13
. hello 300236
statement out of context
r(119);

大师就是这个

作者: sergio 时间: 2013-1-28 21:11
这贴实在！

作者: awer 时间: 2013-2-5 14:40

maodeng 发表于 2013-1-28 13:13 1 R6 r2 Q: T' `$ c5 \6 V
. hello 300236! _$ Y# z/ S7 ]/ Q5 k) c( m
statement out of context
7 H5 Z+ f0 M1 N& S( o5 Z. K& or(119);

抱歉，最近太忙了

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

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

capture program drop hello

program hello

quietly {

clear

cd D:\stock\test\source

insheet using `1'.txt

keep if v5>0

gen vt=_n

tsset vt

gen zf=D1.v5

gsort -v1

gen avol=v7/v5

keep if avol>0

drop vt

gen vt=_n

tsset vt

sum v5 [aw =avol]

gen amean=r(mean) in 1

* 总加权平均值 amean

keep in 1/240

sum v5 [aw =avol]

gen mean240=r(mean) in 1

*240加权平均值 mean240

set obs 500

kdensity v5 [aw =avol],gaussian n(500) generate(dx240 dy240) nograph

*240 kdensity(dx240 dy240)

sum dy240

gen yxis=r(max) in 1

gen yxis2=yxis/2

*yxis

gen dy240max=r(max)

gen dy240minp=r(min)

gen vtt=_n

tsset vtt

gen ddy240=D1.dy240/D1.dx240

sum ddy240

gen ddy240max=r(max)

gen ddy240min=r(min)

gen ddy240g=ddy240*(dy240max-dy240min)/(ddy240max-ddy240min)

sum ddy240g

gen yxis3=r(min) in 1

gen yxis4=0 in 1

  *ddy240g

  preserve

  sum v5 [aw =avol]

  keep if dx240<r(mean)

  sum ddy240

  if r(min)<0 {

keep if ddy240<=0

count

if r(N)>2 {

   pkexamine dx240 dy240

local tempminp=r(tomc)
restore
   }

      else {

   di "cannot find minp"

restore
   sum dx240
   local tempminp=r(mean)

   }
}
else {

   keep if ddy240<=r(min)

   sum dx240

   local tempminp=r(min)
   restore
   }

gen minp=`tempminp' in 1
*240 minp

preserve

sum v5 [aw =avol]

keep if dx240>r(mean)

sum ddy240

if r(max)>0 {

keep if ddy240>=0

count

if r(N)>2 {

   pkexamine dx240 dy240

local tempmaxp=r(tomc)
restore
   }


      else {
restore
   di "cannot find maxp"

   sum dx240
   local tempmaxp=r(mean)

   }

   }

else {

   keep if ddy240>=r(max)

   sum dx240

   local tempmaxp=r(max)
   restore
   }



gen maxp=`tempmaxp' in 1

*240maxp

*-------------ddy240pminp左侧高峰值-----

preserve

keep if dx240<`tempminp'

keep if ddy240>0

count

if r(N)>2 {

pkexamine dx240 ddy240

local tempminphddy240=r(tomc)
restore
   }

      else {
restore
   di "cannot find pminphddy240"


   sum dx240
   local tempminphddy240=r(mean)


   }


gen pminphddy240=`tempminphddy240' in 1

*-------------ddy240pminp右侧低峰值 -----------------
preserve

keep if dx240>`tempminp'

keep if dx240<`tempmaxp'

keep if ddy240<0

gen temddy240=-ddy240

count

if r(N)>2 {

   pkexamine dx240 temddy240

  local tempminplddy240=r(tomc)
restore
   }

      else {

   di "cannot find pminplddy240"
restore

   sum dx240
   local tempminplddy240=r(mean)


   }

gen pminplddy240=`tempminplddy240' in 1

*-------------ddy240pmaxp左侧高峰值 -------------------

   preserve

keep if dx240>`tempminp'

keep if dx240<`tempmaxp'

keep if ddy240>0
count

if r(N)>2 {

   pkexamine dx240 ddy240

local tempmaxphddy240=r(tomc)

restore
   }

      else {

   di "cannot find pmaxphddy240"
   restore

   sum dx240
   local tempmaxphddy240=r(mean)


   }

gen pmaxphddy240=`tempmaxphddy240' in 1

*-------------ddy240pmaxp右侧低峰值 -------------------

   preserve

keep if dx240>`tempmaxp'

gen temddy240=-ddy240

keep if temddy240>0

count

if r(N)>2 {

pkexamine dx240 temddy240

local tempmaxplddy240=r(tomc)
restore
   }

      else {
restore
   di "cannot find pmaxlhddy240"


        sum dx240
   local tempmaxplddy240=r(mean) in 1


   }

gen pmaxplddy240=`tempmaxphddy240' in 1

*-------------------------------------

gen price5d=v5 in 1/20

gen price5yy=_n in 1/20

gen price5y=price5yy*(dy240max-dy240min)/20

*price5y &price5d

gen priceny=price5y in 1

gen pricen=v5 in 1


*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)

* ---------可选项目----------------------
*(spike yxis amean,lc(blue))(scatter yxis2 amean, lc(blue) msymbol(+) mlabel(amean))
}
end

作者: maodeng 时间: 2013-2-7 00:03

awer 发表于 2013-2-5 14:40 / B& \. C7 g0 @" g+ j
抱歉，最近太忙了$ Z3 h5 Q4 j8 X; I! C

+ |! }1 X" F4 I* \& B300236 我用下面的程序试过，没问题

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

作者: 喜欢喝冰茶 时间: 2013-2-7 04:05
stata是不是商用的，R里应该有类似的东西吧

作者: 喜欢喝冰茶 时间: 2013-2-8 07:16

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

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

作者: 汉水云潇 时间: 2013-4-17 17:14
准备试一试
先找Stata软件

作者: 汉水云潇 时间: 2013-4-17 22:21
本帖最后由汉水云潇于 2013-4-17 22:27 编辑

晕呀
论坛复制程序代码
乱码夹杂在其中
很烦
对了两次
以为没问题了

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

作者: 汉水云潇 时间: 2013-4-18 07:50
试验成功，下面是三一重工的图

作者: peterbev 时间: 2013-4-26 12:46
先用用试一试，不知道准确不

作者: 中秋不圆 时间: 2016-3-12 15:46
噢。。。这个帖子好，比较直观。

作者: bknite 时间: 2016-3-21 23:51
这个帖子值得学习。

作者: jeanny 时间: 2016-3-22 10:41
感谢！！！

作者: Hubert2000 时间: 2016-3-23 10:23

明月回春发表于 2011-8-13 20:17( p" ] q$ ` @1 G. K: i Z/ t4 o% M
程序写的不错.
# i/ B1 M, Y3 `. w; \不过有一点要注意.通达信导出的数据,只对价格复权,并没有对成交量复权.所以数据需要手工处理 ...

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

作者: axiutime 时间: 2016-8-3 13:28
提示: 该帖被管理员或版主屏蔽

欢迎光临爱吱声 (http://aswetalk.net/bbs/)