Lightweight and robust wrapper around lm.fit() for use in bootstrap procedures.
Designed to handle possible rank-deficiency gracefully by returning NA-coefficients
for linearly dependent terms. Primarily used internally within simultaneous_ci().
fit_model_q(X_full, y, q_indices)Named numeric vector of estimated coefficients. If fitting fails or coefficients are dropped due to collinearity, NA values are returned with expected names.
Numeric matrix. Full design matrix including intercept if present.
Column names must be unique. Typically derived from X + intercept inside simultaneous_ci().
Numeric vector. Response variable, same length as nrow(X_full).
Integer vector. Column indices (1-based) specifying the submodel to fit.