Check if object is potential impact fraction
is_pif.Rd
Method for checking whether the object is a potential impact fraction
Note
The population attributable fraction is a type of potential impact fraction so
the is_pif
method will return TRUE
even when estimating paf()
.
Examples
#Example 1
data(ensanut)
options(survey.lonely.psu = "adjust")
design <- survey::svydesign(data = ensanut, ids = ~1, weights = ~weight, strata = ~strata)
rr <- function(X, theta) {
exp(-2 +
theta[1] * X[, "age"] + theta[2] * X[, "systolic_blood_pressure"] / 100)
}
cft <- function(X) {
X[, "systolic_blood_pressure"] <- X[, "systolic_blood_pressure"] - 5
return(X)
}
pifsim <- pif(design,
theta = log(c(1.05, 1.38)), rr, cft,
additional_theta_arguments = c(0.01, 0.03), n_bootstrap_samples = 10,
)
is_pif(pifsim)
#> [1] TRUE