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-boun...@r-project.org [mailto:r-help-bounces@r- > project.org] On Behalf Of John Kane > Sent: Friday, August 10, 2012 9:54 AM > To: aprendiz programa; r-help@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: aprendizprog...@hotmail.com > Sent: Fri, 10 Aug 2012 17:08:17 +0300 > To: jrkrid...@inbox.com, r-help@r-project.org > Subject: RE: [R] help error histograma > > is attached data and script > > Date: Fri, 10 Aug 2012 05:53:45 -0800 > > From: jrkrid...@inbox.com > > Subject: Re: [R] help error histograma > > To: aprendizprog...@hotmail.com; kri...@ymail.com; r-help@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: aprendizprog...@hotmail.com > > > Sent: Fri, 10 Aug 2012 16:08:16 +0300 > > > To: kri...@ymail.com, r-help@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: kri...@ymail.com > > >> To: aprendizprog...@hotmail.com > > >> CC: r-help@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@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@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@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. ______________________________________________ R-help@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.
