time_spline(Theoph$Time, df = 3)
# Or, compute the spline basis beforehand, and then pass it to time_spline()
basis <-
splines::bs(Theoph$Time, df = 3, Boundary.knots = c(0, max(Theoph$Time)))
time_spline(Theoph$Time, basis = basis)
Run the code above in your browser using DataLab