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
|
% Generated by roxygen2: do not edit by hand
% Please edit documentation in R/Rcode.r
\name{nessie}
\alias{nessie}
\title{Net Expected Sample Size Is Estimated}
\usage{
nessie(
formula = formula(data),
data = parent.frame(),
ratetable = relsurv::slopop,
times,
rmap
)
}
\arguments{
\item{formula}{a formula object, same as in \code{rs.surv}. The right-hand
side of the formula object includes the variable that defines the subgroups
(a variable of type \code{factor}) by which the expected sample size is to
be calculated.}
\item{data}{a data.frame in which to interpret the variables named in the
\code{formula}.}
\item{ratetable}{a table of event rates, organized as a \code{ratetable}
object, such as \code{slopop}.}
\item{times}{Times at which the calculation should be evaluated - in years!}
\item{rmap}{an optional list to be used if the variables are not organized
and named in the same way as in the \code{ratetable} object. See details of
the \code{rs.surv} function.}
}
\value{
A list of values.
}
\description{
Calculates how the sample size decreases in time due to population mortality
}
\details{
The function calculates the sample size we can expect at a certain time
point if the patients die only due to population causes (population survival
* initial sample size in a certain category), i.e. the number of individuals
that remains at risk at given timepoints after the individuals who die due
to population causes are removed. The result should be used as a guideline
for the sensible length of follow-up interval when calculating the net
survival.
The first column of the output reports the number of individuals at time 0.
The last column of the output reports the conditional expected (population)
survival time for each subgroup.
}
\examples{
data(slopop)
data(rdata)
rdata$agegr <-cut(rdata$age,seq(40,95,by=5))
nessie(Surv(time,cens)~agegr,rmap=list(age=age*365.241),
ratetable=slopop,data=rdata,times=c(1,3,5,10,15))
}
\references{
Pohar Perme, M., Pavlic, K. (2018) "Nonparametric Relative
Survival Analysis with the R Package relsurv". Journal of Statistical
Software. 87(8), 1-27, doi: "10.18637/jss.v087.i08"
}
\seealso{
\code{rs.surv}
}
\keyword{survival}
|
