0

我正在编写一个函数，它打开 dir() 中的所有文件，进行一些计算和重新排列，然后在 dir() 中再次写入（write.table）每个文件

dir()
 [1] "v2p1"  "v2p10" "v2p11" "v2p12" "v2p2"  "v2p3"  "v2p4"  "v2p5"  "v2p6"  "v2p7" 
[11] "v2p8"  "v2p9"  "v2r1"  "v2r10" "v2r11" "v2r12" "v2r2"  "v2r3"  "v2r4"  "v2r5" 
[21] "v2r6"  "v2r7"  "v2r8"  "v2r9"  "v2s1"  "v2s10" "v2s11" "v2s12" "v2s2"  "v2s3" 
[31] "v2s4"  "v2s5"  "v2s6"  "v2s7"  "v2s8"  "v2s9"

代码看起来像这样，除了循环不起作用之外，它做的一切都很好。

mysub <- function(x) {sub(",",".",x)}
abrir<-function(etapa){  
### volnum c(1:10) indica en qué vol quiero,
### etapa c("p","r","s") indica en qué etapa es del experimento

i<-1
for (i in dir()){

##### ABRIR ARCHIVOS DESDE DIR() ######
Archivo<-paste(getwd(),"/",i,sep="")
r1<-read.table(Archivo,fill=T,sep=",")

##### CALCULAR LA VELOCIDAD DE CADA CORRIDA ######
r11 <- (apply(r1, 2, mysub ))
r1<- data.frame(apply(r11, 2, as.numeric)) ## ojo que pueden aparecer NA
R1<-r1[1:3136,]     ### OJO, DEPENDE LA CORRIDA, SACA LOS VALORES DEL FINAL QUE NO SIRVEN
V<-data.frame(rbind(c(0,0,0,0),(abs(diff(as.matrix(R1[,2:5]))))))
Tiempo<-R1[,1]
Velocidad<-data.frame(Tiempo,V)  

##### GUARDAR LAS VELOCIDADES QUE CORRESPONDEN  ######

if (grep("p",i)==TRUE){
salida<-paste("VEL",i,".txt",sep="")
write.table(Velocidad,salida,sep="\t",row.names=FALSE)
} else if (grep("r",i)==TRUE){
salida<-paste("VEL",i,".txt",sep="")
write.table(Velocidad,salida,sep="\t",row.names=FALSE)
}else if(grep("s",i)==TRUE){
salida<-paste("VEL",i,".txt",sep="")
write.table(Velocidad,salida,sep="\t",row.names=FALSE)
}else stop("etapa inválida")
i<-i+1
}
}

abrir() 为我保存了正确的“VELv2p1.txt”并返回此错误 i + 1 中的错误：二进制运算符的非数字参数

我不知道如何解释 R 完成后继续使用第二个元素 i=1 我猜 i<-i+1 这不是这样。我会很感激你的帮助。

4

1 回答 1

1

错误消息与i循环中是一个character对象而不是numeric. 该语法for(i in dir())产生实际的文件i，而不是索引。在for这样的使用中，不需要设置i到下一个项目，for循环将采取。因此，只需删除即可i <- i + 1解决您的问题。

关于您的代码的一些说明：

  • 缩进将提高代码的可读性。就像添加空格一样。
  • 您可以在函数末尾调用 write.table，从而节省三行相同的代码。实际上，这组if语句似乎都执行相同的代码。
  • for循环通常不是 R 中的最佳选择。矢量化或apply基于解决方案产生更短，有时更快的代码。

我大致采取以下方法：

list_of_files = list.files('.')
do_some_calculations = function(file, output_name) {
    output_name = sprintf('VEL%d.txt', list_of_files)
    Archivo<-paste(getwd(), "/", i, sep = "")
    r1<-read.table(Archivo,fill = T, sep = ",")
    r11 <- (apply(r1, 2, mysub ))
    r1<- data.frame(apply(r11, 2, as.numeric)) ## ojo que pueden aparecer NA
    R1<-r1[1:3136,]     ### OJO, DEPENDE LA CORRIDA, SACA LOS VALORES DEL FINAL QUE NO SIRVEN
    V <- data.frame(rbind(c(0,0,0,0), (abs(diff(as.matrix(R1[,2:5]))))))
    Tiempo <- R1[,1]
    Velocidad <- data.frame(Tiempo, V)  
    write.table(Velocidad, output_name,s ep="\t", row.names=FALSE)
}
于 2014-02-15T18:21:56.583 回答