{"id":20232,"date":"2012-07-13T10:07:54","date_gmt":"2012-07-13T16:07:54","guid":{"rendered":"http:\/\/rankexploits.com\/musings\/?p=20232"},"modified":"2012-07-13T14:53:20","modified_gmt":"2012-07-13T20:53:20","slug":"effect-of-covariance-screening-on-method-i","status":"publish","type":"post","link":"https:\/\/rankexploits.com\/musings\/2012\/effect-of-covariance-screening-on-method-i\/","title":{"rendered":"Effect of covariance &#038; Screening on Method I."},"content":{"rendered":"<p>Bo emailed me with questions and comments on my <a href=\"http:\/\/rankexploits.com\/musings\/2012\/methods-suppose-we-had-5000-proxies\/\">previous post<\/a>. He was curious what happens to some of my graphs if I remove the covariance between the target temperature and the proxy response parameter $latex \\lambda_{i} $. I thought you guys might be too. So I created these graphs. Each compares the screened by correlation during the calibration period and unscreened results if I screening using <I>local<\/I> temperatures. Each case is done at a different magntidue of correlation between $latex \\lambda_{i} $ and $latex T_{i} $. (The convention followed in the legend is negative is based on the <i>trend<\/i> during the calibration period which is a sign error. But I started that way and haven&#8217;t fixed the sign error in the legend yet.)  I&#8217;m showing screening with local values because that a <I>lighter<\/I> level of screening (and so better!) <\/p>\n<p>Below, I&#8217;ve compared the synthetic results with  $latex cov [ \\lambda_{i} $ and $latex mean[T_{i}] = 0. ] $.<\/p>\n<p><a href=\"http:\/\/rankexploits.com\/musings\/wp-content\/uploads\/2012\/07\/MethodI_covaraince_none.png\"><img loading=\"lazy\" decoding=\"async\" src=\"http:\/\/rankexploits.com\/musings\/wp-content\/uploads\/2012\/07\/MethodI_covaraince_none-500x500.png\" alt=\"\" title=\"MethodI_covaraince_none\" width=\"500\" height=\"500\" class=\"aligncenter size-medium wp-image-20237\" srcset=\"https:\/\/rankexploits.com\/musings\/wp-content\/uploads\/2012\/07\/MethodI_covaraince_none-500x500.png 500w, https:\/\/rankexploits.com\/musings\/wp-content\/uploads\/2012\/07\/MethodI_covaraince_none-300x300.png 300w, https:\/\/rankexploits.com\/musings\/wp-content\/uploads\/2012\/07\/MethodI_covaraince_none.png 1008w\" sizes=\"auto, (max-width: 500px) 100vw, 500px\" \/>  <\/a><\/p>\n<p>Above you can see that the screening picked out 4123 out of 5000 proxies (false rejecting 877 as &#8216;bad&#8217;).  The consequence is the screened results have a slight bias during the reconstruction period; this bias results in the aqua trace being consistently higher than the black line representing the target. In contrast, the red line is unbiased. <\/p>\n<p><strong>Positive covariance between &#8220;T and gamma&#8221;:<\/strong><br \/>\nI repeated the exact same runs, but this time set the synthetic runs to make proxies were &#8220;warm&#8221; during the pre-calibration periods more reponsive to temperature variations and those that were &#8220;cold&#8221; during the historic period less so.<br \/>\n<a href=\"http:\/\/rankexploits.com\/musings\/wp-content\/uploads\/2012\/07\/MethodI_positive_covariance.png\"><img loading=\"lazy\" decoding=\"async\" src=\"http:\/\/rankexploits.com\/musings\/wp-content\/uploads\/2012\/07\/MethodI_positive_covariance-500x500.png\" alt=\"\" title=\"MethodI_positive_covariance\" width=\"500\" height=\"500\" class=\"aligncenter size-medium wp-image-20233\" srcset=\"https:\/\/rankexploits.com\/musings\/wp-content\/uploads\/2012\/07\/MethodI_positive_covariance-500x500.png 500w, https:\/\/rankexploits.com\/musings\/wp-content\/uploads\/2012\/07\/MethodI_positive_covariance-300x300.png 300w, https:\/\/rankexploits.com\/musings\/wp-content\/uploads\/2012\/07\/MethodI_positive_covariance.png 1008w\" sizes=\"auto, (max-width: 500px) 100vw, 500px\" \/><\/a><br \/>\n several things:<\/p>\n<ol>\n<li>Looking at the red (unscreened) trace, we can see method I is too warm during the reconstruction period. But also, the reconstruction in method I does not follow the target well during the calibration period. What this means is that when the bias exists, someone is likely to notice the issue. (Inspecting the algebra shows why this happens.)  This effect could be called the the fundamental bias in method I. It&#8217;s in the math.  This effect is peculiar to method I (and is a reason one might wish to avoid method I.)<\/li>\n<li>Looking at the aqua trace wee see that method I becomes <em>even warmer<\/em> when the proxies are screened. This extra warming over an above the amount seen in the red line happens because screening  picks out proxies that response better to the local temperature during the calibration period.  Absent covariance local temperature and the proxy, this results in a warming bias seen above.  Since the responsiveness is better at higher temperatures, the 4123 proxies that were retained were on average warmer during the reconstruction period than the original 5000 proxies. This results in a warm bias. And on top of that, we have the warming bias introduced by the existence of the correlation which we see in the red trace. So we have three separate effects all resulting in warming bias when correlation between the average temperature at a proxy during the reconstruction period and it&#8217;s responsiveness to temperature is positive.\n <\/li>\n<\/ol>\n<p>Note that when correlation exists <i>and<\/I> and we screen we will find that the underlying temperature at the selected proxies (i.e. the the 4123 proxies ones retained here) is biased relative to the full 5000 samples. This will happen with any method and is not unique to method A.   So, when this correlation is not zero, some bias can be introduced into many methods even if it&#8217;s not introduced absent screening. (Some methods may, however, be immune to this. We would need to examine each one individually.)<\/p>\n<p>Bo also wants to see how things look if I use fewer proxies&#8211; something closer to the number of proxies one might obtain for a real reconstruction. I&#8217;ll be posting some of those on Monday. I need to figure out how I want to show that because with fewer proxies, I think it&#8217;s best to show several realizations for each method.  <\/p>\n<p>I want to look add method IV for rescaling results&#8211; it will match the final step in Mann08.  I haven&#8217;t yet figured out how I do RegEm with proxies. But I&#8217;ll think about that. I do want to see how each method works on toy problems. \ud83d\ude42<\/p>\n","protected":false},"excerpt":{"rendered":"<p>Bo emailed me with questions and comments on my previous post. He was curious what happens to some of my graphs if I remove the covariance between the target temperature and the proxy response parameter $latex \\lambda_{i} $. I thought you guys might be too. So I created these graphs. Each compares the screened by &hellip; <a href=\"https:\/\/rankexploits.com\/musings\/2012\/effect-of-covariance-screening-on-method-i\/\" class=\"more-link\">Continue reading <span class=\"screen-reader-text\">Effect of covariance &#038; Screening on Method I.<\/span> <span class=\"meta-nav\">&rarr;<\/span><\/a><\/p>\n","protected":false},"author":2,"featured_media":0,"comment_status":"closed","ping_status":"closed","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[3],"tags":[420],"class_list":["post-20232","post","type-post","status-publish","format-standard","hentry","category-statistics","tag-screening"],"_links":{"self":[{"href":"https:\/\/rankexploits.com\/musings\/wp-json\/wp\/v2\/posts\/20232","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/rankexploits.com\/musings\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/rankexploits.com\/musings\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/rankexploits.com\/musings\/wp-json\/wp\/v2\/users\/2"}],"replies":[{"embeddable":true,"href":"https:\/\/rankexploits.com\/musings\/wp-json\/wp\/v2\/comments?post=20232"}],"version-history":[{"count":0,"href":"https:\/\/rankexploits.com\/musings\/wp-json\/wp\/v2\/posts\/20232\/revisions"}],"wp:attachment":[{"href":"https:\/\/rankexploits.com\/musings\/wp-json\/wp\/v2\/media?parent=20232"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/rankexploits.com\/musings\/wp-json\/wp\/v2\/categories?post=20232"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/rankexploits.com\/musings\/wp-json\/wp\/v2\/tags?post=20232"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}