SECR error message

[Cheryl Lohr]

OFFICIAL

Dear Murray,

Hopefully a very simple question: but this is an error message I've not encountered before.

> fit5<-secr.fit(BLwarren, mask = m1, detectfn = 'HR',
+                model = list(sigma~g), groups = 'Sex')
Checking data 
Preparing detection design matrices 
Preparing density design matrix 
Finding initial parameter values... 
Initial values  D = 0.25578, g0 = 0.03345, sigma = 347.85941 
Maximizing likelihood... 
Eval     Loglik              D             g0          sigma       sigma.gM sigma.gUnknown              z 
Error in .localstuff$Eng[sessnum, g] <- N * meanpdot : 
  number of items to replace is not a multiple of replacement length

All previous models in this set running very smoothly. No errors in input data. The only thing that is different is this model includes an individual grouping covariate. I'm using the same mask and detectfn. Never had a problem with this script before today. Is there a recent update issue I need to deal with?

Thanks,

Cheryl

Cheryl Lohr

Senior Research Scientist, Animal Science Program

DBCA Biodiversity and Conservation Science

Location: 37 Wildlife Pl, Woodvale, WA 6026

Mail: Woodvale Wildlife Research Centre

Locked Bag 104 Bentley Delivery Centre

WA 6983

Mob: 0407335004

 

The Department acknowledges all Aboriginal people as the traditional owners of the land and waters it manages throughout Western Australia. We pay our respects to them, their culture and to their Elders past and present. Our office is located on the lands of the Whadjuk Noongar People.

 

 

---

This message is confidential and is intended for the recipient named above. If you are not the intended recipient, you must not disclose, use or copy the message or any part of it. If you received this message in error, please notify the sender immediately by replying to this message, then delete it from your system.

[Murray Efford]

Hello Cheryl

I can't put my finger on the cause right away and I'm away from my desk for a few days. It might be informative to try the alternative conditional likelihood relative density model 

fit6 <- secr.fit(BLwarren, mask = m1, detectfn = 'HR', model = list(sigma~Sex), CL=TRUE)

(not really documented yet). I'll check this out when I get back.

Murray