similar to: image stops detecting wrong sized z matrix (PR#352)

Dear all, I'm new on R and this mailing list. We work on spatial rainfall estimation with R and Grass. We have a problem with the krige function from the sgeostat package. We would like to limit the estimated area with a polygon limit. I use a 50 points polygon to describe my work area. The krige function work quiet well without limit. But if I use this option I have the following error

Hello all! I'm a newcomer to R, so please bear with me. My problem is the following: I'm trying to use R to do some geostatistics stuff (mostly kriging), and am having trouble printing a multiple figure. When I create a single figure (one image/contour command), it works fine. The instant I try to create a multiple figure file, I start getting postscript errors. I'm using the

I?m trying to make an estimation of sea-depth for 2000 animal locations using 2 million depth soundings. The result was: Cannot allocate vector of size 20351 Kb I've tried to overcome this problem by making the following adaptations. -Increased the max-mem-size to 1.5Gb (Thanks to Thomas Lumley for his help) -Use a max.dist (10km) for the function Krige -Run the Krige computation for each

Hi, I am fitting a weighted least square regression and trying to compute SSE,SST and SSReg but I am not getting SST = SSReg + SSE and I dont know what I am coding wrong. Can you help please? xnam <-colnames(X) # colnames Design Matrix fmla1 <- as.formula(paste("Y ~",paste(xnam, collapse=

Dear List I would like to apply kriging (using sgeostat or fields package) The result I got is as follows: >Error: Cannot allocate vector of size 41148 Kb >In addition: Warning Message: >Reached total allocation of 255 Mb. See help(memory) >memory.size(max=TRUE) 250,834,944 >memory.size() 95,052,192 In other words. R should be able to allocate the vector of 41Mb!!!! What is the

Hello, R-users! When I load the sgeostat library I obtain some warnings messages. >library(sgeostat) What kind of problem will I expect to have with these messages? Loading required package: mva Loading required package: tripack Warning messages: 1: `provide' is deprecated. See ?Deprecated. in: .Deprecated() 2: `provide' is deprecated. See ?Deprecated. in: .Deprecated() Thank you

Greetings all, Ran into a strange problem with the krige function from geoR. The problem that I am having is that while the krige function seems to work well, the resulting predicted values are all NAs. Given the size of the datasets I am working with can't attach it, but I can provide snippets of the datasets. > casedata station year month day obs mpe bias type

Hello, I'm interested in doing Bayesian kriging using R. I see that the package geoR has a function that will allow one to do this (krige.bayes). However, my data are not in the form of points, but rather they are blocks that represent spatial averages (i.e., the number of fishing hooks per month in a given lat x long square). I am therefore interested in treating the data as

Hi, all, Have anyone used kringing included in R? How is it? Does it handle anisotropy data well? How does it compare with Kriging in Arc/Info? or other geostatistics software customized to do kriging or other geostatistics functions? I tried Easykriging, a geostatistics tool developed for Matlab. It has very nice GUI, but it does not provide library which i can call in my programs. so it is

Dear all, When a conventional kriging and then a contour plot is limited with a polygon (as possible with krige.conv and contour.kriging), the polygon border is displayed in black by default. > kc<-krige.conv(CZdata,loc=pred.grid,borders=czpoly,krige=krige.control(obj.m=ls)) > plot(CZcoord,xlab="x",ylab="y",type="n",asp=1) >

Hi all, I tried to produce some kriged surfaces with geoR (latest version). The size of the grid should be around 900 x 650 cells (what I find is not a very big grid), and the number of points is around 2500. The command krige.conv stopped after arround 5 min saying it can not allocate a vector with around 1.5 billion units. Sounds reasonable. Is there a workaround? How would I partition the

HI, I wonder does anyone have experience with doing sequential gaussian simulation with krige() function in gstat? I find it VERY slow compared to use krige() to achieve kriging function itself.. I wonder why, is that because it has to model the variogram, and do the kriging separately for each point to be simulated? so it would be N times slower to achieve the simulation than the kriging if

Hi all, I use kriging to interpolate the precipitation from stations, but the map of this results show lots of stripes. (please see the attachment)I think there's something wrong with the setting of the dimension of this matrix, however, I have no idea how to know or test to see if this setting is correct or not.I've tried to switch the latitude and longitude, but still got the same

Hello, Thomas > > I'd suggest to start with filling necessary bits in the > > X86RegisterInfo.cpp. This includes frame moves information and > > description of stack layout. > I looked at X86RegisterInfo.cpp, but I think that it already supports > x86-64. At least there were no obvious places were the code would only > work for 32bit mode. Right. See my last e-mail

ear ALL, The last part of my thesis analysis is the cross validation. Right now I am having difficulty using the cross validation of gstat. Below are my commands with the tsport_ace as the variable: nfold <- 3 part <- sample(1:nfold, 69, replace = TRUE) sel <- (part != 1) m.model <- x2[sel, ] m.valid <- x2[-sel, ] t<- fit.variogram(v,vgm(0.0437, "Exp", 26, 0)) cv69

Software information: R Version 1.5.1 (2002-06-17) ESS 5.1.20 Emacs 21.1.1 Mozilla 1.0 Linux Mandrake 8.2 Colleagues I am having a little trouble installing packages. Having downloaded the gzipped tarball of a package, put it in a /tmp directory, and using this command: [root at localhost /]# R CMD INSTALL -l 'usr/lib/R/library/' '/home/smc/tmp/sgeostat_1.0-18.tar.gz' I

Hi, list, I am using 'sgeostat' package by Albrecht Gebhardt and I am trying to put a correlation coefficient of some kind on the lagplots. Is there a possiblity to do so? Thanks Mahdi -- ----------------------------------- Mahdi Osman (PhD) E-mail: m_osm at gmx.net ----------------------------------- GMX im TV ... Die Gedanken sind frei ... Schon gesehen? Jetzt Spot online ansehen:

Hi R-people: Working in geoRglm, it shows me, according to AIC criterion, that the non-spacial model describes the process in a better way. It's the first time that I'm facing up to. These are my results: OP2003Seppos.AICnonsp-OP2003Seppos.AICsp #[1] -4 (OP2003Seppos.lf0.p<-exp(OP2003Seppos.lf0$beta)/(1+exp(OP2003Seppos.lf0$beta))) #P non spatial #[1] 0.9717596

Hi R-Users Is there any package in R that work Indicator Kriging? Or somebody is working in a link between GSLIB library in FORTRAN for R? Thank you very much for your help. --

Hola, # Hacemos el KED. Ver función "krige()": KED.rad <- krige( formula=pluvPcp~layer, # covariable -> radar locations=lluvia.rad.pluv.spdf, newdata=radarGrid, # podría ser cualquier objeto Spatial model=v.fit, # modelo de semivariograma. maxdist=Inf