round(trend_ar1$model_params[["sigma_hatsq"]], 3) (`actual`) not equal to 0.29 (`expected`).

  `actual`: 0.5
`expected`: 0.3

expect_s3_class(fit_meanshift_norm(CET, tau = NULL), "mod_cpt")
expect_equal(changepoints(fit_meanshift_norm(CET, tau = NULL)), integer())
y <- fit_meanshift_norm_ar1(CET, tau = ids)
expect_true(all(abs(fitted(z) - fitted(x)) < 0.000000001))
expect_true(all(abs(fitted(w) - fitted(y)) < 0.000000001))
trend_wn <- fit_trendshift(CET, tau = ids)
trend_ar1 <- fit_trendshift_ar1(CET, tau = ids)
spline_wn <- fit_lmshift(CET, tau = ids, deg_poly = 4)
trend_wn_trunc <- fit_lmshift(CET_trunc, tau = tau, deg_poly = 2)
MDL(trend_wn_trunc) + 2 * log(nobs(trend_wn_trunc))
trend_ar1_trunc <- fit_lmshift_ar1(CET_trunc, tau = tau, deg_poly = 1)
MDL(trend_ar1_trunc) + 2 * log(nobs(trend_ar1_trunc))

round(as.numeric(logLik(trend_ar1))) (`actual`) not equal to round(-288.8) (`expected`).

  `actual`: -403
`expected`: -289

round(BIC(trend_ar1)) (`actual`) not equal to round(654.19) (`expected`).

  `actual`: 883
`expected`: 654

floor(MDL(trend_ar1)) (`actual`) not equal to floor(656.52) (`expected`).

  `actual`: 885
`expected`: 656

round(trend_ar1$model_params[["phi_hat"]], 3) (`actual`) not equal to 0.058 (`expected`).

  `actual`: 1.0
`expected`: 0.1

round(model_variance(trend_ar1_trunc), 3) (`actual`) not equal to 0.305 (`expected`).

  `actual`: 0.6
`expected`: 0.3

round(trend_ar1_trunc$model_params[["phi_hat"]], 3) (`actual`) not equal to 0.073 (`expected`).

  `actual`: 1.0
`expected`: 0.1