{"id":2206,"date":"2016-09-27T09:01:14","date_gmt":"2016-09-27T08:01:14","guid":{"rendered":"http:\/\/www.lookingatnothing.com\/?p=2206"},"modified":"2016-09-27T09:01:14","modified_gmt":"2016-09-27T08:01:14","slug":"self-correlation-a-unique-measure-of-saxs-information-content","status":"publish","type":"post","link":"https:\/\/lookingatnothing.com\/index.php\/archives\/2206","title":{"rendered":"Self-correlation: a unique measure of SAXS information content?"},"content":{"rendered":"<p>A funny point came up last week when we were trying to assess the Q-uncertainty inherent in our measurement. We were trying a method inspired by one from Christian Gollwitzer (part 3.3.2, paper <a href=\"http:\/\/dx.doi.org\/10.1039\/C6AY00419A\">here<\/a>). What we did is the following:<\/p>\n<p><!--more--><\/p>\n<p>&nbsp;<\/p>\n<figure id=\"attachment_2213\" aria-describedby=\"caption-attachment-2213\" style=\"width: 300px\" class=\"wp-caption alignright\"><a href=\"http:\/\/www.lookingatnothing.com\/wp-content\/uploads\/2016\/09\/sharpShift.png\"><img loading=\"lazy\" decoding=\"async\" class=\"size-medium wp-image-2213\" src=\"http:\/\/www.lookingatnothing.com\/wp-content\/uploads\/2016\/09\/sharpShift-300x200.png\" alt=\"A pattern with sharp features compared to itself shifted by 0.1 1\/nm.\" width=\"300\" height=\"200\" srcset=\"https:\/\/lookingatnothing.com\/wp-content\/uploads\/2016\/09\/sharpShift-300x200.png 300w, https:\/\/lookingatnothing.com\/wp-content\/uploads\/2016\/09\/sharpShift-768x512.png 768w, https:\/\/lookingatnothing.com\/wp-content\/uploads\/2016\/09\/sharpShift-1024x683.png 1024w, https:\/\/lookingatnothing.com\/wp-content\/uploads\/2016\/09\/sharpShift.png 1800w\" sizes=\"auto, (max-width: 300px) 100vw, 300px\" \/><\/a><figcaption id=\"caption-attachment-2213\" class=\"wp-caption-text\">Figure 1: A pattern with sharp features compared to itself shifted by 0.1 1\/nm.<\/figcaption><\/figure>\n<p>We take a given measured pattern (absolute scale, complete with uncertainties), and shift it by a tiny bit in <img decoding=\"async\" src=\"https:\/\/s0.wp.com\/latex.php?latex=q&#038;bg=ffffff&#038;fg=000&#038;s=0&#038;c=20201002\" alt=\"q\" class=\"latex\" \/>. We then compare the fit of the shifted pattern with its original (including scaling the intensity), and calculate the goodness-of-fit parameter <img decoding=\"async\" src=\"https:\/\/s0.wp.com\/latex.php?latex=%5Cchi%5E2_r&#038;bg=ffffff&#038;fg=000&#038;s=0&#038;c=20201002\" alt=\"&#92;chi^2_r\" class=\"latex\" \/> (Figure 1). In order to calculate this self-similar <img decoding=\"async\" src=\"https:\/\/s0.wp.com\/latex.php?latex=%5Cchi%5E2_r&#038;bg=ffffff&#038;fg=000&#038;s=0&#038;c=20201002\" alt=\"&#92;chi^2_r\" class=\"latex\" \/>, we do need to perform a (linear) interpolation on the intensity and the uncertainties.<\/p>\n<figure id=\"attachment_2214\" aria-describedby=\"caption-attachment-2214\" style=\"width: 300px\" class=\"wp-caption alignright\"><a href=\"http:\/\/www.lookingatnothing.com\/wp-content\/uploads\/2016\/09\/smoothShift.png\"><img loading=\"lazy\" decoding=\"async\" class=\"size-medium wp-image-2214\" src=\"http:\/\/www.lookingatnothing.com\/wp-content\/uploads\/2016\/09\/smoothShift-300x200.png\" alt=\"A smoothly decaying pattern compared to itself shifted by 0.1 1\/nm.\" width=\"300\" height=\"200\" srcset=\"https:\/\/lookingatnothing.com\/wp-content\/uploads\/2016\/09\/smoothShift-300x200.png 300w, https:\/\/lookingatnothing.com\/wp-content\/uploads\/2016\/09\/smoothShift-768x512.png 768w, https:\/\/lookingatnothing.com\/wp-content\/uploads\/2016\/09\/smoothShift-1024x683.png 1024w, https:\/\/lookingatnothing.com\/wp-content\/uploads\/2016\/09\/smoothShift.png 1800w\" sizes=\"auto, (max-width: 300px) 100vw, 300px\" \/><\/a><figcaption id=\"caption-attachment-2214\" class=\"wp-caption-text\">Figure 2: A smoothly decaying pattern compared to itself shifted by 0.1 1\/nm.<\/figcaption><\/figure>\n<p>With this method, we can calculate the dependency of <img decoding=\"async\" src=\"https:\/\/s0.wp.com\/latex.php?latex=%5Cchi%5E2_r&#038;bg=ffffff&#038;fg=000&#038;s=0&#038;c=20201002\" alt=\"&#92;chi^2_r\" class=\"latex\" \/> for a given shift in <img decoding=\"async\" src=\"https:\/\/s0.wp.com\/latex.php?latex=q&#038;bg=ffffff&#038;fg=000&#038;s=0&#038;c=20201002\" alt=\"q\" class=\"latex\" \/>. As you may remember, <img decoding=\"async\" src=\"https:\/\/s0.wp.com\/latex.php?latex=%5Cchi%5E2_r+%5Cleq+1&#038;bg=ffffff&#038;fg=000&#038;s=0&#038;c=20201002\" alt=\"&#92;chi^2_r &#92;leq 1\" class=\"latex\" \/> when the two patterns overlap on average to within the uncertainty of the data. The <img decoding=\"async\" src=\"https:\/\/s0.wp.com\/latex.php?latex=q&#038;bg=ffffff&#038;fg=000&#038;s=0&#038;c=20201002\" alt=\"q\" class=\"latex\" \/> shift at which we reach\u00a0<img decoding=\"async\" src=\"https:\/\/s0.wp.com\/latex.php?latex=%5Cchi%5E2_r+%5Capprox+1&#038;bg=ffffff&#038;fg=000&#038;s=0&#038;c=20201002\" alt=\"&#92;chi^2_r &#92;approx 1\" class=\"latex\" \/> is our maximum allowed shift <img decoding=\"async\" src=\"https:\/\/s0.wp.com\/latex.php?latex=q_%7B%5Cmathrm%7Bmaxshift%7D%7D&#038;bg=ffffff&#038;fg=000&#038;s=0&#038;c=20201002\" alt=\"q_{&#92;mathrm{maxshift}}\" class=\"latex\" \/> (Figure 4).<\/p>\n<figure id=\"attachment_2211\" aria-describedby=\"caption-attachment-2211\" style=\"width: 300px\" class=\"wp-caption alignright\"><a href=\"http:\/\/www.lookingatnothing.com\/wp-content\/uploads\/2016\/09\/flatShift.png\"><img loading=\"lazy\" decoding=\"async\" class=\"size-medium wp-image-2211\" src=\"http:\/\/www.lookingatnothing.com\/wp-content\/uploads\/2016\/09\/flatShift-300x200.png\" alt=\"A mostly flat pattern compared to itself shifted by 0.1 1\/nm.\" width=\"300\" height=\"200\" srcset=\"https:\/\/lookingatnothing.com\/wp-content\/uploads\/2016\/09\/flatShift-300x200.png 300w, https:\/\/lookingatnothing.com\/wp-content\/uploads\/2016\/09\/flatShift-768x512.png 768w, https:\/\/lookingatnothing.com\/wp-content\/uploads\/2016\/09\/flatShift-1024x683.png 1024w, https:\/\/lookingatnothing.com\/wp-content\/uploads\/2016\/09\/flatShift.png 1800w\" sizes=\"auto, (max-width: 300px) 100vw, 300px\" \/><\/a><figcaption id=\"caption-attachment-2211\" class=\"wp-caption-text\">Figure 3: A mostly flat pattern compared to itself shifted by 0.1 1\/nm.<\/figcaption><\/figure>\n<p>This maximum allowed shift will vary depending on the shape in the datasets: the allowed <img decoding=\"async\" src=\"https:\/\/s0.wp.com\/latex.php?latex=q_%7B%5Cmathrm%7Bmaxshift%7D%7D&#038;bg=ffffff&#038;fg=000&#038;s=0&#038;c=20201002\" alt=\"q_{&#92;mathrm{maxshift}}\" class=\"latex\" \/> is small (0.005) for patterns with sharp features (Figure 1), broader (0.04) for smooth patterns (Figure 2), and very large (0.1) for flat profiles (Figure 3). These values should be compared with the smallest q value in our data: <img decoding=\"async\" src=\"https:\/\/s0.wp.com\/latex.php?latex=q_%7Bmin%7D+%3D+0.1&#038;bg=ffffff&#038;fg=000&#038;s=0&#038;c=20201002\" alt=\"q_{min} = 0.1\" class=\"latex\" \/>. So what we seem to have, then, is some sort of measure for the <img decoding=\"async\" src=\"https:\/\/s0.wp.com\/latex.php?latex=q&#038;bg=ffffff&#038;fg=000&#038;s=0&#038;c=20201002\" alt=\"q\" class=\"latex\" \/>-dependent information content in the scattering pattern.<\/p>\n<figure id=\"attachment_2212\" aria-describedby=\"caption-attachment-2212\" style=\"width: 300px\" class=\"wp-caption alignright\"><a href=\"http:\/\/www.lookingatnothing.com\/wp-content\/uploads\/2016\/09\/selfcorrelation.png\"><img loading=\"lazy\" decoding=\"async\" class=\"size-medium wp-image-2212\" src=\"http:\/\/www.lookingatnothing.com\/wp-content\/uploads\/2016\/09\/selfcorrelation-300x200.png\" alt=\"The evolution of the reduced Chi-squared as a function of the shift along q.\" width=\"300\" height=\"200\" srcset=\"https:\/\/lookingatnothing.com\/wp-content\/uploads\/2016\/09\/selfcorrelation-300x200.png 300w, https:\/\/lookingatnothing.com\/wp-content\/uploads\/2016\/09\/selfcorrelation-768x512.png 768w, https:\/\/lookingatnothing.com\/wp-content\/uploads\/2016\/09\/selfcorrelation-1024x683.png 1024w, https:\/\/lookingatnothing.com\/wp-content\/uploads\/2016\/09\/selfcorrelation.png 1800w\" sizes=\"auto, (max-width: 300px) 100vw, 300px\" \/><\/a><figcaption id=\"caption-attachment-2212\" class=\"wp-caption-text\">Figure 4: The evolution of the reduced Chi-squared as a function of the shift along q.<\/figcaption><\/figure>\n<p>This could maybe be used to estimate the uncertainty in our <img decoding=\"async\" src=\"https:\/\/s0.wp.com\/latex.php?latex=q&#038;bg=ffffff&#038;fg=000&#038;s=0&#038;c=20201002\" alt=\"q\" class=\"latex\" \/>-vector, provided we take a pattern with sharp features, and preferably know the size of the object that it represents. Perhaps something silver behenate would do. In any case, we know that the uncertainty determined in this manner lies well below the width the datapoints themselves describe.<\/p>\n<p>For now, at least, we have an interesting measure and some food for thought.<\/p>\n","protected":false},"excerpt":{"rendered":"<div class=\"mh-excerpt\"><p>A funny point came up last week when we were trying to assess the Q-uncertainty inherent in our measurement. We were trying a method inspired <a class=\"mh-excerpt-more\" href=\"https:\/\/lookingatnothing.com\/index.php\/archives\/2206\" title=\"Self-correlation: a unique measure of SAXS information content?\">[&#8230;]<\/a><\/p>\n<\/div>","protected":false},"author":2,"featured_media":2212,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"jetpack_post_was_ever_published":false,"_jetpack_newsletter_access":"","_jetpack_dont_email_post_to_subs":false,"_jetpack_newsletter_tier_id":0,"_jetpack_memberships_contains_paywalled_content":false,"_jetpack_memberships_contains_paid_content":false,"activitypub_content_warning":"","activitypub_content_visibility":"","activitypub_max_image_attachments":4,"activitypub_interaction_policy_quote":"anyone","activitypub_status":"","footnotes":"","jetpack_publicize_message":"","jetpack_publicize_feature_enabled":true,"jetpack_social_post_already_shared":true,"jetpack_social_options":{"image_generator_settings":{"template":"highway","default_image_id":0,"font":"","enabled":false},"version":2}},"categories":[1],"tags":[],"class_list":["post-2206","post","type-post","status-publish","format-standard","has-post-thumbnail","hentry","category-uncategorized"],"jetpack_publicize_connections":[],"jetpack_featured_media_url":"https:\/\/lookingatnothing.com\/wp-content\/uploads\/2016\/09\/selfcorrelation.png","jetpack_shortlink":"https:\/\/wp.me\/p1gZ2v-zA","jetpack_sharing_enabled":true,"_links":{"self":[{"href":"https:\/\/lookingatnothing.com\/index.php\/wp-json\/wp\/v2\/posts\/2206","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/lookingatnothing.com\/index.php\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/lookingatnothing.com\/index.php\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/lookingatnothing.com\/index.php\/wp-json\/wp\/v2\/users\/2"}],"replies":[{"embeddable":true,"href":"https:\/\/lookingatnothing.com\/index.php\/wp-json\/wp\/v2\/comments?post=2206"}],"version-history":[{"count":5,"href":"https:\/\/lookingatnothing.com\/index.php\/wp-json\/wp\/v2\/posts\/2206\/revisions"}],"predecessor-version":[{"id":2215,"href":"https:\/\/lookingatnothing.com\/index.php\/wp-json\/wp\/v2\/posts\/2206\/revisions\/2215"}],"wp:featuredmedia":[{"embeddable":true,"href":"https:\/\/lookingatnothing.com\/index.php\/wp-json\/wp\/v2\/media\/2212"}],"wp:attachment":[{"href":"https:\/\/lookingatnothing.com\/index.php\/wp-json\/wp\/v2\/media?parent=2206"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/lookingatnothing.com\/index.php\/wp-json\/wp\/v2\/categories?post=2206"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/lookingatnothing.com\/index.php\/wp-json\/wp\/v2\/tags?post=2206"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}