ZeligGAM: Zelig for Generalized Additive Models

set.seed(0)
n <- 100
sig <- 2
x0 <- runif(n, 0, 1)
x1 <- runif(n, 0, 1)
x2 <- runif(n, 0, 1)
x3 <- runif(n, 0, 1)
f0 <- function(x)
  2 * sin(pi * x)
f1 <- function(x)
  exp(2 * x)
f2 <- function(x)
  0.2 * x ^ 11 * (10 * (1 - x)) ^ 6 + 10 * (10 * x) ^ 3 * (1 - x) ^ 10
f2 <- function(x)
  0 * x
f <- f0(x0) + f1(x1) + f2(x2)

g <- (f - 5) / 3
g <- binomial()$linkinv(g)

y <- rbinom(g,1,g)

my.data <- as.data.frame(cbind(y, x0, x1, x2, x3))


# TEST LOGIT

z5 <- zlogitgam$new()
z5
z5$zelig(y ~ s(x0) + s(x1) + s(x2) + s(x3), data = my.data)

z.out1 <- zelig(y ~ s(x0) + s(x1) + s(x2) + s(x3), data = my.data, model = "normal.gam")
z.out2 <- zelig(y ~ s(x0) + s(x1) + s(x2) + s(x3), data = my.data, model = "logit.gam")
z.out3 <- zelig(y ~ s(x0) + s(x1) + s(x2) + s(x3), data = my.data, model = "probit.gam")
z.out4 <- zelig(y ~ s(x0) + s(x1) + s(x2) + s(x3), data = my.data, model = "poisson.gam")