Tejada, J., Bosco, G. G., Morato, S., & Roque, A. C. (2010). Characterization of the rat exploratory behavior in the elevated plus-maze with Markov chains. Journal of Neuroscience Methods, 193(2), 288–295. http://doi.org/10.1016/j.jneumeth.2010.09.008
Summary
First-order Markov chain describing a rat behavior in the elevated plus-maze. It does well qualitatively, not so great quantitatively. It has a sample size of 63 rodents.
Notes
- Suggested that since both previous models simulated the rat using sequential probabilistic processes, they could also use Markov Chains.
- Markov Chains: the behavior of the system is described with a set of discrete states. Time is also discrete. The state in a given time is dependent on the previous k states and has a specific probability. A first-order Markov Chain depends only on the current state. It is also called “homogeneous” if the transition probability does not depend on time.
- Sought to define the rodent behavior using first-order homogeneous Markov chains.
- Sample of 63 rats. 27 in the control group, 26 on anxiolytic drugs, and 10 on anxiogenic.
- Divided the maze into 11 squares. One in the center, five in each arm. What is considered as two open arms by other researchers, here is a “longer” one open arm. The same applies to the closed arm. The idea is to remove the number of states and complexity from something that’s not the goal of the research.
- Auto-transitions – go to the same state from one time unit to the other – are forbidden, since the representation intends to model the transition between squares.
- Markov Chain states: 1-5 open arms, 7-11 closed arms, and 6 center.
- The transition is considered a Bernoulli trial: only two possible outcomes. The state can change to the state to the right – probability p – or to the left – probability p-1.
- Stochastic matrix (pij):
- pij = nij / nj, where i, j = 1, 2… 11, i ≠ j
- ni = {j=1, 11} ∑ nij (i ≠ j)
- Stationary probabilities (π):
- π = πP
- Stationary probabilities (π) can be calculated in two ways:
- Solve {i=1, 11} ∑ πi = 1; or
- Pn = P × P × P × P × … × P (n times)
- They used both and considered n=30 to give a good-enough estimate. So they used the second option.
- Results show that rats under anxiolytic drug don’t have a preference for a square. Suggests that such drugs smooth out the preference difference between fear and curiosity. “as in a symmetric random walk.”
- Authors affirmed that it was already possible to describe rat’s behavior using only first-order Markov Chains. They say that preference between open and closed arms, and their transitions, seem to be simple enough for first-order. But exploratory patterns may be more complex.
- Claimed that it’s possible to tell if an unknown drug has real anxiolytic or anxiogenic effects using the models.
- Qualitatively good, quantitatively not so much.
Thoughts
- Considering the scope of the paper, explains considerably well Markov Chains and how to interpret the results.
- Why using only first-order Markov Chains if the assumption – included in previous researches – is that the action depends on past actions, like an inertial movement? An experiment could be to try to reproduce it using more previous actions.
- Although much more understandable than last times I read about Markov Chains, it isn’t completely clear how to solve the matrices. Also, if there are other ways of representing or solving the chain problem.
- It could be interesting to also test with a more significant number of squares/states and allow auto-transitions, by measuring a transition in a time unit such as 1 second. It could improve results.
- Would it be possible to design a Markov Chain that also considers behavioral actions such as grooming?
- Based on my current readings, it is the research using the larger sample size so far: n=63. Would it have a specific reason for them to have a larger sample size, while others don’t have?