Difference between revisions of "Precision Models"
(→Inference from Extreme Spreads: Deleted because we now have much better analysis in Range Statistics) |
(→Statistical Analysis of Dispersion) |
||
(14 intermediate revisions by the same user not shown) | |||
Line 3: | Line 3: | ||
= Models of Dispersion = | = Models of Dispersion = | ||
− | We | + | We present five models for measuring and analyzing precision: |
# [[Closed Form Precision]] | # [[Closed Form Precision]] | ||
# [[Circular Error Probable]] | # [[Circular Error Probable]] | ||
# [[Elliptic Error Probable]] | # [[Elliptic Error Probable]] | ||
# [[Range Statistics]] | # [[Range Statistics]] | ||
+ | # [[Order Statistics]] | ||
Before selecting one consider the following background: | Before selecting one consider the following background: | ||
== General Bivariate Normal == | == General Bivariate Normal == | ||
− | [http://en.wikipedia.org/wiki/Normal_distribution The normal, a.k.a. Gaussian, distribution] is the broadly accepted model of a random variable like the dispersion of a physical gunshot from its center point. The normal distribution is parameterized by its mean and standard deviation, or <math>(\mu, \sigma)</math>. As explained in ''[[What is Precision?]]'' we are only interested in the dispersion component, since the center point of impact is controlled by sighting in the gun (i.e., adjusting its aiming device). Therefore we will assume that a gunner can dial <math>\mu | + | [http://en.wikipedia.org/wiki/Normal_distribution The normal, a.k.a. Gaussian, distribution] is the broadly accepted model of a random variable like the dispersion of a physical gunshot from its center point. The normal distribution is parameterized by its mean and standard deviation, or <math>(\mu, \sigma)</math>. As explained in ''[[What is Precision?]]'' we are only interested in the dispersion component, since the center point of impact is controlled by [[FAQ#How_many_shots_do_I_need_to_sight_in.3F|sighting in the gun]] (i.e., adjusting its aiming device). Therefore we will assume that a gunner can dial <math>\mu \approx 0</math> and leave that parameter out of the question in what follows. |
Since we are interested in shot dispersion on a two-dimensional target we will look at a [http://en.wikipedia.org/wiki/Bivariate_normal_distribution bivariate normal distribution], which has separate parameters for the standard deviation in each dimension, <math>\sigma_x, \sigma_y</math>, as well as a correlation parameter ''ρ''. | Since we are interested in shot dispersion on a two-dimensional target we will look at a [http://en.wikipedia.org/wiki/Bivariate_normal_distribution bivariate normal distribution], which has separate parameters for the standard deviation in each dimension, <math>\sigma_x, \sigma_y</math>, as well as a correlation parameter ''ρ''. | ||
== Uncorrelated Bivariate Normal == | == Uncorrelated Bivariate Normal == | ||
− | We don't have any | + | We don't have any evidence that there is, or should be, correlation between the horizontal and vertical dispersion of gunshots. Therefore, throughout our analysis we will assume ''ρ'' = 0. |
− | We do know that targets can often exhibit vertical or horizontal stringing, and therefore <math>\sigma_x \neq \sigma_y</math>. To the extent these parameters are not equal they produce | + | We do know that targets can often exhibit vertical or horizontal stringing, and therefore <math>\sigma_x \neq \sigma_y</math>. To the extent these parameters are not equal they produce [[Elliptic Error Probable|elliptic]] instead of circular shot groups. |
However, we know some of the significant sources of stringing and can potentially factor them out: | However, we know some of the significant sources of stringing and can potentially factor them out: | ||
− | # The primary source of x-specific variance is crosswind. If we measure the wind while shooting we can bound and remove a “wind variance” term from that axis. E.g., | + | # The primary source of x-specific variance is crosswind. If we measure the wind while shooting we can bound and remove a “wind variance” term from that axis. E.g., ''Suppose the orthogonal component of wind is ranging at random from 0-10mph during the shooting. Given lag-time ''t'' this will expand the no-wind horizontal dispersion at the target by <math>\sigma_w</math>.''<ref>Wind deflection is a function of the ballistic curve and distance, but can be expressed as a simple product of the cross-wind velocity and lag time. For more information on the "lag rule" see Bryan Litz, ''Applied Ballistics for Long Range Shooting, 2<sup>nd</sup> Edition'' (2011) A4; or Robert McCoy, ''Modern Exterior Ballistics, 2<sup>nd</sup> Edition'' (2012) 7.27.</ref> Since variances are additive we could adjust <math>\sigma_x</math> via the equation <math>{\sigma'}_x^2 = \sigma_x^2 - \sigma_w^2</math>. |
# The primary source of y-specific variance is muzzle velocity, which we can actually measure with a chronograph (or assert) and then remove from that axis. E.g., "If standard deviation of muzzle velocity is <math>\sigma_{mv}</math> then, given the bullet's ballistic model for the given target distance, the vertical spread attributable to that is some <math>\sigma_v</math>. Here too we can remove this known source of dispersion from our samples via the equation <math>{\sigma'}_y^2 = \sigma_y^2 - \sigma_v^2</math>. This adjustment is shown in several of the examples: | # The primary source of y-specific variance is muzzle velocity, which we can actually measure with a chronograph (or assert) and then remove from that axis. E.g., "If standard deviation of muzzle velocity is <math>\sigma_{mv}</math> then, given the bullet's ballistic model for the given target distance, the vertical spread attributable to that is some <math>\sigma_v</math>. Here too we can remove this known source of dispersion from our samples via the equation <math>{\sigma'}_y^2 = \sigma_y^2 - \sigma_v^2</math>. This adjustment is shown in several of the examples: | ||
#* [[22LR CCI 40gr HV 40-shot 100-yard Example]] | #* [[22LR CCI 40gr HV 40-shot 100-yard Example]] | ||
Line 31: | Line 32: | ||
In view of the preceding: | In view of the preceding: | ||
− | # The [[Closed Form Precision]] model requires that we assume the shot group is, or can be normalized to be, a fairly symmetric bivariate Gaussian process | + | # The [[Closed Form Precision]] model requires that we assume the shot group is, or can be normalized to be, a fairly symmetric bivariate Gaussian process. This assumption is the most amenable to statistical analysis. |
+ | # [[Order Statistics]] are slightly less efficient and amenable to abstract analysis, but are both more robust and easier to apply "in the field." | ||
# [[Circular Error Probable]] disregards any ellipticity in the actual shot process in order to characterize precision using a single parameter. Since most of precision estimation is for the purposes of comparing loads, rifles, and shooters, we need a single number and we don't care if the dispersion is elliptic: tighter is always better. | # [[Circular Error Probable]] disregards any ellipticity in the actual shot process in order to characterize precision using a single parameter. Since most of precision estimation is for the purposes of comparing loads, rifles, and shooters, we need a single number and we don't care if the dispersion is elliptic: tighter is always better. | ||
− | # [[Elliptic Error Probable]] allows for a full characterization of the General Bivariate Normal model. For some applications – e.g., computing hit probabilities on non-circular targets – we want to preserve statistically significant ellipticity | + | # [[Elliptic Error Probable]] allows for a full characterization of the General Bivariate Normal model. For some applications – e.g., computing hit probabilities on non-circular targets – we want to preserve statistically significant ellipticity. |
− | # Extreme Spread and the other [[Range Statistics]], which increase with group | + | # Extreme Spread and the other [[Range Statistics]], which increase with number of shots per group ''n'', do not have any useful functional forms. The characteristics of these measures have to be derived from Monte Carlo simulation. They are the least efficient statistics but are also the most commonly used because they are so easy to measure in the field and so familiar to shooters. |
− | = | + | One practical question that many shooters raise is what to do with outliers, known in the sport as "fliers." We address [[Fliers|fliers here]]. |
+ | |||
+ | = Tools = | ||
See [[Measuring Tools]] for convenient ways of measuring and analyzing precision. | See [[Measuring Tools]] for convenient ways of measuring and analyzing precision. | ||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
= References = | = References = | ||
<references /> | <references /> |
Latest revision as of 14:13, 21 November 2023
Previous: Describing Precision
Contents
Models of Dispersion
We present five models for measuring and analyzing precision:
- Closed Form Precision
- Circular Error Probable
- Elliptic Error Probable
- Range Statistics
- Order Statistics
Before selecting one consider the following background:
General Bivariate Normal
The normal, a.k.a. Gaussian, distribution is the broadly accepted model of a random variable like the dispersion of a physical gunshot from its center point. The normal distribution is parameterized by its mean and standard deviation, or \((\mu, \sigma)\). As explained in What is Precision? we are only interested in the dispersion component, since the center point of impact is controlled by sighting in the gun (i.e., adjusting its aiming device). Therefore we will assume that a gunner can dial \(\mu \approx 0\) and leave that parameter out of the question in what follows.
Since we are interested in shot dispersion on a two-dimensional target we will look at a bivariate normal distribution, which has separate parameters for the standard deviation in each dimension, \(\sigma_x, \sigma_y\), as well as a correlation parameter ρ.
We don't have any evidence that there is, or should be, correlation between the horizontal and vertical dispersion of gunshots. Therefore, throughout our analysis we will assume ρ = 0.
We do know that targets can often exhibit vertical or horizontal stringing, and therefore \(\sigma_x \neq \sigma_y\). To the extent these parameters are not equal they produce elliptic instead of circular shot groups.
However, we know some of the significant sources of stringing and can potentially factor them out:
- The primary source of x-specific variance is crosswind. If we measure the wind while shooting we can bound and remove a “wind variance” term from that axis. E.g., Suppose the orthogonal component of wind is ranging at random from 0-10mph during the shooting. Given lag-time t this will expand the no-wind horizontal dispersion at the target by \(\sigma_w\).[1] Since variances are additive we could adjust \(\sigma_x\) via the equation \({\sigma'}_x^2 = \sigma_x^2 - \sigma_w^2\).
- The primary source of y-specific variance is muzzle velocity, which we can actually measure with a chronograph (or assert) and then remove from that axis. E.g., "If standard deviation of muzzle velocity is \(\sigma_{mv}\) then, given the bullet's ballistic model for the given target distance, the vertical spread attributable to that is some \(\sigma_v\). Here too we can remove this known source of dispersion from our samples via the equation \({\sigma'}_y^2 = \sigma_y^2 - \sigma_v^2\). This adjustment is shown in several of the examples:
Statistical Analysis of Dispersion
In view of the preceding:
- The Closed Form Precision model requires that we assume the shot group is, or can be normalized to be, a fairly symmetric bivariate Gaussian process. This assumption is the most amenable to statistical analysis.
- Order Statistics are slightly less efficient and amenable to abstract analysis, but are both more robust and easier to apply "in the field."
- Circular Error Probable disregards any ellipticity in the actual shot process in order to characterize precision using a single parameter. Since most of precision estimation is for the purposes of comparing loads, rifles, and shooters, we need a single number and we don't care if the dispersion is elliptic: tighter is always better.
- Elliptic Error Probable allows for a full characterization of the General Bivariate Normal model. For some applications – e.g., computing hit probabilities on non-circular targets – we want to preserve statistically significant ellipticity.
- Extreme Spread and the other Range Statistics, which increase with number of shots per group n, do not have any useful functional forms. The characteristics of these measures have to be derived from Monte Carlo simulation. They are the least efficient statistics but are also the most commonly used because they are so easy to measure in the field and so familiar to shooters.
One practical question that many shooters raise is what to do with outliers, known in the sport as "fliers." We address fliers here.
Tools
See Measuring Tools for convenient ways of measuring and analyzing precision.
References
- ↑ Wind deflection is a function of the ballistic curve and distance, but can be expressed as a simple product of the cross-wind velocity and lag time. For more information on the "lag rule" see Bryan Litz, Applied Ballistics for Long Range Shooting, 2nd Edition (2011) A4; or Robert McCoy, Modern Exterior Ballistics, 2nd Edition (2012) 7.27.