Description Usage Arguments Details Value Author(s) References See Also Examples
This function estimates prevalence ratios (PRs)
and bootstrap confidence intervals using logistic models for marginal standardization.
The estimation of standard errors for PRs is obtained through use of bootstrapping.
Confidence intervals of (1alpha)% for PRs are available for standard logistic regression
and for randomeffects logistic models (Santos et al, 2008). The function
prLogisticBootMarg
allows estimation of PRs using marginal standardization procedure
(Wilcosky and Chambless, 1985).
1  prLogisticBootMarg(object, data, conf = 0.95, R = 99, ...)

object 
any fitted model object from which fixed effects estimates can be extracted. The details of model specification are given below. 
data 
a required data frame containing the variables named in 
conf 
scalar or vector specifying confidence level(s) for estimation. The default is

R 
the number of bootstrap replicates. The default is 
... 
optional additional arguments. Currently none are used in any methods. 
The fitted model object can be obtained using glm()
function for binary responses
when unit samples are independent. The glmer()
function should be used
for correlated binary responses. Only binary predictors are allowed. If categorization for predictors
is other than (0,1), factor()
should be considered.
Returns prevalence ratio and its 95% bootstrap confidence intervals for marginal standardization. Both normal and percentile bootstrap confidence intervals are presented.
Raydonal Ospina, Department of Statistics, Federal University of Pernambuco, Brazil
([email protected])
Leila D. Amorim, Department of Statistics, Federal University of Bahia, Brazil
([email protected]).
Localio AR, Margolis DJ, Berlin JA (2007). Relative risks and confidence intervals were easily computed indirectly from multivariate logistic regression. Journal of Clinical Epidemiology, 60, 874882.
Oliveira NF, Santana VS, Lopes AA (1997). Ratio of proportions and the use of the delta method for confidence interval estimation in logistic regression. Journal of Public Health, 31(1), 9099.
Santos CAST et al (2008).
Estimating adjusted prevalence ratio in clustered crosssectional epidemiological data.
BMC Medical Research Methodology, 8 (80). Available from
http://www.biomedcentral.com/14712280/8/80.
Wilcosky TC, Chambless LE (1985). A comparison of direct adjustment and regression adjustment of epidemiologic measures. Journal of Chronic Diseases, 34, 849856.
glm
, glmer
,
prLogisticDelta
,prLogisticBootCond
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29  ### For independent observations:
## Estimates from logistic regression with bootstrap confidence intervals 
## marginal standardization
# Not run:
# data("titanic", package = "prLogistic")
# attach(titanic)
# fit.logistic=glm(survived~ sex + pclass + embarked, family=binomial,
# data = titanic)
# prLogisticBootMarg(fit.logistic, data = titanic)
# End (Not run:)
# Another way for fitting the same model:
# Not run:
# prLogisticBootMarg(glm(survived~ sex + pclass + embarked,
# family=binomial, data = titanic), data=titanic)
# End (Not run:)
### For clustered data
# Estimates from randomeffects logistic regression
## with bootstrap confidence intervals  marginal standardization
# Not run:
# library(lme4)
# data("Thailand", package = "prLogistic")
# attach(Thailand)
# ML = glmer(rgi ~ sex + pped + (1schoolid),
# family = binomial, data = Thailand)
# prLogisticBootMarg(ML, data = Thailand)
# End (Not run:)

Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.