作者eri820503 (snoopyahoo)
看板Statistics
标题[程式] R changepoint package使用
时间Tue Jan 26 22:34:55 2016
[软体程式类别]:
R
[程式问题]:
change point的寻找
[软体熟悉度]:
熟悉
[问题叙述]:
最近在看change point相关的paper,想先用package找出change point
看paper上讲常用的为"changepoint","DNAcopy"这两个
但是"DNAcopy"我在R上面找不到??
(不知道有没有其他人用过,paper作者写说这方法满热门的)
而"changepoint"我是使用cpt.meanvar这个函数,但过程中有些问题
想问一下有没有人有过相关经验
[程式范例]:
x=c(rnorm(50,0,1),rnorm(50,5,3),rnorm(50,10,1),rnorm(50,3,10))
cpt.meanvar(x,penalty="Manual",pen.value="4*log(n)",method="BinSeg",
Q=5,class=FALSE)
上面这一段是package上的范例,
而我将x改成我的资料,後面都不变的情况,出现了:
Error in if (pen.return < 0) { : missing value where TRUE/FALSE needed
而将後面设定改成其它的也无法执行
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 36.236.243.83
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/Statistics/M.1453818899.A.810.html
2F:→ eri820503: 谢谢~ 01/27 20:23