[R-sig-Geo] Error message from likfit.glsm in geoRglm

Ruben Roa RRoa at fisheries.gov.fk
Thu Jun 30 14:14:56 CEST 2005


Hi:

I am trying to fit a spatial generalized linear mixed model
with geoRglm (v.0.8-11, geoR version 1.5-7, R 2.1.0). The 
family is binomial and the link is logit. There are no covariates, 
so beta is only the intercept. In the data i have 4 columns, 
namely easting, northing, number of successes and number of trials.
These are the features of the model for generating the MCMC samples:

$cov.pars     $beta    $cov.model   $nugget $kappa  $family        $link
[1] 0.35 7.00 [1] 0.6  [1] "matern" [1] 0.1 [1] 1.5 [1] "binomial" [1] "logit"

I get the following error message:

Error in rep(no, length.out = length(ans)) : 
        argument "no" is missing, with no default

However, i don't see the "no" argument in the likfit.glsm
function. Neither does it appear as an argument for optim(). 
So I am lost. Can anyone give a clue of which specification i
am missing. Btw, this was my command:

c1f.lf<-likfit.glsm(c1i.pre.lf,ini.phi=7,hessian=TRUE)

Thanks in advance,

Ruben




More information about the R-sig-Geo mailing list