[R] Histogram for Left Censored Data

[Jim Lemon]

Hi Steve,
Maybe something like this:

Sconc<-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)
hist(Sconc[,1],breaks=c(1450,1550,1650,1750,1900))
abline(v=1450,col="red")
text(1550,0.003,"limit of measurement")
arrows(1480,0.003,1450,0.003)

Jim


On Thu, Dec 31, 2015 at 9:20 PM, Steven Stoline <sstoline at gmail.com> 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)
>
>
> *Column 2* is an indicator for censoring "*1*" for left censored
> observations and "*0*" for non-censored (fully measured) observations.
>
> with many thanks
> steve
>
> --
> 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.
>

	[[alternative HTML version deleted]]