[R] help error histograma
David L Carlson
dcarlson at tamu.edu
Fri Aug 10 18:25:24 CEST 2012
Now that John has put your data into a readable format, there are a number of issues with your histogram that don't make much sense. You have enlarged the text of the labels and greatly enlarged the size of the title, but then printed no title (you have cex.main=6 and main="") and you have set probability=TRUE to get a density plot, but then labeled the axis with ylab="frequência" indicating it is a frequency plot. Also at .5 intervals the bars are very small. You might try the following alternative:
hist(dados[[1]], 0:110, xlab="chuva (mm/dia)", ylab="frequência", main="", las=1)
-------
David
> -----Original Message-----
> From: r-help-bounces at r-project.org [mailto:r-help-bounces at r-
> project.org] On Behalf Of John Kane
> Sent: Friday, August 10, 2012 9:54 AM
> To: aprendiz programa; r-help at r-project.org
> Subject: Re: [R] help error histograma
>
>
> Thank you but it is not a good idea to send attachments. I received
> them
> because the email came directly to me butw probablyw others did not
> because
> the R-list usually removes them to protect against viruses. It is
> better to
> put everything in the email.
> I have included your code and data below.
> I do not understand what some of your code is doing but the
> immediate
> problem is that your data has a range of 0 - 109 and you have a seq
> of
> seq(0,80,0.5) in the hist command
> Change to seq(0,110,0.5) and it should work.
> What is this supposed to do?
> is.na(dados) <- dados == -999.9
>
>
>
> John Kane
> Kingston ON Canada
> #=============Original Code=====================#
>
> library(MASS)
> dados<-read.table("/home/john/rdata/pelotas.txt",header=FALSE)
> is.na(dados) <- dados == -999.9 # what the devil? jrk
> png(filename="teste1.png",width=600,height=600)
> par(mfrow=c(2,1))
> #breaks<- seq(0,max(dados[[1]],by=0.5))
> hist(dados[[1]],seq(0,110,0.5), prob=TRUE),
> xlab="chuva
> (mm/dia)",ylab="frequência", main="",
> cex.lab=1.6,
> cex.sub=3,cex.axis=3,cex.main=6)
> # calculo do quantil ou percentil de 90%
> q=quantile(dados[1],0.90,na.rm=T)
> #===============================================
> dados <- dput(dados)
> structure(list(V1 = c(0, 0.4, 0, 0, 0.4, 0, 0, 0, 0, 0, 1, 0,
> 0, 0.4, 5, 16, 0, 6.6, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 9.8,
> 10.6, 5.6, 2.2, 3.2, 0, 0, 0.1, 0.3, 0, 15.8, 0, 0, 0, 0, 0,
> 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 30.5, 0, 0, 2.2, 44.8, 20.6,
> 9, 3.6, 2.2, 0, 0.6, 0, 0, 0, 0, 0, 0, 0, 13.2, 0.6, 0, 0, 0,
> 0, 0, 21.4, 0, 0, 3.4, 6.4, 0, 0, 0, 1, 0, 0, 0, 1.4, 0, 0, 0,
> 0, 0, 7, 0.2, 1.4, 0, 0, 0, 13.4, 18.6, 3.6, 7, 1.7, 0, 0, 0,
> 0, 0, 0, 0.2, 0, 0, 0, 77.6, 0, 22.6, 1.2, 15.4, 0, 0, 0, 0,
> 0, 0, 0, 0, 0.8, 33.2, 7, 6, 1.2, 0.4, 0.4, 0, 0, 0, 0.6, 74.8,
> 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 10.2, 3.6, 0.4, 1.2, 47.4, 0,
> 0, 0, 0, 5.2, 8.2, 4.2, 9.8, 0, 0, 0, 0, 0, 0, 0, 0, 0, 109,
> 0, 0, 75.6, 16, 0, 0, 0, 0, 0, 4.2, 20.6, 19.2, 4.4, 4.4, 0,
> 0, 0, 0, 0, 0, 0, 0, 0.4, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
> 0, 0, 0, 0, 0, 16, 0, 0, 0, 0, 0, 3, 0.5, 0, 0, 0, 0, 0, 0, 1.8,
> 54.6, 54, 3.4, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 27.8,
> 0, 9.8, 25.6, 21.8, 0.8, 0, 9.2, 43.2, 0, 18.6, 6.4, 0, 7, 0.8,
> 0, 0, 0, 0, 0, 3.8, 0, 0, 43.4, 0.6, 3.2, 0, 0, 0.6, 14.8, 0,
> 3.6, 0, 0, 0, 0, 3, 4.4, 0, 9.2, 0, 0, 0, 10.8, 0, 5, 12.6, 2.2,
> 1.2, 0, 0, 0, 16.8, 21.4, 0, 0, 0, 0, 8.6, 0, 0, 0, 52.8, 0,
> 0, 0, 0, 0, 0, 0, 2.2, 0, 0, 0, 14.8, 0.2, 0, 0, 0, 20.6, 0,
> 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 4, 3, 11.8, 7.4, 0, 0,
> 0, 0, 10, 0, 0.4, 0.2, 0, 9.2, 3.2, 0, 3.4, 0, 0, 0, 0, 0, 0,
> 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
> 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
> 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
> 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
> 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
> 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
> 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
> 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
> 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
> 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
> 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
> 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
> 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
> 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
> 0, 0, 0, 0, 0.4, 0.2, 0, 0, 0.4, 19.4, 6.8, 0, 0, 0, 0, 0, 0,
> 0, 0, 0, 0, 0, 0, 0, 1.8, 8.6, 5.4, 8.2, 24.4, 0, 20.6, 0, 0,
> 0, 13.6, 0, 0, 5.8, 0, 0, 0, 0, 8, 77.2, 0, 0, 0, 0.4, 0, 0,
> 0, 9, 0, 0, 38.8, 24.6, 0.4, 0, 8.6, 0.6, 0, 0, 0, 11.6, 2.2,
> 0, 0, 0.2, 0, 0, 0, 0.2, 15, 0, 0, 0, 2, 62.4, 0, 0, 0, 0, 15.4,
> 0, 28.8, 30.4, 0, 0, 0, 0, 0, 0, 0, 0, 11, 13.8, 0.6, 0, 34.6,
> 0, 5, 0.2, 59, 19.6, 0, 3.2, 14, 0, 0, 0, 0, 5, 1, 0, 19.2, 6.2,
> 1.5, 0, 0, 8.8, 14.6, 0, 0.4, 0.2, 0, 0, 0, 4.6, 0, 0, 14.2,
> 0.4, 0, 0, 0, 0, 1.2, 0, 0, 0, 0, 2.2, 8.2, 0, 7.8, 0, 8, 19,
> 0, 0, 0.6, 36.2, 5, 0, 0, 6.8, 4.4, 70.6, 16.8, 0, 4.2, 37.6,
> 81.2, 0, 0, 0, 0, 2, 0, 49.8, 0, 3.8, 2.4, 2.6, 0, 0, 4.8, 0,
> 0, 62.2, 2.6, 0, 0, 0, 0, 0, 0, 0, 0, 0.2, 0, 0, 0, 0, 0, 0,
> 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 18.4, 0, 0, 10.2, 0, 0, 0, 0,
> 0, 0, 0, 0, 0, 0, 0, 20.8, 52.2, 69.4, 0, 22, 0, 4.4, 0.6, 14.8,
> 1, 0, 0, 1.2, 0, 0, 0, 12.2, 0, 0, 0, 16.4, 31.8, 0, 0, 49, 36.2,
> 1, 10.8, 5.4, 0, 0, 0, 0, 7.6, 0, 32, 0.6, 0, 3, 2.4, 0, 0.2,
> 0, 0, 2.4, 21.8, 0, 0, 0, 1, 0, 0, 0, 0, 2, 17.4, 0, 0, 0, 19,
> 24.8, 6.6, 0, 0, 0, 1.4, 0, 0, 0, 0, 0, 0, 0, 0, 7.4, 0, 0, 0,
> 0, 0, 8.6, 9.8, 45, 20.4, 0, 7.4, 4.6, 0.2, 0, 0, 0, 9, 0, 0,
> 0, 0, 0, 0, 0, 0, 0, 0, 0, 18.2, 5.6, 3, 0, 0, 20, 0, 0, 0, 0,
> 0, 0, 0, 0, 3.4, 5.4, 51.6, 0.4, 0, 0, 0, 0, 0, 3, 102, 0, 4.8,
> 0, 0, 1.6, 0, 0, 65.2, 0.6, 0, 0.4, 0, 0, 0, 0.4, 5.4, 2, 0,
> 0, 0, 0, 0, 0, 2.6, 4.2, 0, 0, 14.4, 0, 0, 0, 0, 0, 0, 18.2,
> 0, 0, 0, 0, 0, 0, 5, 3.2, 26.4, 16.4, 0, 0, 0, 0, 0, 0, 0, 0,
> 1.2, 0, 0, 6, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 29.2, 51.6,
> 2.6, 0, 0, 1.8, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 5.8, 8.6, 0, 0,
> 0, 0, 0, 0, 0, 0, 0, 0, 11.8, 0, 21.6, 41.6, 21.4, 0.3, 0, 3.6,
> 0, 0.8, 0, 0, 1.8, 0, 1.6, 6.4, 0.2, 0, 0, 7.6, 0, 0, 0, 0, 0,
> 0, 0, 0, 0, 6.8, 0, 0, 0, 0, 0, 13.4, 8.4, 2, 0, 0, 0, 0, 9.6,
> 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 7.6, 9.8, 2.2, 0, 17.4, 2, 1.4,
> 0.6, 0.2, 0, 0, 0, 14.2, 0, 0, 22, 0, 17.8, 0, 7.4, 0, 0, 0,
> 11.2, 0, 0, 0, 4.2, 0, 39, 13.2, 0.2, 0, 2.2, 0.3, 0, 0, 0, 0,
> 0, 3.2, 0, 0, 0, 0, 0, 0.4, 0.8)), .Names = "V1", class =
> "data.frame",
> row.names = c(NA,
> -1214L))
>
> -----Original Message-----
> From: aprendizprogram at hotmail.com
> Sent: Fri, 10 Aug 2012 17:08:17 +0300
> To: jrkrideau at inbox.com, r-help at r-project.org
> Subject: RE: [R] help error histograma
>
> is attached data and script
> > Date: Fri, 10 Aug 2012 05:53:45 -0800
> > From: jrkrideau at inbox.com
> > Subject: Re: [R] help error histograma
> > To: aprendizprogram at hotmail.com; kridox at ymail.com; r-help at r-
> project.org
> >
> > You probably need to post the data here, or somewhere such
> as
> [url=http://www.mediafire.com/][b]MediaFire[/b][/url] for other
> readers to
> examine.
> >
> > To post it here use dput() and paste the results into the email.
> >
> > John Kane
> > Kingston ON Canada
> >
> >
> > > -----Original Message-----
> > > From: aprendizprogram at hotmail.com
> > > Sent: Fri, 10 Aug 2012 16:08:16 +0300
> > > To: kridox at ymail.com, r-help at r-project.org
> > > Subject: Re: [R] help error histograma
> > >
> > >
> > > Hi,
> > > I checked using the summary, and my data has 1000
> elementsranging from 0
> > > to 30 (values)
> > > but I do not understand the error
> > >
> > >
> > >> Date: Fri, 10 Aug 2012 11:18:01 +0900
> > >> From: kridox at ymail.com
> > >> To: aprendizprogram at hotmail.com
> > >> CC: r-help at r-project.org
> > >> Subject: Re: [R] help error histograma
> > >>
> > >> Hello,
> > >>
> > >> David L Carlson already suggested to you to check the range of
> your
> > >> data. Did you do it?
> > >>
> > >> Regards
> > >>
> > >>
> > >>
> > >> Le 10/08/2012 09:36, aprendiz programa a icrit :
> > >>>
> > >>> Hi,
> > >>> My error isErro em hist.default(dados[[1]], freq = TRUE,
> seq(0, 30,
> > >>> 0.5), prob = FALSE, : some 'x' not counted; maybe 'breaks' do
> not
> > >>> span range of 'x'
> > >>> hist(dados[[1]],seq(0, 30, 0.5), prob=TRUE, xlab="chuva
> > >>> (mm/dia)",ylab="frequjncia", main="", cex.lab=1.6,
> > >>> cex.sub=3,cex.axis=3,cex.main=6)
> > >>> Someone help me?
> > >>> [[alternative HTML version deleted]]
> > >>>
> > >>>
> > >>>
> > >>> ______________________________________________
> > >>> R-help at r-project.org mailing list
> > >>> 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 at r-project.org mailing list
> > > 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.
> >
> > ____________________________________________________________
> > FREE ONLINE PHOTOSHARING - Share your photos online with your
> friends and
> family!
> > Visit http://www.inbox.com/photosharing to find out more!
> >
> >
> _________________________________________________________________
>
> [1]Smileys Preview
> Get Free Smileys for Your IM & Email - Learn more
> at
> [2]www.crawler.com/smileys
> Works with AIM^®, MSN^® Messenger, Yahoo!^® Messenger, ICQ^®,
> Google
> Talk^™ and most webmails
>
> References
>
> 1. http://www.inbox.com/smileys
> 2. http://www.crawler.com/smileys
> ______________________________________________
> R-help at r-project.org mailing list
> 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.
More information about the R-help
mailing list