Roman’s Anomaly Regression and GISTEMP

(Post revised May 09, 2910: the simple linear regression line compared to Roman’s anomaly regression should of course have been fitted to the anomalies, not the raw data. Additional plots added. Must learn not to post in haste!)
At RomanM discusses monthly anomaly regression, and when asked in comments whether he has tried this with the GISS data replied

GISS is left as an exercise for the reader. 😉

Nick Stokes in another comment cautions

Your example relates to taking a set of station values, computing an anomaly presumably by CAM, and then running a regression through the serial monthly anomalies for individual stations. But who does this? Not Gistemp – I don’t believe they compute trends anywhere, and they use a more complicated anomaly procedure.

This post is the first in a series attempting Roman’s exercise, and at the same time providing a detailed example of the GISTEMP anomaly and urban station adjustment procedures which may be useful for anyone who would like to learn more about these procedures, but does not have the time or inclination to wade through the FORTRAN code.

The series will include the following  (and this table of contents will be updated as the relevant posts are added):

(These posts will be based on the March Gistemp analysis. Identification of stations as rural or urban remains an issue while the coordinates in v2.inv remain suspect. An appropriate disclaimer should be assumed where adjusted or gridded results are shown)

First, a couple of comparisons for GISTEMP stations. These two were chosen as examples where the trends found were substantial, and where Roman’s procedure gives a noticeably different trend when compared to simple linear regression. Trends shown here are degrees C per year.

Byrd Station

Byrd Station monthly anomalies. The heavy black line ("ALL") is the simple linear regression line, the other lines are fitted using Roman's procedure. 0.384843 is the trend from the simple linear regression line, 0.018286 the trend from Roman's procedure.

Smith RS, AL

Smith RS, AL. See previous caption for description of lines and trends

So using Roman’s procedure instead of the GISTEMP anomaly procedure (which is not the simple linear regression fit) could potentially produce different results. This will be explored in a later post in this series.

Revision: May 08, 2010. As can be seen  below, where both a simple linear regression line and Roman’s anomaly regression are fitted to from the anomalies, the trend using Roman’s anomaly regression is unchanged, but the trend derived from the simple linear regression line is now quite similar to that  derived by Roman’s anomaly regression.Byrd Station anomalies

Byrd Station anomalies
Smith Rs, AL anomalies

Smith Rs, AL anomalies

This entry was posted in Uncategorized and tagged , , . Bookmark the permalink.

2 Responses to Roman’s Anomaly Regression and GISTEMP

  1. vjones says:

    Hi Peter,
    Glad to see more from you. This is a great visual representation of the comparison. I’ve been hoping someone would have a go at station and box-level GISS data with Roman’s method. I’ve even wondered how long it would take me to learn enough R ‘from the basics up’.

    I have been wading through notes on the GIStemp anomaly proceedure just in the last few days. Every time I revisit this area I become a little more sure of what the effects of the procedure are, and I am now eagerly anticipating what you have to say on it.

  2. oneillp says:

    Peter O’Neill: This comment below is probably spam, and has been deleted as a posted comment, and the related url removed. Mail to fails, but if this is in fact a genuine comment, please post your actual question or provide a working e-mail address. ( appears to be a newly created site without content)

    Hello! Please e-mail me your contacts. I have a question

    Thank you!!!…

Leave a Reply

Fill in your details below or click an icon to log in: Logo

You are commenting using your account. Log Out /  Change )

Google+ photo

You are commenting using your Google+ account. Log Out /  Change )

Twitter picture

You are commenting using your Twitter account. Log Out /  Change )

Facebook photo

You are commenting using your Facebook account. Log Out /  Change )


Connecting to %s

This site uses Akismet to reduce spam. Learn how your comment data is processed.