similar to: polygons on scatterplot3d

I have been attempting to estimate the periodic contribution of an effect to some data but have not been able to fit a sine wave within R. It would be nice to start by being able to fit a sine wave with an amplitude and frequency. x<-seq(0,20,by=0.5) y<-2*sin(2*pi*.5*x) #amplitude =2, frequency=0.5 # This failed to converge r<-nls(y ~ A*sin(2*pi*F*x), start=list(A = 1, F = 1),

Greetings, I am attempting to run a function, which produces a vector and requires two input variables, across two nested factor levels. I can do this using by(X, list(factor1, factor2), function), however I haven't found a simple way to extract the list output into an organized vector form. I can do this using nested loops but it isn't exactly an optimal approach. Thank you

Hi R list, I would like to automate, or speed up the process from which I take several separate datasets, stored in .csv formate, import and merge them by a common variable. So far I have greatly sped up the loading process but cannot think of a way to automate the merger of all datasets into a common data.frame. My apologies if this has been covered, any R search suggestions are

Greetings, I have been working on a script that conducts repeated statistics and plots to my data. In this case it is sub-setting the dataframe by month. The intent is to develop a custom analysis and plotting that I can run on a large number of data sets. Unfortunately, a small portion of my subsets (~1%) cause an error with one of the wrapped subroutines that results in the whole

Hi R-Users, I have a time series of bivariate observations (x,y,t). I've developed a few panel routines to explore certain aspects of the data, and now would like to break the time series up into smaller chunks to explore the temporal trends. I would like to know if anyone has any experience breaking up time series with lattice. Base graphics offers coplot, utilizing the

Hi, I am having an issue with a custom panel for lattice. The problem comes when I try passing a groups argument. Here is the custom panel, a wrapper for smooth spline. I copied panel.loess and replaced the loess arguments with smooth.spline(). [Note: I would like to use the cross-validation fitting properties of smooth.spline.] library(lattice)

Hi list, Well, this time I’ve a doubt with mapping generation. I was already able to read and plot shapefiles, plot point on this map. All this without any problems, but now I want to do something that I think, only Golden Software Surfer is capable of. I would like to plot a kriging result on the existing map (map script below). Well I looked for on the net, but I was not able to be sure of

Deemed colleagues I would appreciate your help with a sentence to calculate the linear distance between 2 geographical points (coordinates in UTM). In advance thnks for your attention, -- Ricardo Bandin Llanos rbandin@udec.cl Estudiante - Magíster Cs. m. Pesquerías Universidad de Concepción, Región del Bio-Bio, Chile Celular: (0056-41) 97949957 [[alternative HTML version deleted]]

I have been crashing R while using lattice. My system consists of a Mac Intel, R 2.5.1, lattice 0.15-5, and I plot through the default quartz. R crashes after plotting several lattice functions. Has anyone else encountered this problem and does anyone know the cause, or a solution? Is this a lattice:mac, lattice:quartz, Lattice:me issue? example... This crashes my system after plotting.

[Env: Win XP] I've just upgraded from R 2.12.2 to R 2.13.1. As part of my upgrade process, I typically install some in-development packages from R-Forge that are not on cran. But for the first time, it doesn't work. e.g., > install.packages("p3d", repos="http://R-Forge.R-project.org") trying URL

Hi, I would like to put the mean and + / - the standard deviation as lines on the x axis of a histogram. My attempts using the histogram function have been unable to do this. My searches are unsucessful on this subject. Any ideas are appreciated. Thanks stemp <- 5 6 5 5 5 5 6 5 m <-5.25 stanD <- 0.46291005 using windows xp R 2.6 --------------------------------------------------

Hi, Can anyone guide me as to how I can add points to a p3d() plot from the onion package?? I want to plot points with different colors on the same 3D plot.? Perhaps I can do this without adding points but somehow directing the 'h' parameter to give different color to points based on a factor I assign to them? FYI, I can do this using using scatterplot3d() and points3d(), but these plots

Hello I've created a 3d scatterplot, and had no problems labeling the points. However, I've been really struggling to change the color of the points based on a factor (see 'group' below). Is such a thing possible? My data look like this: food group x y z apple fruit 0.216 -0.112 -0.893 orange fruit 0.814 0.097 0.460 broccoli veg -0.239 0.240 -0.425 banana

Hi R-help I am trying to mask the ocean from an image plot I have made. Here is some example code: library(mapdata) image(x=110:155, y =-40:-10, z = outer(1:45, 1:30, "+"), xlab = "lon", ylab = "lat") outline <- map("worldHires", plot=FALSE) # returns a list of x/y coords xrange <- range(outline$x, na.rm=TRUE) # get bounding box yrange

Dear Colleagues: I have 3 numeric variables (say x, y, z) in a data.frame and want to draw a 2D-scatterplot with the 3rd dimension as color, if possible without having to convert the whole thing to an appropiate matrix for image() manually. Can e.g. scatterplot3d do something like this? There's a function xyz.convert, but I couldn't find the way to use it properly. Sincerly,

Hi R-Help, I would like to make filled contour maps of ocean data overlaid by costlines from the map package. I can draw the filled contours and the coastlines om the same plot, but the filled contour also covers part of the land. To get rid of that I tried to draw a filled coastline map on top of the filled contour, but the filled map only draws the closed contours - so most of the land is

Hi All, Does anybody know of any existing functions that will calculate distance inland from a coastline? It's possible to test if a lon,lat location is land or sea using map.where(), but I need to add a buffer to this of say 2km, to allow for points that are just on the coast, and below the resolution of the worldHires database. I'm working with a marine mammal satellite telemetry

I am running R 2.5 on a Mac platform and have a difficulty loading the Tcl/Tk package. Specifically: R>library(tcltk) Loading Tcl/Tk interface ... Error in fun(...) : couldn't connect to display ":0" Error : .onLoad failed in 'loadNamespace' for 'tcltk' Error: package/namespace load failed for 'tcltk' Is there an explanation to this problem? Thank you

I have a 3xN matrix of parameters obtained from N regressions where the 3 parameters are jointly statistically significant. I would like to reproduce a 3D confidence ellipsoid projecting 2D ellipses onto the XY plane as in Figure 5.2 in this

Hi I'm using R Version 1.4.0 Patched (2002-01-16) in a SuSE 7.1 box. I've a plot with the european coastline. After I want to plot a single point in the map using points, but it plots this point in the wrong place. Why is this happening ? Is there a known problem of the points function ? Regards EJ > plot(eucoast,type="l") > points(c(-15.0,57.5)) ps: shall I send the