[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