[Gstat-info] Using 'average' data option in 2.5.1

[John A Stevenson]

Hi,

I am using gstat2.5.1 to do some kriging.  When I run gstat, I get a 
singular matrix error (below - 1).  If I change my command file to 
include the 'average' command (below - 2), it gets confused and can't 
read it (below - 3). 

Is there a particular syntax that I am missing?

Cheers

John

p.s. I am running this in the GRASS shell, but I don't think that this 
has anything to do with it.

------------------------------------------------

(1) - Gstat original output

GRASS 6.4.0RC5 (thingy):~/laptop/thingy/chirp > ~/usr/local/bin/gstat 
papmap_krige.gstat
gstat: Linux version 2.5.1 (02 April 2007)
Copyright (C) 1992, 2006 Edzer J. Pebesma
using Marsaglia's random number generator
data(papmap): papmap_combo.xyz     (Ascii table file)
attribute:             col[3]     [x:] x_1        : [    485010,    495190]
n:                     162768     [y:] y_2        : 
[7.10721e+06,7.11859e+06]
sample mean:           279.29     sample std.:           100.74    
[using ordinary kriging]
initializing maps .
Warning: at least some data are outside the map boundaries

"chfactor.c", line 130: singular matrix in function LDLfactor()

gstat caught an error that occurred in the matrix library,
the reason for it was: singular matrix

HINT: Read the manual at http://www.gstat.org/ ;
look for: Trouble shooting -> Error messages -> From meschach

gstat: matrix library error: singular matrix


--------------------------------------------------
(2) Command file

#
# gstat command file, Linux version 2.5.1 (02 April 2007)
# Thu Nov 26 12:19:06 2009
#
data(papmap): 'papmap_combo.xyz', x=1, y=2, v=3, average, min=10, 
max=40, radius=1000;
variogram(papmap): 11376.7 Sph(6410.99);
set fit = 1;
set zero = 0.01;
mask: 'mask.asc';
predictions(papmap): 'papmap_kriged.asc';





--------------------------------------------------
(3) Gstat new output

GRASS 6.4.0RC5 (thingy):~/laptop/thingy/chirp > ~/usr/local/bin/gstat 
papmap_krige.gstat
gstat: Linux version 2.5.1 (02 April 2007)
Copyright (C) 1992, 2006 Edzer J. Pebesma
data(papmap): 'papmap_combo.xyz', x=1, y=2, v=3, average, min=10, max=40
                                                       ^
gstat:papmap_krige.gstat:5: `,' unexpected
gstat: syntax error: papmap_krige.gstat


-- 


Dr John Stevenson
Postdoctoral Research Associate
School of Earth, Atmospheric and Environmental Sciences
Williamson Building (Room 2.42)
University of Manchester
Manchester M13 9PL, UK
tel. +44(0)161 306 6585; fax. +44(0)161 306 9361;
john.stevenson at manchester.ac.uk