| binApply2D {oce} | R Documentation |
The function FUN is applied to f in bins specified by
xbreaks and ybreaks. (If FUN is mean,
consider using binMean2D instead, since it should be faster.)
binApply2D(x, y, f, xbreaks, ybreaks, FUN, ...)
x |
a vector of numerical values. |
y |
a vector of numerical values. |
f |
a vector of data to which the elements of |
xbreaks |
values of x at the boundaries between bins; calculated using
|
ybreaks |
values of y at the boundaries between bins; calculated using
|
FUN |
function to apply to the data |
... |
arguments to pass to the function |
A list with the following elements: the breaks in x and y
(xbreaks and ybreaks), the break mid-points (xmids and
ymids), and a matrix containing the result of applying function
FUN to f subsetted by these breaks.
Dan Kelley
Other bin-related functions: binApply1D,
binAverage, binCount1D,
binCount2D, binMean1D,
binMean2D
library(oce)
## Not run:
## secchi depths in lat and lon bins
if (require(ocedata)) {
data(secchi, package="ocedata")
col <- rev(oce.colorsJet(100))[rescale(secchi$depth,
xlow=0, xhigh=20,
rlow=1, rhigh=100)]
zlim <- c(0, 20)
breaksPalette <- seq(min(zlim), max(zlim), 1)
colPalette <- rev(oce.colorsJet(length(breaksPalette)-1))
drawPalette(zlim, "Secchi Depth", breaksPalette, colPalette)
data(coastlineWorld)
mapPlot(coastlineWorld, longitudelim=c(-5, 20), latitudelim=c(50, 66),
grid=5, fill='gray', projection="+proj=lcc +lat_1=50 +lat_2=65")
bc <- binApply2D(secchi$longitude, secchi$latitude,
pretty(secchi$longitude, 80),
pretty(secchi$latitude, 40),
f=secchi$depth, FUN=mean)
mapImage(bc$xmids, bc$ymids, bc$result, zlim=zlim, col=colPalette)
mapPolygon(coastlineWorld, col='gray')
}
## End(Not run)