## Astrostats

A few weeks ago I posted an item on the theme of how gambling games were good for the development of probability theory. That piece contained a mention of one astronomer (Christiaan Huygens), but I wanted to take the story on a little bit to make the historical connection between astronomy and statistics more explicit.

Once the basics of mathematical probability had been worked out, it became possible to think about applying probabilistic notions to problems in natural philosophy. Not surprisingly, many of these problems were of astronomical origin but, on the way, the astronomers that tackled them also derived some of the basic concepts of statistical theory and practice. Statistics wasn’t just something that astronomers took off the shelf and used; they made fundamental contributions to the development of the subject itself.

The modern subject we now know as physics really began in the 16^{th} and 17^{th} century, although at that time it was usually called Natural Philosophy. The greatest early work in theoretical physics was undoubtedly Newton’s great *Principia*, published in 1687, which presented his idea of universal gravitation which, together with his famous three laws of motion, enabled him to account for the orbits of the planets around the Sun. But majestic though Newton’s achievements undoubtedly were, I think it is fair to say that the originator of modern physics was Galileo Galilei.

Galileo wasn’t as much of a mathematical genius as Newton, but he was highly imaginative, versatile and (very much unlike Newton) had an outgoing personality. He was also an able musician, fine artist and talented writer: in other words a true Renaissance man. His fame as a scientist largely depends on discoveries he made with the telescope. In particular, in 1610 he observed the four largest satellites of Jupiter, the phases of Venus and sunspots. He immediately leapt to the conclusion that not everything in the sky could be orbiting the Earth and openly promoted the Copernican view that the Sun was at the centre of the solar system with the planets orbiting around it. The Catholic Church was resistant to these ideas. He was hauled up in front of the Inquisition and placed under house arrest. He died in the year Newton was born (1642).

These aspects of Galileo’s life are probably familiar to most readers, but hidden away among scientific manuscripts and notebooks is an important first step towards a systematic method of statistical data analysis. Galileo performed numerous experiments, though he certainly carry out the one with which he is most commonly credited. He did establish that the speed at which bodies fall is independent of their weight, not by dropping things off the leaning tower of Pisa but by rolling balls down inclined slopes. In the course of his numerous forays into experimental physics Galileo realised that however careful he was taking measurements, the simplicity of the equipment available to him left him with quite large uncertainties in some of the results. He was able to estimate the accuracy of his measurements using repeated trials and sometimes ended up with a situation in which some measurements had larger estimated errors than others. This is a common occurrence in many kinds of experiment to this day.

Very often the problem we have in front of us is to measure two variables in an experiment, say X and Y. It doesn’t really matter what these two things are, except that X is assumed to be something one can control or measure easily and Y is whatever it is the experiment is supposed to yield information about. In order to establish whether there is a relationship between X and Y one can imagine a series of experiments where X is systematically varied and the resulting Y measured. The pairs of (X,Y) values can then be plotted on a graph like the example shown in the Figure.

** **

** **

In this example on it certainly looks like there is a straight line linking Y and X, but with small deviations above and below the line caused by the errors in measurement of Y. This. You could quite easily take a ruler and draw a line of “best fit” by eye through these measurements. I spent many a tedious afternoon in the physics labs doing this sort of thing when I was at school. Ideally, though, what one wants is some procedure for fitting a mathematical function to a set of data automatically, without requiring any subjective intervention or artistic skill. Galileo found a way to do this. Imagine you have a set of pairs of measurements (x_{i},y_{i}) to which you would like to fit a straight line of the form *y=mx+c*. One way to do it is to find the line that minimizes some measure of the spread of the measured values around the theoretical line. The way Galileo did this was to work out the sum of the differences between the measured y_{i} and the predicted values *mx+c* at the measured values x=x_{i}. He used the absolute difference *|y _{i}-(mx_{i}+c)|* so that the resulting optimal line would, roughly speaking, have as many of the measured points above it as below it. This general idea is now part of the standard practice of data analysis, and as far as I am aware, Galileo was the first scientist to grapple with the problem of dealing properly with experimental error.

The method used by Galileo was not quite the best way to crack the puzzle, but he had it almost right. It was again an astronomer who provided the missing piece and gave us essentially the same method used by statisticians (and astronomy) today.

Karl Friedrich Gauss was undoubtedly one of the greatest mathematicians of all time, so it might be objected that he wasn’t really an astronomer. Nevertheless he was director of the Observatory at Göttingen for most of his working life and was a keen observer and experimentalist. In 1809, he developed Galileo’s ideas into the method of least-squares, which is still used today for curve fitting.

This approach involves basically the same procedure but involves minimizing the sum of [y_{i}-(mx_{i}+c)]^{2} rather than |y_{i}-(mx_{i}+c)|. This leads to a much more elegant mathematical treatment of the resulting deviations – the “residuals”. Gauss also did fundamental work on the mathematical theory of errors in general. The normal distribution is often called the Gaussian curve in his honour.

