[R] help error histograma
John Kane
jrkrideau at inbox.com
Fri Aug 10 20:07:49 CEST 2012
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: aprendizprogram at hotmail.com
> Sent: Fri, 10 Aug 2012 20:51:48 +0300
> To: r-help at 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: dcarlson at tamu.edu
>> To: jrkrideau at inbox.com; aprendizprogram at hotmail.com;
>> r-help at 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-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="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: 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^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 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 3D EARTH SCREENSAVER - Watch the Earth right on your desktop!
More information about the R-help
mailing list