# NOT RUN {
# 90% confidence intervals for \code{p} and \code{k}
mod <- pkm(formula_p = p ~ Visibility * Season, formula_k = k ~ Season,
data = wind_cleared$SE)
qpk(p = c(0.05, 0.95), model = mod)
# }
Run the code above in your browser using DataLab