View source: R/api-prediction.R
| predict_mfrm_population | R Documentation |
Forecast population-level MFRM operating characteristics for one future design
predict_mfrm_population(
fit = NULL,
sim_spec = NULL,
n_person = NULL,
n_rater = NULL,
n_criterion = NULL,
raters_per_person = NULL,
reps = 50,
fit_method = NULL,
model = NULL,
maxit = 25,
quad_points = 7,
residual_pca = c("none", "overall", "facet", "both"),
seed = NULL
)
fit |
Optional output from |
sim_spec |
Optional output from |
n_person |
Number of persons/respondents in the future design. Defaults to the value stored in the base simulation specification. |
n_rater |
Number of rater facet levels in the future design. Defaults to the value stored in the base simulation specification. |
n_criterion |
Number of criterion/item facet levels in the future design. Defaults to the value stored in the base simulation specification. |
raters_per_person |
Number of raters assigned to each person in the future design. Defaults to the value stored in the base simulation specification. |
reps |
Number of replications used in the forecast simulation. |
fit_method |
Estimation method used inside the forecast simulation. When
|
model |
Measurement model used when refitting the forecasted design. Defaults to the model recorded in the base simulation specification. |
maxit |
Maximum iterations passed to |
quad_points |
Quadrature points for |
residual_pca |
Residual PCA mode passed to |
seed |
Optional seed for reproducible replications. |
predict_mfrm_population() is a scenario-level forecasting helper built
on top of evaluate_mfrm_design(). It is intended for questions such as:
what separation/reliability would we expect if the next administration had 60 persons, 4 raters, and 2 ratings per person?
how much Monte Carlo uncertainty remains around those expected summaries?
The function deliberately returns aggregate operating characteristics (for example mean separation, reliability, recovery RMSE, convergence rate) rather than future individual true values for one respondent or one rater.
If fit is supplied, the function first constructs a fit-derived parametric
starting point with extract_mfrm_sim_spec() and then evaluates the
requested future design under that explicit data-generating mechanism. This
should be interpreted as a fit-based forecast under modeling assumptions, not
as a guaranteed out-of-sample prediction.
An object of class mfrm_population_prediction with components:
design: requested future design
forecast: facet-level forecast table
overview: run-level overview
simulation: underlying evaluate_mfrm_design() result
sim_spec: simulation specification used for the forecast
settings: forecasting settings
ademp: simulation-study metadata
notes: interpretation notes
forecast contains facet-level expected summaries for the requested
future design.
Mcse* columns quantify Monte Carlo uncertainty from using a finite number
of replications.
simulation stores the full design-evaluation object in case you want to
inspect replicate-level behavior.
This helper does not produce definitive future person measures or rater severities for one concrete sample. It forecasts design-level behavior under the supplied or derived parametric assumptions.
The forecast is implemented as a one-scenario Monte Carlo / operating-
characteristic study following the general guidance of Morris, White, and
Crowther (2019) and the ADEMP-oriented reporting framework discussed by
Siepe et al. (2024). In mfrmr, this function is a practical wrapper for
future-design planning rather than a direct implementation of a published
many-facet forecasting procedure.
Morris, T. P., White, I. R., & Crowther, M. J. (2019). Using simulation studies to evaluate statistical methods. Statistics in Medicine, 38(11), 2074-2102.
Siepe, B. S., Bartos, F., Morris, T. P., Boulesteix, A.-L., Heck, D. W., & Pawel, S. (2024). Simulation studies for methodological research in psychology: A standardized template for planning, preregistration, and reporting. Psychological Methods.
build_mfrm_sim_spec(), extract_mfrm_sim_spec(),
evaluate_mfrm_design(), summary.mfrm_population_prediction
spec <- build_mfrm_sim_spec(
n_person = 40,
n_rater = 4,
n_criterion = 4,
raters_per_person = 2,
assignment = "rotating"
)
pred <- predict_mfrm_population(
sim_spec = spec,
n_person = 60,
reps = 2,
maxit = 10,
seed = 123
)
s_pred <- summary(pred)
s_pred$forecast[, c("Facet", "MeanSeparation", "McseSeparation")]
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.