
Hi Jim,
thanks for your reply, I will try using this instead!
best,
Laura



On Feb 24, 2012, at 3:09 AM, Jim Lemon wrote:

> On 02/24/2012 09:30 AM, Laura Matrajt wrote:
>> 
>> Hi,
>>  I found some weird behavior of the function multhist in the plotrix package and I would like to ask you about it:
>> so, it works well to put two histograms next to each other, but I don't understand the breaks…
>> so for example,
>> library(plotrix)
>> mh<- list(rnorm(200, mean=100, sd=10), rnorm(200, mean=100, sd=10))
>> multhist(mh)                                           #produces a plot with 12 bins
>> multhist(mh, breaks=2)                        #produces a plot with 2 bins, as expected
>> multhist(mh, breaks=3)                         #produces a plot with 4 bins???
>> multhist(mh, breaks=4)                         #produces a plot with 4 bins
>> multhist(mh, breaks=7)                         #produces a plot with 6 bins???
>> 
>> etc. So odd number of breaks doesn't seem to work, and even more, some times it does the lowest even number of bins (so if I put breaks=7 I will get 6 bins) but other times it does the highest even number of bins (so if I put breaks=3 I get 4 bins)
>> 
>> In addition, things become even worse when you try to add labels to the bins.
>> 
>> 
>> So, all this to say that I am trying to put two histograms side to side, and put some correct number of bins and labels… does anyone have an idea of how to do this?
> 
> Hi Laura,
> This may be due to the way breaks are calculated when you just specify a suggested number of bins. However, when I tried this with barp:
> 
> mhmat<-rbind(table(cut(mh[[1]],breaks=3)),
> table(cut(mh[[2]],breaks=3)))
> library(plotrix)
> barp(mhmat)
> 
> I got what I expected, three groups of two bars. Perhaps this will do what you want.
> 
> Jim


	[[alternative HTML version deleted]]

