[R-sig-Geo] Handling realisations in geoRglm

Paulo Justiniano Ribeiro Jr paulojus at c3sl.ufpr.br
Thu Aug 31 05:53:51 CEST 2006


I'm afraid with the current implementation the answers are no and no ...
geoRglm cannot handle realisations as geoR does (e.g. considering
different years as independent replications of the process).

Furthermore neither geoR nor geoRglm can fit spatial-temporal models
(with both, space and time dependence)
RandomFields can do this but, to my knowledge, only for Gaussian data


On Wed, 30 Aug 2006, Chris Barker wrote:

> Dear R users:
> I want to model mosquito count data based on landcover attributes and
> meteorological variables using a Poisson GLSM in the geoRglm package.  I
> have monthly mosquito counts over more than 20 years with repeated
> observations from individual trap sites over time.  I have used
> as.geodata() to successfully read my dataset into the geodata format
> utilized by geoR and geoRglm, including specification of the
> realisations.  (I know the realisations have been specified correctly
> because I can get data summaries by realisation from the geodata
> object.)  However, after reading the packages' documentation and
> searching the mailing lists and other sources, it seems that the
> functions in geoRglm do not acknowledge the existence of multiple
> realisations per site.  I see that the likfit() function in geoR has an
> argument for realisations, but I cannot find anything similar in
> geoRglm.  Is it possible to model data with repeated realisations from
> the same sites using geoRglm?  If so, how can this be done and is there
> a way to model dependence among the realisations over time?
> Thank you,
> Chris Barker
> Ph.D. Candidate
> University of California, Davis
> cmbarker at ucdavis.edu
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo at stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo

More information about the R-sig-Geo mailing list