Category Archives: Scripts

Replicating McShane and Wyner

R coder mind of a Markov chain has replicated portions of the M&W work. They write: There are a bunch of “hockey sticks” that calculate past global temps. through the use of proxies when instrumental data is absent. There is a new one out there by McShane and Wyner (2010) that’s creating quite a stir […]

More on Retrieving KNMI Data

I’ve done a considerable upgrade to my function for retrieving model data from KNMI within R. This builds on the KNMI webpage but IMO is a considerable enhancement of it. I’ve made the script available here . The function read.knmi.models is built as an emulator of the radio buttons. Geert’s radio buttons (if I’m understanding […]

Downloading UWisc Data

The Jeffs have been exploring the UWisc AVHRR data which is stored as a lot of gzipped ncdf files. I couldn’t figure out how to download and open this data into R. So, as I usually do in these cases, I asked CA reader Nicholas who, as usual, has a solution. I edited this slightly […]

R – the choice for serious analysis

While Steve is a little “under the weather” (it must be all the snow that Al Gore sent him), I thought I’d mention an interesting article in the New York Times which sings the praises of the programming language R. R is similar to other programming languages, like C, Java and Perl, in that it […]

Emulating CRUTem Graphics

Because of the puzzles in trying to replicate the NOAA graphic from archived data, I’ve also tried to replicate the HadCRU plots from archived data. Hu McCulloch has criticized the highly inappropriate use of Mercator projections by climate scientists and, as an exercise, I figured out a fairly simple method of plotting these results in […]

A Collation Utility for GISS dset1 and dset2

GISS has been providing a considerable amount of intermediate information on their results. Unfortunately, it’s been provided in binary format that is presumably suited for people who speak Fortran with a Unix accent. I presume that such people converse with one another in medieval Latin. It’s not very handy for people who use modern languages. […]

An R Package by a CA Reader Solves the Z Problem

CA reader, Nicholas, an extremely able computer analyst, has helped me with a number of problems with downloading data in compressed formats into R. One of the most annoying and heretofore unsolved problems was how to get Z files into R without having to handle them manually – a problem that I revisited recently when […]

Gridding from CRN1-2

In order to assess the 2006 anomaly in John V’s 2006 calculation reported previously, I’ve calculated a gridded anomaly for the U.S. using standard spatial smoothing software in R. This enables the calculation to be done in a few lines as explained below. I’ll also show results which reconcile most of the results to date. […]

Peterson (2003)

Peterson (2003) online here is an influential study cited by IPCC AR4 purporting to show that the urbanization effect is negligible. It concluded: Using satellite night-lights—derived urban/rural metadata, urban and rural temperatures from 289 stations in 40 clusters were compared using data from 1989 to 1991. Contrary to generally accepted wisdom, no statistically significant impact […]

USHCN Trends: Red States and Blue States

I’ve calculated the trends for all 1221 USHCN stations after 1900 and plotted up the contour maps for the raw, time-of-observation adjusted and fully adjusted results. Pretty interesting results. Here’s the results for the “raw” data, the adjusted versions are shown below. Contour map of 20th century trends for 1221 USHCN stations – “raw” data. […]