mi problema és el siguiente: tengo que hacer muchos test (shapiro.test, chisqr.test... ) en varias variables de un dataframe. Diseñé una funcion para ahorrar tiempo, pero no me funciona, el tapply me falla:
getwd()
load("ALA_plu_any.Rdata")
trest <- function(dafra, fact1, fact2, tesT){
tapply(dafra$fact1, dafra$fact2, tesT)
dafra_pvalor <- as.data.frame(sapply(tapply(dafra$fact1, dafra$fact2, tesT), "[[", 'p.value'))
colnames(dafra_pvalor)<-paste(as.character(fact1), as.character(fact2), as.character(tesT), "_", "p-valor", collapse="_", sep="")
dafra_xsquared <- as.data.frame(sapply(tapply(dafra$fact1, dafra$fact2, tesT), "[[", 'statistic'))
colnames(dafra_xsquared)<-paste(as.character(fact1), as.character(fact2), as.character(tesT), "_", "p-valor", collapse="_", sep="")
tttt <- cbind(dafra_pvalor, dafra_xsquared)
return(tttt)
}
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))
TTT <- trest(ALA_plu_any, Pluviometria, Any, chisq.test)
#dafra = dataframe ; fact1 y fact 2= factores ; tesT = test
#Muchas gracias!!