Hi, The problem solved! thank you very much !!! Kane, the command > is.na(dados) <- dados == -999.9, substituting missing values ââ(-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="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. > [[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.
