Hi,
First of all, thank you for this package and your continuous support.
My question is simple, but I'll provide some context in case you have any further advice to give (skip to the red text below if you want to save time).
I have a question relating to a multi-year fish presence/absence dataset. I have 39 years of data collected over the management area of a large city that has ~30 river basins. Think of it a bit like a citizen science database or something like that. The data was collected using every method known to fish ecologists, but the methods were always effective ones. Sampling was not spatially or temporally consistent within river basins, but non the less I've combined all observations within each river basin for each year. In the dynamic occupancy model I am using each year of data as secondary surveys (replicates) and three years are combined to make one primary sampling occasion (13 primary sampling occasions over 39 years).
The species I'm interested in (one species at a time) are migratory and spawn in response to Autumn and Spring freshes (small floods). Flows in most of the river basins are regulated and historically the magnitude of Autumn and Spring freshes has been impacted. Efforts have been made to rectify this situation over the last couple of decades though and the expectation is that these species have benefited.
I'd like to see if the percentage of days where Autumn and Spring freshes have occurred in each primary sampling period relate to changes in each species' occupancy.
It's all a bit rough around the edges, but I think it could work well enough if the patterns are strong and the interpretation is appropriate.
My question relates to the detection probability. As the data was collected with multiple methods, at multiple times of the year, from multiple people, I don't want to relate it to any environmental variables. I just want it to vary between years without relating it to a observation covariate. How do I write that in the equation? If ~1 is constant between years and ~observation covariate has it varying between years and relating the variation to the variable, what is the formula for just varying between years?
I've provided a basic example below..
dynamic_occ_m1 <- colext(
# Psi depends on initial Autumn flow conditions
psiformula = ~AutumnFlow1,
# colonization depends on change in Autumn flow conditions
gammaformula = ~autumflow_time,
# extinction
depends on change in Autumn flow conditions
epsilonformula = ~autumflow_time,
# I don't want detection to depend on anything. I want it to vary through time as per it's base calculation
pformula = ~?,
# data must be a unmarkedMultFrame
data = Tupong_unmarkedMultFrame,
# method is optim method, leave as "BFGS"
method = "BFGS")
Thanks in advance for your time. It's much appreciated.
Regards,
James