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
|
cr.setup <- function(y) {
yname <- as.character(substitute(y))
if(!is.category(y)) y <- factor(y, exclude=NA) # was category 20may02
y <- oldUnclass(y) # in case is.factor
ylevels <- levels(y)
kint <- length(ylevels) - 1
y <- as.integer(y-1)
reps <- ifelse(is.na(y), 1, ifelse(y < kint-1, y+1, kint))
subs <- rep(1:length(y), reps)
cuts <- vector('list',kint+2)
cuts[[1]] <- NA
for(j in 0:kint) cuts[[j+2]] <- if(j < kint-1) 0:j else 0:(kint-1)
cuts <- unlist(cuts[ifelse(is.na(y),1,y+2)])
y <- rep(y, reps)
Y <- 1*(y==cuts)
labels <- c('all', paste(yname,'>=',ylevels[2:kint],sep=''))
cohort <- factor(cuts, levels=0:(kint-1), labels=labels)
list(y=Y, cohort=cohort, subs=subs, reps=reps)
}
|
