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.

Reply via email to