Congratulations on maintaining such a friendly forum, no need to say that your input is much useful to us, the users
I’m running into a frustrating problem myself. I’d like to fit capture-recapture models within a robust design framework. To do so, the ‘robust’ help file of the RMark package (?robust) is a very nice starting point. However, I cannot reproduce the example.
More precisely, when I run this code:
- Code: Select all
library(RMark)
?robust
data(robust)
# Random emigration model, p=c varies by time and session, S by time
time.intervals=c(0, 1, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0)
S.time=list(formula=~time)
p.time.session=list(formula=~-1+session:time,share=TRUE)
GammaDoublePrime.random=list(formula=~time,share=TRUE)
model.1=mark(data = robust, model = "Robust",
time.intervals=time.intervals,
model.parameters=list(S=S.time,
GammaDoublePrime=GammaDoublePrime.random,p=p.time.session),threads=2)
I get the following error and warning messages in the R console (pardon My French...):
MARK did not run properly. If error message was not shown, re-run MARK with invisible=FALSE
Model output is not available
Error in wtable[i, cols] = real[pim[i, cols]] :
le nombre d'objets à remplacer n'est pas multiple de la taille du remplacement
De plus : Warning message:
In grep("Effective sample size ", out, ignore.case = TRUE) :
la chaîne de caractères entrée 101 est incorrecte dans cet environnement linguistique
I also checked the output file, which gave complementary information on my problem:
INPUT --- group=1 f0 Session 1 rows=1 cols=1 Square ;
ERROR -- Type of input matrix: Ì, S, p, r, c, N, Nbar, Nstar, U, F, F', Gamma, Gamma', Gamma'', R, R', f, M, Lambda, Psi, pent, pi, Omega, a'', a', Theta, Theta', or Delta was not specified.’
Interestingly (or not), I get the same error message for all models of the ‘robust’ example, except for the Huggins emigration model (the last one) which works perfectly well.
I found on the archives of the RMark forum a problem similar to mine. It was posted by Sergi in November 2014 with the title ‘Problems with the new version 2.1.9’. The issue was solved by Jeff who found out that it was an issue of version. The results above were obtained with R version 3.2.2 and RMark 2.1.13. I’m using a Mac.
Apologies in advance if this is a naive question, any help would be appreciated.
Cheers,
Olivier
