Read and follow AGENTS.md before changing this repository. The repository-wide contracts in that file are normative. Use DOCUMENTATION.md to locate the canonical PRD/TRD/Architecture/UML/ERD/API/security/privacy/operability/traceability authorities.
- Use test-driven development for every behavior change.
- Keep one branch/PR internally coherent and independently reviewable, but do not stop the invocation merely because one bounded slice completed while another safe action exists.
- After every mutation, merge, proof, or defer decision, re-enumerate the executable queue; waiting on one branch is local, not a repository-wide blocker.
- Prefer explicit types and small modules with stable interfaces.
- Preserve source spans, temporal provenance, uncertainty, purpose/authorization, and model-version metadata end to end.
- Preserve standalone operation and integrate with other CWL services only through versioned contracts; never use hidden cross-service database coupling.
- Do not replace statistical estimation with an LLM judgment.
- Do not convert association, temporal precedence, or document links into causal language without identification evidence.
- Do not remove repeated report language with global stopword lists or use TF-IDF/BM25 as inferential weights. Model template, section, copied-text, style, modality, and corpus-background sources explicitly.
- Do not treat raw topic proportions as ordinary Euclidean indicators. Use logistic-normal coordinates or valid log-ratio coordinates and propagate posterior uncertainty into ESEM/DSEM.
- Do not treat metric/weak invariance as a latent-mean license. Strong (equal loading and intercept) or strict is required;
#84metriclicenses shared metric meaning only. Putnick and Bornstein (2016, PMC5145197 opened 2026-08-19T22:15Z) require scalar invariance before latent-mean comparison; residual invariance is not a prerequisite. Two-observation series have no residual degrees of freedom (ordinary_least_squares_fitreturns residual variance0) and cap at strong/scalar; they still license means. This is two-group OLS, not MGCFA. Meredith (1993) names remain unread labels (Unpaywall/OpenAlex 2026-08-25T11:32Z: closed). - Do not use the difference quotient as a continuous-time rate. The scalar map is
a = ln(φ) / Δton event time. Discrete lags from unequal event intervals are not one coefficient; remap them through that log-rate. Binary64exp(a Δt) = 0is not a discrete lag. A constant predictor's discrete effect is Voelkle et al. (2012, Eq. 12), evaluated asa_yx (expm1(z) / a_xx)withz = a_xx Δtso a finite result is not lost whenzoverflows to-∞or whena_yx Δtoverflows. Whenexpm1(z)overflows at a finitez, rewrite in log space; a zero continuous effect is exactly zero; an overflowinga_yx/a_xxrewrite term fails closed. The first-order product is the underflow limit of that equation, not the general constant-predictor discrete effect. A time-varying predictor whose sampling interval equals its constancy interval uses Voelkle et al. (2012, Eq. 14):b* = a_yx Δt. Unmatched intervals fail closed (Oud & Jansen, 2000, unread). Discrete process noise is Driver et al. (2017, Eq. 3):Q_Δt = 0.5 q (expm1(z) / a)withz = 2 (a Δt)andq = G G⊤ ≥ 0; do not form2 afirst;a = 0andz → 0recoverq Δt; a zero diffusion is exactly zero; an overflowing rewrite scale0.5 q / afails closed; this is not a Kalman filter.Q_Δtiscov(η_t | η_{t-1}), notVar(η_t). The lagged covariance isexp(a Δt) pand the unconditional variance isexp(2 a Δt) p + Q_Δt(Driver et al., 2017, Eq. 3–4, pp. 4–5; JSS has no numbered §2.2). A zero diffusion whose2 (a Δt)overflows to+∞is not a finiteVar(η_t). The stationary within-subject variance is theΔt → ∞limit of Eq. 4:-q / (2 a)for stablea < 0(JSS p. 16asymDIFFUSION; §4.3). When2 ais finite, formq / -(2 a)soq / aoverflow does not lose a finite result (q = MAX,a = -0.75→MAX / 1.5). When2 aoverflows, form(q / a) * -0.5. Do not form0.5 qfirst (q = from_bits(1)underflows).a ≥ 0has no finite stationary variance. Finite-intervalQ_Δtis not that limit. Trait-plus-state variance istrait + stateand lagged covariance istrait + exp(a Δt) p(Driver et al., 2017, §4.3, p. 9). Trait variance is not process noise and notasymDIFFUSION. Evolving the summed variance as if it were all state is not that map. This is not RI-CLPM. Observed-indicator variance isλ² Var(η) + θwhenMANIFESTTRAITVARis zero andλ² Var(η) + θ + ψotherwise (Driver et al., 2017, Eq. 5, p. 5; Table 2, p. 12). Lagged observed covariance isλ² cov(η_t, η_{t-1}) + ψ;MANIFESTVARdoes not enter. Observed-indicator mean isτ + λ μ(Driver et al., 2017, Eq. 5; Table 2, p. 12).MANIFESTMEANSisτ, notE(y).E(η)is notE(y).CINTis notMANIFESTMEANS.T0MEANSis notE(y). The discrete latent mean isμ_t = exp(a Δt) μ_0 + (exp(a Δt) − 1)/a κ(Driver et al., 2017, Eq. 3, p. 4; Table 2, p. 12).T0MEANSis notμ_t.CINTis not that discrete increment. A zero drift isκ Δt. Underflow ofexp(a Δt)to+0drops the carriedT0MEANSand keeps−κ / a. The evolved observed mean isτ + λ μ_t(Driver et al., 2017, Eq. 5 of that Eq. 3 map). The first-occasion mapτ + λ μ_0is notE(y_t).μ_tis notE(y_t). The contemporaneous time-dependent predictor impulse ism x(Driver et al., 2017, Eq. 3 fourth summand; Table 2TDPREDEFFECTisM). Formμ_tfirst, then addm x.TDPREDEFFECTis notCINT.M xis notA^{-1}[e^{A Δt} − I] B zand is not Voelkle et al. (2012, Eq. 14). The §7.2 level-change form is not that impulse. The observed mean of that contemporaneous impulse isτ + λ(μ_t + m x)(Driver et al., 2017, Eq. 5 of the Eq. 3 fourth-summand composition). The evolved mapτ + λ μ_tis not that observed mean. The carry mapτ + λ(μ_t + e^{a(t−u)} m x)is not that observed mean whenu ≠ t. The evolved-plus-impulse latent mean is notE(y_t). The time-independent predictor increment isA^{-1}[e^{A Δt} − I] B z(Driver et al., 2017, Eq. 3 second summand; Table 2TIPREDEFFECTisB). FormB zfirst, then the discrete intercept map. A zero drift isB z Δt.TIPREDEFFECTisB, not that discrete increment.A^{-1}[e^{A Δt} − I] B zis notCINT, notM x, and not Voelkle et al. (2012, Eq. 14). The observed mean of that increment isτ + λ(μ_t + A^{-1}[e^{A Δt} − I] B z)(Driver et al., 2017, Eq. 5 of the Eq. 3 printed addend after theT0MEANScarry and theCINTincrement). The evolved mapτ + λ μ_tis not that observed mean. The contemporaneous mapτ + λ(μ_t + m x)is not that observed mean. The carry mapτ + λ(μ_t + e^{a(t−u)} m x)is not that observed mean whenu ≠ t. The evolved-plus-increment latent mean is notE(y_t). The within-interval time-dependent impulse carry ise^{A(t−u)} M xfort0 < u < t(Driver et al., 2017, Eq. 1–2 Green-function integral of Eq. 2; §7.2 dissipation). Formm xfirst, thene^{a(t−u)} m x. A zero drift ism xwith no dissipation. Underflow ofe^{a(t−u)}to+0is vanishing dissipation and is kept.e^{A(t−u)} M xis not the contemporaneous Dirac, notCINT, notTIPREDEFFECT, and not Voelkle et al. (2012, Eq. 14). An impulse atu = tis the contemporaneous map. An impulse atu ≤ t0is already inη(t0). The observed mean of that carry isτ + λ(μ_t + e^{a(t−u)} m x)(Driver et al., 2017, Eq. 5 of the Eq. 1–2 carried latent mean). The evolved mapτ + λ μ_tis not that observed mean. The contemporaneous mapτ + λ(μ_t + m x)is not that observed mean whenu ≠ t.MANIFESTMEANSis notE(y_t). The carried latent mean is notE(y_t). The first-occasion time-independent predictor shift ist0_b z(Driver et al., 2017, Table 3T0TIPREDEFFECT; Eq. 3 first summand). Formt0_b zfirst, thene^{a Δt} t0_b z. Formμ_tfirst, then add that carry. A zero drift ist0_b z. Underflow ofe^{a Δt}to+0is a vanishing carry of the first-occasion shift and is kept.t0_b zis notA^{-1}[e^{A Δt} − I] B z, notCINT, and notM x.e^{A Δt} t0_b zis nott0_b z.T0TIPREDEFFECTis the coefficient, not the shift. The observed mean of that first-occasion carry isτ + λ(μ_t + e^{a Δt} t0_b z)(Driver et al., 2017, Eq. 5 of the Table 3 / Eq. 3 first-summand composition). The evolved mapτ + λ μ_tis not that observed mean. The process-increment mapτ + λ(μ_t + A^{-1}[e^{A Δt} − I] B z)is not that observed mean. The contemporaneous mapτ + λ(μ_t + m x)is not that observed mean. The impulse-carry mapτ + λ(μ_t + e^{a(t−u)} m x)is not that observed mean whenu ≠ t0. The evolved-plus-carry latent mean is notE(y_t). The first-occasion time-dependent predictor shift ist0_m x0(Driver et al., 2017, Table 3T0TDPREDEFFECT; Eq. 3 first summand; JSS PDF re-opened 2026-08-20T19:10Z). Formt0_m x0first, thene^{a Δt} t0_m x0. Formμ_tfirst, then add that carry. A zero drift ist0_m x0. Underflow ofe^{a Δt}to+0is a vanishing carry of the first-occasion shift and is kept.t0_m x0is notM x, note^{A(t−u)} M xfort0 < u < t, nott0_b z, notA^{-1}[e^{A Δt} − I] B z, and notCINT.e^{A Δt} t0_m x0is nott0_m x0.T0TDPREDEFFECTis the coefficient, not the shift. An impulse atu ≤ t0that usedMis already inη(t0)asTDPREDEFFECT, not asT0TDPREDEFFECT. The observed mean of that first-occasion TD carry isτ + λ(μ_t + e^{a Δt} t0_m x0)(Driver et al., 2017, Eq. 5 of the Table 3 / Eq. 3 first-summand TD composition; JSS PDF re-opened 2026-08-20T19:07Z). The evolved mapτ + λ μ_tis not that observed mean. The process-increment mapτ + λ(μ_t + A^{-1}[e^{A Δt} − I] B z)is not that observed mean. The contemporaneous mapτ + λ(μ_t + m x)is not that observed mean. The impulse-carry mapτ + λ(μ_t + e^{a(t−u)} m x)is not that observed mean whenu ≠ t0. The first-occasion TI mapτ + λ(μ_t + e^{a Δt} t0_b z)is not that observed mean. The evolved-plus-carry latent mean is notE(y_t). The lasting level-changeCINTisκ = −a m x(Driver et al., 2017, §7.2, pp. 20–21; JSS PDF re-opened 2026-08-20T19:45Z). Formm xfirst, then multiply by−a. Stablea < 0is required so−κ / a = m xis an equilibrium offset.a ≥ 0cannot hold a new process mean.−a m xis not the dissipating Diracm x, not a freeCINT, and notA^{-1}[e^{A Δt} − I] B z. The extra near-zero-drift latent process also named in §7.2 is a different specification and is not thisCINTsetting. Equation 3 maps that intercept as(1 − e^{a Δt}) m x(JSS PDF re-opened 2026-08-20T19:50Z). Form the level-changeCINTfirst, then the discrete intercept map. Underflow ofe^{a Δt}to+0keepsm x.(1 − e^{a Δt}) m xis notm x, notκ, and notA^{-1}[e^{A Δt} − I] B z. The printed §7.2 lasting level change is an extra near-zero-drift latent process (Driver et al., 2017, §7.2, pp. 22–23; JSS PDF re-opened 2026-08-20T23:10Z).T0MEANS,CINT,T0VAR,DIFFUSION, andTRAITVARof that process are fixed to 0;TDPREDEFFECTon it is fixed to 1; itsDRIFTdiagonal is very close to 0 (printed example−0.000001; precisely 0 causes computational problems); the original process is driven by theDRIFTcouplinga_{ηξ}. After a unit identification impulse the scalar contribution isa_{ηξ} x (e^{ε Δt} − e^{a Δt}) / (ε − a)(ε = aisa_{ηξ} x Δt e^{a Δt}). Forma_{ηξ} xfirst. A zero coupling or zero predictor is exactly zero.ε ≥ 0fails closed. That contribution is notκ = −a m x, not(1 − e^{a Δt}) m x, and not the dissipating Diracm x. The observed mean of that extra-process contribution isτ + λ(μ_t + a_{ηξ} x (e^{ε Δt} − e^{a Δt}) / (ε − a))(Driver et al., 2017, Eq. 5 of that §7.2 contribution; JSS PDF re-opened 2026-08-21T06:12Z). The extra process hasLAMBDA0 and is not an observed indicator. Original indicators load on the original process after theDRIFTcoupling. The evolved mapτ + λ μ_tis not that observed mean. The contemporaneous mapτ + λ(μ_t + m x)is not that observed mean. The contribution is notE(y_t). The evolved-plus-contribution latent mean is notE(y_t).T0TDPREDEFFECTon the extra process begins att = 0and usesΔt = t − t0for both the original-process evolution and the extra drive.TDPREDEFFECTaftert0usest − uwitht0 < u < tfor the extra drive whileμ_tstill usesΔt. The observed mean of that after-t0 extra-process contribution isτ + λ(μ_t + a_{ηξ} x (e^{ε(t−u)} − e^{a(t−u)}) / (ε − a))(Driver et al., 2017, Eq. 5 of that §7.2 after-t0 contribution; JSS PDF re-opened 2026-08-21T06:32Z). The first-occasion extra-process observed mean is not that observed mean whenu ≠ t0. The impulse-carry mapτ + λ(μ_t + e^{a(t−u)} m x)is a Dirac on the original process and is not thatDRIFTdrive. An impulse atu = t0oru = tis not interior. The asymptotic time-independent predictor effect is-B z / a(Driver et al., 2017, §7.2, pp. 20–21; JSS PDF opened 2026-08-21T13:08Z). FormB zfirst, then divide by-a. Stablea < 0is required.a ≥ 0cannot hold a finite process-mean change.-B z / ais not the coefficientB, notA^{-1}[e^{A Δt} − I] B z, notCINT, and notM x. The asymptotic time-independent predictor variance is(B / a)² v(Driver et al., 2017, §7.2, pp. 20–21addedTIPREDVAR). Form the unit asymptotic effect first, then square, then multiply byv.(B / a)² vis notTRAITVAR, notasymDIFFUSION, and not-B z / a. The asymptotic continuous intercept is-κ / a(Driver et al., 2017, Table 2, p. 12asymCINT; Eq. 3 asΔt → ∞; JSS PDF opened 2026-08-21T16:13Z). Formκfirst, then divide by-a. Stablea < 0is required.-κ / ais notκ, notA^{-1}[e^{A Δt} − I] κ, notT0MEANS, and not-B z / a. The p. 16 stationaryT0MEANSconstraint is-κ / a + −B z / a. Form the intercept contribution first, then include the TI extra effect, then add. That constrained first-occasion mean is not freeT0MEANS, notasymCINTalone, notasymTIPREDEFFECTalone, and not the finite-interval discrete latent mean. Equation 5 of that constrained mean isτ + λ(−κ / a + −B z / a)(Driver et al., 2017, §4.3, pp. 9–10; Eq. 5, p. 5; JSS PDF re-opened 2026-08-21T20:07Z). Form the stationary latent mean first, thenτ + λof that mean.τ + λ μ_0for freeT0MEANSis not that composition.τ + λ(−κ / a)is not that composition whenB z ≠ 0.τ + λ μ_tis not that composition.MANIFESTMEANSis notE(y_0). The constrained latent mean is notE(y_0). The p. 16 constrained first-occasion variancetrait + −q / (2 a) + (B / a)² vis not freeT0VAR, notasymDIFFUSIONalone, notTRAITVARalone, notaddedTIPREDVARalone, and not the finite-interval discrete latent variance. Eq. 5 of that constrained variance isλ²(trait + −q / (2 a) + (B / a)² v) + θ + ψ(JSS PDF re-opened 2026-08-22T03:20Z; form the stationary latent variance first, thenλ² p + θ + ψ;λ² p_0is not that observed variance;λ²(−q / (2 a)) + θis not that observed variance whenTRAITVARoraddedTIPREDVARis nonzero;MANIFESTVARis notVar(y_0); the constrained latent variance is notVar(y_0)). The lagged covariance of that constrained process istrait + e^{a Δt}(−q / (2 a)) + (B / a)² v(Driver et al., 2017, Eq. 3–4 of §4.3 / p. 16T0VAR; JSS PDF re-opened 2026-08-22T19:13Z). Trait andaddedTIPREDVARdo not decay withe^{a Δt}. ContemporaneousT0VARis not that lagged map. Decaying the constrained total as if it were all state is not that lagged map. Equation 5 of that lagged covariance isλ²(trait + e^{a Δt}(−q / (2 a)) + (B / a)² v) + ψ.Θdoes not enter. ContemporaneousVar(y_0)is not that lagged observed covariance. The lagged latent covariance is not that observed covariance. The later-occasion variance of that constrained process istrait + e^{2 a Δt}(−q / (2 a)) + Q_Δt + (B / a)² v(Driver et al., 2017, Eq. 3–4 of §4.3 / p. 16T0VAR; JSS PDF re-opened 2026-08-22T23:12Z). Trait andaddedTIPREDVARdo not enterQ_Δt. Under stationarity that composition equals contemporaneousT0VAR. Evolving the constrained total as if it were all state is not that later map. The lagged covariance omitsQ_Δtand is not that later map.Q_Δtis not that later map. Equation 5 of that later-occasion variance isλ²(trait + e^{2 a Δt}(−q / (2 a)) + Q_Δt + (B / a)² v) + θ + ψ. The lagged observed covariance omitsQ_Δtandθ.MANIFESTVARis notVar(y_t). The later-occasion latent variance is notVar(y_t). The later-occasion variance of §4.3 predeterminedT0VARistrait + e^{2 a Δt} p_0 + Q_Δt + (B / a)² v(Driver et al., 2017, Eq. 3–4 of §4.3 predetermined first occasion; JSS PDF re-opened 2026-08-23T05:12Z). Trait andaddedTIPREDVARdo not enterQ_Δt. FreeT0VARp_0is not that later map. Settingp_0 = −q / (2 a)recovers the stationary later-occasion map. Stationary later variance uses−q / (2 a)in place ofp_0and is not that later map whenp_0is free. Evolvingtrait + p_0 + (B / a)² vas if it were all state is not that later map. AsΔt → ∞with stablea < 0the composition approaches contemporaneous stationaryT0VAR. AsΔt → 0+the composition approachestrait + p_0 + (B / a)² v. Nonzero diffusion witha ≥ 0is a growing process and is kept. Equation 5 of that predetermined later-occasion variance isλ²(trait + e^{2 a Δt} p_0 + Q_Δt + (B / a)² v) + θ + ψ.MANIFESTVARis notVar(y_t). The predetermined later-occasion latent variance is notVar(y_t). Stationary later observed variance is not that observed variance whenp_0is free. The lagged covariance of §4.3 predeterminedT0VARistrait + e^{a Δt} p_0 + (B / a)² v(Driver et al., 2017, Eq. 3–4 of §4.3 predetermined first occasion; JSS PDF re-opened 2026-08-23T09:04Z). Trait andaddedTIPREDVARdo not decay withe^{a Δt}. FreeT0VARp_0is not that lagged map. Settingp_0 = −q / (2 a)recovers the stationary lagged map. Stationary lagged covariance uses−q / (2 a)in place ofp_0and is not that lagged map whenp_0is free. Evolvingtrait + p_0 + (B / a)² vas if it were all state is not that lagged map. Later-occasion variance includesQ_Δtand is not that lagged map. AsΔt → ∞with stablea < 0the state term vanishes. AsΔt → 0+the composition approachestrait + p_0 + (B / a)² v. Equation 5 of that predetermined lagged covariance isλ²(trait + e^{a Δt} p_0 + (B / a)² v) + ψ.MANIFESTVARdoes not enter. The predetermined lagged latent covariance is not that observed covariance. Predetermined later observed variance includesQ_Δtandθand is not that lagged observed covariance. Stationary lagged observed covariance is not that observed covariance whenp_0is free. The predetermined first-occasion variance of §4.3 predeterminedT0VARistrait + p_0 + (B / a)² v. Freep_0is not that map. Stationary first-occasion variance uses−q / (2 a)in place ofp_0and is not that map whenp_0is free. Lagged covariance decays the state and is not that map. Later-occasion variance includesQ_Δtand is not that map. Equation 5 of that predetermined first-occasion variance isλ²(trait + p_0 + (B / a)² v) + θ + ψ.MANIFESTVARis not that first-occasion observed variance. The predetermined first-occasion latent variance is not that observed variance. Stationary first-occasion observed variance is not that observed variance whenp_0is free. Predetermined later observed variance includesQ_Δtand is not that first-occasion observed variance. Later-start lagged covariance of predeterminedT0VARistrait + e^{a s}(e^{2 a u} p_0 + Q_u) + (B / a)² v(Driver et al., 2017, §4.3startoffset; Eq. 4; JSS PDF re-opened 2026-08-23T10:27Z). First-occasion lagged omitse^{a s} Q_u. Later-occasion variance does not lag. Stationary lagged uses−q / (2 a). Decaying the later total is not that map. Equation 5 of that later-start lagged covariance isλ²of it plusψ. Independentε_tdoes not enter. First-occasion lagged observed omitse^{a s} Q_u. Predetermined later observed variance includesQ_uandθand is not that later-start lagged observed covariance. Later-start later-occasion variance of predeterminedT0VARistrait + e^{2 a s}(e^{2 a u} p_0 + Q_u) + Q_s + (B / a)² v(Driver et al., 2017, §4.3startoffset; Eq. 3–4 Chapman–KolmogorovQ_{u+s} = e^{2 a s} Q_u + Q_s; JSS PDF re-opened 2026-08-23T11:05Z). Later-occasion variance atuomitsQ_s. Later-start lagged covariance omitsQ_s. Stationary later uses−q / (2 a). Evolving the later total as if it were all state is not that map. Ignoringstartoffsetomitse^{2 a s} Q_u. Equation 5 of that later-start later-occasion variance isλ²of it plusθ + ψ.MANIFESTVARis not that observed variance. Page 16discreteDRIFTstdise^{a Δt}after strictly positiveasymDIFFUSION-q / (2 a)(Driver et al., 2017, p. 16; footnote 4; §7.1; JSS PDF re-opened 2026-08-23T11:40Z). Unstandardisede^{a Δt}is defined for growinga ≥ 0and for zero diffusion and is notdiscreteDRIFTstd. The §7.1 trait-plus-state autocorrelation(trait + e^{a Δt} p + added) / (trait + p + added)usesTRAITVARand is notdiscreteDRIFTstd.TRAITVARis not the standardisation variance. Page 16discreteDIFFUSIONstdisQ_Δt / (−q / (2 a))after strictly positiveasymDIFFUSION-q / (2 a)(Driver et al., 2017, p. 16; Eq. 4; footnote 4; JSS PDF re-opened 2026-08-23T13:06Z). UnstandardisedQ_Δtis defined for growinga ≥ 0and for zero diffusion and is notdiscreteDIFFUSIONstd. The continuous standardisation−2 ais notdiscreteDIFFUSIONstd.Q_Δt / (trait + p + added)usesTRAITVARand is notdiscreteDIFFUSIONstd.TRAITVARis not the standardisation variance. Page 16DIFFUSIONstdisq / (−q / (2 a)) = −2 aafter strictly positiveasymDIFFUSION-q / (2 a)(Driver et al., 2017, p. 16; Eq. 4; footnote 4; JSS PDF re-opened 2026-08-23T13:20Z). Unstandardisedqis defined for growinga ≥ 0and for zero diffusion and is notDIFFUSIONstd. The discrete standardisationQ_Δt / (−q / (2 a))depends onΔtand is notDIFFUSIONstd.q / (trait + p + added)usesTRAITVARand is notDIFFUSIONstd.TRAITVARis not the standardisation variance. Page 16DRIFTstdis the continuous auto-effect after strictly positiveasymDIFFUSION-q / (2 a)(Driver et al., 2017, p. 16; Eq. 1; footnote 4; JSS PDF re-opened 2026-08-23T13:28Z). Unstandardisedais defined for growinga ≥ 0and for zero diffusion and is notDRIFTstd. The discrete standardisatione^{a Δt}depends on the event interval and is notDRIFTstd.a p / (trait + p + added)usesTRAITVARand is notDRIFTstd.TRAITVARis not the standardisation variance. Page 16asymTIPREDEFFECTstdis(-B / a) · √v / √(-q / (2 a))after strictly positiveasymDIFFUSION-q / (2 a)and strictly positive predictor variancev(Driver et al., 2017, p. 16; §7.2; footnote 4; JSS PDF re-opened 2026-08-23T14:25Z). Unstandardised-B / ais defined for a zero coefficient and for zero predictor variance and is notasymTIPREDEFFECTstd. The finite-interval standardisationA^{-1}[e^{A Δt} − I] B · √v / √pdepends on the event interval and is notasymTIPREDEFFECTstd.(-B / a) · √v / √(trait + p + added)usesTRAITVARand is notasymTIPREDEFFECTstd.TRAITVARis not the standardisation variance. Page 16TIPREDEFFECTstdisB · √v / √(-q / (2 a))after strictly positiveasymDIFFUSION-q / (2 a)and strictly positive predictor variancev(Driver et al., 2017, p. 16; §7.2; footnote 4; JSS PDF re-opened 2026-08-23T16:21Z). UnstandardisedBis defined for a zero coefficient and for zero predictor variance and is notTIPREDEFFECTstd. The asymptotic standardisation(-B / a) · √v / √pis the total change and is notTIPREDEFFECTstd. The finite-interval standardisationA^{-1}[e^{A Δt} − I] B · √v / √pdepends on the event interval and is notTIPREDEFFECTstd.B · √v / √(trait + p + added)usesTRAITVARand is notTIPREDEFFECTstd.TRAITVARis not the standardisation variance. Page 16 / Table 3T0TIPREDEFFECTstdist0_b · √v / √p_0after strictly positive freeT0VARp_0and strictly positive predictor variancev(Driver et al., 2017, Table 3, p. 13; p. 16; footnote 4; 2017-era ctsemsummary.ctsemFit.R; JSS PDF re-opened 2026-08-23T17:20Z). The affected variance is free first-occasionT0VAR, notasymDIFFUSION. Unstandardisedt0_bis defined for a zero coefficient and for zero predictor variance and is notT0TIPREDEFFECTstd.TIPREDEFFECTstdB · √v / √(-q / (2 a))is the continuous coefficient and is notT0TIPREDEFFECTstd.asymTIPREDEFFECTstd(-B / a) · √v / √pis the total change and is notT0TIPREDEFFECTstd.t0_b · √v / √(trait + p_0 + added)usesTRAITVARand is notT0TIPREDEFFECTstd.TRAITVARis not the standardisation variance. 2017-eraaddedT0TIPREDVARist0_b² vafter a first-occasion time-independent predictor (Driver et al., 2017, Table 3, p. 13; p. 16; §7.2; 2017-era ctsemsummary.ctsemFit.R; JSS PDF re-opened 2026-08-23T18:20Z). Formt0_bfirst, then square, then multiply byv. A zero coefficient or zero predictor variance is exactly zero. FreeT0TIPREDEFFECTdoes not requirea < 0.(B / a)² visaddedTIPREDVARand is not this first-occasion map.t0_b · √v / √p_0isT0TIPREDEFFECTstdand is not this variance. FreeT0VARis not this extra TI variance.TRAITVARis not this extra TI variance. Equation 5 of 2017-eraaddedT0TIPREDVARisλ² t0_b² v(Driver et al., 2017, Eq. 5, p. 5; Table 3, p. 13; Table 2, p. 12; 2017-era ctsemsummary.ctsemFit.R; JSS PDF re-opened 2026-08-23T19:10Z). Formt0_b² vfirst, then(λ extra) λwithθ = 0. A zero loading or zero extra is exactly zero.t0_b² vis the latent extra, not the observed extra.λ² p_0 + θis first-occasion observed variance, not this extra.λ² (B / a)² vis Eq. 5 ofaddedTIPREDVAR, not this first-occasion observed extra.MANIFESTVARθis not this extra. Equation 5 of §7.2addedTIPREDVARisλ² (B / a)² v(Driver et al., 2017, Eq. 5, p. 5; Table 2, p. 12; §7.2, pp. 20–21; 2017-era ctsemsummary.ctsemFit.R; JSS PDF re-opened 2026-08-23T19:23Z). Form(B / a)² vfirst, then(λ extra) λwithθ = 0. A zero loading or zero extra is exactly zero. Lasting asymptotic extra requiresa < 0.(B / a)² vis the latent extra, not the observed extra.λ² t0_b² vis first-occasion extra observed TI variance, not this extra.λ² p + θis stationary observed variance, not this extra.MANIFESTVARθis not this extra. Page 16TDPREDEFFECTstdism · √v / √(-q / (2 a))after strictly positiveasymDIFFUSIONand strictly positive time-dependent predictor variance. UnstandardisedMis notTDPREDEFFECTstd.TIPREDEFFECTstdis notTDPREDEFFECTstdeven whenM = B. intercept-styleA^{-1}[e^{A Δt} − I] M · √v / √pis notTDPREDEFFECTstd.m · √v / √(trait + p + added)usesTRAITVARand is notTDPREDEFFECTstd. Table 3 / p. 16T0TDPREDEFFECTstdist0_m · √v / √p_0after strictly positive freeT0VARand strictly positive TD predictor variance. Unstandardisedt0_mis notT0TDPREDEFFECTstd.TDPREDEFFECTstdusesasymDIFFUSIONand is notT0TDPREDEFFECTstd.T0TIPREDEFFECTstdis notT0TDPREDEFFECTstdeven whent0_m = t0_b.t0_m · √v / √(trait + p_0 + added)usesTRAITVARand is notT0TDPREDEFFECTstd. FreeT0VARdoes not requirea < 0. Page 16T0VARstdisp_0 / p_0 = 1after strictly positive freeT0VAR(solve(sqrt(diag(T0VAR))) %&% T0VAR; OpenMx%&%ist(A) %*% B %*% A; the default ridge is 0). UnstandardisedT0VARis notT0VARstd.T0TDPREDEFFECTstdis notT0VARstd.addedT0TIPREDVARis notT0VARstd. Page 16TRAITVARstdistrait / trait = 1after strictly positiveTRAITVAR(solve(sqrt(diag(TRAITVAR))) %&% TRAITVAR; OpenMx%&%ist(A) %*% B %*% A; no ridge addend). UnstandardisedTRAITVARis notTRAITVARstd.T0VARstdis notTRAITVARstdeven when both equal 1.addedT0TIPREDVARis notTRAITVARstd. Page 16MANIFESTTRAITVARstdisψ / ψ = 1after strictly positiveMANIFESTTRAITVAR(solve(sqrt(diag(MANIFESTTRAITVAR))) %&% MANIFESTTRAITVAR; OpenMx%&%ist(A) %*% B %*% A; 2017-era source adds ridging; default ridge is 0). UnstandardisedMANIFESTTRAITVARis notMANIFESTTRAITVARstd.TRAITVARstdis notMANIFESTTRAITVARstdeven when both equal 1.MANIFESTVARis notMANIFESTTRAITVARstd. Page 16MANIFESTVARstdisθ / θ = 1after strictly positiveMANIFESTVAR(solve(sqrt(diag(MANIFESTVAR))) %&% MANIFESTVAR; OpenMx%&%ist(A) %*% B %*% A; 2017-era source adds ridging; default ridge is 0; 2017-eradimnamesassignment tolatentNamesis a source bug). UnstandardisedMANIFESTVARis notMANIFESTVARstd.MANIFESTTRAITVARstdis notMANIFESTVARstdeven when both equal 1. Equation 5Var(y)is notMANIFESTVARstd. Page 16TIPREDVARstdisv / v = 1after strictly positiveTIPREDVAR(solve(sqrt(diag(TIPREDVAR))) %&% TIPREDVAR; OpenMx%&%ist(A) %*% B %*% A; 2017-era source adds ridging; default ridge is 0;dimnamesareTIpredNames). UnstandardisedTIPREDVARis notTIPREDVARstd.MANIFESTVARstdis notTIPREDVARstdeven when both equal 1. Section 7.2addedTIPREDVARis notTIPREDVARstd. Page 16asymDIFFUSIONstdisp / p = 1after strictly positiveasymDIFFUSION(solve(sqrt(diag(asymDIFFUSION))) %&% asymDIFFUSION; OpenMx%&%ist(A) %*% B %*% A; 2017-era source adds ridging; default ridge is 0;dimnamesarelatentNames). UnstandardisedasymDIFFUSIONis notasymDIFFUSIONstd.TIPREDVARstdis notasymDIFFUSIONstdeven when both equal 1.DIFFUSIONstd−2 ais notasymDIFFUSIONstd. Page 16discreteCINTstdisA^{-1}[e^{A Δt} − I] κ / √pafter strictly positiveasymDIFFUSION. UnstandardiseddiscreteCINTis notdiscreteCINTstd.κ / √pis notdiscreteCINTstd.(-κ / a) / √pis notdiscreteCINTstd.asymCINTstdis(-κ / a) / √pafter strictly positiveasymDIFFUSION. UnstandardisedasymCINTis notasymCINTstd.κ / √pis notasymCINTstd.discreteCINTstdis notasymCINTstd.T0MEANSstdisμ_0 / √p_0after strictly positive freeT0VAR. UnstandardisedT0MEANSis notT0MEANSstd.T0VARstdis notT0MEANSstd.μ_0 / √asymDIFFUSIONis notT0MEANSstd. Page 16MANIFESTMEANSstdisτ / √θafter strictly positiveMANIFESTVAR. UnstandardisedMANIFESTMEANSis notMANIFESTMEANSstd.MANIFESTVARstdis notMANIFESTMEANSstd.τ / √(λ² Var(η) + θ)is notMANIFESTMEANSstd. Page 16CINTstdisκ / √pafter strictly positiveasymDIFFUSION. UnstandardisedCINTis notCINTstd.asymCINTstdis notCINTstd.discreteCINTstdis notCINTstd.κ / √(trait + p + added)is notCINTstd. Evolving from that stationary start withCINTandTIPREDEFFECTstays at the stationary mean. Equation 1 is the latent SDE, not the measurement model. Form(λ p) λthen addθ, then addψ.MANIFESTVARisΘ, notVar(y).MANIFESTTRAITVARisΨ_τ, notΘ.TRAITVARis latent and scaled byλ².Var(η)is notVar(y). - Separate cluster means before within-unit lag. CWC plus an event-time lag is not DSEM. Subtracting the person-specific mean from a raw autoregressive series does not isolate the lagged within-person effect (Curran & Bauer, 2011, pp. 607–608); already-centered residuals with irregular event intervals use the exact scalar map.
- Do not treat the CWC cluster-mean coefficient as the between-cluster effect. It is the contextual effect
between − within(Enders & Tofighi, 2007, Table 2, pp. 124–127). - Never use future-available evidence in historical model fits.
- Do not blanket-mask PII when identity/role/linkage is scientifically required. Follow the purpose-bound separation, opaque-ID, encryption, retention, and audit contract in
docs/PRIVACY_DATA_GOVERNANCE.md. - Treat documents and LLM outputs as untrusted. Model routing/orchestration may vary reasoning effort, decomposition, recursion and roles, but deterministic/statistical gates remain authoritative.
- Treat CSAP/SOC 2/ISO/NIST mappings as readiness evidence, never certification or attestation.
Before stating that a task is complete, run the exact focused tests, complete test suite, line/branch coverage gate, docstring gate, formatter, linter, dependency/security checks, build/package checks, documentation contracts, and any required CPU/GPU parity or true-parameter study. Report actual evidence and unresolved external gates. If another safe executable repository action remains, continue rather than using the verification result as a stopping point.