Re: [R] Lattice Histogram Scaling

2017-08-15 Thread Richard M. Heiberger
The original was conditioning the entire histogram, not just the kernel estimate, on bw. I am assuming that the bw in the condition should be bandwidths. What you want is easily done by using the latticeExtra package. This works from a new R session: url <- "http://www.econ.uiuc.edu/~roger/rese

Re: [R] Lattice Histogram Scaling

2017-08-15 Thread Roger Koenker
My apologies, the data can now be found at: url <- "http://www.econ.uiuc.edu/~roger/research/ebayes/velo.d"; x <- scan(url,skip = 1) If I could get each of the histograms to mimic what is produced by hist(x, 100, freq = FALSE) I’ve experimented with xlim, ylim, without success so far... url:

Re: [R] Lattice Histogram Scaling

2017-08-14 Thread Richard M. Heiberger
Your example is not easily reproducible. The REBayes requires Rmosek which requires a system command MOSEK. Please try again with an example using data in base R. Meanwhile, my guess is that you will need to do something like explicitly specifying xlim and ylim so all panels have the same limits.

[R] Lattice Histogram Scaling

2017-08-14 Thread Roger Koenker
I am trying to do some comparisons of density estimators using lattice. The code below attempts to plot the same histogram in each panel and then overplots a kernel estimate with different bandwidths. Finding packet.number() was a bit arduous, but seems to do what I want. My concern now is that cl