Estimada comunidad, escribi el sgte script en R para clasificar las
filas de una planilla y no me funciona:
filenames <- list.files(path =
"/home/neo/Documents/Dropbox2/Tesis/4tesis/2.objesp/experimento/expnov/4fames/")
i <- 1
for (i in 1:length(filenames))
{
tmp <- read.csv(filenames[i], header=FALSE, sep="\t")
j <- 1
while (tmp[j,3] != max(tmp[,3])) { j <- j + 1 }
k <- j
#k <- 1
for (k in j:nrow(tmp))
{
if (tmp[k,6] == 1)
{
tmp[,7] <- 16
tmp[,8] <- 0
tmp[,9] <- "s"
}
if ((1.078 <= tmp[k,6]) & (tmp[k,6] <= 1.081))
{
tmp[,7] <- 18
tmp[,8] <- 0
tmp[,9] <- "s"
}
k <- k + 1
}
write.table(tmp,
paste("/home/neo/Documents/Dropbox2/Tesis/4tesis/2.objesp/experimento/expnov/4fames/",
filenames[i], sep=""), sep="\t", row.names = FALSE, col.names=FALSE)
i <- i+1
}
no entiendo bien xq no funciona ya que en otra parte tengo uno casi
identico que hace lo mismo y funciona bien.
Parece que el error esta en el segundo IF, pero no se bien como
escribirlo en una sola linea, pues no quiero usar IF ELSE ...
A veces me ocurre q el tipo de dato se cambia por alguna razon q
desconozco y entonces la clasificacion no funciona por eso, pero se
arregla usando as.numeric(9 o as.character() segun corresponda, pero en
este caso no me ha funcionado (lo use pero no lo puse en este script).
Adjunto un archivo simple de datos para q se pueda reproducir el script,
a ver si alguien me hecha una mano,
Muchas gracias,
eric.
--
Forest Engineer
Master in Environmental and Natural Resource Economics
Ph.D. student in Sciences of Natural Resources at La Frontera University
Member in AguaDeTemu2030, citizen movement for Temuco with green city
standards for living
Nota: Las tildes se han omitido para asegurar compatibilidad con algunos
lectores de correo.
1 11.4702415466309 0.0418207980692387 0.2470912784338 0.000148864927808787 NA 18 0 "s"
2 13.4479007720947 1.90983533859253 0.470895916223526 0.00679823228948302 NA 18 0 "s"
3 13.575249671936 1.7584000825882 0.462930619716644 0.00625918474630924 NA 18 0 "s"
4 16.8946704864502 2.30076265335083 0.646402716636658 0.00818977356025525 NA 18 0 "s"
5 17.0956058502197 5.80301666259766 0.872533559799194 0.0206563646901374 NA 18 0 "s"
6 17.2105121612549 5.56320810317993 1.41665434837341 0.0198027443496899 NA 18 0 "s"
7 17.3523559570312 4.87942123413086 1.12816715240479 0.0173687429055028 NA 18 0 "s"
8 17.4699363708496 2.37469887733459 0.763153433799744 0.00845295626249767 NA 18 0 "s"
9 17.6279544830322 0.986213862895966 0.349221140146255 0.00351051778736909 NA 18 0 "s"
10 17.7760181427002 2.29407668113708 0.866581737995148 0.00816597423511364 NA 18 0 "s"
11 19.1464252471924 0.932058095932007 0.335037440061569 0.00331774541783726 NA 18 0 "s"
12 19.2906913757324 0.977343797683716 0.354994148015976 0.00347894398489658 NA 18 0 "s"
13 19.6240005493164 0.308019518852234 0.180336236953735 0.00109642344370665 NA 18 0 "s"
14 19.8428802490234 1.14335989952087 0.462329685688019 0.00406989337266693 NA 18 0 "s"
15 20.1201648712158 1.21129238605499 0.435258090496063 0.00431170522635347 NA 18 0 "s"
16 22.3271789550781 2.09697031974792 0.72647362947464 0.00746435624565608 NA 18 0 "s"
17 24.3413333892822 160.024963378906 52.9973831176758 0.569623386468252 1 18 0 "s"
18 26.2848339080811 86.3257064819336 22.9554481506348 0.307284190086465 1.07984363418869 18 0 "s"
_______________________________________________
R-help-es mailing list
[email protected]
https://stat.ethz.ch/mailman/listinfo/r-help-es
