Hi Hp: I'm not sure what you want. Below looks slightly better but if you want the x-axis to have better labels, then there is a names.arg argument to barplot. see ?barplot.
I don't know what each of the bins represents but if you want to include all of those zeros then I don't think there's any way of getting around the fact that there is so much space in the barplot. Maybe someone else has a suggestion. breaks=c(-1.55,-1.50,-1.45,-1.40,-1.35,-1.30,-1.25,-1.20,-1.15,-1.10,-1.05,-1.00,-0.95,-0.90,-0.85,-0.80,-0.75,-0.70, -0.65,-0.60,-0.55,-0.50,-0.45,-0.40,-0.35,-0.30,-0.25-0.20,-0.15,-0.10,-0.05,0.00,0.05,0.10,0.15,0.20,0.25,0.30,0.35,0.40, 0.45,0.50,0.55) counts=c(287,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,212,2624,2918,0,0,0,75,36317,4963,0,0,2462,0,0,0,0,0,142) percentage=counts/sum(counts) barplot(percentage,xlab="breaks",xlim=c(1,42)) On Mon, Jan 21, 2013 at 7:47 PM, hp wan <huaping....@gmail.com> wrote: > Ok, that is no problem. > > > 2013/1/22 Mark Leeds <marklee...@gmail.com> > >> let me look at but it's probably best to send to the whole list because >> there are many >> people on it way more knowledgable than myself. I'm ccing the list and >> hope you >> don't mind. my fault for replying privately initially. >> >> >> >> On Mon, Jan 21, 2013 at 7:36 PM, hp wan <huaping....@gmail.com> wrote: >> >>> Thanks for your reply! >>> >>> >>> breaks=c(-1.55,-1.50,-1.45,-1.40,-1.35,-1.30,-1.25,-1.20,-1.15,-1.10,-1.05,-1.00,-0.95,-0.90,-0.85,-0.80,-0.75,-0.70,-0.65,-0.60,-0.55,-0.50,-0.45,-0.40,-0.35,-0.30,-0.25-0.20,-0.15,-0.10,-0.05,0.00,0.05,0.10,0.15,0.20,0.25,0.30,0.35,0.40,0.45,0.50,0.55) >>> >>> counts=c(287,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,212,2624,2918,0,0,0,75,36317,4963,0,0,2462,0,0,0,0,0,142) >>> percentage=counts/sum(counts) >>> barplot(percentage,xlab=breaks) >>> >>> The horizontal value (that is xlabe) looks very ugly. I hope it looks >>> like the xlab of hist, that is x axis correspond to breaks. >>> >>> After ?barplot, I also have no idea to implement it. >>> >>> 2013/1/22 Mark Leeds <marklee...@gmail.com> >>> >>>> I'm not sure that I understand but can't you just take the data and >>>> divide it by the sum of the data and plot that ? >>>> >>>> >>>> >>>> On Mon, Jan 21, 2013 at 6:36 PM, hp wan <huaping....@gmail.com> wrote: >>>> >>>>> Thanks for your reply. >>>>> >>>>> If I set the probability = FALSE, the column values are corresponding >>>>> to the refrequency (the numbers of values falling in intervals). I want >>>>> the >>>>> coulumn values are percentage, that is x$counts/sum(x$counts). >>>>> >>>>> >>>>> 2013/1/22 Mark Leeds <marklee...@gmail.com> >>>>> >>>>>> Hi: the density integrates to 1 but the actual height of the density >>>>>> at each point is not less necessarily than 1. for what you want, you >>>>>> should >>>>>> be using probability = FALSE. >>>>>> >>>>>> you can do pnorm(x=0,0,1) to see this. >>>>>> >>>>>> >>>>>> >>>>>> On Mon, Jan 21, 2013 at 5:18 PM, hp wan <huaping....@gmail.com>wrote: >>>>>> >>>>>>> Hi All, >>>>>>> >>>>>>> When carrying out hist(samples,breaks=50,probability=TRUE), the >>>>>>> column >>>>>>> values are considerably greater than 1, which seams very >>>>>>> unreasonable. The >>>>>>> plot is attached. >>>>>>> >>>>>>> I think the column value of the hist plot should correspond to >>>>>>> x$counts/sum(x$counts) >>>>>>> (x=hist(samples,breaks=50,probability=TRUE)). The size of data is >>>>>>> a >>>>>>> little bit larger, causing failure of uploading. If you need the >>>>>>> data, I >>>>>>> can email it to you. >>>>>>> >>>>>>> Can anyone help me? >>>>>>> >>>>>>> Thanks! >>>>>>> >>>>>>> Best regares, >>>>>>> >>>>>>> Huaping Wan >>>>>>> >>>>>>> ______________________________________________ >>>>>>> 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<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.
