Association or variance cancellation, not evidence of causal buffering. Variance_ratio is 1 - var(rowSums(X))/sum(var(X_j)); positive values indicate cancellation. mean_neg_cor is the negative mean correlation.
Usage
rri_compensation_index(
res,
per_group = FALSE,
group_cols = NULL,
id = NULL,
method = c("mean_neg_cor", "variance_ratio"),
scale_output = TRUE
)Arguments
- res
RRI result.
- per_group
Compute by trajectory/group.
- group_cols
Group identifiers.
- id
Optional aligned identifiers.
- method
mean_neg_cor or variance_ratio.
- scale_output
For correlations maps
[-1, 1]to[0, 1]; for variance cancellation truncates negative values to zero. A correlation score of 0.5 means zero mean correlation, not moderate biological compensation.
Examples
# \donttest{
sim <- simulate_redox_holobiont(seed = 1)
res <- rri_pipeline(soil = sim$Eh_stability, plant = sim$ROS_flux)
#> Warning: Unanchored latent axes have arbitrary signs; RRI is exploratory, not directionally validated resilience.
#> Warning: Excluding simulator-derived hidden columns from scoring: Cacc_EAC, Cacc_EDC, Cacc_total, Cacc_fraction, net_oxidative_balance, alpha_accept, alpha_donate, k_accept_h, k_donate_h
rri_compensation_index(res)
#> $compensation_index
#> [1] 0.2863531
#>
#> $pairwise_cor
#> Physio_Soil
#> 0.4272939
#>
#> $method_used
#> [1] "mean_neg_cor"
#>
#> $interpretation
#> [1] "Descriptive asynchrony only; does not establish functional compensation."
#>
# }