Construct a list of hero_bsplines using the
default values suggested by Ruppert, Wand, and Carroll
(2003). Specifically, if
r = range(evalargs[[i]]) and
l = length(evalargs[[i]]), then Ruppert, Wand,
and Carroll (2003) suggest
nknots = min(ceiling(l/4), 35) and the function
returns the hero_bspline for that dimension as
bspline(r, nknots = nknots).