search for: maxdistance

Dear all, I'm new on R and this mailing list. We work on spatial rainfall estimation with R and Grass. We have a problem with the krige function from the sgeostat package. We would like to limit the estimated area with a polygon limit. I use a 50 points polygon to describe my work area. The krige function work quiet well without limit. But if I use this option I have the following error

...forth in any number of Clevelend and Devlin articles, ie, nearest neighbor smoothing. However, I'm interested in looking at fixed bandwidth smoothing and I'm wondering if the second option, (alpha>1) is some variation of this theme. I don't understand the "alpha times actual maxdistance" statement. If my estimation point is 10 meters from the farthest data point in physical space and I set my alpha = 2, then according to the statement above the smoothing window would be 20m???? Obviously this can't be right. Where have I gone wrong? Thanks in advance for your comments...

...forth in any number of Clevelend and Devlin articles, ie, nearest neighbor smoothing. However, I'm interested in looking at fixed bandwidth smoothing and I'm wondering if the second option, (alpha>1) is some variation of this theme. I don't understand the "alpha times actual maxdistance" statement. If my estimation point is 10 meters from the farthest data point in physical space and I set my alpha = 2, then according to the statement above the smoothing window would be 20m???? Obviously this can't be right. Where have I gone wrong? Thanks in advance for your comments...

Hi, I would know how to plot two variograms on a same graph. I can plot one by one but I would draw both on the same one. Is it possible? Do i need any special package? Thanks! Cordialement Damien Landais

Hola, # Hacemos el KED. Ver función "krige()": KED.rad <- krige( formula=pluvPcp~layer, # covariable -> radar locations=lluvia.rad.pluv.spdf, newdata=radarGrid, # podría ser cualquier objeto Spatial model=v.fit, # modelo de semivariograma. maxdist=Inf

...oked at the first one, and for the strings "UndeclaredSoFar" and "_Imaginary" and a maxEditDistance of 5 (= the difference in their lengths, if that matters), the algorithm returns an edit distance of 2. One of the other tests (Parser/objc-quirks.m) is actually a crash. + // maxDistance (aka D) directly influences the width of the strip (2*D+1) + // we arbitrarily restrict it in case the caller did not see fit to do it + // (from clang\lib\Sema\SemaLookup.cpp in TypoCorrectionConsumer::FoundName) + unsigned const maxDistance = std::min(unsigned((otherLength+2)/3), +...

Dear all, I got point data of trees. I was wondering if anybody has experience in searching the neighbors within a specified distance efficiently. X Y Z 99 34 65 98 35 29 98 34 28 99 33 33 98 32 23 99 33 21 99 33 22 99 32 24 99 30 23 ... What I want to do is : searching for the neighbors with a distance R for each tree & the neighbor must have a bigger Z.

Dear all, I would like to estimate rainfall with the krige function. First, I produce a polygon of my region of interest (where poly_test.txt is a x,y suite of points defining a polygon, obtained with the grass v.out.ascii command) : user>border_limite<-read.table("/home/lionel/rwork/poly_test.txt",header=FALSE) user>polygone<-list(x=border_limite[,1],y=border_limite[,2])

  Hola,   Por favor, ¿podríais indicarme qué recursos (librerías o ideas) pueden resultar de utilidad para crear un gráfico del estilo del de la figura 3.8 del siguiente link?   http://www.tsc.uvigo.es/BIO/Bioing/ChrLDoc3.html#3.5   Actualmente estoy utilizando funciones muy básicas y la verdad es que no me encuentro muy satisfecha con el resultado.   Muchas gracias.   Eva [[alternative HTML

ear ALL, The last part of my thesis analysis is the cross validation. Right now I am having difficulty using the cross validation of gstat. Below are my commands with the tsport_ace as the variable: nfold <- 3 part <- sample(1:nfold, 69, replace = TRUE) sel <- (part != 1) m.model <- x2[sel, ] m.valid <- x2[-sel, ] t<- fit.variogram(v,vgm(0.0437, "Exp", 26, 0)) cv69

Sale plano sí. Ya se que sin tener los datos y el código es un poco difícil, pero es que mis datos ocupan mucho, es imposible. Seguiré mirando por internet. Muchas gracias Rubén. Un saludo, > To: r-help-es en r-project.org > From: rubenfcasal en gmail.com > Date: Thu, 6 Aug 2015 14:21:47 +0200 > Subject: Re: [R-es] Duda interpolación (package ' gstat ') > > Hola

Hello Doug, *putting llvmdev in copy since they are concerned too* I've finally got around to finish a working implementation of the typical Levenshtein Distance with the diagonal optimization. I've tested it against the original llvm implementation and checked it on a set of ~18k by randomly generating a variation of each word and checking that both implementations would return the

Hi, Well after trying to fix the assign commands in version 3.1.7 of your nlme package I still cannot get the plot functions to work. Basically I have no other packages installed (apart from those provided with R) so nlme is the first. my R version info is; platform sparc-sun-solaris2.7 arch sparc os solaris2.7 system sparc, solaris2.7 status

Hi the list, I am writing a package in S4 and I do not manage to understand a bug. The "R CMD check" and the "R CMD build" both work. Here is links to the package (not on CRAN yet for the raison that I explain bellow): http://christophe.genolini.free.fr/aTelecharger/kml_0.5.zip http://christophe.genolini.free.fr/aTelecharger/kml_0.5.tar.gz Then I install the package and I

Dear Aaron, I am really a tool user and not a tool maker (actually an ecologist doing some biostatistics)... so, I take the liberty of sending a copy of this e-mail to the r-help list where capable computer persons and true statisticians may provide more relevant information and also to Paulo Ribeiro and Peter Diggle, the authors of geoR.. I really feel that your huge matrix cannot be

I am trying to use the nlme package in R following the book "Mixed-Effects Models in S and S-Plus". I have the latest versions of R (1.6.1) and nlme (3.1-33) for Windows (95 or later). I had problems with the instruction Variogram using a gls class object (Section 5.4, pg 261-262). Even when I use the scripts included in the package, I get an error message and not the expected

HI, I wonder does anyone have experience with doing sequential gaussian simulation with krige() function in gstat? I find it VERY slow compared to use krige() to achieve kriging function itself.. I wonder why, is that because it has to model the variogram, and do the kriging separately for each point to be simulated? so it would be N times slower to achieve the simulation than the kriging if

Hi all, >str(df) `data.frame': 31837 obs. of 3 variables: $ x : num 410683 410700 410720 410740 410324 ... $ y : num 43136 43126 43123 43125 42709 ... $ wz: num -101.1 -94.9 -93.3 -94.5 30.8 ... >library(gstat) >g<-gstat(id="rv",form=wz~1,loc=~x+y,data=df,model=mat,nmax=500,set=list(average=1)) >str(g) List of 3 $ data :List of 1 ..$ rv:List of 10 ..

Hi, I am out of town and will get back to you on the 13th of July. Leo >>> "r-help at stat.math.ethz.ch" 06/27/03 00:32 >>> Send R-help mailing list submissions to r-help at stat.math.ethz.ch To subscribe or unsubscribe via the World Wide Web, visit https://www.stat.math.ethz.ch/mailman/listinfo/r-help or, via email, send a message with subject or body