Comments on: ARMA (1,1) Coefficients for CRU Gridcells

By: Steve McIntyre

Steve McIntyre — Sun, 14 Aug 2005 19:24:07 +0000

###CALCULATION OF ARMA COEFFICIENTS load("c:/climate/data/jones/hadcruv.tab") stat<-array(rep(NA,3*2592),dim=c(2592,3)) count<-apply(!is.na(v),2,sum) #this calculaiton had to be re-started on a couple of occasions due to arima calculation problems; I just determined where it stopped ansd started at the next number for (i in 1:2592) { if ( count[i]>50) stat[i,]<- coef( arima(v[,i],order=c(1,0,1)) ) } dummy<-array(c(1:2592,jonesinv(1:2592)),dim=c(2592,3)) stat<-cbind(dummy,stat) trend<-rep(NA,2592) #this is OLS trend year<-(1:1769)/12 for (i in 1:2592) { if (count[i]>50) { fm<-lm(v[,i]~ year) trend[i]<-coef(fm)[2] } } stat<-cbind(stat,trend) arma.coefficients<-stat save(arma.coefficients,file="c:/climate/data/jones/arma.coefficients.tab")##PLOT BY LATITUDE nf <- layout(array(1:2,dim=c(2,1)),heights=c(1.1,1.3)) layout.show(nf) par(mar=c(0,4,2,2)) plot(stat[,1],stat[,3],xlab="",ylab="AR1 Coef",axes=FALSE) axis(side=2) axis(side=1,at=seq(-60,+60,30),labels=FALSE) box() par(mar=c(4,4,0,2)) plot(stat[,1],stat[,4],xlab="",ylab="MA1 Coef",axes=FALSE) axis(side=2) axis(side=1,at=seq(-60,+60,30)) box()

By: Sandy Taylor

Sandy Taylor — Sun, 14 Aug 2005 03:19:05 +0000

You wanted one of the 500,000, well I’m about 200 of those.
Thanks for your efforts to keep the PC world and their true believers honest. I know how hard it can be to be the one eyed man in the world of the blind.

I’ve been reading this blog for a few months. I’m naturally skeptical and found this site through at general search on the topic of global warming. It’s one of the few I come back to and check at least weekly. I don’t trust the “global warming is obvious believers” about how obvious it is, but I am still willing to read their sites. No matter how many times I email them articles to read, they never do (move into the light, Carol Anne). I did have my faith shaken by reading Roy Spencer, a big U Alabama climate scientist, come down on the ID side of the evolution debate, but I’m still skeptical.
Keep on keeping them honest, if you don’t I don’t know who will.

By: JerryB

JerryB — Sat, 13 Aug 2005 23:50:16 +0000

> There’s lots of interesting exercises with the temperature data set.

Yup!

> I wish that I’d left it alone.

Best punchline of the month!

By: Ross McKitrick

Ross McKitrick — Sat, 13 Aug 2005 22:20:19 +0000

This graphic reminds me of a similar one in Karner’s latest paper (Central European Jnl of Physics, see http://www.aai.ee/~olavi/ ffor pdf) where he graphs the Hurst coefficient H by latitude band. People have been looking for evidence of polar greenhouse effect amplification by looking for increased trend coefficients near the poles. But increased trend coefficients is consistent with other hypotheses as well so it’s a confounded test. It seems to me a more precise test is decreased speed of mean-reversion, which is to say that, given a + greenhouse shock there is a slow, and possibly zero-speed reversion to the previous mean (ie nonstationarity). If greenhouse warming is amplified and more dominant in the polar regions then persistency should increase, or in this case the AR1 coefficient should approach 1, though it’s not measuring quite the same thing as H. A lower AR coefficient near the poles implies stronger mean-reversion on a short time scale; the negative AR coefficient just implies an oscillation rather than a decay back to the mean. If polar amplification is really so strong at the surface it seems to me the graphs should be U-shaped.

By: Peter Hartley

Peter Hartley — Sat, 13 Aug 2005 18:30:17 +0000

An interesting research exercise might be to take the satellite temperatures for different latitude bands (if not CRU grid cells), estimate separate ARMA models for them and then compare them to the corresponding ARMA processes for the CRU data as you heve presented them here. Do you know whether it is easy to get the satellite data by latitude and longitude grid cells?

Steve: Yes. It retrieves quite easily. The cell size is 2.5 x 2.5 instead of 5×5 so a little fiddling is needed to make them compatible. There’s lots of interesting exercises with the temperature data set. I wish that I’d left it alone.