similar to: akima interpolation and triangulation question

Hello, dear R-community. This is a question about TukeyHSD between factor combinations of a Three-Way ANOVA, which is - since it is a multi measure ANOVA - not a simple ANOVA but a Generalised Linear Mixed Model (GLMM), calculated with "lmer". > growth <- groupedData(length~meas|box_id,outer=~spec*comp*water,data=all.spec) > model <-

Hi, I have a data set on germination and plant growth with the following variables: dataset=fm mass (response) sub (fixed effect) moist (fixed effect) pop (fixed effect) mum (random effect nested within population) iheight (covariate) plot (random effect- whole plot factor for split-plot design). I want to see if moist or sub interacts with mum for any of the pops, but I am getting an error

Okay, been working with tripack, seems the most mature package for this. Got it to work well with their test data set - data(tritest). When i tried random numbers to explore further, i am getting some results that don't reconcile. example run this: library(tripack) y <- runif(100) x <- runif(100) vm <- voronoi.mosaic(x,y) plot(vm) par(new=T) plot(x,y,col='blue') when

Hi, I'm currently hunting a bug in the akima library, especially in the code behind the interpp.old function (bi-variate linear interpolation). It is based on a triangulation algorithm, interpolation at a given point needs to know the triangle which contains this point, then the interpolation is a straightforward calculation based on the three vertexes. The problem is: Sometimes the triangle

maybe it's in the data? So here it comes. > sv.growth Grouped Data: length ~ meas | box_id meas spec comp water box_id sprouts leaves length long.sprout 1 1 Sv control moist 1 8.800000 37.80 211.2000 60.6 2 1 Sv xfull moist 2 7.000000 8.00 174.8000 62.8 3 1 Sv control moist 3 9.000000

Hi, I'd like to unsubscribe from the list. Thanks Neyra ________________________________ De: "r-help-request@r-project.org" <r-help-request@r-project.org> Para: r-help@r-project.org Enviado: sáb, octubre 30, 2010 5:30:07 AM Asunto: R-help Digest, Vol 92, Issue 31 Send R-help mailing list submissions to r-help@r-project.org To subscribe or unsubscribe via the

Dear R helpers, I try to use the 'delaunayn' function in the 'geometry' package for Delaunay triangulation in 2 dimensions. For the four following points, I get a warning message : > coord=matrix(ncol=2,byrow=TRUE,c(622633,7073452, + 621228,7073517, + 621879,7071762, +

Hello. I'd like to convert voronois object created by tripack to polygons to use them in GIS(Geographic Information Systems) software. I tried to create voronoi objects by using following code. > library(tripack) > x <- rnorm(10) > y <- rnorm(10) > plot(x,y) > v<-voronoi.mosaic(x,y) > plot(v) But from here, I could not create polygons. Of course, I know

Dear R users, I want to do a fitted.contour plot of selected columns of a dataframe M with M$AM and M$Irradiance as x and y axes respectively. The level of the contour shall be determined by M$PR. Some words on my data first. Dataframe M looks like: head(M$Irradiance) [1] 293 350 412 419 477 509 head(M$AM) [1] 2.407 2.161 1.964 1.805 1.673 1.563 head(M$PR) [1] 70.102 72.600 75.097 80.167

Dear, I have XYZ data available in a MySQL database. I get it out, can plot the data with the plot() function, load it into a geoR datastructure. But what I actually would like to do is a simple contouring of the data based on a no Kriging interpolation such as TIN based. I know the first thing I shold do is interpolate a full matrix for the region I have my points for, then contour should

> From: Bill Simpson <wsimpson at uwinnipeg.ca> > Subject: Re: contour labelling [was "Re: [R] image ()"] > MIME-Version: 1.0 > > > We may be missing the crux here: was it to have a means of > > displaying such surfaces? > > Yes, bang on. I knew of the > existence of the Delaunay code in the R libraries. Yes I want > to make a picture that

Just been trying to get the polygon triangulation code from this package: http://www-2.cs.cmu.edu/~quake/triangle.html to dyn.load into R. Uh oh. Floating point exceptions. Track it down to some FPU diddling that the author deems is necessary. Here's my minimal code that breaks: flipme.c: #include <fpu_control.h> void flipme(){ int cword; cword=4210; _FPU_SETCW(cword);

Full_Name: Michael R. Allen Version: 1.01 OS: linux redhat Submission from: (NULL) (149.149.32.46) Hello, Thought you might be interested. I tried to load the package tripack needed by the sgeostat package. After some careful searching, three packages were found which had the same error. Here is the output. Mike > library(acepack) Error in dyn.load(x, as.logical(local),

I am trying to triangulate a point set as follows: > head(cbind(x,y)) x y [1,] -78.1444 -60.4424 [2,] -78.1444 -58.4424 [3,] -78.1444 -56.4424 [4,] -78.1444 -54.4424 [5,] -76.1444 -60.4424 [6,] -76.1444 -58.4424 > length(x) [1] 5000 > tri <- tri.mesh(x, y) Fehler in tri.mesh(x, y) : error in trmesh > tri <- tri.mesh(x, y, "remove") Fehler in tri.mesh(x, y,

Hello, I'd like to make a z(x,y) plot for irregularly spaced x,y. What are routines are available in R for this purpose? Thanks, Mark

Stephen John Smoogen wrote: > On 1 March 2018 at 12:26, hw <hw at gc-24.de> wrote: >> Stephen John Smoogen wrote: >>> >>> On 1 March 2018 at 08:42, hw <hw at gc-24.de> wrote: >>> >>>> >>>> I didn?t say I want that, and I don?t know yet what I want. A captive >>>> portal may >>>> be nice, but I haven?t

This isn''t specific to RSpec, but is hopefully on-topic for this list. I like (especially when "ping pong pairing") to write a spec, then write the smallest amount of code I can to pass it (especially when "ping pong pairing"). Sometimes this means hard-coding a return value, which means another spec is needed to prove that the code is really behaving as it

I have a problem with including extra data in a lattice graphic. I am trying to put some calculated values in an extra column of a boxplot. A minimal example should show what I am trying to do: foo <- data.frame( Treatment=rnorm(1:12,2),

Hi. I've tried to create a polygon with one hole by gpclib using following example script. holepoly <- read.polyfile(system.file("poly-ex/hole-poly.txt", package ="gpclib"), nohole = FALSE) area.poly(holepoly) plot(holepoly,poly.args=list(col="red",border="blue")) And I noticed plot function couldn't draw polygons with holes correctly.

Hi R helpers, Is there a function (package?) in R available which tests "spatial autocorrelation" between points (e.g. vector layer of weather stations)? (e.g. Moran's I...) Via the archives we found out that there is a package 'spdep' which uses grid data for testing spatial autocorrelation. Thanks a lot, Jan