feat(psychometric): restore Driver p.16 asymDIFFUSIONstd p/p=1 on main

CHANGELOG.md

@@ -36,6 +36,7 @@ All notable changes to TEPP are documented here. The format follows Keep a Chang
 
 ## [Unreleased]
 
+- `psychometric_core` recovers the Driver, Oud, and Voelkle (2017, p. 16 `asymDIFFUSIONstd`; footnote 4; Eq. 4, p. 5; Table 2, p. 12; 2017-era ctsem `summary.ctsemFit.R`; JSS PDF re-opened 2026-08-26T17:20Z from https://www.jstatsoft.org/index.php/jss/article/download/v077i05/1104) scalar standardised asymptotic within-subject variance on current main after `0ce16e8` dropped the pre-consolidation code while research notes already named the map (register items 89–90). Page 16 prints standardised matrices with the suffix `std` when appropriate, and names `asymDIFFUSION` the total within-subject variance as `Δt → ∞`. The printed example on p. 16 is `discreteDRIFTstd`, not `asymDIFFUSIONstd`. Footnote 4 standardises using only the relevant variance, not the total. The relevant variance for that named process-dynamics correlation is within-subject `asymDIFFUSION` `p = −q / (2 a)`, not free first-occasion `T0VAR`. The 2017-era source forms `asymDIFFUSIONstd` as `solve(sqrt(diag(asymDIFFUSION))) %&% asymDIFFUSION` when `verbose = TRUE`. OpenMx `%&%` is `t(A) %*% B %*% A`. The default `ridging = FALSE` adds 0, not `0.0001`; that ridge is a numerical hack and is not this exact map. The scalar correlation is `p / p = 1` after strictly positive `p`. Form strictly positive `p` first, then `1 / √p`, then `(1 / √p) p (1 / √p)`. Unstandardised `p` is defined for a zero process; standardised `asymDIFFUSION` is not. Zero `q` has no positive SD and fails closed. Lasting `p` requires stable `a < 0`. A non-event clock fails closed. Distinct positive `p` recover the same 1. `p_0 / p_0 = 1` is `T0VARstd` and recovers the same number and remains a distinct named quantity. `q / p = −2 a` is `DIFFUSIONstd` and is not this correlation. `v / v = 1` is `TIPREDVARstd` and recovers the same number and remains a distinct named quantity. Meredith (1993) remains unread (Unpaywall 2026-08-26T17:20Z: `is_oa: false`; OpenAlex closed; Springer `content/pdf` is an HTML stub). Mislevy (1991, *Psychometrika, 56*, 177–196) remains unread on the same terms (DOI `10.1007/bf02294457`; Unpaywall `is_oa: false`). Still not a Kalman filter, not a matrix `expm`, not ESEM estimation, not DSEM, and not ctsem estimation.
 - `psychometric_core` recovers the Driver, Oud, and Voelkle (2017, p. 16 `T0VARstd`; Table 2, p. 12; footnote 4; 2017-era ctsem `summary.ctsemFit.R`; JSS PDF re-opened 2026-08-26T07:17Z from https://www.jstatsoft.org/index.php/jss/article/download/v077i05/1104) scalar standardised initial latent variance on current main after `0ce16e8` dropped the pre-consolidation code while research notes already named the map (register items 79–80). Page 16 prints standardised matrices with the suffix `std` when appropriate. The printed example on p. 16 is `discreteDRIFTstd`, not `T0VARstd`. Footnote 4 standardises using only the relevant variance, not the total. Table 2 names `T0VAR` the latent process initial variance/covariance. The first-occasion relevant variance is free `T0VAR` `p_0`, not process-dynamics `asymDIFFUSION` `-q / (2 a)`. The 2017-era source forms `T0VARstd` as `solve(sqrt(diag(T0VAR))) %&% T0VAR` when `verbose = TRUE`. OpenMx `%&%` is `t(A) %*% B %*% A`. The default `ridging = FALSE` adds 0, not `0.0001`; that ridge is a numerical hack and is not this exact map. The scalar correlation is `p_0 / p_0 = 1` after strictly positive `p_0`. Form strictly positive `p_0` first, then `1 / √p_0`, then `(1 / √p_0) p_0 (1 / √p_0)`. A zero first-occasion variance has no positive SD and fails closed. `T0` is an event-time occasion, so a non-event clock fails closed. Free `T0VAR` does not require stable `a < 0`. Distinct positive `p_0` recover the same 1. `μ_0 / √p_0` is `T0MEANSstd` and recovers the same number when `μ_0 = √p_0` and remains a distinct named quantity. `p / p = 1` is `asymDIFFUSIONstd` and recovers the same number and remains a distinct named quantity. Meredith (1993) remains unread (Unpaywall 2026-08-26T07:17Z: `is_oa: false`; OpenAlex closed; Springer `content/pdf` is a 3038-byte HTML stub). Mislevy (1991, *Psychometrika, 56*, 177–196) remains unread on the same terms (DOI `10.1007/bf02294457`; Unpaywall `is_oa: false`). Still not a Kalman filter, not a matrix `expm`, not ESEM estimation, not DSEM, and not ctsem estimation.
 - `psychometric_core` recovers the Driver, Oud, and Voelkle (2017, p. 16 `asymCINTstd`; Eq. 3, p. 4; Table 2, p. 12; footnote 4; 2017-era ctsem `summary.ctsemFit.R`; JSS PDF re-opened 2026-08-26T00:20Z from https://www.jstatsoft.org/index.php/jss/article/download/v077i05/1104) scalar standardised asymptotic continuous intercept on current main after consolidation dropped the pre-consolidation `(-κ / a) / √p` slice. Page 16 prints standardised matrices with the suffix `std` when appropriate, and asymptotic values as `Δt → ∞`. Footnote 4 standardises using only the relevant variance, not the total. Table 2 names `κ` `CINT`. The relevant variance for that process intercept is within-subject `asymDIFFUSION` `p = −q / (2 a)`. The 2017-era source forms unstandardised `asymCINT` whenever `verbose = TRUE` as `-solve(DRIFT) %*% CINT` and does not form an `asymCINTstd` matrix. Form strictly positive `p` first, then the asymptotic intercept, then divide by `√p`. A zero intercept is exactly zero after that positive SD. Zero `q` has no positive process SD and fails closed. Lasting `p` requires stable `a < 0`. A non-event clock fails closed. `κ / √p` is `CINTstd` and is not this total-change map. `A^{-1}[e^{A Δt} − I] κ / √p` is `discreteCINTstd` and depends on the event interval. Meredith (1993) remains unread (Unpaywall 2026-08-25T18:22Z: `is_oa: false`; Springer `content/pdf` is a 3038-byte HTML stub). Mislevy (1991) remains unread on the same terms. Still not a Kalman filter, not a matrix `expm`, not ESEM estimation, not DSEM, and not ctsem estimation.
 - `psychometric_core` recovers the Driver, Oud, and Voelkle (2017, p. 16 `T0MEANSstd`; Table 2, p. 12; footnote 4; 2017-era ctsem `summary.ctsemFit.R`; JSS PDF re-opened 2026-08-26T04:09Z from https://www.jstatsoft.org/index.php/jss/article/download/v077i05/1104) scalar standardised initial latent mean on current main after `0ce16e8` dropped the pre-consolidation code while research notes already named the map. Page 16 prints standardised matrices with the suffix `std` when appropriate. Footnote 4 standardises using only the relevant variance, not the total. Table 2 names `T0MEANS` the `n.latent × 1` matrix of latent process means at the first time point `T0` and names `T0VAR` the latent process initial variance/covariance. The first-occasion relevant variance is free `T0VAR` `p_0`, not process-dynamics `asymDIFFUSION` `-q / (2 a)`. The 2017-era source forms unstandardised `T0MEANS` and does not form a `T0MEANSstd` matrix; the scalar map is `μ_0 / √p_0` after strictly positive `p_0`. A zero mean is exactly zero. Zero `p_0` has no positive SD and fails closed. `T0` is an event-time occasion, so a non-event clock fails closed. Free `T0MEANS` does not require stable `a < 0`. `p_0 / p_0 = 1` recovers the same number when `μ_0 = √p_0` and remains a distinct named quantity. `μ_0 / √asymDIFFUSION` uses process-dynamics variance and is not this first-occasion map. Meredith (1993) remains unread (Unpaywall 2026-08-26T00:22Z: `is_oa: false`; OpenAlex closed; Springer `content/pdf` is a 3038-byte HTML stub). Mislevy (1991, *Psychometrika, 56*, 177–196) remains unread on the same terms. Still not a Kalman filter, not a matrix `expm`, not ESEM estimation, not DSEM, and not ctsem estimation.

crates/psychometric_core/src/error.rs

@@ -585,6 +585,33 @@ pub enum PsychometricError {
     /// `asymDIFFUSIONstd` is the correlation form of process-
     /// dynamics `asymDIFFUSION`.
     StandardisedAsymptoticDiffusionIsNotStandardisedInitialLatentVariance,
+    /// Driver p. 16 `asymDIFFUSIONstd` was requested without a
+    /// strictly positive `asymDIFFUSION`. Footnote 4
+    /// standardisation of the 2017-era `asymDIFFUSION` matrix
+    /// requires strictly positive `−q / (2 a)`.
+    StandardisedAsymptoticDiffusionRequiresPositiveStationaryVariance,
+    /// Driver p. 16 unstandardised `asymDIFFUSION` `p` was treated
+    /// as `asymDIFFUSIONstd`. Unstandardised `p` is defined for a
+    /// zero process; standardised `asymDIFFUSION` is not.
+    UnstandardisedAsymptoticDiffusionIsNotStandardisedAsymptoticDiffusion,
+    /// Driver p. 16 `T0VARstd` was treated as p. 16
+    /// `asymDIFFUSIONstd`. Equal numbers of 1 after a strictly
+    /// positive relevant variance are still distinct named
+    /// quantities. `asymDIFFUSIONstd` is the correlation form of
+    /// process-dynamics `asymDIFFUSION`; `T0VARstd` is the
+    /// correlation form of free `T0VAR`.
+    StandardisedInitialLatentVarianceIsNotStandardisedAsymptoticDiffusion,
+    /// Driver p. 16 `DIFFUSIONstd` `q / p = −2 a` was treated as
+    /// `asymDIFFUSIONstd`. The continuous-diffusion ratio is not
+    /// the correlation form of `asymDIFFUSION`.
+    StandardisedContinuousDiffusionIsNotStandardisedAsymptoticDiffusion,
+    /// Driver p. 16 `TIPREDVARstd` was treated as p. 16
+    /// `asymDIFFUSIONstd`. Equal numbers of 1 after a strictly
+    /// positive relevant variance are still distinct named
+    /// quantities. `asymDIFFUSIONstd` is the correlation form of
+    /// process-dynamics `asymDIFFUSION`; `TIPREDVARstd` is the
+    /// correlation form of `TIPREDVAR`.
+    StandardisedTimeIndependentPredictorVarianceIsNotStandardisedAsymptoticDiffusion,
     /// Driver p. 16 `discreteCINTstd` was requested without a strictly
     /// positive `asymDIFFUSION`. Footnote 4 standardises using only the
     /// relevant variance; zero `q` has no positive process SD.
@@ -1063,6 +1090,21 @@ impl fmt::Display for PsychometricError {
             Self::StandardisedAsymptoticDiffusionIsNotStandardisedInitialLatentVariance => {
                 "standardised asymptotic diffusion is not standardised initial latent variance"
             }
+            Self::StandardisedAsymptoticDiffusionRequiresPositiveStationaryVariance => {
+                "standardised asymptotic diffusion requires strictly positive stationary within-subject variance"
+            }
+            Self::UnstandardisedAsymptoticDiffusionIsNotStandardisedAsymptoticDiffusion => {
+                "unstandardised asymptotic diffusion is not standardised asymptotic diffusion"
+            }
+            Self::StandardisedInitialLatentVarianceIsNotStandardisedAsymptoticDiffusion => {
+                "standardised initial latent variance is not standardised asymptotic diffusion"
+            }
+            Self::StandardisedContinuousDiffusionIsNotStandardisedAsymptoticDiffusion => {
+                "standardised continuous diffusion is not standardised asymptotic diffusion"
+            }
+            Self::StandardisedTimeIndependentPredictorVarianceIsNotStandardisedAsymptoticDiffusion => {
+                "standardised time-independent predictor variance is not standardised asymptotic diffusion"
+            }
             Self::StandardisedDiscreteContinuousInterceptRequiresPositiveStationaryVariance => {
                 "standardised discrete continuous intercept requires strictly positive stationary within-subject variance"
             }
@@ -1802,6 +1844,35 @@ mod tests {
         );
     }
 
+    #[test]
+    fn standardised_asymptotic_diffusion_boundary_messages_are_stable() {
+        assert_eq!(
+            PsychometricError::StandardisedAsymptoticDiffusionRequiresPositiveStationaryVariance
+                .to_string(),
+            "standardised asymptotic diffusion requires strictly positive stationary within-subject variance"
+        );
+        assert_eq!(
+            PsychometricError::UnstandardisedAsymptoticDiffusionIsNotStandardisedAsymptoticDiffusion
+                .to_string(),
+            "unstandardised asymptotic diffusion is not standardised asymptotic diffusion"
+        );
+        assert_eq!(
+            PsychometricError::StandardisedInitialLatentVarianceIsNotStandardisedAsymptoticDiffusion
+                .to_string(),
+            "standardised initial latent variance is not standardised asymptotic diffusion"
+        );
+        assert_eq!(
+            PsychometricError::StandardisedContinuousDiffusionIsNotStandardisedAsymptoticDiffusion
+                .to_string(),
+            "standardised continuous diffusion is not standardised asymptotic diffusion"
+        );
+        assert_eq!(
+            PsychometricError::StandardisedTimeIndependentPredictorVarianceIsNotStandardisedAsymptoticDiffusion
+                .to_string(),
+            "standardised time-independent predictor variance is not standardised asymptotic diffusion"
+        );
+    }
+
     #[test]
     fn standardised_discrete_continuous_intercept_boundary_messages_are_stable() {
         assert_eq!(