Hi, I did what you have told me but I can not do it well, NUMERO Provincia CENAES IE.2005 IE.2006 IE.2007 IE.2008 29 183 Alicante 7711 42.000 46.000 55.000 47.00030 479 Alicante 8299 1158.000 476.000 503.000 300.00031 591 Alicante 8129 540.000 764.000 1037.000 938.00032 1429 Alicante 8219 5233.000 6234.000 4659.000 2212.00033 2681 Alicante 8299 17.000 8.000 10.000 9.00034 3374 Alicante 8299 346.000 416.000 343.000 330.00035 5341 Alicante 8299 1198.000 1531.000 235.000 0.00036 5772 Alicante 7732 709.000 689.000 584.000 419.00037 6127 Alicante 7711 225.000 125.000 175.000 71.00038 6132 Alicante 7711 1084.000 634.000 592.000 291.00039 6133 Alicante 7711 5143.000 5838.000 7136.000 7835.00040 6134 Alicante 7711 327.000 394.000 442.000 331.00041 6138 Alicante 7739 301.000 341.000 408.000 470.00042 6140 Alicante 7732 274.000 483.000 468.000 280.000...56 337 Valencia 7911 261.000 623.000 1302.000 1487.00057 772 Valencia 7911 163.000 178.000 192.000 127.00058 1082 Valencia 7739 399.000 528.000 778.000 758.00059 1628 Valencia 8130 269.000 353.000 396.000 486.00060 1645 Valencia 7739 848.000 1015.000 1217.000 876.00061 1960 Valencia 8110 168.000 180.000 205.000 235.00062 2415 Valencia 7911 221.000 139.000 110.000 184.00063 2416 Valencia 7911 262.000 265.000 296.000 207.00064 2446 Valencia 8299 44.000 193.000 337.000 415.00065 3485 Valencia 8122 880.000 726.000 774.000 1081.00066 4053 Valencia 8130 243.000 649.000 880.000 734.00067 5566 Valencia 8299 325.000 767.000 97.000 158.00068 5779 Valencia 7739 298.000 356.000 344.000 348.00069 6167 Valencia 7711 2483.000 2461.000 2363.000 1915.000.... This is the code that you have told meresult <- vector("list", ncol(datos))for( i in 2:ncol(datos) {result [i-1] <- bartlett.test(datos[,i]~datos$Provincia))} This is what I get Error en `[.data.frame`(DATOS, , i) : undefined columns selectedAdemás: Mensajes de aviso perdidos1: In result[i - 1] <- bartlett.test(DATOS[, i] ~ DATOS$Provincia) : número de items para para sustituir no es un múltiplo de la longitud del reemplazo2: In result[i - 1] <- bartlett.test(DATOS[, i] ~ DATOS$Provincia) : número de items para para sustituir no es un múltiplo de la longitud del reemplazo3: In result[i - 1] <- bartlett.test(DATOS[, i] ~ DATOS$Provincia) : número de items para para sustituir no es un múltiplo de la longitud del reemplazo4: In result[i - 1] <- bartlett.test(DATOS[, i] ~ DATOS$Provincia) : número de items para para sustituir no es un múltiplo de la longitud del reemplazo5: In result[i - 1] <- bartlett.test(DATOS[, i] ~ DATOS$Provincia) : número de items para para sustituir no es un múltiplo de la longitud del reemplazo6: In result[i - 1] <- bartlett.test(DATOS[, i] ~ DATOS$Provincia) : número de items para para sustituir no es un múltiplo de la longitud del reemplazo7: In result[i - 1] <- bartlett.test(DATOS[, i] ~ DATOS$Provincia) : número de items para para sustituir no es un múltiplo de la longitud del reemplazo (Sorry, I have R in Spanish, número de items para para sustituir no es un múltiplo de la longitud del reemplazo= number of items to replace is not a multiple of replacement length) SO, I have changed a bit your code : result <- vector("list", ncol(datos))for( i in 2:ncol(datos)+1) {result [i-1] <- bartlett.test(datos[,i]~datos$Provincia))}And now R answer is: Error en model.frame.default(formula = datos[, i] ~ datos ~ Provincia) : object is not a matrix. What can I do? Is it necessary that datos works as a matrix? Thanks in advance [[alternative HTML version deleted]]
______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.