Leave a Reply

4 Comment threads
13 Thread replies
Most reacted comment
Hottest comment thread
5 Comment authors
Notify of

Where Anthony Watts gets it wrong yet again and pastes Australian data into the post.
And you get a hearty endorsment from Ken Ring.

Still, that’s no fault of the paper.


Actually, I quite appreciated Ken’s post, did you read it, Simon? About alleged changes in pressure levels around NZ due to climate change.

Richard C (NZ)

I read it too. Can’t say I know much about Ken Ring, controversial I gather, but his post on historical NZ pressure was interesting.

Richard C (NZ)

Steve McIntyre commented at BH. I don’t think he’d read the paper and he’s asking for code which RS93 never supplied, neither SI. Why then de Freitas et al for which there is not necessarily code as RS93? The method is there to be implemented in whatever way, not necessarily code. I note Spencer and Braswell (2014) used both Fortran and Excel to get the same results.

He’s also asking for data on an ftp server. Given it’s readily available in CliFlo for which I gave him a link, I would have thought that’s unnecessary. Whether he goes from “semi-interested” and “insufficiently interested to try” to interested enough to replicate remains to be seen. I hope i haven’t antagonized him but what he’s asking is superfluous to my mind, although granted he operates in a different and loftier sphere than I do. Perhaps I should stop helping in that respect.


Richard C (NZ)

I’ve thought about coding RS93 and PM-95 (ACORN) but I’m not skilled in any modern scripting and probably too late now to pick up from donkey’s years ago. Don’t think I have the stats nouse anyway which is why I backed out early on when it became obvious that the rigour of RS93 was required. It would be very useful though to have these breakpoint analyses techniques coded in software to run against other series e.g. PM-95 against 7SS and RS93 against ACORN. BEST’s Scalpel is another one. Few of us are developers or fluent with stats but most can be users when an easy tool is at hand for rudimentary demonstration. RS93 looks a bit daunting at first but I’m wondering how the process can be eased by modern IT techniques rather than low level coding e.g diagrams that automatically generate code (think Bizagi for business) or modular pseudo code, pseudo English, or something. Not sure of the appropriate jargon nowadays, let alone higher level tools. I’m thinking in terms of making the stats simpler so that more people understand the breakpoint analysis process. Huub Bakker swears by VB6 but I’m wondering… Read more »


Steve McIntyre, Steven Mosher and all those guys are all “R” programmers, (which I don’t know) so it would actually be really helpful if someone got some R code and data out in the public domain

Richard C (NZ)

Yes, I see this in RHtestsV4 users manual “This simple manual is to provide a quick reference to the usage of the functions included in the RHtestsV4 package (and also to the usage of the equivalent FORTRAN functions, which are available by sending a request in English to Xiaolan.Wang@ec.gc.ca). Users are assumed to have the general knowledge of R (how to start and end an R session and how to call an R function).” Plenty of tutorials e.g Getting Started with the R Data Analysis Package http://heather.cs.ucdavis.edu/~matloff/r.html But I’ve yet to discover BOM’s code i.e. whether it’s R or something else. And if BOM has a bespoke package then is it really necessary to know R anyway except how to start it? The key is the stats though. I’ve never coded anything like that. My experience, years ago now, was information crunching in business situations – RDBMS in COBOL with embedded SQL at most academically, then just ACCESS RDBMS commercially. I did do a lot of systems analysis diagramming manually because of my drafting background before the automated tools came in which I’ve used a bit since but not commercially. I think modular… Read more »


You know Huub Bakker? I worked for him for a while, a good man. Well, two levels under him, to be fair.
I think it may be a different Huub, actually.

Richard C (NZ)

>”You know Huub Bakker?”

I don’t know him but he used to frequent this blog. Commented a bit recently at JoNova in the solar model threads. That was where I saw he uses VB6 for teaching:

Dr Huub Bakker – Senior Lecturer – Massey University

Same Huub?


No, different, sorry.

Richard C (NZ)

Hey Andy, Huub Bakker’s wind turbine/human effects work:


Dickinson, P. (2010). The sounds from wind turbines. Theory, practice, assumptions, and reality.. In B. Rapley, & H. Bakker (Eds.) Sound, Noise, Flicker and the Human Perception of Wind Farm Activity. New Zealand: Atkinson and Rapley Consulting
[Chapter]Edited by: Bakker, H.

Bakker, HH., & Rapley, BI. (2010). Sound Characteristics of Multiple Wind Turbines. In BI. Rapley, & HHC. Bakker (Eds.) Sound, Noise, Flicker and the Human Perception of Wind Farm Activity. : Atkinson and Rapley Consulting
[Chapter]Authored by: Bakker, H.Edited by: Bakker, H.
Read Abstract: abstract icon

