Re: [R] if statement

From: Ken Knoblauch <knoblauch_at_lyon.inserm.fr>
Date: Wed, 13 Jun 2007 07:05:55 +0000 (UTC)

Simon Blomberg <s.blomberg1 <at> uq.edu.au> writes:

>
> My solutions are usually too baroque, but does this do what you want?
>
> x <- rnorm(100)
> quants <- quantile(x, c(.3, .7))
> Case <- rep(2, length(x)) # 2 lies in the middle of the distribution
> Case[x <= quants[1]] <- 0
> Case[x >= quants[2]] <- 1
> Case
>
> On Mon, 2007-06-11 at 15:14 -0700, Jiong Zhang, PhD wrote:
> > Hi all,
> >
> > I have a rather naive question. I have the height of 100 individuals in
> > a table and I want to assign the tallest 30% as Case=1 and the bottom
> > 30% as Case=0. How do I do that?
> >

Or, how about,

x <- rnorm(100)
Case <- cut(x, quantile(x, c(0, 0.3, 0.7, 1)), c(0, 2, 1), TRUE)

ken

-- 
Ken Knoblauch
Inserm U846
Institut Cellule Souche et Cerveau
Département Neurosciences Intégratives
18 avenue du Doyen Lépine
69500 Bron
France
tel: +33 (0)4 72 91 34 77
fax: +33 (0)4 72 91 34 61
portable: +33 (0)6 84 10 64 10
http://www.pizzerialesgemeaux.com/u846/

______________________________________________
R-help_at_stat.math.ethz.ch 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.
Received on Wed 13 Jun 2007 - 07:14:24 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Wed 13 Jun 2007 - 07:31:38 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.