clinDR: fitEmaxB – R documentation

Become an expert in R — Interactive courses, Cheat Sheets, certificates and more!

fitEmaxB

Bayesian fit of hyperbolic or sigmoidal Emax models to continuous/binary dose response data.

Description

Uses Rpackage rstan to fit a Bayesian hyperbolic or sigmoidal Emax model. Different intercepts for multiple protocol-data are supported. For binary data, the Emax model is on the logit scale.

Usage

fitEmaxB(y, dose, prior, modType = 4, prot = rep(1, length(y)), 
	count = rep(1, length(y)), xbase=NULL,
	binary = FALSE, msSat = NULL,
	pboAdj = FALSE, mcmc = mcmc.control(), estan = NULL, 
	diagnostics = FALSE, nproc = getOption("mc.cores", 1L))

Arguments

`y`	Outcome for each patient. Missing `Y` values are are not permitted. Dose/protocol group means for grouped continuous data. For binary data, `y` must be 0/1 and counts must be supplied for each 0/1 value.
`dose`	Dose for each patient.
`prior`	Prior specification through an object of type 'emaxPrior' or 'prior'. See `emaxPrior.control` and `prior.control` for details. The 'emaxPrior' specifies the magnitude of the potential effect for a specified dose (typically the highest anticipated dose and/or the dose in a POC study), while the 'prior' specifies the theoretical maximum effect (the emax parameter). The 'prior' specification is deprecated and will be removed.
`modType`	modType=3 (default) for the 3-parameter hyperbolic Emax model. modType=4 for the 4-parameter sigmoidal Emax model.
`prot`	Protocol (group) membership used to create multiple intercepts. The default is a single protocol. The prior disribution for the placebo response is re-used independently for each intercept.
`count`	Counts for the number of patients when the `Y` are dose continuous group means or binary 0/1 values. Default is 1 (ungrouped data).
`xbase`	A matrix of baseline covariates with rows corresponding to `y` that enter as linear additive predictors. The baseline covariates must be centered about their (protocol-specific) means. `xbase` does not include an intercept or protocol indicators. Covariates cannnot be specified with PBO adjusted or aggregated input.
`binary`	When `TRUE`, the `y` are assumed to be coded `0/1`, and the means reported are proportions. The Emax model is specified on the logit scale, and proportions are estimated from the model by back-transformation.
`msSat`	If continuous `Y` are dose/protocol group means rather than individual measurements, the within group variance, `msSat`, should be supplied. This variance is the mean square from the model saturated in dose and protocol. It is used to improve the residual variance estimate for the Emax model. If it is not supplied, the residual SD (and associated SE) will have low degrees of freedom.
`pboAdj`	For published data with only pbo-adjusted dose group means and SEs, the model is fit without an intercept(s). If initial parameters are supplied, the intercept (E0) should be assigned `0`. A zero for the placebo mean should not be included in `Y`. This option is not available for binary data. Potential correlation between between placebo-adjusted means is ignored.
`mcmc`	Inputs controlling `rstan` execution. See `mcmc.control` for details.
`estan`	The compiled `rstan` Emax model is usually loaded automatically. It can be load to an object using the function `selEstan` and passed to `fitEmaxB` for repeated executions to improve efficiency and stability.
`diagnostics`	Printed output from rstan. See `Details` for more information.
`nproc`	The number of processor requested for `STAN` MCMC computations. Defaults to the value set by the `rstan` installation. When set explicitly, `nproc` is usually 1 or the number of MCMC chains. If greater than the number of chains, it is set to the number of chains.

Details

The function compileStanModels must be executed once to create compiled STAN code before fitEmaxB can be used.

MCMC fit of a Bayesian hyperbolic or sigmoidal Emax model. The prior distributions available are based on the publication Thomas, Sweeney, and Somayaji (2014) and Thomas and Roy (2016).

The posterior distributions are complex because the distributions of the Emax and ED50 parameters change substantially as a function of the lambda, often creating 'funnel' type conditions. Small numbers of divergences are common and do not appear easily avoided. Extensive simulation using evaluations with emaxsimB support the utility of the resulting approximate posterior distributions. The number of divergences can be viewed using diagnostics=TRUE. The usual convergence diagnostics should always be checked.

Value

A list assigned class "fitEmaxB" with:

`estanfit`	The `rstan` object with the model fit.
`y, dose, prot, count, modType, binary, pboAdj, nbase, msSat, prior, mcmc`	Input values.

Note

The default modType was changed from 3 to 4 for clinDR version >2.0

Author(s)

Neal Thomas

References

Thomas, N., Sweeney, K., and Somayaji, V. (2014). Meta-analysis of clinical dose response in a large drug development portfolio, Statistics in Biopharmaceutical Research, Vol. 6, No.4, 302-317. <doi:10.1080/19466315.2014.924876>

Thomas, N., and Roy, D. (2016). Analysis of clinical dose-response in small-molecule drug development: 2009-2014. Statistics in Biopharmaceutical Research, Vol. 6, No.4, 302-317 <doi:10.1080/19466315.2016.1256229>

Examples

## Not run: 

data("metaData")
exdat<-metaData[metaData$taid==1,]

prior<-emaxPrior.control(epmu=0,epsca=4,difTargetmu=0,difTargetsca=4,dTarget=20,
				p50=(2+5)/2,
				sigmalow=0.01,sigmaup=3)
										
mcmc<-mcmc.control(chains=3)

msSat<-sum((exdat$sampsize-1)*(exdat$sd)^2)/(sum(exdat$sampsize)-length(exdat$sampsize))
fitout<-fitEmaxB(exdat$rslt,exdat$dose,prior,modType=4,prot=exdat$protid,
				count=exdat$sampsize,msSat=msSat,mcmc=mcmc)
plot(fitout)

## End(Not run)

clinDR

Simulation and Analysis Tools for Clinical Dose Response Modeling

v2.3.5

GPL (>= 2)

Authors

Neal Thomas [aut, cre] (<https://orcid.org/0000-0002-1915-8487>), Jing Wu [aut], Mike K. Smith [aut]

Initial release

2021-04-11