Christmas gifts from Just Simple Enough!

And, no, not the gift of homework. The gifts of song and movie!

Song. An original song written by a graduate student about graduate student-supervisor meetings. It’s a catchy tune! Click Here.

Movie. This is a movie that I took in 2009 while attending a Mathematical Biology Summer School in Botswana. Click here. (Unfortunately, I encountered some technical difficulties uploading this to YouTube, but it’s still watchable, albeit in micro-mini).

Happy holidays everyone!

Who’s ready to Rock?

Several faculty at MUN have photos of icebergs on their websites. Luckily, I was able to get a jump start on this during my recent visit.

In Atlantic Canada, which university has:

  1. The largest undergraduate enrollment,
  2. The most research funding, and
  3. A non-zero abundance of icebergs?

If you said Memorial University (MUN), then you win. You’ll find it all in St. John’s, Newfoundland (a.k.a. The Rock).

If you, or someone you know, are interested in graduate studies at MUN, please visit my website.

Snails shells: the logarithmic spiral

There’s a great post up on The Atavism by David Winter where he explains why the shape of the snail’s shell is a logarithmic spiral. What I find interesting is that the shell has the logarithmic spiral shape under two assumptions:

  • The radius of the shell increases exponentially as the snail ages (i.e., the rate of radial growth is proportional to the shell radius) and,
  • The angle between the centre of the snail and the end of the shell changes at a constant rate – like the second hand of a clock which moves at a constant rate of 360-degrees per minute.

Does anyone know what kind of snail this is? Photo credit: Robyn Hurford

I was hoping that we could use this simple model to do a couple of quick experiments.

Firstly, can we validate the hypothesis that snail shells are logarithmic spirals using something like a Turing test?

That is, among a set of spirals are blog-readers going to pick out the logarithmic spiral as being snail-like?

Secondly, since the art of model derivation is subjective, I want to solicit opinions on this particular model derivation – to see if everyone has the same instinctive appraisal of model assumptions or if there are a range of different tastes on the matter.

… and after you’ve voted be sure to go check out The Atavism for a nice explanation and some great snail pictures!

Overfitting – the elephant in the Cartesian plane

Just because it’s fun and because I’m still not on the other side of my blog holiday, here’s a nice picture of an elephant from Burnham and Anderson (2002) who got it from Wei (1975).

It’s all in response to John von Neumann who said,

“with four parameters I can fit an elephant, and with five I can make him wiggle his trunk.”

But that was more hyperbole than quantitatively accurate. Further investigations by Wei (1975) revealed that the 30-term elephant

“may not satisfy the third-grade art teacher, but would carry most chemical engineers into preliminary design.”


Burnham K, Anderson D. 2002. Model selection and multimodel inference. Springer, USA.

Wei J. 1975. Least square fitting of an elephant. Chemtech 5: 128–129.

Blog challenge: Sheep cyclone

In relation to multiscale modelling, a category of problem that receives a good amount of attention is understanding how collective (group) movement arises from movement rules defined at the individual level. For example, the direction of each individual’s next move might depend on where that individual’s nearest neighbours are. What type of rules are necessary for the group to stay together? If the group stays together, what does the collective movement look like (directed or meandering)?

Some nice examples are:

and more generally see here, here, here and here. So this is my question:

What type of individual movement rules are needed to produce a Sheep cyclone?

Bonus question: How close together do the parallel walls (relative to the car width) need to be?

Hit me with your best shot

Choosing and designating some models as Great has been causing me anxiety and so here’s a model that I’m going to write about just for fun. Thanks to Titia Praamsma for sending me this mathematical model of shot selection in basketball by Dr. Brian Skinner of the University of Minnesota. Skinner even acknowledges the parsimony issue:

While the complex nature of decision-making in basketball makes such a description seem prohibitively difficult, it is nonetheless natural as belonging to the class of “optimal stopping problems”…

Fig. 1. Occassionally, I have wanted to derive a mathematical model for basketball, however, the complex nature of the decision-making process in basketball makes it prohibitively difficult for me.

Now, I follow basketball fairly closely, yet I’ve never been able to come up with any good ideas for applicable mathematical models (see Fig. 1). After reading Skinner’s work I realize that this is because I was framing the question so that I would have to tackle it using complicated and possibly uninformative methodology.*

Let me explain. When I think about basketball, I think that as a coach your job is (partly) to invent set plays that will increase your team’s chance of getting a high quality shot. This problem is spatial, game theoretical (i.e., it depends on whether the other team is playing a zone or a man-to-man defense) and it’s probably necessary to coordinate the movement of all five players (i.e., n=1 or 2 is fundamentally a different problem and so studying these simple cases is likely uninstructive). I know these types of problems are hard and so it seemed like too much work for me to derive a mathematical model for basketball.

The above discussion is meant to illustrate that anything can seem complicated if you choose to look at it that way.**  After acknowledging the complexity of basketball, Skinner goes on to come up with a simple model of shot selection. I think that his conceptual approach is quite clever and I hope this is underscored by my admission that I had (naively) written basketball off as not something that was amenable to simple models.

Fig. 2. Tiffany Hayes (with the ball) of the University of Connecticut: dropping 35 points, shooting 11-for-15, and lending support to the controversial "hot hand" phenomenon (see Skinner (2012)).

Certainly, part of what Skinner does that’s clever is come up with a good question. The question is this:

Given a shot opportunity of quality, q, should a player take the shot?

The answer is that it’s a good idea to shoot when q > f where f is a threshold quality value that depends on n, the number of shot opportunities remaining in a possession and the article explains how n might depend on the team’s turnover rate, the existence of a shot clock and how fast the team moves the ball.

Interestingly, the article reports that a team that moves the ball well has a higher n, and therefore a higher f, and would want to execute more passes (until the shot opportunity q > f arises) than a team with poor ball movement (i.e., lower n and lower f). This is counter-intuitive because you’d think that if a team passes the ball quickly, then they can shoot sooner. This is true, but the result isn’t about whether a team with good ball movement can shoot earlier in the possession and win, it’s about what’s optimal: it’s optimal to expect a better shot to arise during a possession if shots are created at a faster rate – and so the team with good ball movement makes more passes and they might even wait longer to shoot.

Let’s revisit the model construction. You might say that if Brittany Griner gets the ball in the post then she should dunk it. That’s a good decision because this would be a high quality shot where q > f, but what I like about the model is that it abstracts away the name of the player, where the player is on the court, and what type of offense was run to generate the shot opportunity, and simply summarizes all this complexity into the one variable, q, shot quality.

Another aspect of this paper that I like is the comparison with data from the NBA. In fact, section 4. of Skinner’s paper is solely dedicated to recasting  f in terms of a shooting rate under the assumption that shot opportunities arise at a constant rate τ, so that the number of remaining shot opportunities n at time t is Poisson distributed. This is a nice final step because how would one ever know if q > f without these addition assumptions? My point is that q (shot quality) is not an especially useful quantity because how would one measure that? On the other hand, the shot rate can be estimated from a play-by-play box score which reports when shots were taken during the course of a game.

So that’s basketball. Some time in the future we might talk about soccer or coffee, but I have some other posts to get to before that. There’s also a homework problem if you continue to scroll down.

* I think I’m suggesting if the only ideas that you have involve deriving a complicated model, the solution might be to refine the question.

** I stated that there’s always a complicated way of looking at a problem. What I mean is that you can take a complex phenomenon at face value and then obviously it will appear complex. Today’s philosophical question is:

Given a complex process does there necessarily exist a simple way of looking at it that will yield productive insight?