Hi Edzer, - number of variables = 1 (might be 15 in the future) - number of observations per variable = 1,400 - number of cells in mask maps = 9,262,500 (but ~6 million missing values) - number of mask (input, output) maps = 1 - number of simulations = 500 - (max.) size of kriging neighbourhood = 20 points
Maybe there is a more efficient way of providing input maps? In a different more compressed format and without missing values? > BTW did you resolve the matrix singularity problem? If so, how? Yes. It was caused by duplicates. The average option did not work for simulations but did work for kriging. Cheers, Ulrich -- __________________________________________________ Ulrich Leopold MSc. Department of Physical Geography Institute for Biodiversity and Ecosystem Dynamics Faculty of Science University of Amsterdam Nieuwe Achtergracht 166 NL-1018WV Amsterdam Phone: +31-(0)20-525-7456 (7451 Secretary) Fax: +31-(0)20-525-7431 Email: [EMAIL PROTECTED] http://www.frw.uva.nl/soil/Welcome.html Check us also out at: Netherlands Centre for Geo-ecological Research http://www.frw.uva.nl/icg