There's an interesting contrast between the laws of nature and the laws which constitute our legal system. The laws of nature are compact and precise; written in standard notation without accompanying explanation, the fundamental laws fit on a few pages. The laws of the legal system span thousands of volumes and are frequently ambiguous and ever-changing. On the other hand, we know what the laws of the land actually are. The laws of nature are not completely explored; there's large regions of the parameter space where we just don't know the laws at all. Still, in that sense physicists and lawyers could very very roughly be said to be in the same sort of business.

My sister happens to be in the latter class, finishing up her first semester in law school. When the semester ended she had a math question for me: grades were posted (hers were very good) but class percentile rankings were not. Is it possible to estimate the latter given the former?

Normally the answer is no, since only knowing your grades is not helpful without knowing how other students did on average. But law school - or at least this one - has a presumably anti-grade-inflation policy wherein each class must have its grades assigned in such a way as to result in an overall mean GPA of 2.8 for the students. Therefore your grade can at least establish where you are with respect to that average. I presume it's the mean anyway, I suppose they might require it to be the median.

But unfortunately that's not enough to allow a good estimate of class rank either. We would need to know how the students were distributed about the mean. We don't, so I had to admit it was impossible to say anything meaningful. Even if we assumed the grades were normally distributed, we don't know the standard deviation. But since this is Sunday Function. let's pretend we did know what the grade distribution was. How can we compute percentile rank?

Assume the grades are normally distributed with a mean of 2.8 and a standard deviation of 0.5. This isn't quite possible since grades can't go above 4.0, but all things considered this isn't a terribly implausible way for grades to be distributed. Given these parameters, the function describing this distribution is just the gaussian function scaled and shifted:

Here mu is the mean and sigma is the standard deviation. Plot it:

The area under the curve in a particular range of grades represents the fraction of the students with grades in that range. The complete area under the curve is exactly equal to 1, representing all of the students. The area between 0 (well, technically negative infinity) and 2.8 is 0.5, meaning half of the students are below 2.8 in their grades.

So we need a function to make this procedure systematic. Given a particular GPA, calculate the area under the curve between zero and that given GPA, which we'll call x (again, technically between negative infinity and x but that's pretty irrelevant here since the normal distribution is negligibly small below x = 0).

This is just an integral. Integrate the function f(t) over all t between minus infinity and x, and call that a new function g(x) - our Sunday Function. The jazz with f(t) instead of f(x)is because you're not technically supposed to have the same variable names in the integral and its limits. They're the same function though, we're just giving the argument a different name.

This Sunday Function is so important it has its own name. It's called the error function. Plug in a GPA, and it tells you how much of the area under the normal distribution is lower than that value. Let's plot it, with parameters appropriately adjusted to match our particular grade distribution:

So if you had a 3.5 GPA, you'd plug in and see that the error function yields a value of 0.919, leaving you higher than just about 92% of the class.

You do have to do a little work to compute this, your average calculator doesn't have a button that'll just do this for you. Generally the procedure is to take the value of the GPA (or whatever else you happen to be calculating), changing variables such that it's scaled to the normal distribution of mean 0 and standard deviation 1, and plugging that value into the usual error function that I've linked above. The error function is tabulated in books, and most scientific and graphing calculators can do it from scratch.

It sounds a little involved, but compared to the rules of civil procedure I imagine it's pretty much trivial.

- Log in to post comments

A good analysis, though I think one weakness is that sometimes grade distributions can be bimodal (they certainly were for the class I graded) which would make the analysis harder (and would probably lower your percentile, in this example).

As a bad lower bound you can use Chebyshev's inequality for percentile. There's a (1/1.4)^2 probability of being outside the GPA (2.1, 3.5) (about 50%), so if we assume the distribution is symmetric around the mean then you KNOW you're at least at the 75th percentile, and very probably doing much better.

And immediately after posting and looking through Wikipedia I found there's a one-sided Chebyshev inequality, Cantelli's inequality. So no matter how pathological the distribution is, a 3.5 GPA in this example will at least be at the 1-(1+1.4^2)^-1 = 66th percentile.

In the many things I've graded, I have noticed that distributions tend to be multi-modal, especially when the average is high relative to a 100. In large exams where the average is low (say, the 50s) things tend to look a lot more Gaussian. The last large-class final whose grades I logged was a rather depressing affair, with an average right at 50 if I remember properly. But it was also a gorgeous Gaussian, the bell curve fit was just astonishing. I still have that Excel data (properly FERPA-sanitized) sitting around somewhere; I should do a post on it one day.

Interesting---any reason you can think of that higher average tests have less Gaussian distributions than lower average tests (usually---I recall getting to see the histogram for a test that had an average of 40 something that was pretty bimodal, actually)?

Well, the easier and higher average tests tend to run into the "can't get higher than a 100%" wall, and so the distribution fails to be continuous there. There the studied/slacked distinction tends to be especially dominant since the bright but not stellar students are still nearly maxing out. Once things get hard the "many tiny influences" criteria of the central limit theorem kicks in and things start looking more Gaussian. That's my guess, anyway.

It's also a function of time too; early in the semester the lowest-modal students haven't dropped yet. By finals they're mostly gone.

A typical upper level physics exam will have three or four questions. Often these questions are not of equal difficulty: many professors aim to have one or two questions that everybody should get, one that requires some in-depth understanding of the material, and one that you really have to know your stuff to get. With an exam like this I would not be surprised to see a multi-modal distribution: that is the intent of the exam designer. But if you have an exam with more true/false, multiple choice, or short answer questions, it's more likely that your curve will resemble a Gaussian (but watch for power law tails). For me, even the freshman level exams were of the former kind, but many professors give the latter sort of exam in introductory courses.

You can get pretty much the distribution you want by varying the point values. When I taught an elective class in high school, I didn't want the weak students to fail and didn't want to frustrate the diligent, so I included 10-15 higher level questions with very low point values so that the bright/diligent could get A's that actually reflected their proficiency while not allowing the weaker/lazier students to unfairly get A's or be punished with F's in what should be an enjoyable elective. As high school teachers we were largely hostage to the arbitrary 100 point system. So, we had to be, uh, flexible to maintain the confidence of students, parents and admin in our grading practices.

I assure you, physicists and lawyers are most definitely not in the same sort of business, no matter how roughly put.

More's the pity for us lawyer-types, I'd say.

Honestly, this is nice and all (especially as a spur to brush up on my skills, god I'm getting rusty), but it's still wholly and painfully inappropriate to do this at all. Since the Gaussian admits a negative score, this needs to be done with some other distribution. If only those others ::cough::lognormal::cough:: were so easy to teach and use.