# Conditional Simulation question Classic List Threaded 2 messages Reply | Threaded
Open this post in threaded view
|

## Conditional Simulation question

 Hello, 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: > dat<-read.table("examplepb.dat",header=TRUE) > geodat<-as.geodata(dat) > plot(geodat) > vario<-variog(geodat) > lines.variomodel(cov.model="gaussian",nugget=1800,cov.pars=c(1400,60),max.dist=150) 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: > simdat<-grf(625,grid="reg",nugget=1800,cov.pars=c(1400,60)) 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 following properties: 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: Loci<-expand.grid(seq(250,500,l=250),seq(115,320,l=205)) So now I run: > simdat<-grf(625,grid=loci,nugget=1800,cov.pars=c(1400,60)) Which produces so much data that I can not R console will not display it. So I run the following R commands: geosim<-as.geodata(simdat) > variosim<-variog(geosim) 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 simulation? 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. -Matthew Sokol
Reply | Threaded
Open this post in threaded view
|

## Conditional Simulation question

 Matthew 1. the "grid" argument can also taka a matrix with coordinates values. Terefore if you want N locations in coordinates which are not in the 1x1 square just generate a Nx2 matrix witht the coordinates and pass it to "grid" 2. grf produces "uncoditional" simulations. To obtain conditional simulations you will need to use either krige.conv (for "tradidional" kriging algorithms)  or krige.bayes(for "Bayesian" kriging.) See for instance help(krige.conv) hope this helps P.J. On Tue, 1 Feb 2005, Matthew Sokol wrote: > Hello, > > 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: > > dat<-read.table("examplepb.dat",header=TRUE) > > geodat<-as.geodata(dat) > > plot(geodat) > > vario<-variog(geodat) > > lines.variomodel(cov.model="gaussian",nugget=1800,cov.pars=c(1400,60),max.dist=150) > > 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: > > > simdat<-grf(625,grid="reg",nugget=1800,cov.pars=c(1400,60)) > > 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 > following properties: > > 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: > Loci<-expand.grid(seq(250,500,l=250),seq(115,320,l=205)) > So now I run: > > simdat<-grf(625,grid=loci,nugget=1800,cov.pars=c(1400,60)) > > Which produces so much data that I can not R console will not display it. > > So I run the following R commands: > geosim<-as.geodata(simdat) > > variosim<-variog(geosim) > 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 > simulation? > > 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. > > > -Matthew Sokol > > _______________________________________________ > R-sig-Geo mailing list > R-sig-Geo at stat.math.ethz.ch > https://stat.ethz.ch/mailman/listinfo/r-sig-geo> > Paulo Justiniano Ribeiro Jr LEG (Laborat?rio de Estat?stica e Geoinforma??o) Departamento de Estat?stica Universidade Federal do Paran? Caixa Postal 19.081 CEP 81.531-990 Curitiba, PR  -  Brasil Tel: (+55) 41 361 3573 Fax: (+55) 41 361 3141 e-mail: paulojus at est.ufpr.br http://www.est.ufpr.br/~paulojus