similar to: k-d tree for loess

Dear R listers, When using smooth.spline to interpolate data, results are generally good. However, some cases produce totally unreasonable results. The data are values of pressure, temperature, and salinity from a probe that is lowered into the ocean, and the objective is to interpolate temperature and salinity to specified pressures. While smooth.spline provides excellent values at the

The University of Miami's Rosenstiel School of Marine and Atmospheric Sciences (RSMAS) is seeking a research associate to work with hydrographic data. The objective of the project is to develop statistical models for estimating salinity and density from observations of temperature together with knowledge of location and time. Because the relationship between salinity and temperature can be

The University of Miami's Rosenstiel School of Marine and Atmospheric Sciences (RSMAS) is seeking a research associate to work with hydrographic data. The objective of the project is to develop statistical models for estimating salinity and density from observations of temperature together with knowledge of location and time. Because the relationship between salinity and temperature can be

Can someone point me to the appropriate functions for fitting multiple statistical models that are coupled to each other. The data are measurements of salinity s and temperature t at stations id and pressures at p as well as surface elevations h at stations id. The problem is, for any new station, to estimate s at all p, given t at all p and given h. If h is ignored, the for each p, there

R 1.8.0 on Windows XP Professional. A huge THANK YOU to the R Team for this marvelous software. I am making lattice plots of oceanographic data. The usual layout does not conform to plotting conventions that marine scientists use when depth is the independent variable. Under those conventions, plots are made with the origin at the upper left, depth on the vertical axis (increasing as it

Dear R experts, When using R 2.1 under Red Hat Enterprise Linux 4 with Mozilla 1.7.6 as browser, help.start() almost works correctly. When I follow one of the search result links everything is fine, but when I go back to follow another link, I find them all inactive. So I have to redo the search. Any suggestions? Thanks, Carlisle Thacker > version _ platform

Dear R experts, I would greatly appreciate your help with R's hypertext online documentation, which I can pass on to our system manager. We are running version 1.7 on sparc-sun-solaris2.9 using netscape7 browser for displaying the help. The system is working, but there are a few frustrating details that I hope might be fixed. Here are three problems we are having: 1. When the R session

Colleagues, I am analyzing data collected during oceanographic cruises. We have conducted many cruises over the last decade. On each cruise we visit ~50 stations. At each station (termed EventNum)we lower an instrument that measures depth, temperature, salinity and oxygen every few seconds as it is lowered through the water. Data from all EventNums on all cruises are stacked, generating a data

Hi all, I am a beginner trying to use R to work with large amounts of oceanographic data, and I find that computations can be VERY slow. In particular, computational speed seems to depend strongly on the number and size of the objects that are loaded (when R starts up). The same computations are significantly faster when all but the essential objects are removed. I am running R on a machine

Dear R-listers, It seems that predict() behaves differently within panel.xyplot. Am I doing something stupid? Thanks, Carlisle First, without xyplot(): > lmtest <- lm(t~s,data=subset(P100,whichLon100==1 & whichLat100==1)) > lmtest Call: lm(formula = s ~ t, data = subset(P100, whichLon100 == 1 & whichLat100 == 1)) Coefficients: (Intercept) t 33.3307

How to combine three plots so that they partially overlap? As the data on all three are near the left side and top, it would be nice to draw the first, shift axes down and to the right and draw the second, and shift again to draw the third. I could shift the data for the 2nd and 3rd plots and construct their axes by hand, but it would be great if R has an easier way.

Dear Rlisters, When trying to indicate which data belong to which of 7 groups, the following plot shows only 5 groups: ss <- trellis.par.get("superpose.symbol") ss$col <- c("red","blue","green","yellow","orange","magenta","cyan") ss$pch <- format(1:7) ss$cex <- 1.2 ss <-

How to add units (e.g. "cm") to the color key of a lattice levelplot? The plots looks fantastic, but it would be nice to indicate somewhere near the end of the color key that the values associated with its colors are in centimeters or some other physical units. The only thing I find is the possibility to specify the labels so that one explicitly includes the units. That leaves

Hi there, I'm using the lattice package to create an xy plot of abundance vs. depth for 5 stages of barnacle larvae from 5 species. Each panel of the plot represents a different stage, while different loess smoothers within each panel should represent different species. However, I would like depth to be on the y-axis and abundance to be on the x-axis, because this is more intuitive as an

Hello, I'm doing the following: library(ncdf) library(fields) library(animation) saline <- open.ncdf("salinity_1990.nc") salt = get.var.ncdf(nc=saline, varid="Salinity") # create an animation of the complete temporal domain in the ncdf file. saveHTML({ for (i in 1:364) { image.plot(salt[, , i]) } }, img.name = "salinity.img",

The NOAA National Centers for Coastal Ocean Science is hiring a Computational Ecologist, a statistical/computational ecologist with experience fitting advanced spatial models to marine wildlife survey data (e.g., seabirds and marine mammal transects, fisheries trawl surveys) in R and other statistical languages. This is a full-time, long-term stable contract position. We are looking for an

My x axis is a series of daily dates (e.g., 01/01/2000, 01/02/2000, etc.) from 2000 to end of 2008. The default only gives me 4 ticks. I want more. Why doesn't this work? sdate<-as.POSIXct(strptime(date,format="%m/%d/%Y")) plot(ppt~sdate,type="l",ylim=c(0,47),col=1,lwd=1,pch=16,ylab="Salinity, psu",xlab="Year",las=1,main="Duck Key

Hi list. I'm plotting pluviometric (Rain) data as a barplot, and then adding the salinity variable to this plot as lines. Obviously as these Y scales are completely different the salinity appears at the lower part of the graph extremely compacted. I need to plot the line at the exactly same area of the barplot but with its own Y axis (at the right), so the salinity can use the plot area

Hi, I have used R for some time, but managed to avoid writing loops. But this time I am afraid there is no way around it. I have a dataframe with time and salinity (see below). I would like to extract the time intervals where salinity changes by less than 0.05. So using the values below this would mean that a subset was made from 2003-07-19 to 2003-07-24, where the change in salinity exceeded

Dear R Wizards, After much frustration and days of confusion I have finally broken down and am asking for help, which I don’t like doing, but I just can’t figure this one out on my own. I’ve conducted a laboratory experiment testing the effects of temperature and salinity on whether or not a biological event will occur (Go or NoGo). I’ve coded the factors temperature and salinity as factors for