search for: totr

...| D_2), family = tweedie(var.power=1,link.power=0), offset=log(Dt), data=data) ] data=read.csv(file="GLMM_data.csv") > models<-read.table("GLMM_tweedie.txt", sep="\t") > data$Ggrec_Dtlog = log(data$Ggrec_Dt+1) > models<-as.vector(models[,1]) > totres=c() > for (i in 1:79) {model=models[i] + res=eval(parse(text=model)) + res=AIC(logLik(x)) + res=as.vector(res) + totres=rbind(totres,res)} The output would then be just a list of all the AIC of each model. For 1 of the models (the one in the [] above) I'm getting the following error me...

...y, my problem seems to be solved by writing a FORTRAN wrapper for the Fortran code! (As long as the check doesn't get smarter...). This is the relevant part of my Fortran code: ----------------------------------------------------------- subroutine gmlfun(what, & totevent, totrs, ns, & antrs, antevents, size, & totsize, eventset, riskset, & nn, antcov, covar, offset, & beta, gamma, & loglik, h1, h2, h11, h21, h22, & score) ...... call gmlfun1(what, & totevent, totrs,...

Sorry for cross-posting, but I realized my question might be more appropriate for r-devel... Thank you, Giovanni ________________________________________ From: R-help <r-help-bounces at r-project.org> on behalf of Giovanni Petris <gpetris at uark.edu> Sent: Tuesday, September 10, 2019 16:44 To: r-help at r-project.org Subject: [R] Calling a LAPACK subroutine from R Hello R-helpers!