Hi All, I have tried every fix on my try or tryCatch that I have found on the internet, but so far have not been able to get my R code to continue with the "for loop" after the lmer model results in an error.
Here is two attemps of my code, the input is a 3D array file, but really any function would do.... metatrialstry<-function(mydata){ a<-matrix(data=NA, nrow=dim(mydata)[3], ncol=5) #colnames(a)=c("sens", "spec", "corr", "sens_se", "spec_se", "counter")#colnames(a)=c("sens", "spec", "corr", "sens_se", "spec_se") k=1 for(ii in 1:dim(mydata)[3]){ tmp<-mydata[,,ii] tmp1<-as.data.frame(tmp) names(tmp1)=c("persons", "d1", "tp", "fn", "fp", "fn", "detect", "d0", "outcome") lm1<-try(lmer(outcome~0+d1+d0+(0+d1+d0 | persons), family=binomial, data=tmp1, nAGQ=3), silent=T) if(class(lm1)[1]!='try-error'){ a[ii,1]=lm1@fixef[1] a[ii,2]=lm1@fixef[2] a[ii,3]=vcov(lm1)[1,2]/prod(sqrt(diag(vcov(lm1)))) a[ii,4:5]=sqrt(diag(vcov(lm1))) } } #k=k+1 #a[ii,6]=k return(a) } ######################################### ######### try / try catch ############### ######################################### metatrialstry<-function(mydata){ a<-matrix(data=NA, nrow=dim(mydata)[3], ncol=5) #colnames(a)=c("sens", "spec", "corr", "sens_se", "spec_se", "counter") colnames(a)=c("sens", "spec", "corr", "sens_se", "spec_se") #a[,6]=rep(0, length(a[,6])) for(ii in 1:dim(mydata)[3]){ tmp<-mydata[,,ii] tmp1<-as.data.frame(tmp) names(tmp1)=c("persons", "d1", "tp", "fn", "fp", "fn", "detect", "d0", "outcome") lm1<-tryCatch(lmer(outcome~0+d1+d0+(0+d1+d0 | persons), family=binomial, data=tmp1, nAGQ=3), error=function(e) e) a[ii,1]=lm1@fixef[1] a[ii,2]=lm1@fixef[2] a[ii,3]=vcov(lm1)[1,2]/prod(sqrt(diag(vcov(lm1)))) a[ii,4:5]=sqrt(diag(vcov(lm1))) } return(a) } Any guidance would be greatly appreciated... thanks! Lisa -- Lisa Wang email: wang.li...@gmail.com cell: +49 -0176-87786557 Tübingen, Germany, 72070 [[alternative HTML version deleted]]
______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.