54
edits
No edit summary |
No edit summary |
||
Line 25: | Line 25: | ||
The average of all sample variances of the mean across all unique reflections of a resolution shell is obtained by calculating the sample variance of the mean for every unique reflection i by: | The average of all sample variances of the mean across all unique reflections of a resolution shell is obtained by calculating the sample variance of the mean for every unique reflection i by: | ||
<math>\sigma^2_{\epsilon i} = \frac{1}{n-1} \cdot \left ( \sum^n_{j} x^2_{j} - \frac{\left ( \sum^n_{j}x_{j} \right )^2}{ n} \right ) | <math>\sigma^2_{\epsilon i} = \frac{1}{n-1} \cdot \left ( \sum^n_{j} x^2_{j} - \frac{\left ( \sum^n_{j}x_{j} \right )^2}{ n} \right ) / \frac{n}{2} </math> | ||
With <math>x_{j} </math> , a single observation j of all observations n of one reflection i. <math>\sigma^2_{\epsilon i} </math> is then divided by the factor <math>\frac{n}{2} </math>, because the variance of the sample mean (the merged observations) is the quantity of interest. As we are considering CC12, the variance <math>\sigma^2_{\epsilon i} </math> is divided by <math>\frac{n}{2} </math> and not only by '''n''' as described in [https://en.wikipedia.org/wiki/Sample_mean_and_covariance#Variance_of_the_sample_mean ], because we are calculating the random errors of the merged intensities of a half-dataset. The single variance terms are then summed up for all reflections i in a resolution shell and divided by N, the total number of unique reflections. | With <math>x_{j} </math> , a single observation j of all observations n of one reflection i. <math>\sigma^2_{\epsilon i} </math> is then divided by the factor <math>\frac{n}{2} </math>, because the variance of the sample mean (the merged observations) is the quantity of interest. As we are considering CC12, the variance <math>\sigma^2_{\epsilon i} </math> is divided by <math>\frac{n}{2} </math> and not only by '''n''' as described in [https://en.wikipedia.org/wiki/Sample_mean_and_covariance#Variance_of_the_sample_mean ], because we are calculating the random errors of the merged intensities of a half-dataset. The single variance terms are then summed up for all reflections i in a resolution shell and divided by N, the total number of unique reflections. | ||
<math>\sum^N_{i} \sigma^2_{\epsilon i} | <math>\sum^N_{i} \sigma^2_{\epsilon i} / N </math> | ||
== Example == | == Example == | ||
Line 74: | Line 68: | ||
<math>\sigma^2_{y} </math>, the variance of all the averaged intensities = 190458.6533 | <math>\sigma^2_{y} </math>, the variance of all the averaged intensities = 190458.6533 | ||
As a result of these calculations CC12 = | |||
== Program == | == Program == |
edits