[程式] R changepoint package使用

[eri820503 (snoopyahoo)]

[軟體程式類別]:
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://www.ptt.cc/bbs/Statistics/M.1453818899.A.810.html
→ Godkin: DNAcopy: https://goo.gl/RZtPEb  (from Bioconductor)
→ eri820503: 謝謝~