Rapley, BI., & Bakker, HH. (Eds.) (2010). Sound, Noise, Flicker and the Human Perception of Wind Farm Activity. Palmerston North: Atkinson & Rapley Consulting Ltd (Palmerston North, New Zealand) in association with Noise Measurement Services Pty Ltd (NMS) (Brisbane, Australia)
[Edited Book]Authored by: Bakker, H.Edited by: Bakker, H.
Read Abstract: abstract icon

Bakker, HH., Bennett, D., Rapley, BI., & Thorne, R.Seismic Effects on Residents from Wind Turbines. In BI. Rapley, & HHC. Bakker (Eds.) Sound, Noise, Flicker and the Human Perception of Wind Farm Activity.
[Chapter]Authored by: Bakker, H.Edited by: Bakker, H.
Read Abstract: abstract icon


Richard C (NZ)

‘Evaluating the impact of wind turbine noise on health related quality of life’ Daniel Shepherd, David McBride1, David Welch2, Kim N. Dirks2, Erin M. Hill (2011) Department of Psychology, School of Public Health, Auckland University of Technology, Auckland, 1Department of Preventive and Social Medicine, University of Otago, Dunedin, 2School of Population Health, The University of Auckland, Auckland, New Zealand Abstract We report a cross-sectional study comparing the health-related quality of life (HRQOL) of individuals residing in the proximity of a wind farm to those residing in a demographically matched area sufficiently displaced from wind turbines. The study employed a nonequivalent comparison group posttest-only design. Self-administered questionnaires, which included the brief version of the World Health Organization quality of life scale, were delivered to residents in two adjacent areas in semirural New Zealand. Participants were also asked to identify annoying noises, indicate their degree of noise sensitivity, and rate amenity. Statistically significant differences were noted in some HRQOL domain scores, with residents living within 2 km of a turbine installation reporting lower overall quality of life, physical quality of life, and environmental quality of life. Those exposed to turbine noise also reported significantly lower… Read more »

Richard C (NZ)

While I’m at it, I’ve laid out the ACORN-SAT methods at JN here (#38.2):


Following on in the next comment (#36.2.1) is a percentile-matching (PM-95 for breakpoint analysis) definition and tutorial. We are dealing with “parameter estimates”, with no guarantee that an estimate will be unique.

Also discussed here:

Richard C (NZ)

I can’t see much to be gained by applying the same breakpoint analysis technique to a series by multiple parties for replication. It’s been done by de Freitas et al and someone can replicate it if they want to verify. OK it needs to be done but what does that really prove?.

It’s the different methods that would really be revealing on the same series – RS93 vs PM-95 vs Scalpel for example. NIWA vs NZCSC vs BOM vs BEST. We’ve got NIWA vs NZCSC now and that’s a real eyeopener but I think there’s many more surprises in store.

Then there’s GISS. What’s their method?

Richard C (NZ)


Apart from breakpoint techniques, the respective homogenization methods are vastly different too. But there are common breakpoints. Then there’s non-site change break adjustments by the others extra to those in NIWA/NZCSC – messy.

In Anthony’s post he featured Rutherglen Research (for some reason) but BEST has a much shorter series for Rutherglen Research. Their earlier Rutherglen data is not homogenized to Rutherglen Research, they go to Rutherglen Post Office.

BEST Mean http://berkeleyearth.lbl.gov/stations/151882
BOM Mincomment image?w=720

Earliest Observations: 1965 vs 1912ish

BEST makes an adjustment at the 1980 Record Gap but as far as I know BOM doesn’t. In the immediate vicinity of Rutherglen, BOM has a continuous series, BEST has a gap 1921 – 1965.


Earliest Observation: Jan 1903
Most Recent Observation: Nov 1921

The NZ 7SS is clean by comparison to Rutherglen, BOM vs BEST.

Richard C (NZ)

BEST runs Albert Park and Auckland Aero in parallel:

AUCKLAND, ALBERT PARK http://berkeleyearth.lbl.gov/stations/157062
AUCKLAND AERO AWS http://berkeleyearth.lbl.gov/stations/157061

NIWA/NZCSC drops Albert Park when Aero starts, homogenizing Albert Park to Aero.

BEST Albert Park is above NZCSC Auckland but below NIWA Auckland. If BEST were to adjust for UHI/sheltering, their Albert Park series would corroborate de Freitas et al (2014)/’Statistical Audit’ Auckland but nowhere near NIWA Auckland.

Post Navigation