Todays post is mainly just a summing up of model results following some recent analysis that has shook the foundations of my trust in my model.
As it turns out, when you have a complex model, analysing results can be complex and assigning blame to any particular parameter/process can be even more challenging.
I've been happily plodding along in the knowledge that I finally have my model churning out some decent results, with the right things happening at the right time in the right places. One final analysis has turned this on its head with the finding that territory sizes are not properly linked to food resources. I have suspected this has been the case for some time, but following the reduction of food resources outside of the protected area I find that male territories in the matrix are about 3/4 the size of territories inside the CBWS forest. Females surprisingly are the reverse, but I think this is mainly due to most females being driven out of CBWS due to avoidance of males. In the matrix it seems they are able to expand their territories and are not so constrained due to the presence of fewer males. The average size of their territories is also unfeasibly small inside CBWS and nowhere near the size that would be needed in the real world to maintain a single individual....
crap.
So it seems a pretty fundamental feature of my model is not working as it should and I am at a pretty big loss as to why this might be the case.
Looking at males, they are attracted to trails, and attracted to male pheromones on trails. This causes larger territory sizes as they move along trails, pop out further into the forest and create new areas to call their own. Individuals are seemingly not linked to any particular area and territories are completely transitory. I wonder whether this is the problem? Males, due to this feature, also exhibit larger territory sizes when there are more males around, presumably because they get drawn onto trails. However, having a process that allows males to move further along trails causes a population crash in pretty short time (i.e. around year 20 of the simulation). I'm not quite sure why, but I'm doing some tests to figure out if the number of moves makes any difference to the fate of the population.
Looks like I might have to go back to the drawing board.
This is not the news I wanted with 6 months left on my PhD!
What do I want:
Males move along trails where possible
Males follow other males along trails, for some suitable distance.
All individuals have larger territories outside of the protected area as they need to move more to consume sufficient food resources
Females avoid males but generally their territories are overlapped to a large extend by males
Males move much faster along trails
Males move much faster when following a female to mate
All individuals would prefer to be inside the protected area
Hi Angela,
ReplyDeleteIn some of your previous blog posts you have noted that it takes quite some time to run a simulation. In my experience, when simulations take a long time to run, bugs tend to go unnoticed much, much longer. Would it be possible for you to run a scaled down version of your model (less agents, smaller, simplified area)? Could you get the run time down to a couple minutes or less, and run simulations in simplified environments? This could help in three ways. First just by being able to get instant feedback, rather than having to wait several hours/days to get simulation results, you run way more diagnostics. Secondly, just knowing that you won’t have to wait days for the results will give you the freedom to follow a random hunch, knowing that if you are wrong you’ve only wasted a few minutes. Finally, by running your simulations in simplified environments, you will probably have a better intuition of how your model should behave, and therefore be more likely to identify when something is amiss. For example, make all that patches the same habitat type. What happens if you increase/decrease the food availability of that habitat type? Is it what you expected?
On a more specific note, I noticed you are working with a Type II functional response with a logistically growing resource. These types of systems are inherently unstable, especially if that half-saturation coefficient is much lower than the maximum prey density (paradox of enrichment). I haven’t gone through how the spatial or territorial dynamics of your model work, but in the absence of other mechanisms your territorial dynamics are going to have to stabilize the inherently unstable Rosenzweig-MacArthur model. One quick way to check if this is responsible for the instability would be to make your resource dynamics semi-hemostat rather than logistic dR/dt = r*(1-R/K).
My apologies if you have already tried theses things. Good luck!
popecology - thanks so much for the comment. I hadn't actually thought of implementing either of these ideas, but they have given me food for thought so I hugely appreciate the input. Its great that someone can sympathize with these problems and give some constructive and useful advice!..
Delete