Oh, okay. I just missed it completely since your data didn't have any missing data so I could not visualise why it was there. I assume -999.9 was in older data.
John Kane Kingston ON Canada > -----Original Message----- > From: aprendizprog...@hotmail.com > Sent: Fri, 10 Aug 2012 20:51:48 +0300 > To: r-help@r-project.org > Subject: Re: [R] help error histograma > > > Hi, > The problem solved! > thank you very much !!! > Kane, > the command > is.na(dados) <- dados == -999.9, substituting missing > values b??b??(-999.9) for NA > > >> From: dcarl...@tamu.edu >> To: jrkrid...@inbox.com; aprendizprog...@hotmail.com; >> r-help@r-project.org >> Subject: RE: [R] help error histograma >> Date: Fri, 10 Aug 2012 11:25:24 -0500 >> >> 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="frequC*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="frequC*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="frequC*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^B., MSN^B. Messenger, Yahoo!^B. Messenger, ICQ^B., >>> Google >>> Talk^b?" 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. >> > > [[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 3D EARTH SCREENSAVER - Watch the Earth right on your desktop! ______________________________________________ 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.
