RDX2 X  relNum source 8function(points=nsw.sam$points, trt=nsw$trt, col="red"){  x <- points[,1]  y <- points[,2]  p <- getPath(col=col) + nn <- which(pointInPoly(x,y, p$x, p$y))  print(table(trt[nn]))  invisible(nn)  }þý points $ nsw.samÿþ trtÿ nswÿþ$ col redþ srcref  srcfileò encoding native.enc timestampAÑ™¶RÀ class POSIXt POSIXctþ filename /tmp/johnm.relNum.R wd /Users/johnm/r/lundþþÿ srcfileþÿ srcrefþ  ÿ ÿÿ srcrefþ  ÿ ÿÿ srcrefþ + ÿ ÿÿ srcrefþ  ÿ ÿÿ srcrefþ  ÿ ÿÿ srcrefþþ { <- x [ÿû?ðþþÿ yÿÿû@þþÿ p getPathÿÿþþÿ nn which pointInPolyÿÿÿÿÿþÿÿÿþþþþ print tableÿÿÿþþþ invisibleÿþþÿÿ function(x0, y0, x, y)  {  cn = rep(0, length(x0))  xp = x[length(x)]  yp = y[length(y)]  for(i in 1:length(x)) { 7 cn = cn + ifelse(((yp <= y0 & y[i] > y0) | 8 (yp > y0 & y[i] <= y0)) & 4 (x0 < xp + ((y0 - yp) / E (y[i] - yp)) * (x[i] - xp)), " 1, 0)  xp = x[i]  yp = y[i] }  cn %% 2 == 1  }þý x0û y0ûÿûÿûþÿ = cn rep lengthÿþþþ ÿ xpÿÿ#ÿÿþþþ ÿ ypÿÿ#ÿÿþþþ for i :?ð#ÿÿþþÿ ÿ!ÿ +!ÿ ifelse & ( |,ÿ+ÿ <=%ÿÿþ >ÿÿ'ÿþÿþþþ,ÿ+ÿ/ÿ%ÿÿþ.ÿÿÿ'ÿþÿþþþþþ,ÿ <ÿ)ÿ$ÿ *,ÿ /,ÿ -ÿ%ÿþþ,ÿ3ÿÿÿ'ÿþ%ÿþþþþ,ÿ3ÿÿÿ'ÿþ$ÿþþþþþþþ?ðþþþ ÿ$ÿÿÿ'ÿþþ ÿ%ÿÿÿ'ÿþþþþ == %%!ÿ@þ?ðþþÿÿ %function(lty = "dotted", col = "red")  {  xp = numeric()  yp = numeric() n = 0 + while(!is.null((l = locator(1)))) {  xp = c(xp, l$x)  yp = c(yp, l$y)  n = n + 1  if(n > 1) : segments(xp[n-1], yp[n - 1], xp[n], yp[n], . lty = lty, col = col) } B segments(xp[n], yp[n], xp[1], yp[1], lty = lty, col = col)  list(x = xp, y = yp)  }þý lty dotted$ÿ redþÿ ÿ$ÿ numericþþ ÿ%ÿ7ÿþþ ÿ nþ while ! is.null,ÿ ÿ l locator?ðþþþþþÿ ÿ$ÿ c$ÿÿ<ÿÿþþþ ÿ%ÿ>ÿ%ÿÿ<ÿÿþþþ ÿ8ÿ)ÿ8ÿ?ðþþ if/ÿ8ÿ?ðþ segmentsÿ$ÿ3ÿ8ÿ?ðþþÿ%ÿ3ÿ8ÿ?ðþþÿ$ÿ8ÿþÿ%ÿ8ÿþ6ÿ6ÿÿÿþþþþ@ÿÿ$ÿ8ÿþÿ%ÿ8ÿþÿ$ÿ?ðþÿ%ÿ?ðþ6ÿ6ÿÿÿþ listÿ$ÿÿ%ÿþþ showLineÿ function(col="cyan"){ $ print("Click at any two points")  xy <- locator(n=2)  b <- diff(xy$y)/diff(xy$x)  a <- xy$y[1] - b * xy$x[1]  abline(a,b, col="cyan") & cat("\ny =", a, "+", b, "x", "\n")  below <- function(x,y){  y - b*x <= a  } below  }þý$ÿ cyanþ ÿ $ ÿò ÿ native.enc ÿAÑš<#ÿ POSIXt POSIXctþÿ /tmp/johnm.showLine.Rÿ /Users/johnm/r/lundþþÿ srcfileþÿ srcrefþ  ÿCÿÿ srcrefþ  ÿCÿÿ srcrefþ  ÿCÿÿ srcrefþ  ÿCÿÿ srcrefþ & ÿCÿÿ srcrefþ    ÿCÿÿ srcrefþ    ÿCÿÿ srcrefþþÿÿ Click at any two pointsþÿ xy=ÿ8ÿ@þþÿ b2ÿ diffÿDÿÿþþFÿÿDÿÿþþþþÿ a3ÿÿÿDÿÿþ?ðþ1ÿEÿÿÿDÿÿþ?ðþþþþ ablineGÿEÿÿ cyanþ cat  y =Gÿ +Eÿ x  þÿ below function$ÿû$ÿûþ ÿ    ÿCÿÿ srcrefþþÿ.ÿ3ÿÿ1ÿEÿÿþþGÿþþ function(x,y){  y - b*x <= a  }þþJÿþþ