[R-br] TABELAS POR GRUPO

[Edson Lira]

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
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120820/a9afb289/attachment.html>