h_jags_get_data: Getting Data for 'JAGS'

In Roche/crmPack: Object-Oriented Implementation of CRM Designs

Getting Data for JAGS

Description

A simple helper function that prepares an object for data argument of rjags::jags.model(), which is invoked by mcmc() method.

Usage

h_jags_get_data(model, data, from_prior)

Arguments

model	(GeneralModel) an input model.
data	(GeneralData) an input data.
from_prior	(flag) sample from the prior only? In this case data will not be appended to the output, i.e. only the variables required by the model@priormodel model will be returned in data.

Examples

# Create some data from the class `Data`.
my_data <- Data(
  x = c(0.1, 0.5, 1.5, 3, 6, 10, 10, 10),
  y = c(0, 0, 0, 0, 0, 0, 1, 0),
  doseGrid = c(0.1, 0.5, 1.5, 3, 6, seq(from = 10, to = 80, by = 2))
)

# Initialize the CRM model.
my_model <- LogisticLogNormal(
  mean = c(-0.85, 1),
  cov = matrix(c(1, -0.5, -0.5, 1), nrow = 2),
  ref_dose = 56
)

jags_data <- h_jags_get_data(my_model, my_data, from_prior = FALSE)
jags_data

Roche/crmPack documentation built on April 30, 2024, 3:19 p.m.