# The data sets were generated as follows
set.seed(1782)
# Simulate training data
nTrain = 1000
B = c(0.1, 0.5, 1.2, -0.75, 0.8)
X = replicate(4, rnorm(nTrain))
p0true = binomial()$linkinv(cbind(1, X) %*% B)
y = rbinom(nTrain, 1, p0true)
colnames(X) = paste0("x", seq_len(ncol(X)))
traindata = data.frame(y, X)
# Simulate validation data
nTest = 500
X = replicate(4, rnorm(nTest))
p0true = binomial()$linkinv(cbind(1, X) %*% B)
y = rbinom(nTest, 1, p0true)
colnames(X) = paste0("x", seq_len(ncol(X)))
testdata = data.frame(y, X)
Run the code above in your browser using DataLab