R help - Oct 1999 - intervals on plot

Dear all R users

I have been trying to plot a group of intervals by using the bottom
code:

> breaks <- factor(cut(col3,breaks=360+5*(0:155)))
> par(tck=0,lty="solid")
> plot(breaks,ann=FALSE,axes=FALSE)
col3 is just a column of 1500 floating point data.

Unfortunately when i do the plot, the x-axis has some of the intervals
as labels dispite me having both 'ann' and 'axes' both as FALSE.

Could anyone tell me how to remove these from my plot or how to bypass
this altogether, maybe by changing the names of the intervals to
something else.

Thanks to all of you for any help.



-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help 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-help-request at
stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._

>>>>> "Nikita" == Nikita Atkins <nda at
mailcall.com.au> writes:
    Nikita> Dear all R users I have been trying to plot a group of
    Nikita> intervals by using the bottom code:

    Nikita> col3 is just a column of 1500 floating point data.
(added a definition for 'col3' (and dropped ">")
  --> completely executable code, please use this in reports like this! --
MM)

    col3 <- 750+ 100* rt(1500, df = 3)

    breaks <- factor(cut(col3,breaks=360+5*(0:155)))
    par(tck=0,lty="solid")
    plot(breaks,ann=FALSE,axes=FALSE)

    Nikita> Unfortunately when i do the plot, the x-axis has some of the
    Nikita> intervals as labels dispite me having both 'ann' and
'axes'
    Nikita> both as FALSE.

The reason is an imperfect implementation of  plot.factor(),
or barplot.default(), respectively,
which are called by the above  plot() command.

I'll have a look at fixing this. 

Note however that you are really just (trying) to plot a histogram of col3
and I wonder why you are *not* just using
    hist(col3)
or
    hist(col3, breaks =360+5*(0:155)), ...)

    Nikita> Could anyone tell me how to remove these from my plot or how to
    Nikita> bypass this altogether, maybe by changing the names of the
    Nikita> intervals to something else.

A workaround:

  tb <- table(breaks); dimnames(tb) <- NULL
  barplot(tb, axes=FALSE)


Apropos, we use locally the function

  unname <- function (obj) 
  {
      if (length(names(obj))) 
	  names(obj) <- NULL
      if (length(dimnames(obj)) && !is.data.frame(obj)) 
	  dimnames(obj) <- NULL
      obj
  }

with which you could abbreviate the above to

     barplot(unname(table(breaks)), axes= FALSE)

Martin Maechler <maechler at stat.math.ethz.ch>
http://stat.ethz.ch/~maechler/
Seminar fuer Statistik, ETH-Zentrum  LEO D10	Leonhardstr. 27
ETH (Federal Inst. Technology)	8092 Zurich	SWITZERLAND
phone: x-41-1-632-3408		fax: ...-1228			<><
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help 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-help-request at
stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._