search for: grf

Hello everyone. I am trying to build up understanding in R by trying to develop just some simple scenarios. I would like to explain you what I am trying to do and what I did so far. I would like to put inside a RasterLayer (raster package) a Gaussian field (for given covariance) using grf function (geoR package) 1. First I created a Raster Layer object r <- raster() # Default values are ok 2. Then I set some values to test how if setValues worked r <- setValues(r,1:ncell(r)/100) # Every cell of the RasterLayer takes as data its cell number/100. 3. Tested that 2 works with...

...ng 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 ------------------------------------------------------------- grf: generating grid 300 * 300 with 90000 points grf: process with 1 covariance structure(s) grf: nugget effect is: tausq= 0 grf: covariance model 1 is: spherical(sigmasq=1, phi=0.1) grf: simulation using the function GaussRF from package RandomFields Error in DoSimulateRF(n = n, reg = register...

hello all, Quick question, how do I get the p value out of the anova? Thanks, Paul > pb<-aov(as.numeric(diff[5,16:33]) ~ grF) > summary(pb) Df Sum Sq Mean Sq F value Pr(>F) grF 3 2.7860e+10 9.2867e+09 4.2236 0.02534 * Residuals 14 3.0783e+10 2.1988e+09 --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 > str(summary(pb)) List...

...> path [1] "D:/rw0641/own/own/procedures/GFR_TIME" > path1 <- paste(path,"gfr.txt",sep="/") > path1 [1] "D:/rw0641/own/own/procedures/GFR_TIME/gfr.txt" -- this is good enough > data <- read.table("D:/rw0641/own/own/procedures/GFR_TIME/grf.txt",h=T) > dat1 <- read.table(path1,h=T) Error: "scan" can't open file but doesn't work > path1 [1] "D:/rw0641/own/own/procedures/GFR_TIME/gfr.txt" although the path is OK - the unpasted variant is of course ok. What happened ? > data day month year...

Been on linux (ubuntu 9.10) for about a month now. Starting to get the hang of it. The only thing from windows I am missing is being able to play this game. I've read a lot of stuff on the net and tried different things. Some say i need to trick the ip but I am unable to figure out the server IP address. A lot of private servers have certain things in place to prevent hacking, so I think

Hi all, I would like to test the relevance of a vector field (i.e. if the vectors are organized or not). To do so, I would like to use an autocorrelation test, so that I have two questions: - is the Watson test applicable to that perpuse ? - is the kuiper test applicable to that purpuse ? Regards ------------------------------------------------ Emmanuel Poizot Cnam/Intechmer B.P. 324 50103

Dear R People: I know that this is COMPLETELY off topic. Does anyone know how to put colors into a LaTex document, please? This is LaTex via WinEdt. I can get colors via PcTex but not in WinEdt. Thanks in advance! Sincerely, Erin Hodgess Associate Professor Department of Computer and Mathematical Sciences University of Houston - Downtown mailto: hodgess at gator.uhd.edu

Hello all, I am trying to create a grid of large number of points 4096*4096, on processing the data I am writing it into a file. phi <- 0.5 N <- 4096 mu <- 90 sim<-grf(N*N,grid="reg",cov.model="spherical",cov.pars=c(1,phi),method="RF") sim$data <- (sim$data - mean(sim$data)) / sd(sim$data) local <- array(sim$data, dim=c(N,N)) sysvar <- local * sigma P <- mu+sysvar write.table(P,"sysv",row.names=F,col.names=F)...

...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 generate a random field, but I need the generated field to have a truncated marginal distributions and a correlation structure with a prescribed range. Is there an R package or functions which can do this? If there is no availabe read-to-use functions or packages,is it poss...

...s, and call the funtion "likfit" (in package geoR) from within a for-loop. "likfit" does again call "optim". After ssuccessfully fitting variograms to some of the data sets, the procedure crashes and I get the error message: Error in optim(par = ini, fn = negloglik.GRF, method = "L-BFGS-B", lower = lower.optim, : non-finite value supplied by optim When I restart the procedure with the data set that first caused the error, it is usually possible to fit the variogram successfully for this data set, but the procedure will crash when fitting t...

dear list I am student M.S. statistics in department statistics . I am working in the function "nls" in the [R 2.3.1] with 246 data and want to fit the "exp" model to vectors( v and u ) but I have a problem to use it u 5.000000e-13 2.179057e+03 6.537171e+03 1.089529e+04 1.525340e+04 1.961151e+04 2.396963e+04 2.832774e+04 3.268586e+04 3.704397e+04 4.140209e+04

Hiya, I've posted this 'enigma' on the newsgroup : http://groups-beta.google.com/group/comp.emulators.ms-windows.wine/browse_thread/thread/57dcb7531ebfb529/65018c77a5e58651#65018c77a5e58651 The answers were pretty rare so I post the problem here, too. " I'm getting connection problems with an application run under Wine which isn't acting like that under a native Windows

...E3Igl1wILPvq lOQR1wDEnmcy15I1IQy1j+jxZUaq3T1dN0um36/4igGL2Uo8FQdDijlITaUDP2sDFwNDVYJvfDvL dC5QLnURagOpHQhiotBEDAOxxm3DEivB/A9cyLu3Qk7CE8vrrCgFaBxMZe7gIqC3C5Jyax/hfCOd fRxoFCYfiBMj1ABYrFPrOoG6uQYJBkZT8MaaoWG0xbU04N0BeRYi3IBc/xB0NASOJO2w/wkhGJsv V7Ms8CBAIInlU1kQoToPyQzc/MZy8C1RHVcgxVMpIFmL8EvmneiYDUmhgrFqEoLTGf397BCzw+u0 UBa3owVOFjjZHww/b31zA3ocMIMT4HQrahlkCrJALtwI2UBsG/DeStwiheMaa6/tByD0k95LADsx 7DHp9HUYunWdwES5qtPkeAq4Noxo8CcYEXLUKAduTnhIS5MQztt/ip2kRzCBULbzVybpRcz12R4+ JE2zJPXJdLhWQG1J8CBxeVsZebJdxX4I5vR0IvBe1tAN4/JkU+R0GLkEfYAwyVMBiJn3J0GtWQR5 6/f/JfxXxQHrNgUEO8wAUT37fiscm6YIchSB6QstBKG3P/aFARd...

...E9vj 4uIlmjjwlXCczp4bp4Eb4m5Vcbcp0+LqSQ1xOr44Tm2FZRE0PmBKKRuuzxKy6HUApx5gxo+ThTG8 DT8lZJIM4AKGqmn8gm2gwFu296ggEhO4wHB7dSnDAutPSSurN3GFX0YMXfx9PhK38JGB1LSyvfVf 5udpP+AoQODnUd7h79AWfg4BtwLXAfYCjgJOAS4D9PwcylmUM/wMWfjbVAqoAywF7AUcBVwGGPjb wFb+J3X5NKzydQDO/wRs5X/EtP4IbOFvgXuLv6UM89/0l1eWHdQYX2mCcRUmmGxHgrFllcX4G/2f THTF+LsDks/1VP1k/ibJAKw4sBUgAVoAbYAegB7caXCnKQJ4GPAUQAboUec06pxGnROAVwGnaTLA D2gBGPnr/egmxk/1extc9Vm46v0Cn11c/CRXP9W4+Kv8ZY2+wn+u0V+B5kF/gr/cn+ei+mTYCXWs oFbQUth1/MWBAptLqU/jRxEkF3ApoA4wB7AUsBOg50d5fv9ylw2NHKYT2I0u3k/vafRH9H0j+Ve6 /N5G5JikIm/Vl8AB7ZX2ernfu2s3RBV5H3oEnIq89+wApyLvXVvBqci7agM...