by monkr » Fri Mar 29, 2019 1:51 pm
I seem to be having problems with the newest update.
1. Although the model from the multi-scale exercise runs without error, the results are different from what I previously had (= those quoted on the exercise). Psi is the same (0.415), but now theta & p are different. They should be :
## est se lower upper
## unit1_1 0.1772 0.0576 0.0903 0.3185
## est se lower upper
## unit1 _ 1-1 0.3683 0.0567 0.2656 0.4845
Instead I get:
> unique(fitted(m1,"theta"))
est se lower upper
unit1_1 0.1515 0.0473 0.0798 0.2686
> unique(fitted(m1,"p"))
est se lower upper
p1(1-1)_1_unit1 0.5 0 0.5 0.5
2. The following error occurs with models in the same exercise with p~week, it doesn't occur with p~survey, e.g,
m2<-occMod(model=list(psi~1,theta~1,p~week),data=camera,type="so.mm")
Error in VC[p.coeff[, "index"], p.coeff[, "index"]] :
only 0's may be mixed with negative subscripts
> traceback()
2: occMod_mm(psi = e$psi, psi.cov = e$psi.cov, theta = e$theta,
theta.cov = e$theta.cov, p = e$p, p.cov = e$p.cov, gamma = e$gamma,
gamma.cov = e$gamma.cov, epsilon = e$epsilon, epsilon.cov = e$epsilon.cov,
modname = modname, paoname = paoname, outfile = outfile,
fixed = fixed, initvals = initvals, data = data, VCoutopt = VCoutopt,
noDerived = noDerived, randinit = randinit)
1: occMod(model = list(psi ~ 1, theta ~ 1, p ~ week), data = camera,
type = "so.mm")
I've uninstalled and reinstalled 2.12.27, restarted R and get the same result.
One other minor issue: The output to the console states the version as 2.12.26, but I'm running 2.12.27.