[R] multi-histogram plotting

[David Winsemius]

On Mar 13, 2012, at 5:33 PM, Sam Steingold wrote:

> I have a vector x:
> table(x)
>
>    2     3     4     5     6     7     8     9    10    11    12     
> 13    14
> 45547 11835  4692  2241  1386   820   593   425   298   239   176    
> 158   115
>   15    16    17    18    19    20    21    22    23    24    25     
> 26    27
>   94    88    76    67    47    46    40    20    30    22    20     
> 33    14
>   28    29    30    31    32    33    34    35    36    37    38     
> 39    40
>   20    10    12    10    11     8     9     8     9     9     8      
> 7     4
>   41    42    43    44    45    46    47    48    49    50    51     
> 52    53
>    5     4     6     5     2     5     4     4     3     1     6      
> 4     3
>   54    55    56    57    58    59    60    61    63    64    65     
> 66    67
>    2     2     1     4     5     2     5     1     3     2     1      
> 1     4
>   71    72    75    78    79    82    83    84    86    88    90     
> 92    93
>    2     3     1     2     2     2     2     1     2     2     1      
> 1     1
>   94    95    96    97    98    99   100   106   109   110   111    
> 112   119
>    1     2     1     1     3     1     1     3     1     2     2      
> 1     1
>  122   125   126   128   132   133   135   140   143   147   148    
> 157   162
>    1     1     1     1     1     1     1     1     1     1     1      
> 1     1
>  165   166   167   169   174   176   193   197   201   208   224    
> 236   339
>    1     1     1     1     1     1     1     1     1     1     1      
> 1     1
>  350   390   391   410   421   447   450   453   479   512   608    
> 679   754
>    1     1     1     1     1     1     1     1     1     1     1      
> 1     1
>  774   788   956   961  9597 14821
>    1     1     1     1     1     1
>
> I want to plot its histogram; moreover,
> I want to plot histograms of its several subsets on the same device
> (plot+lines+lines+... - right?)
> This means that I need a freq=FALSE plots.
> Now, a simple hist(x) produces a useless plot (vertical bar +  
> horizontal
> bar, which could be inferred from the table above).
> So, I want
> 1. a logarithmic vertical scale
> 2. an "interesting" horizontal scale
>
> so I do
>
> h <- hist(x, freq=FALSE, breaks=c(1:20,25,30,40,60,100,100000),  
> plot=FALSE)
>
> now I plot h$density vs 1:length(h$mids) or against
> h$mids[1:(length(h$mids)-1)] if I want to be honest,
> but I still don't see how to effect a log vertical scale.
> I can, of course, plot log(h$density), but then the number labels  
> will be wrong.

You could try apply a log transform to the appropriate component of  
the "h" object and using barplot to display the results.

-- 
David Winsemius, MD
West Hartford, CT