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
|
options(na.action=na.exclude) # preserve missings
options(contrasts=c('contr.treatment', 'contr.poly')) #ensure constrast type
library(survival)
{if (is.R()) mdy.date <- function(m, d, y) {
y <- ifelse(y<100, y+1900, y)
as.Date(paste(m,d,y, sep='/'), "%m/%d/%Y")
}
else mdy.date <- function(m,d,y) {
y <- ifelse(y<100, y+1900, y)
timeDate(paste(y, m, d, sep='/'), in.format="%Y/%m/%d")
}
}
#
# A test of the match.ratetable function, specifically the
# change to allow partial matching of strings
# Note that 10,000 days old is 27.4 years
#
aeq <- function(x,y, ...) all.equal(as.vector(x), as.vector(y), ...)
temp1 <- data.frame(year=mdy.date(2,2,1960:1964),
age = 10000 + 1:5,
sex = c('M', 'fema', 'f', 'ma', 'F'))
temp2 <- ratetable(year=temp1$year, age=temp1$age, sex=temp1$sex)
temp3 <- match.ratetable(temp2, survexp.us)
aeq(temp3$R[,2], c(1,2,2,1,2))
|
