# [R] [FORGED] Histogram for Left Censored Data

David Winsemius dwinsemius at comcast.net
Fri Jan 1 21:42:37 CET 2016

```> On Jan 1, 2016, at 3:45 AM, Steven Stoline <sstoline at gmail.com> wrote:
>
> Dear Rolf:
>
>
> The histogram should contain a bar(s) for the censored data values replaced
> by their detection limit(s) with different color than other bars for the
> noncensored values . In this example there are only 3 censored values with
> only one detection limit of DL = 1450.
>
>
> with many thanks
> steve
>
>
>
> On Thu, Dec 31, 2015 at 4:16 PM, Rolf Turner <r.turner at auckland.ac.nz>
> wrote:
>
>> On 31/12/15 23:20, Steven Stoline wrote:
>>
>>> Dear All:
>>>
>>> I need helps with creating histograms for data that include left
>>> censored observations.
>>>
>>> Here is an example of left censored data
>>>
>>>
>>>
>>> *Sulfate.Concentration*
>>> <-matrix(c(1450,1800,1840,1820,1860,1780,1760,1800,1900,1770,1790,
>>> 1780,1850,1760,1450,1710,1575,1475,1780,1790,1780,1450,1790,1800,
>>> 1,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,1,0,0),24,2)
>>>

myhist <- hist(sulfate[,1], breaks=c(1400,1451,1500,1600,1700,1800,1900), col=c(1,rep(2,5)), xaxt="n")
#  plots with no x axis labeling
myhist
#------------------
\$breaks
[1] 1400 1451 1500 1600 1700 1800 1900

\$counts
[1]  3  1  1  0 14  5

\$density
[1] 0.0024509804 0.0008503401 0.0004166667 0.0000000000 0.0058333333 0.0020833333

\$mids
[1] 1425.5 1475.5 1550.0 1650.0 1750.0 1850.0

\$xname
[1] "sulfate[, 1]"

\$equidist
[1] FALSE

attr(,"class")
[1] "histogram"
#---rebuild the x-axis ----------------
axis(1, at=c(myhist\$mids[1],myhist\$breaks[-(1:2)]), labels=c("<1450", myhist\$breaks[-(1:2)]))

--
David.

>>>
>>> *Column 2* is an indicator for censoring "*1*" for left censored
>>> observations and "*0*" for non-censored (fully measured)
>>> observations.
>>>
>>
>> And what, pray tell, do you want the resulting histogram to look like?
>> See e.g. fortune("mind_read").
>>
>> cheers,
>>
>> Rolf Turner
>>
>> --
>> Technical Editor ANZJS
>> Department of Statistics
>> University of Auckland
>> Phone: +64-9-373-7599 ext. 88276
>>
>
>
>
> --
> Steven M. Stoline
> 1123 Forest Avenue
> Portland, ME 04112
> sstoline at gmail.com
>
> 	[[alternative HTML version deleted]]
>
> ______________________________________________
> R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.

David Winsemius
Alameda, CA, USA

```