I am sorry for my careless. I forgot implementing the last code
corresponding to setting y-axis. Thanks Mark Leeds and all guys.
2013/1/22 hp wan <huaping....@gmail.com>
> 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.
