powered by
Apply Spline Constraint to Treatment Effect Coefficients
apply_spline_constraint(b0, spline_var, knot, zeta, log_hrs, k_treat, verbose)