Hello,

Arun's code is much better.

Rui Barradas
Em 30-11-2012 18:50, arun escreveu:
Hi,
You could also try ?ave()
dat$Age_group<-ave(dat$Age,dat$GENDER,FUN=function(x){cut(x,labels=1:10,breaks=quantile(x,seq(0,1,.1)),include.lowest=TRUE)})
  dat
#  GENDER Age Age_group
#1      2  45         1
#2      1  58        10
#3      1  54         1
#4      2  71        10
#5      2  64         5
A.K.



----- Original Message -----
From: Rui Barradas <ruipbarra...@sapo.pt>
To: R Kozarski <r.kozar...@gmail.com>
Cc: r-help@r-project.org
Sent: Friday, November 30, 2012 1:07 PM
Subject: Re: [R] subgroup-based quantiles

Hello,

If you want Age quantiles by gender, you have to split the data by
gender, apply the same code then recombine the result.


fun <- function(x){
      Age_group <- cut(x[, "Age"], labels=c(1:10),
          breaks=quantile(x[, "Age"], seq(0,1,.1)),
          include.lowest = TRUE)
      cbind(x, Age_group)
}

result <- do.call(rbind, lapply(split(dat, dat[, "GENDER"]), fun))
rownames(result) <- seq_len(nrow(result))
result


Hope this helps,

Rui Barradas
Em 30-11-2012 12:18, R Kozarski escreveu:
Dear R users,
given the patient sample with their Gender and Age
           GENDER Age
      [1,]      2  45
      [2,]      1  58
      [3,]      1  54
      [4,]      2  71
      [5,]      2  64
    ...
I would like to create an another column which groups the patients wrt
Gender specific Age quantiles, following methodology similar to:

Age_group <- cut(Age, labels=c(1:10), breaks=quantile(Age,
seq(0,1,.1)),include.lowest = TRUE)

The function above allows me to group only wrt Age quantiles.

Best, Robert




--
View this message in context: 
http://r.789695.n4.nabble.com/subgroup-based-quantiles-tp4651415.html
Sent from the R help mailing list archive at Nabble.com.

______________________________________________
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.
______________________________________________
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.


______________________________________________
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.

Reply via email to