1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
|
% Generated by roxygen2: do not edit by hand
% Please edit documentation in R/sir.R
\name{sirspline}
\alias{sirspline}
\title{Estimate splines for SIR or SMR}
\usage{
sirspline(
coh.data,
coh.obs,
coh.pyrs,
ref.data = NULL,
ref.obs = NULL,
ref.pyrs = NULL,
ref.rate = NULL,
subset = NULL,
print = NULL,
adjust = NULL,
mstate = NULL,
spline,
knots = NULL,
reference.points = NULL,
dependent.splines = TRUE
)
}
\arguments{
\item{coh.data}{cohort data with observations and at risk time variables}
\item{coh.obs}{variable name for observed cases}
\item{coh.pyrs}{variable name for person-years in cohort data}
\item{ref.data}{aggregated population data}
\item{ref.obs}{variable name for observed cases}
\item{ref.pyrs}{variable name for person-years in population data}
\item{ref.rate}{population rate observed/expected. This overwrites the parameters
\code{ref.pyrs} and \code{ref.obs}.}
\item{subset}{logical condition to subset \code{coh.data} before any computations}
\item{print}{variable names for which to estimate SIRs/SMRs and
associated splines separately}
\item{adjust}{variable names for adjusting the expected cases}
\item{mstate}{set column names for cause specific observations. Relevant only
when coh.obs length is two or more. See help for \code{sir}.}
\item{spline}{variable name(s) for the splines}
\item{knots}{number knots (vector), pre-defined knots (list of vectors) or for optimal number of knots left NULL}
\item{reference.points}{fixed reference values for rate ratios. If left \code{NULL}
the smallest value is the reference point (where SIR = 1).
Ignored if \code{dependent.splines = FALSE}}
\item{dependent.splines}{logical; if TRUE, all splines are fitted in same model.}
}
\value{
A list of data.frames and vectors.
Three spline estimates are named as \code{spline.est.A/B/C} and the corresponding values
in \code{spline.seq.A/B/C} for manual plotting
}
\description{
Splines for standardised incidence or mortality ratio. A useful
tool to e.g. check whether a constant SIR can be assumed for all calendar periods,
age groups or follow-up intervals. Splines can be fitted for these time dimensions
separately or in the same model.
}
\details{
See \verb{[sir]} for help on SIR/SMR estimation in general; usage of splines
is discussed below.
\strong{The spline variables}
The model can include one, two or three splines variables.
Variables can be included in the same model selecting \code{dependent.splines = TRUE}
and SIR ratios are calculated (first one is the SIR, others SIR ratios).
Reference points vector can be set via \code{reference.points}
where first element of the vector is the reference point for first ratio.
Variable(s) to fit splines are given as a vector in argument \code{spline}.
Order will affect the results.
\strong{dependent.splines}
By default dependent.splines is FALSE and all splines are fitted in separate models.
If TRUE, the first variable in \code{spline} is a function of a SIR and other(s) are ratios.
\strong{knots}
There are three options to set knots to splines:
Set the number of knots for each spline variable with a \strong{vector}.
The knots are automatically placed to the quantiles of observed cases in cohort data.
The first and last knots are always the maximum and minimum values, so knot
value needs to be at least two.
Predefined knot places can be set with a \strong{list} of vectors.
The vector for each spline in the list specifies the knot places. The lowest
and the largest values are the boundary knots and these should be checked beforehand.
If \code{knots} is left \strong{NULL}, the model searches the optimal number
of knots by model AIC by fitting models iteratively from 2 to 15 knots and
the one with smallest AIC is selected.
If \code{dependent.splines = TRUE}, the number of knots is searched by fitting each spline
variable separately.
\strong{print}
Splines can be stratified by the levels of variable given in \code{print}. If
\code{print} is a vector, only the first variable is accounted for. The knots
are placed globally for all levels of \code{print}. This also ensures that the likelihood
ratio test is valid.
Splines are also fitted independently for each level of \code{print}.
This allows for searching interactions, e.g. by fitting spline for period
(\code{splines='period'}) for each age group (\code{print = 'agegroup'}).
\strong{p-values}
The output p-value is a test of whether the splines are equal (homogenous)
at different levels of \code{print}.
The test is based on the likelihood ratio test, where the full model
includes \code{print} and is
compared to a null model without it.
When \code{(dependent.splines = TRUE)} the p-value returned is a global p-value.
Otherwise the p-value is spline-specific.
}
\examples{
\donttest{
## for examples see: vignette('sir')
}
}
\seealso{
\verb{[splitMulti]}
\href{../doc/sir.html}{A SIR calculation vignette}
Other sir functions:
\code{\link{lines.sirspline}()},
\code{\link{plot.sirspline}()},
\code{\link{sir}()},
\code{\link{sir_exp}()},
\code{\link{sir_ratio}()}
Other main functions:
\code{\link{Surv}()},
\code{\link{rate}()},
\code{\link{relpois}()},
\code{\link{relpois_ag}()},
\code{\link{sir}()},
\code{\link{survmean}()},
\code{\link{survtab}()},
\code{\link{survtab_ag}()}
}
\author{
Matti Rantanen, Joonas Miettinen
}
\concept{main functions}
\concept{sir functions}
|
