similar to: RandomFields: Controling seed with GaussRF

Not really a bug (I eventually discovered), but ... When running GaussRF()/DoSimulateRF() with n>1, the default pch="!" prints backspaces to the terminal, which causes problems (LaTeX barfs) if you're running this within an Sweave() job. Setting RFparameters(pch="") fixes the problem. Could put this in the documentation, but I don't know where you could

Hi, Please see the email exchanges below. I am having trouble generating output that is large enough for our needs, specifically when using the GaussRF function. However, when I wrote Dr. Schlather (the author of the GaussRF function), he indicated that there is also a limit imposed by R itself. Is this something that we can overcome? Thank you very much for any assistance you may provde.

Hi R-Users, I am ussing GaussRf in the Random Fields Package to generate spatial-temporal fields. I am having a hard time understanding how to use the NSST model. Does anybody have some experience with this model or any of the non-separable space time models available? I got all the papers cited in the manual but I am still lost. Thanks a lot Ele

Dear all, In my 'simple' computer I was running some experiments to help me understand how faster a multicore lapply will be. I thought it might be interesting for some people to look at the results. Even though are not accurate, still might be a good indicator how much improvement there can be. A.Case. The classic: for 1:100 for (i in c(1:dimz)){ print(sprintf('Creating the %d

Dear all, I would like to ask your help understand the subsequent steps for making my program faster. The following code: Gauslist<-array(data=NA,dim=c(dimx,dimy,dimz)) for (i in c(1:dimz)){ print(sprintf('Creating the %d map',i)); Gauslist[,,i]<-f <- GaussRF(x=x, y=y, model=model, grid=TRUE,param=c(mean,variance,nugget,scale,Whit.alpha)) } creates 100 GaussMaps (each

I have asked the same question on stackoverflow but did not get a satisfying answer. I am trying to simulate a lognormal spatial random field but I need the simulated value in a certain range. So I need some easy to use functions to generate a truncated Gaussian field to start with. To be specific, I need a function like GaussRF from the RandomFields package or grf from the geoR package to

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

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

Dear all. I am using the mclapply function to split my code to the many cores my system has. It seems that is working fine. This is the parallel version of lcapply. The only problem that I seem to have is that the printf cannot print messages. The ideal to me is to have fro my function an output of the form Shadowlist<-mclapply(1:dimz, function(i) { print(sprintf('Creating the

I don't understand what is happening! I have a (large) object sim1, an matrix list with dim c(101,101) where each element is an 3*3 matrix. I am subsetting that with a matrix coo, of dim c(100,2), of unique indices, but the resulting object has length 99, not 100 as expected. Code reproducing the problem follows: library(RandomFields) set.seed(123) sim0 <- GaussRF(x=seq(0, 100, by=1),

Hi R people: I am new to R. I am writing a function to (1) produce a sparse stochastic Gaussian 2D field and (2) perform a systematic transect sampling on this field, this carried out many times in a simulation framework. My function does a good job at producing the random field (a matrix of zeros and some manifestations of the stochastic process, depending on a parameter of the function

Hello everyone. I would kindly request your help concerning how R converts data between different structrures. In the following example please keep attention on the following two 1) I create f <- GaussRF(x=x, y=y, model=model, grid=TRUE,param=c(mean, variance, nugget, scale, alpha)) with image(x,y,f) and image(f) I get exactly the same image. then 2)I try to store f inside a raster

Hi All, I am new in R and trying to simulate random normal 2D field with mean trend say north-south. My domain is 10x10 grid and I am trying to use mvnorm but do not know how to specify the domain and the mean field. I would appreciate any help. Cheers, SK --------------------------------- [[alternative HTML version deleted]]

Hello, I am a graduate student of University of Florida. I am trying to run a process variation model and I am facing the following error. ------------------------------------------------------------- Analysis of geostatistical data For an Introduction to geoR go to http://www.leg.ufpr.br/geoR geoR version 1.6-27 (built on 2009-10-15) is now loaded

Hi, I'm having some difficulty installing the RandomFields package in R 1.7.0 under FreeBSD. I've installed the geoR package with no difficulties. Following are the errors that I get: * Installing *source* package 'RandomFields' ... ** libs c++ -I/usr/local/lib/R/include -I/usr/local/include -mieee-fp -fPIC -O -pipe -c MPP.cc -o MPP.o c++ -I/usr/local/lib/R/include

The package list was accessed October 28, 10:13 am MST using Windows XP and mirror USA (CA 1). Bill Morphet, Ph.D.

I have been unable to install the package RandomFields. I am using R 2.9.0-4 on Ubuntu 9.04. To install, I use the command: sudo R CMD INSTALL RandomFields_1.3.37.tar.gz The output follows below. Any help Would be appreciated. D. Hoysak Brandon University * Installing to library ?/usr/local/lib/R/site-library? * Installing *source* package ?RandomFields? ... ** libs g++

I have been unable to install the package RandomFields. I am using R 2.9.0-4 on Ubuntu 9.04. To install, I use the command: sudo R CMD INSTALL RandomFields_1.3.37.tar.gz The output follows below. Any help Would be appreciated. D. Hoysak Brandon University * Installing to library ?/usr/local/lib/R/site-library? * Installing *source* package ?RandomFields? ... ** libs g++

Hi All, Here is my sample data set.. y x 7/4/2009 -0.2368 -1.2727 7/11/2009 -0.5039 -5.2805 7/18/2009 -0.6655 -6.9641 7/25/2009 -0.3936 -3.6937 8/1/2009 -0.3463 -5.6457 8/8/2009 -0.3000 -1.7368 8/15/2009 0.2378 6.4600 8/22/2009 -0.2962 -3.1113 8/29/2009 -0.4346 -4.2039 9/5/2009 -0.6971 -7.8216 9/12/2009 -0.1217 5.1446 9/19/2009 -0.3107 2.0862 9/26/2009 -0.1797 -3.6055 10/3/2009 0.2299 -0.7373

Hello, I have some trouble using step() and stepAIC() functions. I'm predicting recruitment against several factors for different plant species using a negative binomial glm. Sometimes, summary(step(model)) or summary(stepAIC(model) does not select the best model (lowest AIC) but just stops before. For some species, step() works and stepAIC don't and in others, it's the opposite.