next up previous contents index
Next: Producing the diagram Up: Belief comparison diagrams Previous: Belief comparison diagrams

Overview

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, tex2html_wrap_inline36880 and tex2html_wrap_inline36878 , with possibly different specifications about and between tex2html_wrap_inline34520 and tex2html_wrap_inline36880 on the one hand, and about and between tex2html_wrap_inline34520 and tex2html_wrap_inline36878 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 tex2html_wrap_inline38574 and tex2html_wrap_inline38576 are respectively tex2html_wrap_inline38578 and tex2html_wrap_inline38580 . Initially we will be concerned with tex2html_wrap_inline38582 , but we allow the possibility tex2html_wrap_inline38584 . We also restrict ourselves to the case where data is available on tex2html_wrap_inline36880 and tex2html_wrap_inline36878 , so that we can calculate both possible adjusted expectations for tex2html_wrap_inline34520 , and that we also have available the observed value b of tex2html_wrap_inline34520 .

Under model 1, suppose now that the collection tex2html_wrap_inline34520 is adjusted by collection tex2html_wrap_inline36880 . This gives rise to an adjusted version tex2html_wrap_inline38598 which is some vector of linear combinations of the elements in tex2html_wrap_inline34520 . This quantity is observed to be tex2html_wrap_inline38602 , the vector of prediction errors for model 1. We now calculate tex2html_wrap_inline38604 and tex2html_wrap_inline38606 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 tex2html_wrap_inline34520 . Let tex2html_wrap_inline38610 and tex2html_wrap_inline38612 be the smallest and largest canonical quantities (generalised eigenvalues) for this comparison, corresponding to canonical quantities (generalised eigenvectors) tex2html_wrap_inline38614 whose observed magnitudes are tex2html_wrap_inline38616 and tex2html_wrap_inline38618 respectively. The eigenvectors are uncorrelated under both models. The eigenvalues tex2html_wrap_inline38620 , whereas the tex2html_wrap_inline38622 's are normed to have variance unity (or, pathologically, zero) under model 1: tex2html_wrap_inline38624 =1. If model 1 is correct, we should expect to see tex2html_wrap_inline38626 about equal to 1. Observing tex2html_wrap_inline38626 close to tex2html_wrap_inline38630 suggests that the data give more support to model 2. Observing tex2html_wrap_inline38626 distant from both unity and tex2html_wrap_inline38630 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 tex2html_wrap_inline38636 . The observed value tex2html_wrap_inline38638 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 tex2html_wrap_inline38618 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 tex2html_wrap_inline38642 . A large amount of red indicates that model 2 has much more variance for some linear combinations of tex2html_wrap_inline34520 than does model 1.

The top right quadrant is, from the centre outwards, shaded (white,blue) in the proportions tex2html_wrap_inline38646 . The observed value tex2html_wrap_inline38648 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 tex2html_wrap_inline38616 larger than three are plotted on the quadrant boundary to indicate support for neither model. The quadrant is plotted only if tex2html_wrap_inline38652 . A large amount of blue indicates that model 2 has much less variance for some linear combinations of tex2html_wrap_inline34520 than does model 1.

We repeat the adjustment process under model 2, arriving at the vector tex2html_wrap_inline38656 , the observed vector tex2html_wrap_inline38658 , and variance matrices tex2html_wrap_inline38660 and tex2html_wrap_inline38662 .

We now carry out a comparison of these two alternative variance specifications for this second adjusted version of tex2html_wrap_inline34520 . Let tex2html_wrap_inline38666 and tex2html_wrap_inline38668 be the smallest and largest canonical quantities (generalised eigenvalues) for this comparison, corresponding to canonical quantities (generalised eigenvectors) tex2html_wrap_inline38670 whose observed magnitudes are tex2html_wrap_inline38672 and tex2html_wrap_inline38674 respectively. The eigenvectors are uncorrelated under both models. The eigenvalues tex2html_wrap_inline38676 , whereas the tex2html_wrap_inline38678 's are normed to have variance unity (or, pathologically, zero) under model 2: tex2html_wrap_inline38680 =1. If model 2 is correct, we should expect to see tex2html_wrap_inline38682 about equal to 1. Observing tex2html_wrap_inline38682 close to tex2html_wrap_inline38686 suggests that the data give more support to model 1. Observing tex2html_wrap_inline38682 distant from both unity and tex2html_wrap_inline38686 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 tex2html_wrap_inline38692 . The observed value tex2html_wrap_inline38694 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 tex2html_wrap_inline38674 larger than three standard deviations are plotted on the quadrant boundary to indicate support for neither model. This quadrant is plotted only if tex2html_wrap_inline38698 . A large amount of blue indicates that model 1 has much more variance for some linear combinations of tex2html_wrap_inline34520 than does model 2.

The bottom right quadrant is, from the centre outwards, shaded (white,red) in the proportions tex2html_wrap_inline38702 . The observed value tex2html_wrap_inline38704 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 tex2html_wrap_inline38672 larger than three are plotted on the quadrant boundary to indicate support for neither model. The quadrant is plotted only if tex2html_wrap_inline38708 . A large amount of red indicates that model 1 has much less variance for some linear combinations of tex2html_wrap_inline34520 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.


next up previous contents index
Next: Producing the diagram Up: Belief comparison diagrams Previous: Belief comparison diagrams

David Wooff
Wed Oct 21 15:14:31 BST 1998