[R] Ask for histogram

[Steve Lianoglou]

Hi,

On Wed, Jan 13, 2010 at 12:58 PM, Yi Du <abraham.du at gmail.com> wrote:
> Hi,
>
>
> I use a vector of data to draw the histogram, but it is different from the
> graph by SAS. Can you check it for me please?

How are we supposed to check something without data, pictures, etc?
What do you want checking, exactly?

> b is a column vector of 4332
>
> hist(b,probability=T,breaks=30,col='lightblue',ylim=c(0,1))
> rug(b)
>
> When I used rug, I find the records are smaller than 4332. I don't know
> where I did wrong.

What do you mean? Is the histogram that you're getting surprising? Is
the result of adding a "rug" surprising?

Are you actually trying to count 4332 tick marks at the bottom of your
plot? What records are smaller than 4332?

Try to see if what rug returns, eg:

r <- rug(b)

length(r) should be as long as your `b` vector

I'm not sure what you're asking, but hopefully some of the info I
threw at you is helpful. Please be a bit more specific with any follow
up if you still find anything confusing.

-steve

-- 
Steve Lianoglou
Graduate Student: Computational Systems Biology
 | Memorial Sloan-Kettering Cancer Center
 | Weill Medical College of Cornell University
Contact Info: http://cbio.mskcc.org/~lianos/contact