Thursday, December 01, 2011

NRC shelled out real doug for bogus statistics

The National Research Council did so in the 1970s. Professor Dr Michel David was awarded Grant NRC7035 to advance geostatistics. So he plodded away and got his work printed in 1977. But he came up with a peculiar caution. He pointed out that “professional statisticians would find unqualified statements”. The author mentioned it on page VII of what he had come to call Geostatistical Ore Reserve Estimation. But why didn’t he ask a genuine statistician to read his draft? I found out that the author was right after having read my own copy of his book.

David was as smitten with geostatistics in the 1970s as young Matheron was with applied statistics in the 1950s. Alas, Matheron’s pursuit of applied statistics was not to last. On the contrary, Professor Dr Georges Matheron in the 1970s praised geostatistics as a new science. He did so because he had failed to grasp that functions do have variances. Thus, each and every distance-weighted average has its own variance in applied statistics. Just the same, Matheron got hooked on his variance-deprived distance-weighted average. So much so that he got into calling it a kriged estimate. He did it to honor D Krige who worked with distance-weighted averages at Witwatersrand gold mines in South Africa. What went wrong with Matheron’s new science is that the variance of the kriged estimate has gone missing!

A strong case can be made that eulogies be written long before one’s time on this planet comes to an end. One might ask a lawyer to assist in assuring the veracity of one’s credentials. A case in point is Professor Dr George Matheron’s 2000 eulogy. Dr F P Agterberg was his eulogist. He remembered him as the founder of spatial statistics. As luck would have it, Matheron never tested for spatial dependence by applying Fisher’s F-test to the variance of a set of measured values and the first variance term of the ordered set. Given that Matheron’s magnum opus is posted on the web, one might spend a life time to study his work.

Professor Dr Michel David passed away on May 10, 2000. His obituary was put together by Professor Dr Roussos Dimitrakopoulos and Michel Dagbert. The International Association for Mathematical Geology awarded David the W C Krumbein medal in 1988. David became a Fellow of the Royal Society of Canada in the same year. Moreover, CIMMP recognized David’s worldwide achievements in 1989 with the award of the Selwyn G Blaylock medal.

My son and I knew precisely what had gone wrong with geostatistics when I was face-to-face with David on Saturday, March 23, 1991. It was at a seminar on Sampling and Ore Reserves at the Royal York Hotel, Toronto, Ontario. CIM Bulletin in 1990 rejected Precision Estimates for Ore Reserves. We had shown how to test for spatial dependence and how to derive unbiased confidence limits for gold content and grade. David saw fit to nitpick that twenty years of geostatistical literature went missing. He did not ask me a single question. Take a look at what is wrong in Matheron’s new science of geostatistics.

Marechal & Serra, Random Kriging, 1970, Figure 10
David, Geostatistical Ore Reserve Estimation, 1977, Figure 203

David wrote: Writing all the necessary covariances for that system of equations might be a good test to find out whether one really understands geostatistics! Merks and Merks claim: Deriving the variance of each distance-weighted average AKA kriged estimate, and counting degrees of freedom for that his system of equations might be a good test to find out whether one really grasps applied statistics.

The National Research Council is the Government of Canada’s institute for research and development. As such it has been active since 1916. NRC’s task is to stand on guard for ethics and integrity. It did not know in the 1970s that Matheron’s new science of geostatistics is an invalid variant of applied statistics. Neither did Dr Roger A Blais, a Professor in Economic Geology and a Fellow of the Royal Society of Canada. He made it possible for David to write so much about so little. David also confessed to be indebted to Matheron. David’s writing added up to a batch of bogus statistics. NRC7035 was in place not only for his 1977 Geostatistical Ore Reserve Estimation but also for his 1988 Handbook of Applied Advanced Geostatistical Ore Reserve Estimation.

Dr Isobel Clark in her 1979 Practical Geostatistics derived the variance of the distance-weighted average AKA kriged estimate. Professor Dr Michel David wrote about the famous Central Limit Theorem but did not apply it. Professor Dr Roussos Dimitrakopoulos wrote a touching farewell to Michel David (1945-2000). David is no longer listed under Obituaries of Deceased Fellows. Sic Transit Gloria Mundi!

No comments: