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
|
% Generated by roxygen2: do not edit by hand
% Please edit documentation in R/incidence_rates.R
\name{rate}
\alias{rate}
\title{Direct-Standardised Incidence/Mortality Rates}
\usage{
rate(
data,
obs = NULL,
pyrs = NULL,
print = NULL,
adjust = NULL,
weights = NULL,
subset = NULL
)
}
\arguments{
\item{data}{aggregated data (see e.g. \verb{[lexpand]},
\verb{[aggre]} if you have subject-level data)}
\item{obs}{observations variable name in data.
\link[=flexible_argument]{Flexible input}, typically e.g.
\code{obs = obs}.}
\item{pyrs}{person-years variable name in data.
\link[=flexible_argument]{Flexible input}, typically e.g.
\code{pyrs = pyrs}.}
\item{print}{variable name to stratify the rates.
\link[=flexible_argument]{Flexible input}, typically e.g.
\code{print = sex} or \code{print = list(sex, area)}.}
\item{adjust}{variable for adjusting the rates.
\link[=flexible_argument]{Flexible input}, typically e.g.
\code{adjust = agegroup}.}
\item{weights}{typically a list of weights or a \code{character} string
specifying an age group standardization scheme; see
the \link[=direct_standardization]{dedicated help page}
and examples.}
\item{subset}{a logical expression to subset data.}
}
\value{
Returns a \code{data.table} with observations, person-years, rates and
adjusted rates, if available. Results are stratified by \code{print}.
Adjusted rates are identified with suffix \code{.adj} and
\code{.lo} and \code{.hi} are for confidence intervals lower and upper
95\% bounds, respectively.
The prefix \code{SE.} stands for standard error.
}
\description{
\code{rate} calculates adjusted rates using
preloaded weights data or user specified weights.
}
\details{
Input data needs to be in aggregated format with observations
and person-years. For individual data use \verb{[lexpand]}, or
\verb{[ltable]} and merge person-years manually.
The confidence intervals are based on the normal approximation of the logarithm of the rate.
The variance of the log rate that is used to derive the confidence intervals
is derived using the delta method.
}
\examples{
## Prepare data with lexpand and then reformat agegroup.
data(sibr)
x <- lexpand(sibr, birth = bi_date, entry = dg_date, exit = ex_date,
breaks = list(per = c(1990,2000,2010,2020), age = c(0:17*5,Inf)),
aggre = list(agegroup = age, year.cat = per),
status = status != 0)
x$agegroup <- cut(x$agegroup, c(0:17*5,Inf), right = FALSE)
## calculate rates for selected periods with Nordic 2000 weights:
r1 <- rate( data = x, obs = from0to1, pyrs = pyrs, print = year.cat,
adjust = agegroup, weights = 'nordic')
r1
## use total person-years by stratum as weights (some have zero)
w <- ltable(x, by.vars = "agegroup", expr = sum(pyrs))
w[is.na(w$V1),]$V1 <- 0
r2 <- rate( data = x, obs = from0to1, pyrs = pyrs, print = year.cat,
adjust = agegroup,
weights = w$V1)
r2
## use data.frame of weights:
names(w) <- c("agegroup", "weights")
r2 <- rate( data = x, obs = from0to1, pyrs = pyrs, print = year.cat,
adjust = agegroup,
weights = w)
r2
## internal weights (same result as above)
r3 <- rate( data = x, obs = from0to1, pyrs = pyrs, print = year.cat,
adjust = agegroup,
weights = "internal")
r3
}
\seealso{
\verb{[lexpand]}, \verb{[ltable]}
Other main functions:
\code{\link{Surv}()},
\code{\link{relpois}()},
\code{\link{relpois_ag}()},
\code{\link{sir}()},
\code{\link{sirspline}()},
\code{\link{survmean}()},
\code{\link{survtab}()},
\code{\link{survtab_ag}()}
Other rate functions:
\code{\link{rate_ratio}()}
}
\author{
Matti Rantanen, Joonas Miettinen
}
\concept{main functions}
\concept{rate functions}
|
