Thanks. This time the x-axis works well but the y-axis is missing.
2013/1/22 Mark Leeds <marklee...@gmail.com>
> 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.
