Caros,

rotineiramente me deparo com um problema que não sei como resolver. Considerem o seguinte exemplo:

(c1 <- trunc(runif(1000,1,10)))

(c2 <- trunc(runif(1000,5,70)))

(c3 <- trunc(runif(1000,1,4)))

(var <- rnorm(1000,10,2))

(peso <- round(rnorm(1000,1,.2)^2,1))

df <- data.frame(c1,c2,c3,var,peso)

Gostaria de uma função que pudesse me fornecer em uma única tabela do tipo

cols = c1,c2,c3,count,medvar

e, mais genericamente,

cols = c1,c2,c3,somapeso,medpondvar

com todas as combinações existentes (as vezes todas as teoricamente possiveis) entre c1,c2,c3

c1

c2

c3

c1*c2

c1*c3

c2*c3

c1*c2*c3

De forma mais específica, gostaria:

t1 <- aggregate(list(count = df$var), list(c1=rep('*tot',1000) ,c2=rep('*tot',1000),c3=rep('*tot',1000)),length)

tc1 <- aggregate(list(count = df$var), list(c1=df$c1 ,c2=rep('*tot',1000),c3=rep('*tot',1000)),length)

tc2 <- aggregate(list(count = df$var), list(c1=rep('*tot',1000),c2=df$c2 ,c3=rep('*tot',1000)),length)

tc3 <- aggregate(list(count = df$var), list(c1=rep('*tot',1000),c2=rep('*tot',1000),c3=df$c3 ),length)

tc1c2 <- aggregate(list(count = df$var), list(c1=df$c1 ,c2=df$c2 ,c3=rep('*tot',1000)),length)

tc1c3 <- aggregate(list(count = df$var), list(c1=df$c1 ,c2=rep('*tot',1000),c3=df$c3 ),length)

tc2c3 <- aggregate(list(count = df$var), list(c1=rep('*tot',1000),c2=df$c2 ,c3=df$c3 ),length)

tc1c2c3 <- aggregate(list(count = df$var), list(c1=df$c1 ,c2=df$c2 ,c3=df$c3),length)

comp1 <- rbind(t1,tc1,tc2,tc3,tc1c2,tc1c3,tc2c3,tc1c2c3)

rm(t1,tc1,tc2,tc3,tc1c2,tc1c3,tc2c3,tc1c2c3)

t1 <- aggregate(list(medvar = df$var), list(c1=rep('*tot',1000) ,c2=rep('*tot',1000),c3=rep('*tot',1000)),mean)

tc1 <- aggregate(list(medvar = df$var), list(c1=df$c1 ,c2=rep('*tot',1000),c3=rep('*tot',1000)),mean)

tc2 <- aggregate(list(medvar = df$var), list(c1=rep('*tot',1000),c2=df$c2 ,c3=rep('*tot',1000)),mean)

tc3 <- aggregate(list(medvar = df$var), list(c1=rep('*tot',1000),c2=rep('*tot',1000),c3=df$c3 ),mean)

tc1c2 <- aggregate(list(medvar = df$var), list(c1=df$c1 ,c2=df$c2 ,c3=rep('*tot',1000)),mean)

tc1c3 <- aggregate(list(medvar = df$var), list(c1=df$c1 ,c2=rep('*tot',1000),c3=df$c3 ),mean)

tc2c3 <- aggregate(list(medvar = df$var), list(c1=rep('*tot',1000),c2=df$c2 ,c3=df$c3 ),mean)

tc1c2c3 <- aggregate(list(medvar = df$var), list(c1=df$c1 ,c2=df$c2 ,c3=df$c3),mean)

comp2 <- rbind(t1,tc1,tc2,tc3,tc1c2,tc1c3,tc2c3,tc1c2c3)

rm(tc1,tc2,tc3,tc1c2,tc1c3,tc2c3,tc1c2c3)

comp <- merge(comp1,comp2,by=c('c1','c2','c3'),all=T)

rm(comp1,comp2)

Exite algum pacote que faca isso?

_______________________________________________
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.