After Galileo, the development of statistics as a means of data analysis in natural philosophy was dominated by astronomers. I can’t possibly go systematically through all the significant contributors, but I think it is worth devoting a paragraph or two to a few famous names.

I’ve already mentioned Jakob Bernoulli, whose famous book on probability was probably written during the 1690s. But Jakob was just one member of an extraordinary Swiss family that produced at least 11 important figures in the history of mathematics. Among them was Daniel Bernoulli who was born in 1700. Along with the other members of his famous family, he had interests that ranged from astronomy to zoology. He is perhaps most famous for his work on fluid flows which forms the basis of much of modern hydrodynamics, especially Bernouilli’s principle, which accounts for changes in pressure as a gas or liquid flows along a pipe of varying width.

But the elder Jakob’s work on gambling clearly also had some effect on Daniel, as in 1735 the younger Bernoulli published an exceptionally clever study involving the application of probability theory to astronomy. It had been known for centuries that the orbits of the planets are confined to the same part in the sky as seen from Earth, a narrow band called the Zodiac. This is because the Earth and the planets orbit in approximately the same plane around the Sun. The Sun’s path in the sky as the Earth revolves also follows the Zodiac. We now know that the flattened shape of the Solar System holds clues to the processes by which it formed from a rotating cloud of cosmic debris that formed a disk from which the planets eventually condensed, but this idea was not well established in the time of Daniel Bernouilli. He set himself the challenge of figuring out what the chance was that the planets were orbiting in the same plane simply by chance, rather than because some physical processes confined them to the plane of a protoplanetary disk. His conclusion? The odds against the inclinations of the planetary orbits being aligned by chance were, well, astronomical.

The next “famous” figure I want to mention is not at all as famous as he should be. John Michell was a Cambridge graduate in divinity who became a village rector near Leeds. His most important idea was the suggestion he made in 1783 that sufficiently massive stars could generate such a strong gravitational pull that light would be unable to escape from them. These objects are now known as black holes (although the name was coined much later by John Archibald Wheeler). In the context of this story, however, he deserves recognition for his use of a statistical argument that the number of close pairs of stars seen in the sky could not arise by chance. He argued that they had to be physically associated, not fortuitous alignments. Michell is therefore credited with the discovery of double stars (or binaries), although compelling observational confirmation had to wait until William Herschel’s work of 1803.

It is impossible to overestimate the importance of the role played by Pierre Simon, Marquis de Laplace in the development of statistical theory. His book *A Philosophical Essay on Probabilities*, which began as an introduction to a much longer and more mathematical work, is probably the first time that a complete framework for the calculation and interpretation of probabilities ever appeared in print. First published in 1814, it is astonishingly modern in outlook.

Laplace began his scientific career as an assistant to Antoine Laurent Lavoiser, one of the founding fathers of chemistry. Laplace’s most important work was in astronomy, specifically in celestial mechanics, which involves explaining the motions of the heavenly bodies using the mathematical theory of dynamics. In 1796 he proposed the theory that the planets were formed from a rotating disk of gas and dust, which is in accord with the earlier assertion by Daniel Bernouilli that the planetary orbits could not be randomly oriented. In 1776 Laplace had also figured out a way of determining the average inclination of the planetary orbits.

A clutch of astronomers, including Laplace, also played important roles in the establishment of the Gaussian or normal distribution. I have also mentioned Gauss’s own part in this story, but other famous astronomers played their part. The importance of the Gaussian distribution owes a great deal to a mathematical property called the Central Limit Theorem: the distribution of the sum of a large number of independent variables tends to have the Gaussian form. Laplace in 1810 proved a special case of this theorem, and Gauss himself also discussed it at length.

A general proof of the Central Limit Theorem was finally furnished in 1838 by another astronomer, Friedrich Wilhelm Bessel– best known to physicists for the functions named after him – who in the same year was also the first man to measure a star’s distance using the method of parallax. Finally, the name “normal” distribution was coined in 1850 by another astronomer, John Herschel, son of William Herschel.

I hope this gets the message across that the histories of statistics and astronomy are very much linked. Aspiring young astronomers are often dismayed when they enter research by the fact that they need to do a lot of statistical things. I’ve often complained that physics and astronomy education at universities usually includes almost nothing about statistics, because that is the one thing you can guarantee to use as a researcher in practically any branch of the subject.

Over the years, statistics has become regarded as slightly disreputable by many physicists, perhaps echoing Rutherford’s comment along the lines of “If your experiment needs statistics, you ought to have done a better experiment”. That’s a silly statement anyway because all experiments have some form of error that must be treated statistically, but it is particularly inapplicable to astronomy which is not experimental but observational. Astronomers need to do statistics, and we owe it to the memory of all the great scientists I mentioned above to do our statistics properly.

September 20, 2009 at 10:41 pm

I didn’t realise that Galileo did straight-line fitting, and I’ll bet that minimum sum of absolute deviations is not much different in practice from least squares. Impressive.

