data-raw/get_data.R

In rgcstats/loglogitlink: A New Link Function for Binary Regression that Smoothly Combines the Log and Logit Functions

# create an example dataset

set.seed(30035)
N <- 500
x1 <- 2*runif(N)-1
x2 <- 2*runif(N)-1
#lp <- -0.5+1.5*x1+1.5*x2
#p <- exp(lp)/(1+exp(lp))
lp <- -1+1.5*x1/2+1.5*x2/2
p <- loglogit.linkinv(eta=lp,cutover=0.8)
y <- 1*(runif(N)<=p)
loglogit.example <- data.frame(y=y,x1=x1,x2=x2)
save(loglogit.example,file="data/loglogit.example.rdata")