Decomposes the marginal trait variance into contributions from each
active component of the model. Column names follow the fitted report
blocks: rr_B, diag_B, rr_W, diag_W, diag_species, phylo_rr,
phylo_diag, propto, and residual (Gaussian/lognormal observation
variance where present). Non-Gaussian link-implicit residual shares are
handled by extract_proportions(), not by VP().
Arguments
- fit
A fitted multivariate model returned by
gllvmTMB().
Value
An n_traits × n_components matrix of variance shares (rows
sum to 1). Columns are only those active in fit$use.
Details
Mirrors gllvm::VP() / gllvm::plotVP().
Structured trait-intercept rho fits are excluded; use extract_Sigma() and
its source_strength metadata.
See also
extract_proportions() for the canonical per-trait
variance-share decomposition (shared_unit, unique_unit_obs, phy
tiers, link-residual, etc.).