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]]
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) > > > *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
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]]
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) >> >> >> *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]]