<?xml version="1.0" encoding="UTF-8"?><rss version="2.0"
	xmlns:content="http://purl.org/rss/1.0/modules/content/"
	xmlns:dc="http://purl.org/dc/elements/1.1/"
	xmlns:atom="http://www.w3.org/2005/Atom"
	xmlns:sy="http://purl.org/rss/1.0/modules/syndication/"
	xmlns:georss="http://www.georss.org/georss" xmlns:geo="http://www.w3.org/2003/01/geo/wgs84_pos#" xmlns:media="http://search.yahoo.com/mrss/"
		>
<channel>
	<title>Comments on: This Gets Even More Amusing</title>
	<atom:link href="http://climateaudit.org/2008/10/25/this-gets-even-more-amusing/feed/" rel="self" type="application/rss+xml" />
	<link>http://climateaudit.org/2008/10/25/this-gets-even-more-amusing/</link>
	<description>by Steve McIntyre</description>
	<lastBuildDate>Sat, 25 May 2013 07:32:56 +0000</lastBuildDate>
	<sy:updatePeriod>hourly</sy:updatePeriod>
	<sy:updateFrequency>1</sy:updateFrequency>
	<generator>http://wordpress.com/</generator>
	<item>
		<title>By: Is the Santer 17 &#8216;d&#8217; test too liberal? Part 2. &#124; The Blackboard</title>
		<link>http://climateaudit.org/2008/10/25/this-gets-even-more-amusing/#comment-166817</link>
		<dc:creator><![CDATA[Is the Santer 17 &#8216;d&#8217; test too liberal? Part 2. &#124; The Blackboard]]></dc:creator>
		<pubDate>Mon, 03 Nov 2008 16:01:03 +0000</pubDate>
		<guid isPermaLink="false">http://www.climateaudit.org/?p=4195#comment-166817</guid>
		<description><![CDATA[[...] the paper has its merits, it also has some curious features some of which have been discussed at Climate Audit. I&#8217;ve mostly focused examining the conclusions we would draw if the method of Santer17 were [...]]]></description>
		<content:encoded><![CDATA[<p>[...] the paper has its merits, it also has some curious features some of which have been discussed at Climate Audit. I&#8217;ve mostly focused examining the conclusions we would draw if the method of Santer17 were [...]</p>
]]></content:encoded>
	</item>
	<item>
		<title>By: Kenneth Fritsch</title>
		<link>http://climateaudit.org/2008/10/25/this-gets-even-more-amusing/#comment-166816</link>
		<dc:creator><![CDATA[Kenneth Fritsch]]></dc:creator>
		<pubDate>Fri, 31 Oct 2008 15:13:09 +0000</pubDate>
		<guid isPermaLink="false">http://www.climateaudit.org/?p=4195#comment-166816</guid>
		<description><![CDATA[Re: &lt;a href=&quot;#comment-309021&quot; rel=&quot;nofollow&quot;&gt;Hu McCulloch (#64)&lt;/a&gt;,

Hu, I appreciate your professional thoughts on these analyses.  I do not have a comprehensive background on the theoretical underpinnings of some of these statistical issues and thus your explanations help this (very)layperson considerably.

In my view, I think that using the annual data to avoid the serial correlation that results from using the monthly data should reduce the uncertainty in the value of the adjusted trend standard deviation used in Santer et al.  I use an example below, that might be all wet from a theoretical standpoint (and I hope someone here evaluates that aspect and comments), but it makes my point.

Using the RSS T2LT time series from Santer et al. (2008) for monthly data from 1979-1999, one obtains an unadjusted trend standard deviation of 0.0307 degrees C per decade.  With a residual lag 1 versus residual correlation of  r = 0.886, an adjusted trend standard deviation of 0.132 is obtained.

Now, I go to the potentially illegitimate part.  If the lag1 residual versus residual regression is carried out, the correlation of r = 0.887 is obtained and with a range covering the 5-95% CI of 0.827 to 0.944.  Using those limits one can calculate the 5-95% limits on the adjusted trend standard deviation and obtain a range from 0.104 to 0.212 degrees C per decade.

If annual data is used, these uncertainties, if legitimate, are reduced to very small and unimportant values due to the sharply reduced serial correlations.  This approach goes along with UC&#039;s reminders that I have read at CA that in effect states that avoidance of autocorrelation is preferable to using corrections for it.]]></description>
		<content:encoded><![CDATA[<p>Re: <a href="#comment-309021" rel="nofollow">Hu McCulloch (#64)</a>,</p>
<p>Hu, I appreciate your professional thoughts on these analyses.  I do not have a comprehensive background on the theoretical underpinnings of some of these statistical issues and thus your explanations help this (very)layperson considerably.</p>
<p>In my view, I think that using the annual data to avoid the serial correlation that results from using the monthly data should reduce the uncertainty in the value of the adjusted trend standard deviation used in Santer et al.  I use an example below, that might be all wet from a theoretical standpoint (and I hope someone here evaluates that aspect and comments), but it makes my point.</p>
<p>Using the RSS T2LT time series from Santer et al. (2008) for monthly data from 1979-1999, one obtains an unadjusted trend standard deviation of 0.0307 degrees C per decade.  With a residual lag 1 versus residual correlation of  r = 0.886, an adjusted trend standard deviation of 0.132 is obtained.</p>
<p>Now, I go to the potentially illegitimate part.  If the lag1 residual versus residual regression is carried out, the correlation of r = 0.887 is obtained and with a range covering the 5-95% CI of 0.827 to 0.944.  Using those limits one can calculate the 5-95% limits on the adjusted trend standard deviation and obtain a range from 0.104 to 0.212 degrees C per decade.</p>
<p>If annual data is used, these uncertainties, if legitimate, are reduced to very small and unimportant values due to the sharply reduced serial correlations.  This approach goes along with UC&#8217;s reminders that I have read at CA that in effect states that avoidance of autocorrelation is preferable to using corrections for it.</p>
]]></content:encoded>
	</item>
	<item>
		<title>By: Pat Keating</title>
		<link>http://climateaudit.org/2008/10/25/this-gets-even-more-amusing/#comment-166815</link>
		<dc:creator><![CDATA[Pat Keating]]></dc:creator>
		<pubDate>Fri, 31 Oct 2008 02:49:29 +0000</pubDate>
		<guid isPermaLink="false">http://www.climateaudit.org/?p=4195#comment-166815</guid>
		<description><![CDATA[Re: &lt;a href=&quot;#comment-308995&quot; rel=&quot;nofollow&quot;&gt;Stan Palmer (#50)&lt;/a&gt;,
I was once co-author on a conference paper, but could not even attend the session because my security-clearance level was not high enough.]]></description>
		<content:encoded><![CDATA[<p>Re: <a href="#comment-308995" rel="nofollow">Stan Palmer (#50)</a>,<br />
I was once co-author on a conference paper, but could not even attend the session because my security-clearance level was not high enough.</p>
]]></content:encoded>
	</item>
	<item>
		<title>By: MC</title>
		<link>http://climateaudit.org/2008/10/25/this-gets-even-more-amusing/#comment-166814</link>
		<dc:creator><![CDATA[MC]]></dc:creator>
		<pubDate>Thu, 30 Oct 2008 23:02:10 +0000</pubDate>
		<guid isPermaLink="false">http://www.climateaudit.org/?p=4195#comment-166814</guid>
		<description><![CDATA[Re: &lt;a href=&quot;#comment-309021&quot; rel=&quot;nofollow&quot;&gt;Hu McCulloch (#64)&lt;/a&gt;, As I understand it the assumption is that the forcing of the atmosphere can be approximated as a linear trend over time. This of course assumes ENSO/higher order oscillations to be noise and where the signal + noise equation for both papers comes in. Interestingly the post I made earlier for MSU RSS data Re: &lt;a href=&quot;#comment-308727&quot; rel=&quot;nofollow&quot;&gt;MC (#34)&lt;/a&gt;, appears to show that there is a general troposphere ENSO signature independent of height with an additional linear trend
The equation could then be Y0 = ENSO function + linear trend + noise. The errors in the linear trend would be less as the autocorrelation effects would be much less provided the ENSO function was modelled. But then again that seems to be somewhat way off in the future by the look of things
Re: &lt;a href=&quot;#comment-309012&quot; rel=&quot;nofollow&quot;&gt;UC (#58)&lt;/a&gt;, Is that the new Eminem album cover?]]></description>
		<content:encoded><![CDATA[<p>Re: <a href="#comment-309021" rel="nofollow">Hu McCulloch (#64)</a>, As I understand it the assumption is that the forcing of the atmosphere can be approximated as a linear trend over time. This of course assumes ENSO/higher order oscillations to be noise and where the signal + noise equation for both papers comes in. Interestingly the post I made earlier for MSU RSS data Re: <a href="#comment-308727" rel="nofollow">MC (#34)</a>, appears to show that there is a general troposphere ENSO signature independent of height with an additional linear trend<br />
The equation could then be Y0 = ENSO function + linear trend + noise. The errors in the linear trend would be less as the autocorrelation effects would be much less provided the ENSO function was modelled. But then again that seems to be somewhat way off in the future by the look of things<br />
Re: <a href="#comment-309012" rel="nofollow">UC (#58)</a>, Is that the new Eminem album cover?</p>
]]></content:encoded>
	</item>
	<item>
		<title>By: Craig Loehle</title>
		<link>http://climateaudit.org/2008/10/25/this-gets-even-more-amusing/#comment-166813</link>
		<dc:creator><![CDATA[Craig Loehle]]></dc:creator>
		<pubDate>Thu, 30 Oct 2008 21:51:43 +0000</pubDate>
		<guid isPermaLink="false">http://www.climateaudit.org/?p=4195#comment-166813</guid>
		<description><![CDATA[Re: &lt;a href=&quot;#comment-309021&quot; rel=&quot;nofollow&quot;&gt;Hu McCulloch (#64)&lt;/a&gt;, In support of the authorship comment, a colleague was told by dept head that a book he wrote (alone) counted as 1 publication.]]></description>
		<content:encoded><![CDATA[<p>Re: <a href="#comment-309021" rel="nofollow">Hu McCulloch (#64)</a>, In support of the authorship comment, a colleague was told by dept head that a book he wrote (alone) counted as 1 publication.</p>
]]></content:encoded>
	</item>
	<item>
		<title>By: Hu McCulloch</title>
		<link>http://climateaudit.org/2008/10/25/this-gets-even-more-amusing/#comment-166812</link>
		<dc:creator><![CDATA[Hu McCulloch]]></dc:creator>
		<pubDate>Thu, 30 Oct 2008 20:36:00 +0000</pubDate>
		<guid isPermaLink="false">http://www.climateaudit.org/?p=4195#comment-166812</guid>
		<description><![CDATA[Ken Fritsch&#039;s numbers in #30 and 32 above reveal an interesting and important property of trend-coefficient standard errors:  Abstracting from serial correlation, the standard error of an estimate of a mean or the central intercept of a trend line will decrease roughly as 1/n^(1/2) as the sample size n increases, since every observation is equally informative.  However, the standard error of the estimate of the slope will decrease much faster as the observation period is increased, since the end points of a time trend are much more informative about its slope than are the more interior points.  Thus, adding points at either end increases the precision of the slope at a much faster rate than just adding new observations along the original line portion.  In fact, the standard deviation of the slope will decreae roughly as 1/n^(3/2).

This must be the primary reason Ken is getting much smaller standard errors using the full sample (to 2007 or so) than when he stops, as Santer, Nychka et al (2008) did, in 1999.  The autoregression adjustment is important, but since the autoregressive coefficient isn&#039;t much different for the two samples, it is not causing much change in the se.

As I noted in comment 64 of the Oct 22 thread &lt;a href=&quot;http://www.climateaudit.org/?p=4163&quot; rel=&quot;nofollow&quot;&gt;&quot;Replicating Santer Tables 1 and 3&quot;&lt;/a&gt;, the excellent Nychka, Santer, et al (2000) unpublished working paper on serial correlation found by Jean S demonstrates that the adjustment Santer, Nychka et al use in 2008 is in fact inadequate, so that the true standard errors are somewhat larger than are obtained using the 2008 adjustment.  It is unfortunate that Santer, Nychka et al did not heed Nychka, Santer et al!

As Jeff Id points out above, it should be borne in mind that linear trends are highly suspect as literal models of climatic data.  It would be sufficient for the purposes of global warming advocates to make a case that temperature has merely drifted up to a level that is significantly higher than it used to be, without requiring it to be a linear  trend.  Nevertheless, for the data at hand, the significance or insignificance of the &quot;trend&quot; (and differences between measures of the trend) is probably a descriptively adequate way of characterizing such a warming.

As for multiple co-authorship, an unfortunate incentive for co-authors to pile on is that universities (and research centers, I assume) often count raw citations as a measure of performance for salary purposes etc.  If they counted co-author-adjusted citations instead, eg giving each of the 17 authors of the 2008 paper 1/17 citation credit apiece for having co-authored it, gratuitious piling on would quickly come to an end!]]></description>
		<content:encoded><![CDATA[<p>Ken Fritsch&#8217;s numbers in #30 and 32 above reveal an interesting and important property of trend-coefficient standard errors:  Abstracting from serial correlation, the standard error of an estimate of a mean or the central intercept of a trend line will decrease roughly as 1/n^(1/2) as the sample size n increases, since every observation is equally informative.  However, the standard error of the estimate of the slope will decrease much faster as the observation period is increased, since the end points of a time trend are much more informative about its slope than are the more interior points.  Thus, adding points at either end increases the precision of the slope at a much faster rate than just adding new observations along the original line portion.  In fact, the standard deviation of the slope will decreae roughly as 1/n^(3/2).</p>
<p>This must be the primary reason Ken is getting much smaller standard errors using the full sample (to 2007 or so) than when he stops, as Santer, Nychka et al (2008) did, in 1999.  The autoregression adjustment is important, but since the autoregressive coefficient isn&#8217;t much different for the two samples, it is not causing much change in the se.</p>
<p>As I noted in comment 64 of the Oct 22 thread <a href="http://www.climateaudit.org/?p=4163" rel="nofollow">&#8220;Replicating Santer Tables 1 and 3&#8243;</a>, the excellent Nychka, Santer, et al (2000) unpublished working paper on serial correlation found by Jean S demonstrates that the adjustment Santer, Nychka et al use in 2008 is in fact inadequate, so that the true standard errors are somewhat larger than are obtained using the 2008 adjustment.  It is unfortunate that Santer, Nychka et al did not heed Nychka, Santer et al!</p>
<p>As Jeff Id points out above, it should be borne in mind that linear trends are highly suspect as literal models of climatic data.  It would be sufficient for the purposes of global warming advocates to make a case that temperature has merely drifted up to a level that is significantly higher than it used to be, without requiring it to be a linear  trend.  Nevertheless, for the data at hand, the significance or insignificance of the &#8220;trend&#8221; (and differences between measures of the trend) is probably a descriptively adequate way of characterizing such a warming.</p>
<p>As for multiple co-authorship, an unfortunate incentive for co-authors to pile on is that universities (and research centers, I assume) often count raw citations as a measure of performance for salary purposes etc.  If they counted co-author-adjusted citations instead, eg giving each of the 17 authors of the 2008 paper 1/17 citation credit apiece for having co-authored it, gratuitious piling on would quickly come to an end!</p>
]]></content:encoded>
	</item>
	<item>
		<title>By: PhilH</title>
		<link>http://climateaudit.org/2008/10/25/this-gets-even-more-amusing/#comment-166811</link>
		<dc:creator><![CDATA[PhilH]]></dc:creator>
		<pubDate>Thu, 30 Oct 2008 20:33:18 +0000</pubDate>
		<guid isPermaLink="false">http://www.climateaudit.org/?p=4195#comment-166811</guid>
		<description><![CDATA[Re: &lt;a href=&quot;#comment-309016&quot; rel=&quot;nofollow&quot;&gt;Sam Urbinto (#62)&lt;/a&gt;, Oh yes.]]></description>
		<content:encoded><![CDATA[<p>Re: <a href="#comment-309016" rel="nofollow">Sam Urbinto (#62)</a>, Oh yes.</p>
]]></content:encoded>
	</item>
	<item>
		<title>By: Sam Urbinto</title>
		<link>http://climateaudit.org/2008/10/25/this-gets-even-more-amusing/#comment-166810</link>
		<dc:creator><![CDATA[Sam Urbinto]]></dc:creator>
		<pubDate>Thu, 30 Oct 2008 19:23:24 +0000</pubDate>
		<guid isPermaLink="false">http://www.climateaudit.org/?p=4195#comment-166810</guid>
		<description><![CDATA[Re: &lt;a href=&quot;#comment-309014&quot; rel=&quot;nofollow&quot;&gt;PhilH (#60)&lt;/a&gt;,

But does anyone ever get buttered?]]></description>
		<content:encoded><![CDATA[<p>Re: <a href="#comment-309014" rel="nofollow">PhilH (#60)</a>,</p>
<p>But does anyone ever get buttered?</p>
]]></content:encoded>
	</item>
	<item>
		<title>By: lucia</title>
		<link>http://climateaudit.org/2008/10/25/this-gets-even-more-amusing/#comment-166809</link>
		<dc:creator><![CDATA[lucia]]></dc:creator>
		<pubDate>Thu, 30 Oct 2008 19:18:37 +0000</pubDate>
		<guid isPermaLink="false">http://www.climateaudit.org/?p=4195#comment-166809</guid>
		<description><![CDATA[Re: &lt;a href=&quot;#comment-309012&quot; rel=&quot;nofollow&quot;&gt;UC (#58)&lt;/a&gt;,
UC wins the thread!]]></description>
		<content:encoded><![CDATA[<p>Re: <a href="#comment-309012" rel="nofollow">UC (#58)</a>,<br />
UC wins the thread!</p>
]]></content:encoded>
	</item>
	<item>
		<title>By: PhilH</title>
		<link>http://climateaudit.org/2008/10/25/this-gets-even-more-amusing/#comment-166808</link>
		<dc:creator><![CDATA[PhilH]]></dc:creator>
		<pubDate>Thu, 30 Oct 2008 19:13:11 +0000</pubDate>
		<guid isPermaLink="false">http://www.climateaudit.org/?p=4195#comment-166808</guid>
		<description><![CDATA[Re: &lt;a href=&quot;#comment-308999&quot; rel=&quot;nofollow&quot;&gt;Allen63 (#52)&lt;/a&gt;, One caveat:  no one on the Team ever gets &quot;toasted.&quot;]]></description>
		<content:encoded><![CDATA[<p>Re: <a href="#comment-308999" rel="nofollow">Allen63 (#52)</a>, One caveat:  no one on the Team ever gets &#8220;toasted.&#8221;</p>
]]></content:encoded>
	</item>
</channel>
</rss>
