# generate data
n <- 500
x <- data.frame(x1=rnorm(n))
pi.true <- 1/(1+exp(.3+1.5*x$x1))
dat <- RRgen(n, pi.true=pi.true, model="Warner", p=.1)
x$response <- dat$response
# fit and plot model
mod <- RRlog(response ~ x1, data=x, model="Warner", p=.1)
plot(mod, "x1" ,ci=.95)
Run the code above in your browser using DataLab