similar to: Need a hint

I am modeling a trend surface using trmat and want to trim the resulting matrix to the area enclosed by my real data (i.e., remove all the extrapolated areas). I was using chull and in.chull to calculate the convex hull and change all the other values created by trmat to NA. However, my real data has portions that are slightly concave so chull would give me slivers that are extrapolations from

I am using package ks() to build 3D representations of bird territories and calculate territory volume from spatial data (simply x, y, and z coordinates). What I want to do is determine at what sample size (# locations collected) does the territory volume stop increasing. This should give me an idea of the number of points needed for future seasons. So I have a couple of birds each with 200

Hi, Is there a way to plot more than one trend surface using the functions contourplot() and wireframe(). I have found an add=T in contour(), but no equivalent argument in contourplot() and wireframe()? I have taken the example 11-2 (pages 441-451) from Design and analysis of experiments (Montgomery 2001, 5th edition) to see if this could be done in R. I have managed to plot individual

Dear all, I am facing a strange problem with text files formats. I am currently using a C++ code named voro++ (http://math.lbl.gov/voro++/). This code accepts text files with four columns: <id> <x> <y> <z> where id is an identification number and x,y,z are the coordinates of a point in a 3D space. The input file, name it myfile_cpp.txt, is generated by another C++ code

Dear R developers Using the 'foreach' package I encounter warnings like Warning message: In identical(.combine, cbind) : ignoring non-pairlist attributes Warnings appear once in a new R-session when I run a particular piece of code - and not again. I understand from the source code of 'identical' that the function checks and believes that .combine and cbind have non-pairlist

Yesterday I downloaded R-2.2.0.tar.gz, gunzipped and untarred, and did the usual ./configure and make. Everything seemed to go smoothly until it got to the bit about installing recommended packages. It got past ``spatial'' but turned up a fatal error in respect of the ``boot'' package. Here is some of what appeared on the screen:

Dear list, since I have installed GRASS_0.1-8.tar.gz (devel) with R 1.4.1 on RedHat7.1 Linux I get the following message when I want to save the workspace image: -------------------------------------------------- Save workspace image? [y/n/c]: y Error in file(file, "wb") : cannot open file `.RDataTmp' ------------------------------------------------- Where can I find the

>>>>> "Llu?s" == <Lluis.Hurtado at uv.es> writes: Llu?s> Dear all, Llu?s> I am facing a strange problem with text files formats. Llu?s> I am currently using a C++ code named voro++ (http://math.lbl.gov/voro++/). This code accepts text files with four columns: Llu?s> <id> <x> <y> <z> Llu?s> where id is

Dear R users -- I think this question was asked before but there was no reply to it. I would appreciate any suggestion any of you might have. I am interested in plotting several "implicit functions" (F(x,y,z)=0) on the same fig. Is there anyone who has an example code of how to do this? Thank you Yihsu

Dear list, I get the following error message when installing the package GRASS. missing link(s): interp.new missing link(s): vgram Krig predict.Krig predict.se.Krig See also below for the full installation message. What does this mean exactly? Regards, Ulrich ------------------------------------------------------------------ > R CMD INSTALL

Hi All, I have a netcdf gridded file with LCC projection. I can easily use image.plot to visualize it. However, as the axises are in X,Y, not Lat and Lon, I could not add state or country maps onto it (or lat lon information). I do have a grid2d file that describes the lat and lon for each (X,Y) grid, but the lat and lon are not regularly spaced, so I could not use image.plot. Does anyone know

Hi, You might get more help from the R-sig-geo list, which is devoted to spatial topics. However. The *.asc file is an ArcGIS raster export format. You should use whatever the appropriate import commands are for your own gridded rainfall data. If you have a different format, you might or might not be able to import it directly with raster. ?raster will tell you more about the kinds of formats

Hello all, I am trying to use "clim.pact" package for my work, but since this is the beginning for me to use gridded datasets in "R", I am having some trouble. I want to do seasonal analyses like trends, anomalies, variograms, EOF and probably kriging too to downscale my 1 degree gridded data to 0.5.  So, as a first step, I compiled my entire dataset (with 25

Apologies in advance for the question. I am trying to draw a map of the US as a surface plot so that I would be able to drop bars on the different states (something like Uwe Ligges' scatterplot3d example 4). I am not sure where to start looking for such a beast. If anyone has any pointers, ideas, I will be grateful. TIA, Partha

I have been using the sp package to deal with gridded data. I would like to know how to extract the vertices (x,y) of polygons that outline areas of like-valued cells in a grid. Here is a simple 3x3 grid: 2 2 1 1 2 1 1 1 3 x <- c(1,1,1,2,2,2,3,3,3) # define a 3 x 3 array of points y <- c(1,2,3,1,2,3,1,2,3) h <- SpatialPoints(cbind(x,y)) # make these an sp object of the points

Hi! How does one change the background color of the map-panel in spplot()? Example: library(sp) data(meuse.grid) gridded(meuse.grid) = ~x+y spplot(meuse.grid, "part.a") How would I get another background-color for the map-panel (but not for the whole plot) here? Thank you! Marcel

Greetings, I am attempting to build a R package, however the build fails and the following message is ouput: C:\ConstellaGroup\EPA\RAGG\package>R CMD BUILD --binary --force RAGG * checking for file 'RAGG/DESCRIPTION' ... OK * preparing 'RAGG': * checking DESCRIPTION meta-information ... ERROR Error in if (regexpr(dep_regexp, dep) == -1) { : missing value where

Dear All I hope that someone can help. I am working with sp pakage and akima library("akima") library(sp) imagine lots of different dataframes, of row = 100 columns = 3 of x and y coordinates with z values I will call these data frames for the sake of this example akima akima<-as.list(1:100) producing 100 dataframes dataframes of the form akima[[i]] I then wish to interp this

Hi, I am trying to recreate the right graph on page 524 of Gelman's 2006 paper "Prior distributions for variance parameters in hierarchical models" in Bayesian Analysis, 3, 515-533. I am only interested, however, in recreating the portion of the graph for the overlain prior density for the half-Cauchy with scale 25 and not the posterior distribution. However, when I try:

Friends: I wish to overlay data points on a contour graph. The following example produces a nice contour plot, but I have not mastered the concept of using panel functions to modify plots. Can someone show me how to overlay the data points (given after contour plot statement) on the contourplot? --Seth model <- function(a,b,c,X1,X2) # provide model function for contour