Skip to contents

Draws nsim replicates of Y from the fitted ALS mean structure: mu, beta, a, b, U, V, and the family-appropriate noise distribution. For directed unipartite normal and binary fits, reciprocal dyads use the fitted residual rho. The returned object has class "ame.sim" so it is compatible with plot_ppc_*-style consumers (where applicable).

Usage

# S3 method for class 'ame_als'
simulate(object, nsim = 1, seed = NULL, ...)

Arguments

object

an ame_als fit.

nsim

integer; number of replicates to simulate (default 1).

seed

optional RNG seed.

...

ignored.

Value

An object of class "ame.sim" with element Y – a list of nsim simulated outcome arrays (one matrix for a cross-section, one list-of-matrices per slice for a longitudinal fit).

Details

ALS is a point estimator. simulate.ame_als holds mu, beta, a, b, U, V at the point estimate and resamples the noise. For uncertainty over the parameters themselves, use ame_als_bootstrap (whose replicates each carry their own resampled Y) and combine those simulations.