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)
.