[R-sig-Geo] Handling realisations in geoR

[Anders Malmberg]

Hi,

I have replicates of a spatial field to which I would like fit a spatial 
covariance function
using geoR and its function likfit. My problem is that I do not know how 
to specify this in geoR.

I tried to simply rowstack my fields such that I have

       x   y   data
"1" x1 y1 d11
"2" x2 y2 d12
"3" x3 y3 d13
"4" x4 y4 d14
"1" x1 y1 d21
"2" x2 y2 d22
"3" etc...
"4"
"1"
"2"
"3"
"4" x4 y4 d34

stored in an matrix "obj" let's say. How shall I use 
as.geodata(obj,realisations=XXX)?
And how should I use likfit??? It is no problem to rearrange my matrix.

Best regards,
Anders Malmberg
andersm at maths.lth.se