루프에 문제가 있습니다

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") 
} 

이 루프의 수정을 시도하지만 여전히 작동하지 않습니다 (Rstudio에서 작동 중입니다) - 어느 부분이 잘못 되었습니까?

더 복잡한 패턴 (이름 끝에 두 개의 숫자가 들어있는 파일 만 나열)을 사용하고 싶습니다. 그러나 다음과 같이 추가하십시오. pattern="_??.asc"이 작동하지 않는 것 같습니다.

나는 하나의 루프를 추가하여 (my_asc 변수에 수동으로 디렉토리를 삽입하는 대신) 폴더 목록을 얻고 싶지만 아이디어를 얻지 못했습니까? 평균, 최대, 최소 및 표준 편차 값에 대한 벡터를 만드는 내 방식이 작동하지 않는 이유는 무엇입니까 ...

Answer 1

시작 위치. . .

평균, 최대 및 최소 각각에 대해 na.rm = TRUE를 지정하고 sd에 대해 TRUE를 올바르게 입력해야합니다.

hist (z)는 png (파일 이름, ...) 다음에 와야하며 그 뒤에 dev.off() (적어도)가 있어야합니다.

R에 "+"가 붙은 문자열을 붙여 넣을 수 없으면 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() 
}