search

beanumber / tidychangepoint

Changepoint detection with a tidy interface
https://beanumber.github.io/tidychangepoint/
GNU General Public License v3.0
1 stars 1 forks source link

BMDL scores don't match in glance() #21

Closed beanumber closed 9 months ago

beanumber commented 9 months ago 
> x <- segment(DataCPSim, method = "cpt-gbmdl")
method: cpt-gbmdl
  |============================================================| 100%
There were 50 or more warnings (use warnings() to see the first 50)
> glance(x)
# A tibble: 1 × 13
  pkg             version    algorithm   test_stat  BMDL logLik     AIC   BIC num_cpts s3_logLik s3_AIC s3_BIC s3_BMDL
  <chr>           <pckg_vrs> <chr>           <dbl> <dbl> <logLik> <dbl> <dbl>    <int> <logLik>   <dbl>  <dbl>   <dbl>
1 tidychangepoint 0.0.0.9014 GeneticBMDL      23.3  23.3 -3248     6504 6524.        4 -3248       6504  6524.    712.
> bmdl(as.ts(x), changepoints(x))
[1] 712.1098