[R] lattice and histogram

[Deepayan Sarkar]

--- Ted.Harding at nessie.mcc.ac.uk wrote:
> Hello again folks,
> 
> Thanks to all for the advice on getting hold of
> "grid"
> and "lattice". I think I'm getting the hang of it.
> 
> This produces very satisfying arrays of histograms.
> The typical command I am using is
> 
>   histogram( ~ DATA.df$X | DATA.df$F, \
>              type=c("count"), layout=c(5,5) )


A more typical and elegant usage would be

histogram(~X|F, data = DATA, [...])



> For some variables, say again typically "X", for one
> (or sometimes more) level of "F" all values of "X"
> are
> missing (coded as "NA" in the data file).
> 
> If I run the above command with such a variable "X",
> I get
> 
> Error in if (lim[1] > lim[2]) stop("Improper value
> of limit") ;
>         missing value where logical needed
> 
> There is no problem if, at given level of "F", "X"
> is a mixture
> of good values and "NA"s; it seems to do the right
> thing by
> simply ignoring the latter. It's the case of "data =
> empty set"
> that seems to cause the trouble.


This is definitely a bug, which should be fixed in the
next release of lattice. For now, what I can suggest
is to change a couple of lines in the histogram()
function. Replace


xlim.l <- min(xlim.l, tem$xlim[1])
xlim.u <- max(xlim.u, tem$xlim[2])

by


xlim.l <- min(xlim.l, tem$xlim[1], na.rm = TRUE)
xlim.u <- max(xlim.u, tem$xlim[2], na.rm = TRUE)

This should take care of your problem, though there
still might be some other similar loopholes.

PS: variable names like F might cause problems, since
a lot of the lattice code use F/T instead of
FALSE/TRUE. Again, this should not be so, and will be
fixed in the next release.


__________________________________________________
Do You Yahoo!?
Make international calls for as low as $.04/minute with Yahoo! Messenger
http://phonecard.yahoo.com/
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
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
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._