[R] Problem with Zelig - gam.poisson

Miluji Sb milujisb at gmail.com
Fri Dec 23 12:05:52 CET 2016


I am getting the following error when trying to run a gam.possion model
with Zelig.

Error in eval(expr, envir, enclos) : attempt to apply non-function

Even the example won't run.

library(mgcv)
library(Zelig)

n <- 400
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
f3 <- function(x) 0 * x
f <- f0(x0) + f1(x1) + f2(x2)
g <- exp(f/4)
y <- rpois(rep(1, n), g)
my.data <- as.data.frame(cbind(y, x0, x1, x2, x3))

z.out <- zelig(y ~ s(x0) + s(x1) + s(x2) + s(x3), model = "gam.poisson",
data = my.data)

What am I doing wrong? Session info is below. Thanks a lot!

Sincerely,

Milu

####
sessionInfo()
R version 3.3.2 (2016-10-31)
Platform: x86_64-w64-mingw32/x64 (64-bit)
Running under: Windows 7 x64 (build 7601) Service Pack 1

locale:
[1] LC_COLLATE=English_United States.1252  LC_CTYPE=English_United
States.1252    LC_MONETARY=English_United States.1252
[4] LC_NUMERIC=C                           LC_TIME=English_United
States.1252

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base

other attached packages:
[1] Zelig_5.0-13 mgcv_1.8-15  nlme_3.1-128

loaded via a namespace (and not attached):
 [1] Rcpp_0.12.8        nloptr_1.0.4       plyr_1.8.4
miscTools_0.6-22   tools_3.3.2        lme4_1.1-12
 [7] MatchIt_2.4-21     jsonlite_1.1       tibble_1.2
lattice_0.20-34    Matrix_1.2-7.1     DBI_0.5-1
[13] maxLik_1.3-4       parallel_3.3.2     SparseM_1.74       coda_0.18-1
     AER_1.2-4          dplyr_0.5.0
[19] MatrixModels_0.4-1 nnet_7.3-12        stats4_3.3.2       lmtest_0.9-34
     grid_3.3.2         R6_2.2.0
[25] geepack_1.2-1      survival_2.40-1    VGAM_1.0-2
foreign_0.8-67     minqa_1.2.4        Amelia_1.7.4
[31] Formula_1.2-1      car_2.1-3          magrittr_1.5       mcmc_0.9-4
      MASS_7.3-45        splines_3.3.2
[37] pbkrtest_0.4-6     assertthat_0.1     quantreg_5.29
 sandwich_2.3-4     survey_3.31-5      MCMCpack_1.3-8
[43] lazyeval_0.2.0     zoo_1.7-13

	[[alternative HTML version deleted]]



More information about the R-help mailing list