similar to: please help me

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

Dear list - Does anyone have any ideas / comments about why I am receiving the following warning when I run lsoda: 1: lsoda-- at t (=r1), too much accuracy requested in: lsoda(start, times, model, parms) 2: for precision of machine.. see tolsf (=r2) in: lsoda(start, times, model, parms) I have tried changing both rtol and atol but without success. I saw the thread in the

All, I am new to using R and know some basics. I wish to use kriging in R to do the following: given data Y =f(X1,X2,X3,.....,Xn) --1000+ irregular measured data set. I would like to be able to get a single value y given sinle input set (x1,x2,x3,...xn) A google search on this takes me lierally to the same example on involving analysis with soil sampling and I cannot figure out how to

Hello, I need to get a point estimate and SD for a proportion, but the subjects' data are not binary---they are proportions (of doses received). That is, I have a proportion for each subject. In the past I have analysed these data as a continuous (normal) variable, but I really don't want CIs over 100%. This seems like basic stuff, but I don't remember learning it and it's

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

Dear Aaron, I am really a tool user and not a tool maker (actually an ecologist doing some biostatistics)... so, I take the liberty of sending a copy of this e-mail to the r-help list where capable computer persons and true statisticians may provide more relevant information and also to Paulo Ribeiro and Peter Diggle, the authors of geoR.. I really feel that your huge matrix cannot be

Hello again, There is any package which does Z aware (real 3D) interpolations? It can be any method (IDW, kriging or spline) but it should take into consideration not only x and y coordinates for interpolation, but also z coordinate. I looked into different packages but it seems i didn't find the right one. The ultimate goal is to import the output into a GIS (Geographical Information

Hello All! I calculate a variogram using the function variog (package geoR) afterwards I use variofit to fit a spherical model (see code below). Now I just changed the units of the variable (in this case MPa to kPa just a factor of 1000). If I do so, I get a different fit and therefore different ranges etc. Why? The semi-variance is of course 6 orders of magnitude higher but the values

A majority of my data makes a kriged map perfectly using an exponential model for the semivariogram to fit my data and then going through the commands variofit() to define the model and then krige.conv() to use the model to predict values in a grid. But?one set of my data appears to be linearly correlated for the first 5000 meters and not correlated beyond that. I have been having problems

Dear all, I'm trying to fit a variogram model using variofit function in geoR package. There is an option to fix the nugget but is there anyway to force the variogram sill to equal some defined value? I'm working with standardized longitudinal data so the process variance is 1 and I'd like to force the sill of the variogram to reach this value (and prevent it from going

I'm performing ordinary and simple kriging from a set of non-negative values: >krige.control.sk<-krige.control(type.krige="sk",obj.model=my.variogram.model,beta=my.variogram.model$beta) >sk<-krige.conv(cadiz.geo,data=amostragem.cadiz$dia1.std,locations=as.matrix(cadiz.polygrid[,c(1,2)]),krige=krige.control.sk) krige.conv: model with constant mean krige.conv: Kriging

Juan Zwolinski wrote: > > I'm performing ordinary and simple kriging from a set of non-negative > values: > > > >krige.control.sk<-krige.control(type.krige="sk",obj.model=my.variogram.model,beta=my.variogram.model$beta) > >

Dear List, I have fitted a spherical function to my variogram using "variofit(...)" from GeoR. Now I would like to predict some data with the function "predict(object,...)" from package stats. Does anyone know wether this works and if it does how to do it? Thanks a lot! Anna