[程式] R changepoint package使用
[軟體程式類別]:
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
→
01/27 20:08, , 1F
01/27 20:08, 1F
→
01/27 20:23, , 2F
01/27 20:23, 2F