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 <sstol...@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 > sstol...@gmail.com > > [[alternative HTML version deleted]] > > ______________________________________________ > R-help@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]] ______________________________________________ R-help@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.
