Bom dia,estou usando esta rotina para analisar o banco postado abaixo. (Veja no link) g<-read.csv2("C:\\Consulta\\HEMOAM\\Elene_fisio\\gest.csv") g[g==""]<-NA names(g) tabelas<-function(x){         qtde<-table(x)         perc<-100*prop.table(table(x))         cbind(qtde,perc) } options(digits=3) #imc g$IMC_f<-rep(NA,length(g$IMC)) summary(g$IMC) g$IMC_f[g$IMC<18.5]<-'Baixo peso' g$IMC_f[18.5<=g$IMC &g$IMC<25]<-'Normal' g$IMC_f[25<=g$IMC &g$IMC<30]<-'Pré-obesidade' g$IMC_f[30<=g$IMC &g$IMC<35]<-'Obesidade - I' g$IMC_f[35<=g$IMC &g$IMC<40]<-'Obesidade - II' g$IMC_f[g$IMC>=40]<-'Obesidade - III' #Idade g$ida<-rep(NA,length(g$Idade)) summary(g$Idade) g$ida[18<=g$Idade &g$Idade<25]<-'1. 18 a 25' g$ida[25<=g$Idade &g$Idade<32]<-'2. 25 a 32' g$ida[32<=g$Idade &g$Idade<39]<-'3. 32 a 39' apply(g[,c(-1,-2,-3,-7,-11,-14,-15,-16)],2,tabelas) http://www.datafilehost.com/download-64b7d14d.html Gostaria de obter um resulado semelhante a rotina tabelas, só que estratificado por grupo(Controle, Tratado).  Alguém tem uma sugestão? Edson Lira Estatístico Manaus-Amazonas