similar to: Conditionals for Errors or error flags

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),

Does anyone have any experience plotting complex polygons on a scatterplot3d display? I would like to present some spatially binned data in a 3-d type plot with background polygons (i.e.. a filled coastline map). I can get the coastline onto the plot window as a line type but cannot figure out how to fill the shape. I'm certain it may be possible but might require a bit of digging

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

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)

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]]

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

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.

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 --------------------------------------------------

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 am 100% positive there is a better fix for this problem and it has been discussed. Thank you in advance for bringing it to my attention. When upgrading to newer versions of R the old packages are not brought along in the process. I solved that problem by creating a list of the packages in the old version of R, importing them to the new version via a CSV, and re-installing into the newer R

Is it possible to place maps onto lattice plots? With basic plotting you can add a map to a plot library(lattice) long<-c(-69.2, -69.5, -70.1, -70.3) lat<-c(41, 41.5, 43.2, 42.8) plot(long, lat) map('state', c("massachusetts"),add=TRUE) but is it possible with lattice? library(lattice) factor<-c(1,1,2,2) xyplot(lat~long|fact) ...now what? I have looked at panel and

Hi -I am using R v 2.13.0. I am trying to use the aggregate function to calculate the percent at length for each Trip_id and CommonName. Here is a small subset of the data. Trip_id Vessel CommonName Length Count 1 230 Sunlight Shad,American 19 1 2 230 Sunlight Shad,American 20 1 3 230 Sunlight Shad,American 21

Hi - I have a dataframe with 4 variables - long, lat, temp and year. I would like to plot temperature by lat and long in contour plots by year. There are 391 data points. When I try contourplot(temp~long*lat|year,data=hake,cuts=10,labels = TRUE,contour = TRUE,pretty = TRUE,region=TRUE) I get plots with dots around by data, but not contour lines or filled in regions. I have also used the

Hi This is Gaurav Khanna, from MacResearch.org - a community for scientists that use the Mac and related hardware for research. One service we perform for the community is supply g77 & gfortran compilers for Mac OS X (which are not supplied by Apple). We are interested in the progress made by LLVM and have some quick questions. We'd appreciate any comments you can make. 1) What is the

Sorry the first email had no subject. I've searched around and havent had much luck trying to explain why my HTML search no longer works. If I use ?cor.test it is unable to display the webpage. Is this a java problem? I am running windows xp, and use Internet explorer for my default browser. Thanks -------------------------------------------------- David Kowalske Jr Graduate Research

Hi - I am using R version 2.13.0. I have run several GLMMs using the glmmPQL function to model the proportion of fish caught in one net to the total caught in both nets by length. I started with a polynomial regression full model with three length terms: l, l^2, and l^3 (l=length). The length terms and intercept were the fixed effects and the random effect was a paired haul (n=18).

Dear all,    I am having a recurring problem when I attempt to conduct a GLM.  Here is what I am attempting (with fake data): First, I created a txt file, changed the directory in R (to the proper folder containing the file) and loaded the file: #avoid<-read.table("avoid.txt",header=TRUE);avoid #  treatment feeding avoid noavoid #1   control  nofeed     1     357 #2  

in the Rails Recipes book there is a chapter on role-based authentication. it basically set up priviliges based on what page the user was able to access. i went through the chapter and felt like i got a pretty good understanding of how to do it. now i am needing something a little more advanced though. instead of a per page (or action) i would like to have certain sections on a page either