mk2D: Make 2D-spline fits of incidences
Description
Produces two dimensional (2D) spline
fits of cancer incidence versus age and calendar year, with interactions. In conjunction with person years
(PY) at risk, this
is used in csd() to produce expected numbers of cases under a null hypothesis
that prior cancers do not impact subsequent cancer risks.
Usage
mk2D(seerSet, knots=5, write=FALSE, outDir="~/Results",txt=NULL,secondS=NULL)
Arguments
seerSet
Object of class seerSet, i.e. output list of seerSet().
knots
Base number of knots; overrides are in place for some cancers.
write
TRUE = write 2D fits to files. The fits can be >300 MB and take >60 seconds to write,
so leave FALSE unless you need it.
outDir
Folder that will hold the output files.
txt
Additional text to distinguish files with different cancer lists. This may be
useful during spline fit development.
secondS
Charcter vector of second cancers of interest (note: I often capitalize
the final S of vectors of Strings).
Value
The input seerSet with an additional data frame D added to this list.
D holds background/expected incidences over a 1-year resolution age-year grid.
Examples
Run this code# NOT RUN {
library(SEERaBomb)
(pm=simSeerSet())
(pm=mk2D(pm))
names(pm)
head(pm$D)
tail(pm$D)
# }
Run the code above in your browser using DataLab