##Kline (2004) Chapter 5
x1 <- c("relapsed", "not relapsed")
y1 <- c("control", "treatment")
dat <- data.frame(y =
factor(c(rep(x1, c(60, 40)), rep(x1, c(40, 60))), levels=x1),
x = factor(rep(y1, each=100), levels=y1)
)
ind.prop(y~x, data=dat) #count relapse
ind.prop(y~x, data=dat, lev.prop=2) #count not relapse
Run the code above in your browser using DataLab