Page 1 of 1

Closed Models, error from RMark example

PostPosted: Thu Jun 23, 2016 2:27 am
by ctlamb
Hello,

I tried running the script from the
Code: Select all
?edwards.eberhardt
example, but the code fails, returning:

STOP ERROR
MARK did not run properly. If error message was not shown, re-run MARK with invisible=FALSE

Model output is not available

Show Traceback

Rerun with Debug
Error in `colnames<-`(`*tmp*`, value = 1:18) :
attempt to set 'colnames' on an object with less than two dimensions




I tried to reduce the code down and found that the huggins models work fine, but the models with N in the likelihood failed:

Code: Select all


#
data(edwards.eberhardt)
#

  #
  #  Define parameter models
  #
  pdotshared=list(formula=~1,share=TRUE)
  #
  # Run assortment of models
  #
  #
  #   Capture Closed models
  #
  #  constant p=c
 ee.closed.m0=mark(edwards.eberhardt,model="Closed",
                 model.parameters=list(p=pdotshared), invisible=FALSE)


STOP ERROR
MARK did not run properly.  If error message was not shown, re-run MARK with invisible=FALSE

Model output is not available

Error in `colnames<-`(`*tmp*`, value = 1:18) :
  attempt to set 'colnames' on an object with less than two dimensions




Any help on how to get past the errors?

Re: Closed Models, error from RMark example

PostPosted: Thu Jun 23, 2016 11:10 am
by jlaake
MARK changed over from N to f0 awhile back. It seems like you have incompatible versions of MARK and RMark. I suggest you update. Look at the file parameters.txt in the RMark subdirectory of your R library. If it lists f0 for the closed model then it is fine and Mark needs to be updated. If it says N then RMark needs to be updated. The examples ran fine on my machine.

Re: Closed Models, error from RMark example

PostPosted: Thu Jun 23, 2016 8:27 pm
by ctlamb
Ah, yes, that is the fix. Thanks, Jeff!