Hi, you define `knots` but then don't use it in the subsequent `bs` call knots <- unname(quantile(X,probs=seq(from=0, to=1, length.out = num_knots)))