Lo resolví!
trest <- function(fact1, fact2, tesT, chr_var1, chr_var2, chr_test) {
tttt<-tapply(fact1, fact2, tesT)
pval <- as.data.frame(sapply(tttt, "[[", 'p.value'))
stati <- as.data.frame(sapply(tttt, "[[", 'statistic'))
TUT<-cbind(pval, stati)
cc1<-paste(as.character(chr_test),"_", as.character(chr_var1),"_",as.character(chr_var2),"_pvalor", collapse="_", sep="")
cc2<-paste(as.character(chr_test),"_", as.character(chr_var1),"_",as.character(chr_var2), "_statistic", collapse="_", sep="")
colnames(TUT)<-c(cc1, cc2)
return(TUT)
}
ALA_plu_any<-data.frame( Any=c("X1995", "X1996", "X1997", "X1998", "X1999", "X2000", "X2002", "X2003", "X2004", "X2005", "X2006", "X2007", "X2008", "X2009", "X2011", "X1995", "X1996", "X1997", "X1998", "X1999", "X2000", "X2002", "X2003", "X2004", "X2005", "X2006", "X2007", "X2008", "X2009", "X2010", "X2011", "X2012", "X2013", "X2014", "X2015", "X2006", "X2007", "X2009", "X2010", "X2011", "X2012", "X2013", "X2015", "X1995", "X1996", "X1997", "X1998", "X1999", "X2000", "X2002", "X2003", "X2004", "X2005", "X2006", "X2007", "X2008", "X2009", "X2010", "X2011", "X2012", "X2013", "X2014", "X2015", "X2002", "X2003", "X2004", "X2005", "X2006", "X2007", "X2008", "X2009", "X2010", "X2011", "X2012", "X2013", "X2014", "X2015", "X1996", "X1997", "X1998", "X1999", "X2000", "X2002", "X2003", "X2004", "X2005", "X2006", "X2007", "X2008", "X2009", "X2010", "X2011", "X2012", "X2013", "X2014", "X2015", "X1996", "X1997", "X1998", "X1999"), Pluviometria=c(999, 1167, 746, 476, 862.5, 604, 1094.5, 845, 844, 750, 560, 437.7, 1027, 969, 1020.6, 570, 1027, 586.2, 513.4, 649.4, 522.1, 768.6, 624, 551.2, 519.8, 505.6, 440.6, 787.1, 596.6, 780.9, 825.2, 467.8, 640.5, 787.3, 503, 581, 478, 636, 662, 905.3, 466.6, 20.4, 476.4, 618.9, 1183, 627.9, 580, 691.6, 557.1, 869.2, 716.7, 587.8, 592.9, 528.6, 465.3, 874.5, 634.2, 801, 846.3, 487.6, 575.4, 787.7, 411.6, 849.5, 838.5, 632.5, 582.5, 550, 516, 900.5, 643.5, 734.7, 807.9, 495.6, 738.3, 860.3, 410.8, 1268, 778, 781.5, 781.5, 600, 1015, 895.5, 762, 783, 704, 547, 1091.5, 743.5, 916.5, 1035.1, 614.3, 998.7, 1066.9, 696.7, 1558, 787.5, 836.5, 820))
EE<-trest(ALA_plu_any$Pluviometria, ALA_plu_any$Any, chisq.test, "plu", "any","chisqr")