<?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: The Santer &quot;S.D.&quot;</title>
	<atom:link href="http://climateaudit.org/2008/10/24/the-santer-sd/feed/" rel="self" type="application/rss+xml" />
	<link>http://climateaudit.org/2008/10/24/the-santer-sd/</link>
	<description>by Steve McIntyre</description>
	<lastBuildDate>Sun, 26 May 2013 07:21:57 +0000</lastBuildDate>
	<sy:updatePeriod>hourly</sy:updatePeriod>
	<sy:updateFrequency>1</sy:updateFrequency>
	<generator>http://wordpress.com/</generator>
	<item>
		<title>By: Synthetic MSU Temperatures &#8211; Weighting Functions &#171; Trees for the Forest</title>
		<link>http://climateaudit.org/2008/10/24/the-santer-sd/#comment-166748</link>
		<dc:creator><![CDATA[Synthetic MSU Temperatures &#8211; Weighting Functions &#171; Trees for the Forest]]></dc:creator>
		<pubDate>Tue, 25 Aug 2009 23:44:35 +0000</pubDate>
		<guid isPermaLink="false">http://www.climateaudit.org/?p=4185#comment-166748</guid>
		<description><![CDATA[[...] M reported T2LT  weights from John Christy that are entirely pressure level specific.  The weights however correspond only [...]]]></description>
		<content:encoded><![CDATA[<p>[...] M reported T2LT  weights from John Christy that are entirely pressure level specific.  The weights however correspond only [...]</p>
]]></content:encoded>
	</item>
	<item>
		<title>By: The Mystery of CNRM 3.0 &#171; Scientific Prospective</title>
		<link>http://climateaudit.org/2008/10/24/the-santer-sd/#comment-166747</link>
		<dc:creator><![CDATA[The Mystery of CNRM 3.0 &#171; Scientific Prospective]]></dc:creator>
		<pubDate>Sat, 16 May 2009 21:35:32 +0000</pubDate>
		<guid isPermaLink="false">http://www.climateaudit.org/?p=4185#comment-166747</guid>
		<description><![CDATA[[...] into Matlab and computed area-weighted averages for the first 12 pressure levels. Steve M. reports here the T2LT weights by altitude via Christy. I multiplied the 12 time series by their respective [...]]]></description>
		<content:encoded><![CDATA[<p>[...] into Matlab and computed area-weighted averages for the first 12 pressure levels. Steve M. reports here the T2LT weights by altitude via Christy. I multiplied the 12 time series by their respective [...]</p>
]]></content:encoded>
	</item>
	<item>
		<title>By: Ryan O</title>
		<link>http://climateaudit.org/2008/10/24/the-santer-sd/#comment-166746</link>
		<dc:creator><![CDATA[Ryan O]]></dc:creator>
		<pubDate>Thu, 08 Jan 2009 18:51:08 +0000</pubDate>
		<guid isPermaLink="false">http://www.climateaudit.org/?p=4185#comment-166746</guid>
		<description><![CDATA[Something that I don&#039;t understand at all is why models that clearly do not reproduce the temperature profile in the troposphere/stratosphere are included in long-term forecasts?  Even if you use RSS data instead of UAH, many individual model runs fall well outside the RSS data.  Shouldn&#039;t the observational data provide a constraint on which model runs can be included in the ensemble for forecasting?  And if those model runs that fall outside the RSS data are excluded, what does that do to the 100-yr forecast and uncertainty?  How can you justify including model runs that produce unphysical results simply because the resulting ensemble mean uncertainty happens to lie within the observational uncertainty?]]></description>
		<content:encoded><![CDATA[<p>Something that I don&#8217;t understand at all is why models that clearly do not reproduce the temperature profile in the troposphere/stratosphere are included in long-term forecasts?  Even if you use RSS data instead of UAH, many individual model runs fall well outside the RSS data.  Shouldn&#8217;t the observational data provide a constraint on which model runs can be included in the ensemble for forecasting?  And if those model runs that fall outside the RSS data are excluded, what does that do to the 100-yr forecast and uncertainty?  How can you justify including model runs that produce unphysical results simply because the resulting ensemble mean uncertainty happens to lie within the observational uncertainty?</p>
]]></content:encoded>
	</item>
	<item>
		<title>By: UC</title>
		<link>http://climateaudit.org/2008/10/24/the-santer-sd/#comment-166745</link>
		<dc:creator><![CDATA[UC]]></dc:creator>
		<pubDate>Mon, 27 Oct 2008 12:56:59 +0000</pubDate>
		<guid isPermaLink="false">http://www.climateaudit.org/?p=4185#comment-166745</guid>
		<description><![CDATA[No one objected, so we can conclude that this two sample t-test with SEs turns to prediction interval SD test in the case of equal variances (one future observation ). Douglass and Santer are testing trends instead of means, so degrees of freedom need to be corrected accordingly. In addition, variances are not assumed equal, so dof question  gets even more complicated (see for example http://www.itl.nist.gov/div898/handbook/eda/section3/eda353.htm )

Douglass et al assume that observed b_o is the true trend, and in this case one-sample t-test can be performed.  With large n_2, t-distribution approaches rapidly normal d., and the 95 % confidence interval for the trend would be

$latex \bar{x}_2 - \frac{1.96\sigma}{\sqrt{n_2}}   \leq b_o \leq \bar{x}_2 + \frac{1.96\sigma}{\sqrt{n_2}}  $

Length of this interval approaches zero as n_2 gets larger, but there is nothing special about that. That&#039;s what should happen under this assumption.]]></description>
		<content:encoded><![CDATA[<p>No one objected, so we can conclude that this two sample t-test with SEs turns to prediction interval SD test in the case of equal variances (one future observation ). Douglass and Santer are testing trends instead of means, so degrees of freedom need to be corrected accordingly. In addition, variances are not assumed equal, so dof question  gets even more complicated (see for example <a href="http://www.itl.nist.gov/div898/handbook/eda/section3/eda353.htm" rel="nofollow">http://www.itl.nist.gov/div898/handbook/eda/section3/eda353.htm</a> )</p>
<p>Douglass et al assume that observed b_o is the true trend, and in this case one-sample t-test can be performed.  With large n_2, t-distribution approaches rapidly normal d., and the 95 % confidence interval for the trend would be</p>
<p><img src='http://s0.wp.com/latex.php?latex=%5Cbar%7Bx%7D_2+-+%5Cfrac%7B1.96%5Csigma%7D%7B%5Csqrt%7Bn_2%7D%7D+++%5Cleq+b_o+%5Cleq+%5Cbar%7Bx%7D_2+%2B+%5Cfrac%7B1.96%5Csigma%7D%7B%5Csqrt%7Bn_2%7D%7D++&amp;bg=ffffff&amp;fg=000&amp;s=0' alt='&#92;bar{x}_2 - &#92;frac{1.96&#92;sigma}{&#92;sqrt{n_2}}   &#92;leq b_o &#92;leq &#92;bar{x}_2 + &#92;frac{1.96&#92;sigma}{&#92;sqrt{n_2}}  ' title='&#92;bar{x}_2 - &#92;frac{1.96&#92;sigma}{&#92;sqrt{n_2}}   &#92;leq b_o &#92;leq &#92;bar{x}_2 + &#92;frac{1.96&#92;sigma}{&#92;sqrt{n_2}}  ' class='latex' /></p>
<p>Length of this interval approaches zero as n_2 gets larger, but there is nothing special about that. That&#8217;s what should happen under this assumption.</p>
]]></content:encoded>
	</item>
	<item>
		<title>By: UC</title>
		<link>http://climateaudit.org/2008/10/24/the-santer-sd/#comment-166744</link>
		<dc:creator><![CDATA[UC]]></dc:creator>
		<pubDate>Sun, 26 Oct 2008 18:51:54 +0000</pubDate>
		<guid isPermaLink="false">http://www.climateaudit.org/?p=4185#comment-166744</guid>
		<description><![CDATA[I have to back to basics to get an idea what is going on in Santer17, and SD vs. SE debate.

Let&#039;s take two samples of size $latex n_1 $ and $latex n_2 $, compute sample means and variances, $latex \bar{x}_1 $, $latex s_1^2 $ and $latex \bar{x}_2 $, $latex s_2^2 $. If these are from same distribution, the statistic

$latex y=\frac{\bar{x}_1-\bar{x}_2}{(s^2(\frac{1}{n_1}+\frac{1}{n_2}))^{1/2}} $


(where $latex s^2=(n_1s_1^2+n_2s_2^2)/(n_1+n_2-2) $ ), follows a Student&#039;s t with $latex n_1+n_2-2 $ dof. If $latex n_1=1 $, this simplifies to

$latex y=\frac{\bar{x}_1-\bar{x}_2}{(s^2(\frac{1}{1}+\frac{1}{n_2}))^{1/2}} $

This test is more SE like than SD. But when you turn this to as prediction interval for one future sample ($latex x_1 $), you&#039;ll get

$latex \bar{x}_2-t_{\alpha/2}(n_2-1)s\sqrt{1+\frac{1}{n_2}}\leq x_1 \leq \bar{x}_2+t_{\alpha/2}(n_2-1)s\sqrt{1+\frac{1}{n_2}} $

and I think test based on this or y in the above is more or less the same thing. Proof is left to the reader, though ;) This is more SD-like test, specially if $latex n_2 $ is large.

In Santer17, the variances seem to be different for those two samples, so the situation gets more complicated. And it is a bit unclear to me how sample deviations are obtained in H2 case.]]></description>
		<content:encoded><![CDATA[<p>I have to back to basics to get an idea what is going on in Santer17, and SD vs. SE debate.</p>
<p>Let&#8217;s take two samples of size <img src='http://s0.wp.com/latex.php?latex=n_1+&amp;bg=ffffff&amp;fg=000&amp;s=0' alt='n_1 ' title='n_1 ' class='latex' /> and <img src='http://s0.wp.com/latex.php?latex=n_2+&amp;bg=ffffff&amp;fg=000&amp;s=0' alt='n_2 ' title='n_2 ' class='latex' />, compute sample means and variances, <img src='http://s0.wp.com/latex.php?latex=%5Cbar%7Bx%7D_1+&amp;bg=ffffff&amp;fg=000&amp;s=0' alt='&#92;bar{x}_1 ' title='&#92;bar{x}_1 ' class='latex' />, <img src='http://s0.wp.com/latex.php?latex=s_1%5E2+&amp;bg=ffffff&amp;fg=000&amp;s=0' alt='s_1^2 ' title='s_1^2 ' class='latex' /> and <img src='http://s0.wp.com/latex.php?latex=%5Cbar%7Bx%7D_2+&amp;bg=ffffff&amp;fg=000&amp;s=0' alt='&#92;bar{x}_2 ' title='&#92;bar{x}_2 ' class='latex' />, <img src='http://s0.wp.com/latex.php?latex=s_2%5E2+&amp;bg=ffffff&amp;fg=000&amp;s=0' alt='s_2^2 ' title='s_2^2 ' class='latex' />. If these are from same distribution, the statistic</p>
<p><img src='http://s0.wp.com/latex.php?latex=y%3D%5Cfrac%7B%5Cbar%7Bx%7D_1-%5Cbar%7Bx%7D_2%7D%7B%28s%5E2%28%5Cfrac%7B1%7D%7Bn_1%7D%2B%5Cfrac%7B1%7D%7Bn_2%7D%29%29%5E%7B1%2F2%7D%7D+&amp;bg=ffffff&amp;fg=000&amp;s=0' alt='y=&#92;frac{&#92;bar{x}_1-&#92;bar{x}_2}{(s^2(&#92;frac{1}{n_1}+&#92;frac{1}{n_2}))^{1/2}} ' title='y=&#92;frac{&#92;bar{x}_1-&#92;bar{x}_2}{(s^2(&#92;frac{1}{n_1}+&#92;frac{1}{n_2}))^{1/2}} ' class='latex' /></p>
<p>(where <img src='http://s0.wp.com/latex.php?latex=s%5E2%3D%28n_1s_1%5E2%2Bn_2s_2%5E2%29%2F%28n_1%2Bn_2-2%29+&amp;bg=ffffff&amp;fg=000&amp;s=0' alt='s^2=(n_1s_1^2+n_2s_2^2)/(n_1+n_2-2) ' title='s^2=(n_1s_1^2+n_2s_2^2)/(n_1+n_2-2) ' class='latex' /> ), follows a Student&#8217;s t with <img src='http://s0.wp.com/latex.php?latex=n_1%2Bn_2-2+&amp;bg=ffffff&amp;fg=000&amp;s=0' alt='n_1+n_2-2 ' title='n_1+n_2-2 ' class='latex' /> dof. If <img src='http://s0.wp.com/latex.php?latex=n_1%3D1+&amp;bg=ffffff&amp;fg=000&amp;s=0' alt='n_1=1 ' title='n_1=1 ' class='latex' />, this simplifies to</p>
<p><img src='http://s0.wp.com/latex.php?latex=y%3D%5Cfrac%7B%5Cbar%7Bx%7D_1-%5Cbar%7Bx%7D_2%7D%7B%28s%5E2%28%5Cfrac%7B1%7D%7B1%7D%2B%5Cfrac%7B1%7D%7Bn_2%7D%29%29%5E%7B1%2F2%7D%7D+&amp;bg=ffffff&amp;fg=000&amp;s=0' alt='y=&#92;frac{&#92;bar{x}_1-&#92;bar{x}_2}{(s^2(&#92;frac{1}{1}+&#92;frac{1}{n_2}))^{1/2}} ' title='y=&#92;frac{&#92;bar{x}_1-&#92;bar{x}_2}{(s^2(&#92;frac{1}{1}+&#92;frac{1}{n_2}))^{1/2}} ' class='latex' /></p>
<p>This test is more SE like than SD. But when you turn this to as prediction interval for one future sample (<img src='http://s0.wp.com/latex.php?latex=x_1+&amp;bg=ffffff&amp;fg=000&amp;s=0' alt='x_1 ' title='x_1 ' class='latex' />), you&#8217;ll get</p>
<p><img src='http://s0.wp.com/latex.php?latex=%5Cbar%7Bx%7D_2-t_%7B%5Calpha%2F2%7D%28n_2-1%29s%5Csqrt%7B1%2B%5Cfrac%7B1%7D%7Bn_2%7D%7D%5Cleq+x_1+%5Cleq+%5Cbar%7Bx%7D_2%2Bt_%7B%5Calpha%2F2%7D%28n_2-1%29s%5Csqrt%7B1%2B%5Cfrac%7B1%7D%7Bn_2%7D%7D+&amp;bg=ffffff&amp;fg=000&amp;s=0' alt='&#92;bar{x}_2-t_{&#92;alpha/2}(n_2-1)s&#92;sqrt{1+&#92;frac{1}{n_2}}&#92;leq x_1 &#92;leq &#92;bar{x}_2+t_{&#92;alpha/2}(n_2-1)s&#92;sqrt{1+&#92;frac{1}{n_2}} ' title='&#92;bar{x}_2-t_{&#92;alpha/2}(n_2-1)s&#92;sqrt{1+&#92;frac{1}{n_2}}&#92;leq x_1 &#92;leq &#92;bar{x}_2+t_{&#92;alpha/2}(n_2-1)s&#92;sqrt{1+&#92;frac{1}{n_2}} ' class='latex' /></p>
<p>and I think test based on this or y in the above is more or less the same thing. Proof is left to the reader, though <img src='http://s1.wp.com/wp-includes/images/smilies/icon_wink.gif' alt=';)' class='wp-smiley' />  This is more SD-like test, specially if <img src='http://s0.wp.com/latex.php?latex=n_2+&amp;bg=ffffff&amp;fg=000&amp;s=0' alt='n_2 ' title='n_2 ' class='latex' /> is large.</p>
<p>In Santer17, the variances seem to be different for those two samples, so the situation gets more complicated. And it is a bit unclear to me how sample deviations are obtained in H2 case.</p>
]]></content:encoded>
	</item>
	<item>
		<title>By: UC</title>
		<link>http://climateaudit.org/2008/10/24/the-santer-sd/#comment-166743</link>
		<dc:creator><![CDATA[UC]]></dc:creator>
		<pubDate>Sun, 26 Oct 2008 16:05:36 +0000</pubDate>
		<guid isPermaLink="false">http://www.climateaudit.org/?p=4185#comment-166743</guid>
		<description><![CDATA[Re: &lt;a href=&quot;#comment-308248&quot; rel=&quot;nofollow&quot;&gt;lucia (#18)&lt;/a&gt;,



&lt;blockquote&gt;Did someone actually say that weather with AR(1) and lag 1 autocorrelation of 0.9 is reasonable?
&lt;/blockquote&gt;

I get that impression from Santer17. Confusing stuff, just some time ago someone told us that &lt;em&gt;The conclusion is inescapable, that global temperature cannot be adequately modeled as a linear trend plus AR(1) process.&lt;/em&gt;]]></description>
		<content:encoded><![CDATA[<p>Re: <a href="#comment-308248" rel="nofollow">lucia (#18)</a>,</p>
<blockquote><p>Did someone actually say that weather with AR(1) and lag 1 autocorrelation of 0.9 is reasonable?
</p></blockquote>
<p>I get that impression from Santer17. Confusing stuff, just some time ago someone told us that <em>The conclusion is inescapable, that global temperature cannot be adequately modeled as a linear trend plus AR(1) process.</em></p>
]]></content:encoded>
	</item>
	<item>
		<title>By: lucia</title>
		<link>http://climateaudit.org/2008/10/24/the-santer-sd/#comment-166742</link>
		<dc:creator><![CDATA[lucia]]></dc:creator>
		<pubDate>Sun, 26 Oct 2008 15:26:48 +0000</pubDate>
		<guid isPermaLink="false">http://www.climateaudit.org/?p=4185#comment-166742</guid>
		<description><![CDATA[Re: &lt;a href=&quot;#comment-308143&quot; rel=&quot;nofollow&quot;&gt;UC (#15)&lt;/a&gt;,

&lt;blockquote&gt;The interesting point is that AR(1) with p of about 0.9 is suddenly accepted for &#039;climate noise&#039;.&lt;/blockquote&gt;

Did someone actually say that weather with AR(1) and lag 1 autocorrelation of 0.9 is reasonable?

If you were to run AR(1) simulations with this lag 1 autocorrelation, the autocorrelation for observataions of GMST since GMST would be highly unlikely. &lt;em&gt;Highly.&lt;/em&gt;.  Heck, if &quot;weahter noise&quot; is AR(1) with autocorrelation of 0.728, there&#039;s only a 1.7% chance of getting lag 1 autocorrelations as low as we&#039;ve gotten since 2001.

You can read a bit about the analysis &lt;a href=&quot;http://rankexploits.com/musings/2008/result-of-boring-series-gavins-closer-process-falsifies/&quot; rel=&quot;nofollow&quot;&gt;here&lt;/a&gt;.

The autocorrealtions are higher during periods when volcanos like Pinatubo and El Chicon are going off. Otherwise... well... The observational evidence suggests it&#039;s lower.]]></description>
		<content:encoded><![CDATA[<p>Re: <a href="#comment-308143" rel="nofollow">UC (#15)</a>,</p>
<blockquote><p>The interesting point is that AR(1) with p of about 0.9 is suddenly accepted for &#8216;climate noise&#8217;.</p></blockquote>
<p>Did someone actually say that weather with AR(1) and lag 1 autocorrelation of 0.9 is reasonable?</p>
<p>If you were to run AR(1) simulations with this lag 1 autocorrelation, the autocorrelation for observataions of GMST since GMST would be highly unlikely. <em>Highly.</em>.  Heck, if &#8220;weahter noise&#8221; is AR(1) with autocorrelation of 0.728, there&#8217;s only a 1.7% chance of getting lag 1 autocorrelations as low as we&#8217;ve gotten since 2001.</p>
<p>You can read a bit about the analysis <a href="http://rankexploits.com/musings/2008/result-of-boring-series-gavins-closer-process-falsifies/" rel="nofollow">here</a>.</p>
<p>The autocorrealtions are higher during periods when volcanos like Pinatubo and El Chicon are going off. Otherwise&#8230; well&#8230; The observational evidence suggests it&#8217;s lower.</p>
]]></content:encoded>
	</item>
	<item>
		<title>By: Vincent Guerrini Jr</title>
		<link>http://climateaudit.org/2008/10/24/the-santer-sd/#comment-166741</link>
		<dc:creator><![CDATA[Vincent Guerrini Jr]]></dc:creator>
		<pubDate>Sun, 26 Oct 2008 06:41:58 +0000</pubDate>
		<guid isPermaLink="false">http://www.climateaudit.org/?p=4185#comment-166741</guid>
		<description><![CDATA[Steve RE:Sea ice posting above not correct thread here etc... tried your sea ice stretch # 3 to put comment there are 1034 comments ice version #3 and it ends... so not possible unless sea ice stretch version 4? thanks for advice]]></description>
		<content:encoded><![CDATA[<p>Steve RE:Sea ice posting above not correct thread here etc&#8230; tried your sea ice stretch # 3 to put comment there are 1034 comments ice version #3 and it ends&#8230; so not possible unless sea ice stretch version 4? thanks for advice</p>
]]></content:encoded>
	</item>
	<item>
		<title>By: Mike B</title>
		<link>http://climateaudit.org/2008/10/24/the-santer-sd/#comment-166740</link>
		<dc:creator><![CDATA[Mike B]]></dc:creator>
		<pubDate>Sun, 26 Oct 2008 01:44:34 +0000</pubDate>
		<guid isPermaLink="false">http://www.climateaudit.org/?p=4185#comment-166740</guid>
		<description><![CDATA[Re: &lt;a href=&quot;#comment-308137&quot; rel=&quot;nofollow&quot;&gt;craig loehle (#13)&lt;/a&gt;,

I think beaker was sucked through a wormhole into a parallel universe around the time he posited that Santer et. al. contained a typo.  He might be having some trouble posting from there.]]></description>
		<content:encoded><![CDATA[<p>Re: <a href="#comment-308137" rel="nofollow">craig loehle (#13)</a>,</p>
<p>I think beaker was sucked through a wormhole into a parallel universe around the time he posited that Santer et. al. contained a typo.  He might be having some trouble posting from there.</p>
]]></content:encoded>
	</item>
	<item>
		<title>By: UC</title>
		<link>http://climateaudit.org/2008/10/24/the-santer-sd/#comment-166739</link>
		<dc:creator><![CDATA[UC]]></dc:creator>
		<pubDate>Sat, 25 Oct 2008 21:11:41 +0000</pubDate>
		<guid isPermaLink="false">http://www.climateaudit.org/?p=4185#comment-166739</guid>
		<description><![CDATA[Hmmm, this problem should be translated to a test about prediction interval for a single future observation. Or tolerance interval, for Bayesians it doesn&#039;t make a difference ;)

The interesting point is that AR(1) with p of about 0.9 is suddenly accepted for &#039;climate noise&#039;. Mann &amp; Lees 96 paper told us that such value would be unphysical. The problem is, you need high p to keep AGW running even if the temperatures go down. On the other hand, you need low p to be able to say &lt;em&gt;&#039;one cannot simulate the evolution of the climate over last 30 years without including in the simulations mankind&#039;s influence on sulfate aerosols and greenhouse gases.&#039;&lt;/em&gt;]]></description>
		<content:encoded><![CDATA[<p>Hmmm, this problem should be translated to a test about prediction interval for a single future observation. Or tolerance interval, for Bayesians it doesn&#8217;t make a difference <img src='http://s1.wp.com/wp-includes/images/smilies/icon_wink.gif' alt=';)' class='wp-smiley' /> </p>
<p>The interesting point is that AR(1) with p of about 0.9 is suddenly accepted for &#8216;climate noise&#8217;. Mann &amp; Lees 96 paper told us that such value would be unphysical. The problem is, you need high p to keep AGW running even if the temperatures go down. On the other hand, you need low p to be able to say <em>&#8216;one cannot simulate the evolution of the climate over last 30 years without including in the simulations mankind&#8217;s influence on sulfate aerosols and greenhouse gases.&#8217;</em></p>
]]></content:encoded>
	</item>
</channel>
</rss>
