glmm 0.1.3

Standalone f64 GLMM fit kernels (OLS, GLM, LMM, GLMM) in pure Rust on faer — the validation-pinned numerics from the MCPower engine.
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27

//! `StartValues` — the warm-start primitive.
//!
//! Raw optimizer state, *not* high-level variances: the dominant use is threading
//! the previous fit's state into the next refit (loop tier / future inference
//! layer), so it mirrors what the kernels already hold (`ws.params`/`theta`, the β
//! estimate). The stable `fit_warm` takes this as an optional warm start; `fit_cold`
//! (≡ `fit_warm(.., None, ..)`) omits it and cold-starts the kernels with
//! `theta_start = None`. Shared internal primitive: `pub(crate)`
//! always, and now unconditionally `pub` (a stable input to `fit_warm`).
//!
//! Carries only `beta`/`theta`: Gamma φ is profiled (`family::gamma_aic`, no
//! warm-startable state) and the GLMM neg-binomial θ search is a global bracket
//! (`fit::golden_max_ln_theta`), so neither warm-starts anything reachable here.

/// Raw optimizer warm-start state for one model fit.
///
/// - `beta`: fixed-effect start, length `p` (cold = all-zero).
/// - `theta`: RE Cholesky parameters (the kernel's `θ`, column-major vech),
///   length `n_theta` for the model's RE structure (cold = the kernel's flat
///   `THETA0` blind start). Empty for fixed-only (OLS/GLM) models.
#[derive(Debug, Clone, PartialEq)]
pub struct StartValues {
    /// Warm-start β, length `p`.
    pub beta: Vec<f64>,
    /// Warm-start θ (RE Cholesky vech), length `n_theta`; empty for fixed-only models.
    pub theta: Vec<f64>,
}