acp <- cp(log10(pdg) ~ bsplines(age, df = 12, bknots = c(45, 53))
, data = spdg
, keep_fit = TRUE)
acp_pred0 <- predict(acp$fit, se.fit = TRUE)
acp_pred <- predict(acp, se.fit = TRUE)
all.equal(acp_pred0, acp_pred)
Run the code above in your browser using DataLab