problème avec des boucles dans R

Question

library(rgdal) 
my_asc = dir("~/Pulpit/dods/karol/TVDI113_121", 
      pattern=".asc", recursive=TRUE, full.names=TRUE) 
for (i in 1:length(my_asc)) { 
    r <- readGDAL(my_asc[i]) 
    z <- as.matrix(r) 
    vectordata[i] <- mean(z) 
    vectordatamax[i] <- max(z) 
    vectordatamin[i] <- min(z) 
    vectordev[i] <- sd(z, na.rm=TRUE) 
    hist(z) 
    png(filename="hist"+tostring(i)+".png") 
} 

Je tente d'apporter quelques modifications à cette boucle, mais cela ne fonctionne toujours pas (je travaille sous Rstudio) - quel fragment est incorrect?

Je voudrais aussi utiliser un modèle plus compliqué (pour lister uniquement les fichiers qui contiennent à la fin de ce nom deux nombres), mais ajouter quelque chose comme: pattern="_??.asc" ne semble pas fonctionner.

Je voudrais ajouter une boucle supplémentaire pour obtenir la liste des dossiers (au lieu d'insérer manuellement des répertoires dans la variable my_asc), mais je n'ai pas Idea comment je peux le faire? Je ne sais pas, pourquoi ma façon de créer des vecteurs pour les valeurs moyennes, max, min et écart-type ne fonctionne pas ...

Answer 1

Par où commencer. . .

Vous voulez probablement na.rm = TRUE pour chacun des paramètres mean, max et min, et vous devrez entrer TRUE correctement pour sd. L'hist (z) devrait venir après le png (nom de fichier, ...) et cela devra être suivi d'un dev.off() (au moins).

Vous ne pouvez pas coller ensemble des chaînes avec "+" dans R, utilisez paste().

vectordata <- vectordatamax <- vectordatamin <- vectordev <- numeric(length(my_asc)) 
for (i in seq_along(my_asc)) { 
    r <- readGDAL(my_asc[i]) 
    ## as.matrix is not necessary, as the band values are accessible directly 
    ##z <- as.matrix(r) 
    z <- r[[1]] 
    vectordata[i] <- mean(z, na.rm=TRUE) 
    vectordatamax[i] <- max(z) 
    vectordatamin[i] <- min(z) 
    vectordev[i] <- sd(z, na.rm=TRUE) 
    png(filename=paste("hist", i, ".png", sep="")) 
    hist(z) 
    dev.off() 
}