Dear help list,

I am using R version 3.0.1 (2013-05-16) Platform: i386-w64-mingw32/i386 
(32-bit) and mgcv_1.7-24.

I am estimating a geo-additive gam model using the "bam" command in the 
mgcv package. I do have a large dataset consisting of 115866 
observations with 14 potential covariates and an additional spatial effect:

my_k <- 40                # specifies the dimension of the basis
my_k_low <-  5

bam_mod  <- bam(y ~ s(x1, k = my_k, bs = "cr")   +  s(x2, k = my_k, bs = 
"cr") + s(x3, k = my_k_low, bs = "cr") +...+ s(x14, k = my_k, bs = "cr") 
+ s(district, bs = "mrf", xt = list(polys = my_boundary)),
                       data =my_data, gamma = 1.4, gc.level = 2)


where I have reduced my_k to my_k_low for "x3" since for that covariate 
there are not enough unique values to support the dimension of the basis 
set to 40. The covariate "district" carries the discrete spatial 
information (which district does each observations belong to) and 
"my_boundary" is a list defining the neighbourhood structure of the 
districts. The model is estimated without any error or warning messages. 
However, calling the summary gives the following error:

*summary(bam_mod)****Error in eigen(V, symmetric = TRUE) : infinite or 
missing values in 'x'*

I have removed all NA values from the dataset and do not know why the 
error message appears.

Any ideas?

Best
Alexander

	[[alternative HTML version deleted]]

