Hello,
I am planning to use occupancy robust design models in RMark to estimate colonisation and extinction of patches within an experiment, but have some limitations with missing data, so I would appreciate some advice before starting the analysis. Briefly, the experimental design includes 4 experimental treatments with 24 replicate points per treatment, sampled every spring and autumn (with 5 visits in every sampling season) during 8 years. I plan to use treatment as grouping design covariate, and code year and season in the design matrix as covariates because they are of interest.
Now, I see two problems with the database that could affect estimation: i) the location of the points sampled for every treatment changed after 4 years for logistic reasons, so that we have data for one set of 24 points per treat for the first 3 years and another different set of 24 points per treat for the last 5 years; ii) we missed one season where no points at all were sampled in the 3rd year of sampling. So an ultra-summarized version of the dataset could look like:
10110 11011 11111 01001 ..... 11011 ..... ..... ..... ..... ..... ..... ..... ..... ..... .....
..... ..... ..... ..... ..... ..... 10101 11011 11111 01001 10110 11011 11111 01001 11111 01001
(line above first set of points, line below second set of points)
My question is whether RMark (or MARK) is able to deal with that kind of data, provided I am interested in treatment (and its interaction with covariates) effects, not in actual point estimates. As a less preferred alternative, I could analyse the first 3 years on one side, and the last 5 years indepedently (but still will have one missing occasion).
I would appreciate your views thanks. Also I hope this is the right forum to ask (a little bit confused since strong alternatives were MARK and PRESENCE forums!).
Thanks,
N