Tuesday, February 24, 2009

APCOM calling

One more acronym to remember! APCOM stands for Application of Computers and Operations Research in the Mineral Industry. It adds up and it stuck. My son knows a lot about EMF and I know a little about metrology in the mineral industry. We added it up and posted an abstract for our APCOM paper. We did so before the first deadline in 2008 came and went. APCOM is still calling for papers. Another deadline passed on February 15, 2009. Where have all those abstracts gone? Popular myth has it that abstracts always swamp program chairs. APCOM is to be held at Vancouver, British Columbia, Canada, from October 6 to 9, 2009. It 's high time to complete the program.

APCOM’s preliminary program brings up The Future of Mining and Technology’s Role. It doesn’t add up to another acronym but does have a futuristic ring to it. The first day is set aside for Resource Identification Estimation and Planning. That's a load of cool stuff. And it’s comes to my home turf in Vancouver, British Columbia. I would like to show resource planners how to estimate metal contents and grades of reserves and resources. It fits APCOM’s preliminary program like a silk glove. We wrote about it in 1991. We want to talk a bit more about it in 2009.

But wait a pixel-picking second! Where would our paper fit in? Geostatistics 1 on the first day? And Geostatistics 2 on the same day? Looks like Geostatistics for the Next Century all over again. That was quite a krige-and-smooth-fest. It took place at Montreal in 1993 with Chairperson Dimitrakopoulos in charge. At that time he was all wrapped up at the cutting edge of conditional simulation with voodoo variances. What I wanted to do was talk about The Properties of Variances. Nowadays, Dimitrakopoulos is Editor-in-Chief, Journal for Mathematical Geosciences. He is still teaching McGill’s students all he knows about conditional simulation. Despite the fact that signs of change to come are posted on IAMG’s website.

In fact, IAMG in 2009 is promoting mathematics, statistics and informatics. I couldn’t have made it up. No longer does it stand for International Association for Mathematical Geology but for International Association for Mathematical Geosciences. Its mission is “…to promote, worldwide, the advancement of mathematics, statistics and informatics in the Geosciences”. How about that? Agterberg, IAMG’s Past President, has yet to explain why his distance-weighted average point grade doesn’t have a variance. He might one day be asked under oath to explain whether or not his distance-weighted average point grade has a variance.

The following formula gives Agterberg's distance-weighted average point grade of his set of five (5) point grades. Agterberg's real problem is that this formula converges on the Central Limit Theorem if all weighting factors converge on 1/n.

It’s a bit of a stretch for Agterberg and his associates “…to promote…” statistics but declare null and void David’s “…famous Central Limit Theorem”. It doesn't take a genius to prove that the next formula converges on David's famous one when each wi converges on 1/n. The question is then what Agterberg was thinking in 1970 and in 1974. Surely, it's about time he rights his wrong!

I’m a stickler for unbiased sampling practices and sound statistics. So, I’m not at all surprised that scores of geoscientists would rather work with real statistics than with surreal geostatistics. After all, statisticians test for spatial dependence, chart sampling variograms, and count degrees of freedom. No hanky panky with the kriging game. No ifs and buts about degrees of freedom! No functions without variances. All scientists and engineers on our little planet should make sound statistics a way of life.

Sunday, February 08, 2009

Junk statistics on Wikipedia

Wiki’s Kriging bugs me just as much today as did textbooks on geostatistics in the 1990s. What bugs me most of all is that Wiki’s keepers of Kriging didn’t give the set of measured values. That brought back many bad memories of Matheron’s magnum opus. His 1954 Formule des Minerais Connexes is peppered with formulas and symbols, gives but few a derived statistics, and no sets of measured values. This paper is posted as Note géostatistique No 1 but is itself marked Note Statistique No 1. Somebody played a silly game in predating the birthday of Matheron’s new science of geostatistics. Wiki's link to Matheron's seminal work went dead on December 12, 2008. Click Centre de Geosciences, and go to Ressources Documentaires & Logiciels. Next, click Bibliothèque Géostat (en ligne) and take a long look at Matheron’s past. This link is still hot, and I’m tickled pink. I’ve got to get this new link on some of my old blogs.

Wiki's Kriging keepers made me think of Matheron's statistically challenged disciples, and of all their tangled thoughts. Why do formulas and symbols run rampant where sets of measured values are as scarce as hen’s teeth? Little odds and ends of geostat speak such as “…a system of linear equations which is obtained by assuming that ƒ is a sample path of a random process F(x)…” make me cringe. Sounds a bit like Matheron’s take on Brownian motion. Why did degrees of freedom fail to inspire the wardens of Wiki’s Kriging? Of course, it would explain why they just keep on kriging for life!

I had asked for but never got the set of measured values that underpins Figure 1. So, I derived the same measured values in scale units. Fisher’s F-test proved that the ordered set of scale units does not display a significant degree of spatial dependence. The guardians of Wiki’s Kriging pointed out, “From the geological point of view, the practice of kriging is based on assuming continued mineralization between measured values”. What a way to practice kriging! Stanford’s Journel espoused the same sort of assumed nonsense in 1992. I never took him serious but he may well have thought he was. I did what Journel didn't do in 1978. Several years before the Bre-X fraud I derived variances of density- and length-weighted average lead and silver grades of core samples. I worked with weighting factors since the set of measured values in Figure 1 is unevenly spaced.

I’m caught between real statistics and hardcore kriging. I interpolated by kriging between each pair of measured Y- and X-values. The spreadsheet template shows that the first pair gives a Y-value of 103.0 scale units and an X-value of 25.8 scale units, the second pair gives a Y-vale of 96.0 scale units and an X-value of 45.5 scale units, and so on for a set of seventeen (17) pairs. The following chart shows why interpolation by kriging does so much more with less. All it takes is to rig the rules of statistics.

False 95% confidence intervals

Now here’s the clincher. Fisher’s F-test cannot be applied to an ordered set of seventeen (17) values, each of which is either measured or kriged. The problem is sets of measured values do give degrees of freedom whereas kriged values give none. A simple rule of thumb is that measured values do give degrees of freedom whereas kriged values give nothing but headaches. Unless, of course, one grasps the irrefutable fact that each kriged estimate does have its own variance just as much as do central values such as arithmetic means and all sorts of weighted averages.


Listed above are 95% confidence limits for central values of nine (9) measured values only, and of seventeen (17) measured and kriged values. Interpolation by kriging between measured values seem to give a higher degree of precision do than measured values alone. It’s not so much that Krige knew how to work miracles with a few measured values but that Matheron’s disciples have rigged the rules of real statistics. To put it simply, a kriged estimate has its own variance since all functionally dependent values do. A reliable rule of thumb is kriged estimates give big problems whereas measured values give degrees of freedom.

One would expect Wiki’s Kriging squad to show how “95% confidence intervals” in Figure 1 were derived. Surely, the squad was joking when it put forward, “Assuming prior knowledge encapsulates how minerals co-occur as a function of space. Then, given an ordered set of measured grades, interpolation by kriging predicts mineral concentrations at unobserved points”. How about that? Sounds like Wikipedians live in Wonderland. Krige himself couldn’t have cooked up such drivel.