From Martin Maechler <maechler@stat.math.ethz.ch>  Mon Jun  2 08:31:42 1997
From: Martin Maechler <maechler@stat.math.ethz.ch> (Martin Maechler)
Date: Mon, 2 Jun 97 09:31:42 +0200
Subject: R-alpha: plot(1); axis(1, at=10) >>> Seg.fault  -- and a patch
In-Reply-To: <199706020519.HAA01989@aragorn.ci.tuwien.ac.at> (message from
 Kurt Hornik on Mon, 2 Jun 1997 07:19:10 +0200)
Message-ID: <9706020731.AA00382@>

The patched version of  axis(..) which I "seemed" to post here  a while ago,
accidentally was the 0.49 original -- oh dear ! -- S_O_R_R_Y!

Here is the 'real' one:

M> This problem has been mentioned  by Arne Kovac <maak@stats.bris.ac.uk>
M> on May 9-10. 
M> The patch is really a replacement of 'axis' in  RHOME/src/library/base/funs/

M> It contains the fixes which where proposed Arne Kovac
M> and fixes all the errors he reported on May 9.

"axis" <-
function (which, at, labels = TRUE, ...) 
{
	if (which%%2 == 1) {
		axp <- par("xaxp")
		usr <- par("usr")[1:2]
		log <- par("xlog")
	} else {
		axp <- par("yaxp")
		usr <- par("usr")[3:4]
		log <- par("ylog")
	}
	ind <-
	  if (missing(at)) {
	    at <-
	      if (log) {
		if (usr[2] < usr[1] + 1) {
		  seq(axp[1], axp[2], length = axp[3] + 1)
		} else {
		  p1 <- ceiling(min(usr))
		  p2 <- floor(max(usr))
		  if (p2 - p1 < 2) {
		    at <- c(1, 2, 5) * 10^rep(p1:p2, rep(3, p2 - p1 + 1))
		    at[10^usr[1] <= at & at < 10^usr[2]]
		  }
		  else 10^seq(p1, p2, by = 1)
		}
	      } else seq(axp[1], axp[2], length = axp[3] + 1)
	    rep(TRUE, length(at))
	  } else {
	    at <- sort(at)
	    if (log) (10^usr[1] <= at & at <= 10^usr[2])
	    else     (	 usr[1] <= at & at <=	 usr[2])
	  }
	if(any(ind)) {
	  labels <-
	    if (is.logical(labels)) {
	      if (labels) {
		if (!log) 
		  at[abs(at/(max(at) - min(at))) < 0.001] <- 0
		format(at, trim = T)
	      }
	      else rep("", length(at))
	    } else format(labels, trim = T)
	  .Internal(axis(which, as.double(at[ind]), labels[ind], ...))
	} else warning("axis: Couldn't draw anything; maybe 'at' out of range.")
}

=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-
r-devel mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-devel-request@stat.math.ethz.ch
=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-


R-alpha: plot(1); axis(1, at=10) >>> Seg.fault -- and a patch
Martin Maechler