search

rkillick / changepoint

A place for the development version of the changepoint package on CRAN.
127 stars 33 forks source link

This code causes R to crash #83

Open Jiayou-Chao opened 2 days ago

Jiayou-Chao commented 2 days ago

Is this the place for reporting bugs in changepoint.np? Running the following the code will cause the R to crash.

library(changepoint.np)
values <- c(5.22928571428571, 5.19801887517172, 5.0708391923999, 4.52608333333333,                                  
  4.50230323874841, 4.5346807724015, 4.36199173874068, 4.02545652239261, 
  3.69373020979002, 3.51172563631839, 3.58037666802371, 3.61083333333333, 
  3.63166666666667, 3.64308333333333, 3.60344317121439, 3.62616666666667, 
  3.65857132450011, 3.64744365601782, 3.67034674500229, 3.65158333333333, 
  3.66133333333333, 3.64316666666667, 3.64408333333333, 3.61608333333333, 
  3.63991666666667, 3.64866666666667, 3.64966666666667, 3.65391666666667, 
  3.65558333333333, 3.64075, 3.68130934657094, 3.66516666666667, 
  3.7165, 3.67241666666667, 3.693, 3.69258333333333, 3.68628915246162, 
  3.72283333333333)
cpt.np(values, penalty = "MBIC", method = "PELT", test.stat = "empirical_distribution", minseglen = 20)
rkillick commented 1 day ago

This is the place for issues for the changepoint package rather than changepoint.np but I am maintainer for both so here is fine too. Thanks for reporting! I will take a look at this and get back to you. I ran your code the first time and it completed successfully. Then I went to break it down and the next time I ran it, R crashed. So I need to look into it more and try to pin down what is going on and when the problem is triggered and why.

Jiayou-Chao commented 1 day ago

Thank you!