Package 'Romeb'

Bayesian Growth Curve Model for Complete Data, MCAR (Missing Completely at Random, and MAR (Missing at Random

Description

A character string containing the JAGS model specification for complete data (no missing values), MCAR, and MAR.

Usage

model

Format

A character string.

---

Bayesian Growth Curve Model for MNAR (Missing Not At Random)

Description

JAGS model definition for data with MNAR mechanism.

Usage

model_MNAR

Format

A character string.

---

Bayesian Growth-Curve Model for MNAR with Auxiliary Variable (k)

Description

JAGS model definition for MNAR mechanism with auxiliary variable k.

Usage

model_MNAR_k

Format

A character string.

---

Youth Attitudes toward Deviance (NYS, 1976–1980)

Description

A real data set from the 1976–1980 National Youth Survey of U.S. youth.

Usage

NYS

Format

A data frame with 1,725 rows and 7 variables:

age

Participant age (years)

gender

Gender (0 = female, 1 = male)

Atd1

Attitude toward social deviance, wave 1

Atd2

Attitude toward social deviance, wave 2

Atd3

Attitude toward social deviance, wave 3

Atd4

Attitude toward social deviance, wave 4

Atd5

Attitude toward social deviance, wave 5

Source

National Youth Survey, waves 1976–1980 (downloadable at https://www.icpsr.umich.edu/icpsrweb/ICPSR/series/88)

---

Romeb: Robust Median-Based Bayesian Growth Curve Modeling

Description

Implements robust median-based Bayesian growth-curve models that handle MCAR/MAR/MNAR missing-data mechanisms and complete data. Models are fitted via rjags/JAGS and summarized with coda.

Fits a median-based Bayesian growth curve model under MCAR, MAR, MNAR or complete-data assumptions. If K > 0 the first K columns in data are treated as auxiliary variables.

Usage

Romeb(
  Missing_Type,
  data,
  time,
  seed,
  K = 0,
  chain = 1,
  Niter = 6000,
  burnIn = 3000
)

Arguments

Missing_Type	Character; one of MNAR, MAR, MCAR, no missing.
data	Matrix or data frame containing outcome columns (and optionally auxiliary variables).
time	Numeric vector of measurement times (e.g., c(0,1,2,3)).
seed	Integer seed for reproducibility.
K	Integer; number of auxiliary variables (default 0).
chain	Integer; number of MCMC chains (default 1).
Niter	Integer; iterations per chain (default 6000).
burnIn	Integer; burn-in iterations (default 3000).

Value

An object of class RomebResult containing

quantiles

posterior means, SDs and quantiles

geweke

Geweke z-scores

credible_intervals

95% equal-tail credible intervals

hpd_intervals

95% highest posterior density intervals

samps_full

full coda::mcmc.list (including burn-in)

Author(s)

Maintainer: Dandan Tang [email protected] (ORCID)

Authors:

• Xin Tong

See Also

Useful links:

• https://github.com/DandanTang0/Romeb

• Report bugs at https://github.com/DandanTang0/Romeb/issues

Examples

set.seed(123)
  Y <- matrix(rnorm(300), 100, 3)
  fit <- Romeb("no missing", data = Y, time = c(0,1,2), seed = 123, K = 0,
               Niter = 6000, burnIn = 3000)
  print(fit)

---