9.1.3 Visualizing Numerator Relationship Matrices

Two routines are provided for visualization of numerator relationship matrices (NRM), rmuller_pcolor_matrix_pil() and rmuller_spy_matrix_pil().

As an example, we will consider the NRM for the pedigree in Figure 9.1. The matrix is square and symmetric; the diagonal values correspond to $ 1 + f_a$, where $ f_a$ is an animal's coefficient of inbreeding; animals with a diagonal element $ >1$ are inbred.

$\displaystyle \scriptsize
\left[ \begin{array}{llllllllllllllllllll}
1. & 0. & ...
.... & 0. & 0.5 & 0.5 & 0.75 & 0.75 & 0.75 & 1.25
\end{array} \right]
\normalsize
$

Note that the array only contains six distinct values: 0., 0.25, 0.5, 0.75, 1.0, and 1.25. These six values will be used to create the color map used by rmuller_pcolor_matrix_pil().

rmuller_pcolor_matrix_pil() produces pseudocolor plots from NRM. A pseudocolor plot is an array of cells that are colored based on the values the corresponding cells in the NRM. The minimum and maximum values in the NRM are assigned the first and last colors in the colormap; other cells are colored by mapping their values to colormap elements. In the example above, the minimum value is 0.0 and the maximum value is 1.0 (Figure 9.6). The two inbred animals in the population are easily identified as the yellow diagonal elements in the bottom-left corner of the matrix.

Figure 9.6: Pseudocolored NRM from the Boichard et al. (1997) pedigree
rmuller_spy_matrix_pil() is similar to rmuller_pcolor_matrix_pil(), but it is used to visualize the sparsity of a matrix. Cells are either filled, indicating that the value is non-zero, or not filled, indicating that the cell's value is zero. In Figure 9.7 it is easy to see the two separate families in the pedigree.
Figure 9.7: Sparsity of the NRM from the Boichard et al. (1997) pedigree
See About this document... for information on suggesting changes.