MATA: Error in vcv.real[result$par.index[i], result$par.inde

posts related to the RMark library, which may not be of general interest to users of 'classic' MARK

MATA: Error in vcv.real[result$par.index[i], result$par.inde

Postby Eurycea » Tue Sep 01, 2020 3:25 pm

I was wondering if you could help decipher an error message. I am attempting to model average using MATA with the following code:

Code: Select all
model.average(final_models,parameter="GammaPrime",vcv=T,mata=T)


And it throws the error:

Code: Select all
Error in vcv.real[result$par.index[i], result$par.index[i]] :
  subscript out of bounds


The error only occurs when both vcv and mata are TRUE. If vcv=F it works fine, and if mata=F and vcv=T it works fine.

I get the same error for GammaDoublePrime. However, this works fine for S.

Thanks!
Eurycea
 
Posts: 90
Joined: Thu Feb 25, 2010 11:21 am

Re: MATA: Error in vcv.real[result$par.index[i], result$par.

Postby jlaake » Tue Sep 01, 2020 3:29 pm

Obviously a bug. Please use save function in R to save final_models into a .rda file and send to me and I'll have to debug it. Have you been using the argument external=TRUE? If so, I need all the .rda model files as well. My email address is jefflaake@gmail.com.

--jeff
jlaake
 
Posts: 1197
Joined: Fri May 12, 2006 12:50 pm
Location: Escondido, CA

Re: MATA: Error in vcv.real[result$par.index[i], result$par.

Postby Eurycea » Tue Sep 22, 2020 2:31 pm

Upon further exploration, I have also tried the covariate.predictions function. With mata=TRUE, I will sometimes get upper confidence limits above 1.00. Was wondering why my error bars weren't plotting correctly... :D Not sure if this is a feature of MATA (doesn't seem like it after reading the paper) or a bug. Thanks!
Eurycea
 
Posts: 90
Joined: Thu Feb 25, 2010 11:21 am


Return to RMark

Who is online

Users browsing this forum: No registered users and 1 guest