Message-ID: <bbdc7ed01001131010k53ba1d9ch28540ba8eca03594@mail.gmail.com>
Date: 2010-01-13T18:10:10Z
From: Steve Lianoglou
Subject: Ask for histogram
In-Reply-To: <165fd86c1001130958y3c9f6ac5p90a82de0ed588bcd@mail.gmail.com>
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