Description Usage Arguments Details Value Examples
Bayesian Markov Chain Monte Carlo model building for household-community tuberculosis transmission models
1 2 3 4 |
design.matrix |
A data frame type design matrix of observations for rows, variables on each observation for columns. If there are natural clusters (households), then a variable titled 'cluster' must be included. The design matrix must have exactly one column titled 'y' for outcome on each observation. |
categorical.columns |
A scalar (vector) specification of any column(s) of type "factor" to be coded as a dummy variable regression. Default is 0: no columns are factor variables. |
prior.beta.means |
A scalar (vector) specification of Gaussian prior distirbution mean(s) for the coefficient(s) in the household component of the logistic regression model. Default is noninformative on the inverse logit transformed (odds ratio) scale. |
prior.beta.sd |
A scalar (vector) specification of Gaussian prior distirbution standard deviation(s) for the coefficient(s) in the household component of the logistic regression model. Default is noninformative on the inverse logit transformed (odds ratio) scale. |
prior.alpha |
A vector of prior probability parameters for the alpha term in the community component of model, a logistic regression with intercept only. Default is noninformative, with p^C being essentially uniform on the unit interval. |
The package contains long-form documentation in the form of a vignette that cover the use of the main fucntions. Use browseVignettes(package="upmfit") to access them.
A list of: a JAGS model script for specification of the unified probability model to be used; the (possibly recoded) design matrix; a vector of the beta coefficient names to be used in the household component of the UPM.
1 2 3 4 5 | upmbuilder(design.matrix=upmdata)
upmbuilder(design.matrix=upmdata, prior.alpha=c(-2.2,round(1/sqrt(3),3)))
suppressWarnings(cat(upmbuilder(design.matrix=upmdata,
prior.alpha=c(-2.2,round(1/sqrt(3),3)))[[1]]))
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.