polish example

Author: edzer

R/aggregate.R

@@ -133,7 +133,8 @@ aggregate.sf = function(x, by, FUN, ..., do_union = TRUE, simplify = TRUE,
 #' # dasymetric mapping
 #' bir.grd = st_interpolate_aw(nc[c("BIR74","BIR79")], extensive = TRUE, grd.addr, weights = "ones")
 #' xbir.grd = bir.grd # copy for plotting
-#' xbir.grd$BIR74[xbir.grd$BIR74 == 0] = 1 # so that logz shows finite values
+#' smallest_pos = function(x) min(x[x > 0])
+#' xbir.grd$BIR74[xbir.grd$BIR74 == 0] = smallest_pos(xbir.grd$BIR74) # so that logz shows finite values
 #' \donttest{plot(xbir.grd["BIR74"], logz = TRUE, main = "dasym. redistributed birth counts, 1974-")}
 #' # verify sums:
 #' apply(as.data.frame(bir.grd)[1:2], 2, sum)