Skip to contents

Add one batch or node statistics to a running accumulator

Source: R/sufficient_stats.R
accumulate_stats.Rd

Performs the additive aggregation acc <- acc + stats component by component. After processing all batches, the accumulator holds the full-data sufficient summary.

Usage

accumulate_stats(acc, stats)

Arguments

acc

A vcmm_accumulator object from init_accumulator().

stats

A vcmm_ss object from compute_sufficient_stats().

Value

The updated accumulator (class "vcmm_accumulator").

See also

Other sufficient statistics: compute_sufficient_stats(), init_accumulator()

Examples

set.seed(1)
n_batch <- 50; p <- 3; q <- 2
acc <- init_accumulator(p, q)

for (b in 1:3) {
  X <- cbind(1, matrix(rnorm(n_batch * (p - 1)), n_batch, p - 1))
  Z <- matrix(rnorm(n_batch * q), n_batch, q)
  y <- rnorm(n_batch)
  ss <- compute_sufficient_stats(y, X, Z)
  acc <- accumulate_stats(acc, ss)
}

acc$n_obs  # 150
#> [1] 150