Displaying 20 results from an estimated 600 matches similar to: "question about fitted values from geoR - results 'too good'"
2003 Jan 17
3
kriging in R
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
2010 Aug 19
0
2d kriging with anisotropy on an irregular network (RandomFields Package)
Dear List
I am using the RandomFields package, and I have a problem when 2d-kriging, with an anisotropy, some values from an irregular network.
It works well when :
- 2d-kriging, without any anisotropy, some data from an irregular network
- 2d-kriging, with and without anisotropy, some data from a regular network
- 3d-kriging, with and without anisotropy, some data from a regular network
Here is
2012 Feb 10
0
range and anisotropy with RandomFields
Hello,
I am presently trying to get a feel for the various packages out there that
allow me to both analyze and simulate random fields. The package
RandomFields is nice, but there are still a few aspects of its
implementation that are confusing to me and I was hoping someone could help
clarify things for me. It could also be that my questions reflect a lack of
knowledge pertaining to random
2003 Apr 17
0
kriging in R
Hi
If you read the description of "varcov.spatial" you'll see that it is
used to *predict* a covariance matrix, based on the parameters of the
covariance function. So you don't need the oebserved data, you need
parameters for the covariance function.
Regards
EJ
On Thu, 2003-04-17 at 05:43, Yan Yu wrote:
> THanks for the suggestion.
> I have a Q when trying to use it..
2005 Dec 13
1
bug in geoR (?)
I've enconuntered this problem with the last cran version of geoR:
> library(geoR)
> day <- rep(1:2, each=5)
> coords <- matrix(rep(runif(10),2), 10, 2)
> data <- rnorm(10)
> data[1] <- NA
> as.geodata(cbind(coords, data, day), realisations=4)
as.geodata: 1 points removed due to NA in the data
Errore in as.geodata(cbind(coords, data, day), realisations = 4) :
2006 Mar 17
3
can't load "geoR"
Dear R list,
I've just installed R and then followed the instructions to install the
package "geoR" from within R by using the following line of code:
install.packages("geoR", contriburl = "http://www.est.ufpr.br/geoR/windows")
This installed okay under the folder C:\documents and settings\ ...
But when I follow the instructions in the illustrative session
2003 Sep 10
1
geoR variogram problem
Dear GeoR-er,
If I use the variog function in the latest release of geoR, the first lag is
always ignored.
For instance, if you read in geodata, calculates the variogram using the
variog function and give in a uvec like uvec=seq(0,max,by=2.44), it only
starts giving results from distance=4.88 and ignores 2.44!
This wasn't the case in former versions of geoR. Is this done on purpose?
Why?
2000 Oct 20
1
Covariace or variogram surface
Hi
I'm interested in doing a anisotropy analysis of some spatial data and I
want to do a contour plot of a covariance or variogram surface to
visualize the directions of maximal and minimal anisotropy.
I know how to do a contour plot :-) but how can I obtain the variogram
surface having a geoeas format dataset (X, Y, Z) ?
Thanks
EJ
2010 Aug 04
1
geoR - likfit failure
Hi
I'm using geoR package to perform linear spatial interpolation(OK).
The function likfit() fails to compute REML.
The error meassage is : Error in solve.default(v$varcov, xmat);
How I can find out that likfit() is failed to process and retrieving the error message ?
Thank you so much for your help.
2008 Aug 18
1
GeoR model.control - defining covariates at prediction locations
Hi,
Im using geoR and I'm trying to do some predictions, based on an external
trend.
I'm having some problems specifying my model.control, specifically how do I
define my model, and also the source of the covariate data at the prediction
locations?
I am assuming that the covariate data at the prediction locations should be
imported to a geodata object along with the prediction location
2010 Jan 07
1
kriging with geoR package
Hi all,
I have to draw a map of microorganisms repartition in a sample. I use the
"image" function of geoR package to build the map but there are three
problems:
1) the legend is Horizontal (with or without vert.leg=TRUE)
2) I want to plot the position of organisms and level curves of density on
the map but I don't find the script to insert these data on the map.
3)the grey range
2003 Feb 11
2
geoR question from new R user
Hi,
I'm a new R user. My goal is to do a variogram using geoR.
I started by trying to do the example in the geoR Illustrative Session
using my own data.
I am able to read in my Ascii data using: D <- matrix(scan("file.dat",
n=530*3), 530,3, byrow=TRUE).
Then I use: as.geodata(D, coords.col=1:2, data.col=3) to make the
object D geodata.
I check the descriptive statistics,
2010 Jun 09
0
Fixed sill in variogram fitting (geoR)
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
2003 Mar 25
1
geoR vector length error
Dear,
we have a problem with geoR. We try to read an ASCII table (x,y,z) with 40000
lines. With read.geodata we get the error:
Error in vector("double", length) : cannot allocate vector of length 799980000
We can read the file without any problem with read.table, but trying to
convert it to the geodata class gets the same error.
Any help/tips would be appreciated.
Kris, Dave
--
2004 Jan 08
1
Using geoR krige.conv
Hi,
I'm using the function krige.conv from package geoR but I'm getting a
warning that I do not understand:
Warning messages:
1: NaNs produced in: sqrt(variance)
2: NaNs produced in: sqrt(variance)
I'm performimg an ordinary kriging with a log transform (lambda=0).
Does anyone knows where this comes from ?
Thanks and regards
EJ
2003 Jul 24
1
geoR size limit problem
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
2011 Mar 25
1
spatial stats - geoR - variogram - standard deviation
Hello,
I am attempting to get the standard deviation in multiple distance bins in
my spatial data.
It appears as though the 'variog' command in the geoR package will do the
trick, as one of the outputs from 'variog' is 'variog$sd', which,
according to the manual, is the "standard deviation of the values in each
bin".
However, when I run this command, the
2004 Sep 16
1
geoR/variog4() not returning all directions
Mac OS 10.3.5, R 2.0.0
latest version of geoR
I have an incomplete 5 x 20 spatial array of samples (60 out of 100
possible locations) for which I would like to calculate directional
variograms using variog4(). Unfortunately, I can't get it to return all
4 directions. It returns variograms for 45, 90, and 135 degrees,
omitting 0 degrees (pi/4, pi/2, 3pi/4, omitting 0). If I specify 0
2008 Aug 15
1
Strange error message from geoR´s likfit () lik. max. func.
ComRades:
I am geeting the error message
Error in ldots[[which(MET)]] : attempt to select less than one element
when I try to fit the geostatistical model with the likfit() function of
geoR.
I have tried with old data for which likfit() successfully maximised the
likelihood in previous versions of geoR, and yet the current version
fails.
I have tried in Windows Vista and Windows XP (I haven't
2007 Jul 24
0
geoR/likfit: variance explained by model
I have been modelling spatial data using function likfit in package geoR. Now that I have the spatial regression models, I would like to calculate the amount of variance explained by both the trend and the spatial parts of the model. Any advice on how to do this would be much appreciated.
Cheers,
Ailsa