initial_cn <- cn(log10(pdg) ~ btensor(list(day, age)
, df = list(10, 8)
, bknots = list(c(-1, 1), c(44, 53))
)
, data = spdg)
cnr0 <- cnr(initial_cn)
plot(cnr0)
Run the code above in your browser using DataLab