Note that this command is under construction. The purpose of belief comparison diagrams is to summarise relationships between alternative variance specifications and prediction error for a collection of quantities. We have in mind that there are two alternative information sources, and , with possibly different specifications about and between and on the one hand, and about and between and on the other. We will treat the two scenarios as model 1 and model 2 respectively.
Assume that the two alternative variance specifications given for collections and are respectively and . Initially we will be concerned with , but we allow the possibility . We also restrict ourselves to the case where data is available on and , so that we can calculate both possible adjusted expectations for , and that we also have available the observed value b of .
Under model 1, suppose now that the collection is adjusted by collection . This gives rise to an adjusted version which is some vector of linear combinations of the elements in . This quantity is observed to be , the vector of prediction errors for model 1. We now calculate and to be its variance matrix under model 1 and model 2 respectively.
We now carry out a comparison of these two alternative variance specifications for the adjusted version of . Let and be the smallest and largest canonical quantities (generalised eigenvalues) for this comparison, corresponding to canonical quantities (generalised eigenvectors) whose observed magnitudes are and respectively. The eigenvectors are uncorrelated under both models. The eigenvalues , whereas the 's are normed to have variance unity (or, pathologically, zero) under model 1: =1. If model 1 is correct, we should expect to see about equal to 1. Observing close to suggests that the data give more support to model 2. Observing distant from both unity and suggests that the data support neither model.
The upper semicircle of each node on the comparison diagram summarises these features. The top left quadrant is, from the centre outwards, shaded (white,red) in the proportions . The observed value is plotted according to the current value of the cdscaler control, which treats the radius in the top left quadrant as being equal to (by default) three standard deviations relative to model 2. Values of larger than three standard deviations are plotted on the quadrant boundary to indicate support for neither model. The top left quadrant is plotted only if . A large amount of red indicates that model 2 has much more variance for some linear combinations of than does model 1.
The top right quadrant is, from the centre outwards, shaded (white,blue) in the proportions . The observed value is plotted according to the current value of the cdscaler control, which treats the radius in this quadrant as being equal to (by default) three standard deviations relative to model 1 (i.e. 3, as the standard deviation is unity under model 1). Values of larger than three are plotted on the quadrant boundary to indicate support for neither model. The quadrant is plotted only if . A large amount of blue indicates that model 2 has much less variance for some linear combinations of than does model 1.
We repeat the adjustment process under model 2, arriving at the vector , the observed vector , and variance matrices and .
We now carry out a comparison of these two alternative variance specifications for this second adjusted version of . Let and be the smallest and largest canonical quantities (generalised eigenvalues) for this comparison, corresponding to canonical quantities (generalised eigenvectors) whose observed magnitudes are and respectively. The eigenvectors are uncorrelated under both models. The eigenvalues , whereas the 's are normed to have variance unity (or, pathologically, zero) under model 2: =1. If model 2 is correct, we should expect to see about equal to 1. Observing close to suggests that the data give more support to model 1. Observing distant from both unity and suggests that the data support neither model.
The lower semicircle is now drawn similarly to the upper semicircle. The bottom left quadrant is, from the centre outwards, shaded (white,blue) in the proportions . The observed value is plotted according to the current value of the cdscaler control, which treats the radius in this quadrant as being equal to (by default) three standard deviations relative to model 1. Values of larger than three standard deviations are plotted on the quadrant boundary to indicate support for neither model. This quadrant is plotted only if . A large amount of blue indicates that model 1 has much more variance for some linear combinations of than does model 2.
The bottom right quadrant is, from the centre outwards, shaded (white,red) in the proportions . The observed value is plotted according to the current value of the cdscaler control, which treats the radius in this quadrant as being equal to (by default) three standard deviations relative to model 2. Values of larger than three are plotted on the quadrant boundary to indicate support for neither model. The quadrant is plotted only if . A large amount of red indicates that model 1 has much less variance for some linear combinations of than does model 2.
Lines are drawn to connect the observations in diagonally opposite quadrants to indicate that the quadrants convey similar kinds of information. Observations so large as to imply plotting beyond the boundary of a node are plotted on the boundary and shaded a different colour. The radius of the small circle representing the observation can be changed via the cdradius control.