I’d be reluctant to call Galileo the “originator” of modern physics, though. Scholars have increasingly found that the so-called Renaissance was little other than a flowering of the visual arts and an increased intensity of looking at classical Greece for inspiration; and come to understand that thought had its real renaissance in Western Europe in the hundred years from approximately the Norman conquest, leading to the so-called “high Middle Ages” that took a body blow with the terrible plague of the mid-14th century. Physics fits this pattern – Galileo got the ball over the line alright, but it had been passed to him through a long and fascinating development running back to the 11th century. An outstanding book on this for a scientist without extensive background in mediaeval history is “God’s Philosophers” by James Hannam. It explains the place of prominent mediaeval thinkers in the parts of philosophy and theology relevant to the development of science, and does so sympathetically yet without presuming the reader necessarily shares the theology. Galileo and Newton represented *both* a culmination and a start of something new. Hannam also points out that the “revolution in science” was actually a revolution in physics, for biology, or at least medicine, was not revolutionised until several centuries later, in the lifetime of the grandfathers of people alive today.

Anton

September 20, 2009 at 10:48 pm

PS There are at least two English translations of Laplace’s Philosophical Essay on Probabilities, but when O when will somebody translate his great Theorie Analytique des Probabilites, to which the Esssay was merely a preamble? One part of it is merely on generating functions but the other is the greatest work on probability up to Laplace’s time and well beyond. I have been told that ET Jaynes used to translate it word by word with the help of a French dictionary in order to defend Laplace from the attacks of “frequentist” statisticians. On which subject, why do we have statisticians at all? You don’t see people calling themselves “Fourier analysts” today. Probability theory was invented by working scientists…

Anton

September 20, 2009 at 10:56 pm

Anton,

I think least-squares and least-absolute-deviations are indeed very similar unless you have a very strange distribution of errors. Since Galileo didn’t have a calculator his approach was very sensible I think.

Peter

September 20, 2009 at 11:26 pm

Yes, Galileo’s way is what you’d do using trial and error and a ruler.

September 21, 2009 at 11:54 am

I wish I had a hat like the one Gauss is wearing.

September 21, 2009 at 2:11 pm

“A clutch of astronomers”

One of the interesting things about English are words for collections: a gaggle of geese, a pride of lions, a surfeit of lampreys—whatever. What would be a good collective term for astronomers? A galaxy of astronomers is probably too obvious.

Need a hat? There are several vendors of mediaeval and renaissance clothing, bought primarily by those who dance the dances of those times or those into early music etc (or even role-playing types). I’m sure you could buy just the hat, though some tights might be a nice addition.

Why don’t more people spend their normal life in renaissance garb? Strangely, no-one would as much as bat an eye today if a punk in full kit rides on the tube, but dress as Henry VIII or Galileo and you’ll be laughed at. The world isn’t fair.

While Galileo might be the first modern physicist, I think a case could be made that Leondardo was the first modern scientist. (Of course, like Galileo he tapped into an older tradition of the likes of Roger Bacon.)

September 21, 2009 at 2:30 pm

I’m not sure about astronomers, but I once heard at a staff meeting a suggestion that the correct collective noun for professors was an

absence…September 21, 2009 at 3:00 pm

Philip is right Peter, as the scion of a hatting family I can assure you that specialist kit *is* available. (I am from Stockport which used to be a hatting town – Stockport County are known as the Hatters; you don’t have to be mad to support them but it helps.) My friend who ran the Sherlock Holmes Society of London – and also does courtly dancing – would know where such stuff can be had in the south. I think it’s more fun to wear 18th century kit than Renaissance garb. As to whether (Philip) you’d be laughed at on the tube for wearing it: are you sure? Have you tried it?

The physicist CA Truesdell used to give a talk about Euler in which he was dressed for the part, incidentally.

Hannam’s book traces the line from Roger Bacon to Leonardo and Galileo very nicely.

Anton

PS RIP Bobby Robson.

September 21, 2009 at 3:11 pm

I used to have a small black hat similar to a Gaussian pill-box shape which I wore when working at home. It began life as a chef’s hat, I think. I found I could concentrate better with it on, although I don’t know why or what happened to it. I never wore it outside, though, as it probably looked a bit daft.

September 21, 2009 at 3:17 pm

ps. Since, according to Sainsbury’s, Christmas is imminent, I should also mention that I would like a Victorian-style smoking jacket…

September 21, 2009 at 5:34 pm

[…] See the rest here: Astrostats « In the Dark […]

October 11, 2009 at 5:45 pm

[…] the Normal (Gaussian) distribution, with its characteristic Bell-shaped curve. I’ve already mentioned the role that various astronomers played in the development of this bit of mathematics, so I […]

March 19, 2010 at 5:51 pm

[…] bad start. Science didn’t get “seduced” by statistics. As I’ve already blogged about, scientists of the calibre of Gauss and Laplace – and even Galileo – were instrumental […]