时间戳

首页 > 解决方案 > 如何在 R 中使用“'try”来跳过 for 循环中的错误

r - 如何在 R 中使用“'try”来跳过 for 循环中的错误

问题描述

每个人。

我想问一些有关如何避免 for 循环中可能出现的错误消息的问题。我正在运行模拟，在某些重复中可能会出现错误，导致崩溃并停止整个过程，我想要的是简单地自动进入下一次迭代并从我的结果中丢弃“坏”重复。

The code is as below.

    >    for(j in 1:100){ 
    > p<-YHatMulti(theta,X)   
    >  z=dim(p)[1] 
    >  A[[1,1]]<-matrix(apply(apply(cbind(p,X,((1+m-1):(z+m-1))),1,function(x){-x[1]*(1-x[1])*x[3:6]%*%t(x[3:6])*lambda^(((z+m-1)-x[7])/(1000-1)^r)}),1,sum),4,4)
    > A[[2,2]]<-matrix(apply(apply(cbind(p,X,((1+m-1):(z+m-1))),1,function(x){-x[2]*(1-x[2])*x[3:6]%*%t(x[3:6])*lambda^(((z+m-1)-x[7])/(1000-1)^r)}),1,sum),4,4)
    > A[[1,2]]<-matrix(apply(apply(cbind(p,X,((1+m-1):(z+m-1))),1,function(x){x[1]*(x[2])*x[3:6]%*%t(x[3:6])*lambda^(((z+m-1)-x[7])/(1000-1)^r)}),1,sum),4,4)
    > A[[2,1]]<-A[[1,2]]                                                    
    > A_rows<-lapply(1:2,bind_A_row)   Amatrix<-do.call(cbind,A_rows)  
    > Amatrix[is.na(Amatrix)]<-0  
    > B<--apply(cbind(t(sweep(t(X),2,YMat[,1]-p[,1],'*')),t(sweep(t(X),2,YMat[,2]-p[,2],'*')))*(lambda^(((z+m-1)-(1+m-1):(z+m-1))/1000^r)),2,sum)
    > B<-as.matrix(B)   B[is.na(B)]<-0   
    > #B<-B+t(Amatrix%*%(as.matrix(as.vector(t(theta)))))
    >  theta<-as.vector(t(theta))+t(ginv(Amatrix)%*%as.matrix(B))  
    > theta<-matrix(data = theta, ncol= ncol(X), nrow  = typeY-1, byrow =
    > TRUE) 
     }

最后第三行使用“ginv”。我有时发现，矩阵 B 会有 NA 元素。这样循环就结束了。

有人可以帮我使用“尝试”来跳过这个错误吗？因为我在做模拟，所以我可以跳过这个错误并进行下一个模拟。非常感谢。

标签： rloops

解决方案

如果您收到错误消息，我认为您可以使用该tryCatch功能。您可以尝试以下方法：

theta <- tryCatch({ as.vector(t(theta))+t(ginv(Amatrix)%*%as.matrix(B)) },
                    error = function(e){ next })

推荐阅读