[R-sig-Geo] Conditional Simulation question
sokolmatt at gmail.com
Tue Feb 1 21:11:48 CET 2005
My name is Matthew Sokol and I am writing this list serve to ask some
help. Currently I am a master's student at the University of Nevada
and have come across a problem using R that I can not solve.
Here is my problem:
I am using the R package geoR, specifically the function grf to simulate data.
Basically I have a data set that contains an Easting and Northing
coordinate system and a Z value. From this system I wish to perform
Ordinary Kriging. So I run the following command lines in R:
Now since I know the Nugget, Sill, and Range for the variogram model,
I wish to simulate more data that follows the same variogram model. So
I run the following command line:
Now, the problem with this command is that it produces data on grid
that is 1 x 1 (all the Northing and easting coordinates are now a
number 0-1). But my easting and northing coordinates have the
Easting: Min 254.4 Max: 492
Northing: Min 118 Max 315
Which is not a 1 x 1 grid.
I know my problem is in this statement: grid="reg"
I tried the following commands:
So now I run:
Which produces so much data that I can not R console will not display it.
So I run the following R commands:
variog: computing omnidirectional variogram
Error in vector("double", length) : cannot allocate vector of length 1444504375
So now I get an error. I think i may be producing to much data with my
I am lost. I've looked through the manuals but have not found much on
conditional simulation. I was wondering if anyone knew of a problem
like this or where I can look to find more information.
I am sorry if the question is trivial as I have been reading the posts
on this list serve and everyone here has a command of the language far
surpass mine. Thank you for your time.
More information about the R-sig-Geo