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