modelbound.boot: Estimating qgcomp regression line confidence bounds

View source: R/base_bounds.R

modelbound.bootR Documentation

Estimating qgcomp regression line confidence bounds


Calculates: expected outcome (on the link scale), and upper and lower confidence intervals (both pointwise and simultaneous)


modelbound.boot(x, alpha = 0.05, pwonly = FALSE)



"qgcompfit" object from qgcomp.glm.boot,


alpha level for confidence intervals


logical: return only pointwise estimates (suppress simultaneous estimates)


This method leverages the bootstrap distribution of qgcomp model coefficients to estimate pointwise regression line confidence bounds. These are defined as the bounds that, for each value of the independent variable X (here, X is the joint exposure quantiles) the 95% bounds (for example) for the model estimate of the regression line E(Y|X) are expected to include the true value of E(Y|X) in 95% of studies. The "simultaneous" bounds are also calculated, and the 95% simultaneous bounds contain the true value of E(Y|X) for all values of X in 95% of studies. The latter are more conservative and account for the multiple testing implied by the former. Pointwise bounds are calculated via the standard error for the estimates of E(Y|X), while the simultaneous bounds are estimated using the bootstrap method of Cheng (reference below). All bounds are large sample bounds that assume normality and thus will be underconservative in small samples. These bounds may also inclue illogical values (e.g. values less than 0 for a dichotomous outcome) and should be interpreted cautiously in small samples.


Cheng, Russell CH. "Bootstrapping simultaneous confidence bands." Proceedings of the Winter Simulation Conference, 2005.. IEEE, 2005.


A data frame containing


The linear predictor from the marginal structural model


The canonical measure (risk/odds/mean) for the marginal structural model link


the stndard error of linpred


Confidence bounds for the effect measure, and bounds centered at the canonical measure (for plotting purposes)

The confidence bounds are either "pointwise" (pw) and "simultaneous" (simul) confidence intervals at each each quantized value of all exposures.

See Also



## Not run: 
dat <- data.frame(x1=(x1 <- runif(50)), x2=runif(50), x3=runif(50), z=runif(50),
ft <- qgcomp.glm.boot(y ~ z + x1 + x2 + x3, expnms=c('x1','x2','x3'), data=dat, q=5)
modelbound.boot(ft, 0.05)

## End(Not run)

qgcomp documentation built on Aug. 10, 2023, 5:07 p.m.