{"id":20249,"date":"2012-07-13T16:21:18","date_gmt":"2012-07-13T22:21:18","guid":{"rendered":"http:\/\/rankexploits.com\/musings\/?p=20249"},"modified":"2012-07-13T16:21:18","modified_gmt":"2012-07-13T22:21:18","slug":"cps-with-variance-matching-will-be-method-iv","status":"publish","type":"post","link":"https:\/\/rankexploits.com\/musings\/2012\/cps-with-variance-matching-will-be-method-iv\/","title":{"rendered":"CPS with &#8216;variance matching&#8217;: (Will be method IV)"},"content":{"rendered":"<p>When reading Mann08, I realized I need a &#8220;method IV&#8221; in my previous graphs. In method IV, I will rescale the mean proxy value so the variance of the proxy values match the variance of the temperature during the reconstruction period. This will (to some extent) mimic the portion highlighted below. <\/p>\n<blockquote><p>All proxies available within a given instrumental surface temperature grid box were then averaged and scaled to the same mean and decadal standard deviation as the nearest available 5\u00c2\u00b0 latitude by longitude<br \/>\ninstrumental surface temperature grid box series over the calibration period. [&#8230;] The gridded proxy data were then areally weighted, spatially averaged over the target hemisphere, and <em>scaled to have the same mean and decadal standard deviation<\/em> (we refer to the latter as \u00e2\u20ac\u02dc\u00e2\u20ac\u02dcvariance matching\u00e2\u20ac\u2122\u00e2\u20ac\u2122) as the target hemispheric (NH or SH) mean temperature series.\n<\/p><\/blockquote>\n<p>I say to some extent because for the time being I am not adding the decadal smoothing (which will matter for one of the features in the figures below.)<\/p>\n<p>I quickly slapped together a script to illustrate that this method is also biased under screening. Because it can be shown to have a second type of bias that exists only when there are a small number of proxies, I am showing one graph created starting 50 proxies and another with 484. Blue is used for all proxies; green for the reconstructions when we screen those proxies.  Screening is done against global temperatures. The calibration period is 80 years.<\/p>\n<p>In the graph below, red is the target value the reconstruction should match:<br \/>\n<a href=\"http:\/\/rankexploits.com\/musings\/wp-content\/uploads\/2012\/07\/CPS_biasScreened_50.png\"><img loading=\"lazy\" decoding=\"async\" src=\"http:\/\/rankexploits.com\/musings\/wp-content\/uploads\/2012\/07\/CPS_biasScreened_50-500x500.png\" alt=\"\" title=\"CPS_biasScreened_50\" width=\"500\" height=\"500\" class=\"aligncenter size-medium wp-image-20255\" srcset=\"https:\/\/rankexploits.com\/musings\/wp-content\/uploads\/2012\/07\/CPS_biasScreened_50-500x500.png 500w, https:\/\/rankexploits.com\/musings\/wp-content\/uploads\/2012\/07\/CPS_biasScreened_50-300x300.png 300w, https:\/\/rankexploits.com\/musings\/wp-content\/uploads\/2012\/07\/CPS_biasScreened_50.png 1008w\" sizes=\"auto, (max-width: 500px) 100vw, 500px\" \/><\/a><br \/>\nIf you examine the &#8220;blue&#8221; (unscreened) trace, you can see that it is noisy because there are only 50 proxies. More importantly, the mean value over the reconstruction period is biased high- and the bias is statistically significant.  The green trace illustrates what happens if we bias selecting only those proxies that exhibited a statistically significant correlation with the target temperature during the calibration period.  The bias increases.<\/p>\n<p>In the following graph, I increased the numbers of proxies.<\/p>\n<p><a href=\"http:\/\/rankexploits.com\/musings\/wp-content\/uploads\/2012\/07\/BiasCPS_screened484.png\"><img loading=\"lazy\" decoding=\"async\" src=\"http:\/\/rankexploits.com\/musings\/wp-content\/uploads\/2012\/07\/BiasCPS_screened484-500x500.png\" alt=\"\" title=\"BiasCPS_screened484\" width=\"500\" height=\"500\" class=\"aligncenter size-medium wp-image-20257\" srcset=\"https:\/\/rankexploits.com\/musings\/wp-content\/uploads\/2012\/07\/BiasCPS_screened484-500x500.png 500w, https:\/\/rankexploits.com\/musings\/wp-content\/uploads\/2012\/07\/BiasCPS_screened484-300x300.png 300w, https:\/\/rankexploits.com\/musings\/wp-content\/uploads\/2012\/07\/BiasCPS_screened484.png 1008w\" sizes=\"auto, (max-width: 500px) 100vw, 500px\" \/><\/a><\/p>\n<p>Notice with the bias for the blue (unscreened) result is now non-zero but small.  This is CPS that scales using \u00e2\u20ac\u02dc\u00e2\u20ac\u02dcvariance matching\u00e2\u20ac\u2122\u00e2\u20ac\u2122 is biased when few proxies are used whether or not we screen. This bias vanishes as we increase the number of proxies used.  In contrast, the green (screened) trace remains biased. That is because this bias does not vanish when we increase the number of proxies.<\/p>\n<p>Because this method is used, I&#8217;ll be adding this to my other tests using &#8220;toy&#8221; temperatures calling it &#8220;method IV&#8221;.  That way we can compare the relative bias for different numbers of proxies, calibration periods, screening and other analytical choices a reconstruct-o-logist might make.<\/p>\n","protected":false},"excerpt":{"rendered":"<p>When reading Mann08, I realized I need a &#8220;method IV&#8221; in my previous graphs. In method IV, I will rescale the mean proxy value so the variance of the proxy values match the variance of the temperature during the reconstruction period. This will (to some extent) mimic the portion highlighted below. All proxies available within &hellip; <a href=\"https:\/\/rankexploits.com\/musings\/2012\/cps-with-variance-matching-will-be-method-iv\/\" class=\"more-link\">Continue reading <span class=\"screen-reader-text\">CPS with &#8216;variance matching&#8217;: (Will be method IV)<\/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":[15],"tags":[422,151,420],"class_list":["post-20249","post","type-post","status-publish","format-standard","hentry","category-data-comparisons","tag-cps","tag-mann","tag-screening"],"_links":{"self":[{"href":"https:\/\/rankexploits.com\/musings\/wp-json\/wp\/v2\/posts\/20249","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=20249"}],"version-history":[{"count":0,"href":"https:\/\/rankexploits.com\/musings\/wp-json\/wp\/v2\/posts\/20249\/revisions"}],"wp:attachment":[{"href":"https:\/\/rankexploits.com\/musings\/wp-json\/wp\/v2\/media?parent=20249"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/rankexploits.com\/musings\/wp-json\/wp\/v2\/categories?post=20249"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/rankexploits.com\/musings\/wp-json\/wp\/v2\/tags?post=20249"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}