Friday, April 01, 2011

ISO to tackle trueness

ISO has set the stage to tackle trueness. It did so by issuing an NWIP ( New Work In Progress). I am pleased that ISO/TC69/SC6/WG1 has been entrusted with the task. Accuracy (trueness and precision) is way wide of the mark. Trueness (accuracy and precision) reads but a bit better. I would rather work with either Precision and accuracy or Precision and bias. Trueness makes more sense in a court of law. I have worked on a number of ISO Committees since 1974. I would find it a senseless task to keep track of trueness. Testing for bias always makes sense. Student’s t-test not only shows what’s biased and what’s not but also gives intuitive measures for statistical risks. I am pleased that ISO never took to kriging and smoothing. Unlike ASTM which went along with geostatistical thinking in 1994. ASTM was set up in 1898 by chemists and engineers. A precursor of CIMMP was set up in 1898. CIMMP's geostatistical peer review has been biased since the 1990s.

ISO/TC183 has set up a standard method to derive metal contents and grades of mineral concentrates and ores. It would have been just as simple to set up an ISO Technical Committee to derive metal contents and grades of mineral reserves and resources. But the COSMO brains behind the mining industry were too keen to assume, krige, smooth, and rig the rules of applied statistics. That’s why Bre-X’s bogus grades and Busang’s barren rock morphed so smoothly into a massive phantom gold resource. And that’s when Barrick Gold asked me to assist in testing for bias between paired gold assays determined by cyanide leaching and by fire assaying.

Student's t-test for bias

The observed t-value of 11.258 is statistically significant at 99.9% probability. It’s impossible to salt drill core with placer gold. In contrast, salting crushed core with placer gold is a cinch. Early in 1997 I derived confidence limits for the mass of gold for one of Barrick’s many deposits. My report was deemed worth its weight in gold. Here’s what Richard Rohmer wrote in Golden Phoenix, “For Munk and others affected by Bre-X and Busang, the strange news about de Guzman was the first hint that the find might well not be what Walsh and Felderhof were claiming”. He also pointed out, “Peter Munk was appalled when he read the Strathcona report. The damage inflicted on Canada’s national and international mining industry was beyond what anyone would have thought possible a few short weeks before”. Now that’s what crying out late is all about! I took a look in my crystal ball. It seems as if a Munk Debates for and against geostatistics has been put on ice!

Here's what Peter Munk himself pontificated on the jacket of his Golden Phoenix:
"You have to be courageous; you have to learn to take advantages of change. Be non-conventional; don't fritter your energies - be focused; remember to share. Most important, use the biggest weapon of all weapons, the least appreciated yet the most important tool for success, and this is moral integrity; and don't be afraid to dream and don't be afraid to dream big."


Never mind Munk's moral integrity! I dream about scientific integrity in mineral exploration and mining. It’s a piece of cake to derive borehole statistics with spreadsheet software. I like to call such stats the fingerprint of a borehole. Those who have tried like it a lot! Following is a synopsis derived of the fingerprint of one mind-boggling borehole.

Fingerprint of monster borehole

It shows some 5.5 million ounces of gold in 78 million metric tons of ore. SME published in Transactions 2000 a paper on Borehole Statistics with Spreadsheet Software. One reviewer called it “an excellent paper” and the other thought “it would stir up a hornet’s nest”. But the hornets stayed put!

It makes no sense to assume spatial dependence in sample spaces. Fisher’s F-test ought to be applied to the variance of the set of measured values and the first variance term of the ordered set. Interpolation between measured values gives a false positive for spatial dependence. And that’s where geostatistocrats are taking not only mineral reserve and resource estimation but also the study of climate change. Or climate dynamics as I like to call it! I found Loehle’s 2000-year global temperature reconstruction by accident. Loehle and McCulloch pointed out in 2008 that data sets were smoothed with a 30-year running mean. It does create a pretty smooth picture. But here’s what happens with variances when measured values are enriched with variance-deprived kriged estimates.

Test for homogeneity of variances

I partitioned Loehle’s 2000-year set into twenty (20) sets of no more than 100-year each, derived the first variance term of each set, and applied the chi-square test to this set of twenty (20) variances. These variances constitute too homogeneous a set. The probability that this inference is true exceeds 99.9%. The probability that it is false is less than 0.1%.

I trust that ISO/TC69/SC6/WG1 will never assume spatial dependence between measured values in ordered sets. ISO/TC17 on Steel and ISO/TC34 on Food Products are interested in this New Work in Progress. I shall continue to write whatever needs to be put in writing to ensure that applied statistics prevails not only in ISO standards but also in mineral reserve and resource estimation.

No comments: