antedep 0.2.0

New features

• New S3 methods deviance(), confint(), and vcov() for gau_fit, cat_fit, and inad_fit objects, improving compatibility with standard R model-fitting workflows.
• New summary.partial_corr() method prints a lag-by-lag table of intervenor-adjusted partial correlations (mean absolute value, range, and number of significant pairs at each lag).

Improvements

• Matrix inversion in ci_inad(), fit_gau_em(), and cat_test_stats.R now uses chol2inv(chol(.)) instead of solve() for symmetric positive-definite matrices, improving numerical stability and efficiency.

antedep 0.1.0

New features

• fit_inad() gains a nb_inno_size_ub argument (default 50) that caps the upper bound of the negative-binomial innovation size parameter during optimization, improving numerical stability for near-Poisson data.
• test_order_gau() accepts order_null and order_alt as convenience aliases for p and the absolute alternative order; both are also returned in the result object.

Bug fixes

• ci_inad(): fixed a sign error in the observed Fisher information for the negative-binomial innovation size parameter; the Hessian term (r + u) / (r + λ)² was added instead of subtracted, producing confidence intervals that were too wide.
• ci_inad(): the numerical second derivative for nb_inno_size CIs now retries with progressively smaller step sizes (×0.1, ×0.01) before falling back to NA, avoiding spurious failures when the default step lands in a non-finite region.
• test_homogeneity_inad(): degrees of freedom for LRT tests involving innovation = "nbinom" are now computed from the actual number of NB size parameters in the fitted models rather than assuming a fixed count of 1. This corrects LRT statistics and p-values whenever nb_inno_size is fitted as a time-varying vector.
• ci_inad() tau profile CI: nb_inno_size (negative-binomial innovation dispersion) is now held fixed at its full-model MLE during profile refits, consistent with the constrained-fit paradigm used throughout the package. Previously it was re-optimised as a nuisance parameter, which could widen the interval to the point of crossing zero even when the LRT clearly rejects the null (Variant 1 vs Variant 2 fix).
• ci_inad() tau profile CI: the bracket search in .ci_tau_profile_inad no longer imposes an artificial upper cap (max(|tau_mle| + 1, 1)) on the search range. The maximum bracket iterations are increased from 20 to 50 and the initial step size is set to max(0.1, |tau_mle| * 0.2), preventing the search from stalling for large or near-zero MLEs.

Initial release

• Initial CRAN submission candidate for Gaussian AD, categorical AD, and INAD workflows.
• Added/expanded examples for key user-facing modeling functions (fit_*, em_*, simulate_*, logL_*).
• Refreshed missing-data notes and harmonized documentation links and metadata.
• logL_gau() default missing-data behavior is now na_action = "fail" (previously marginalization-first in earlier drafts). For missing inputs, pass na_action = "marginalize" or na_action = "complete" explicitly.
• Added packaged datasets labor_force_cat (categorical labor-force sequences) and race_100km (continuous 100km race split times).