Hi Hp: I decided to use axis because names.arg wasn't working very well.
Below is
an improvement and I'm sure, given enough effort, it can be improved much
more. But
atleast below will get you started.
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)
binstrings <- as.character(breaks)
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)
temp <- barplot(percentage,xlab="",space=1.0, axes=FALSE)
print(temp)
axis(1,labels=binstrings,at=seq(1.5,83.5,2),cex.axis=0.8)
axis(2,labels=seq(0,1,0.1),at=seq(0,1,0.1),cex.axis=0.8)
On Mon, Jan 21, 2013 at 8:26 PM, Mark Leeds <marklee...@gmail.com> wrote:
> Hi: if you mean that you want the labels to be the strings "-1.55",
> "-1.50" etc then you
> can use the names.arg argument for that. You can create a vector of
> character strings
>
> binstrings = c("-1.55","-1.50", .... )
>
> or just convert your data to strings using as.character. Then, in the call
> to barplot, include names.arg = binstrings. I think that should work.
>
>
>
>
>
>
> On Mon, Jan 21, 2013 at 8:16 PM, hp wan <huaping....@gmail.com> wrote:
>
>> But the x-axis of barplot is still not what I want. The xlab is breaks,
>> not -1.55,-1.50,....,0.55.
>>
>>
>> 2013/1/22 Mark Leeds <marklee...@gmail.com>
>>
>>> 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.
