Information Technology Reference
In-Depth Information
0.1
1
0.15
4
-0.02
0.1
2
RAC
0.05
Gaut
AtMr
-0.005
-0.02
-0.05
-0.005
0.01
CrJk
0.05
CmRb
-0.01
0.02
KZN
InAs
DrgR
CmRb
0.02
WCp DrgR
0.2
AtMr
-0.02
InAs
0.01
0
0.1
PubV
Mrd
CmAs
CmAs
Arsn
0
BRs
BNRs
-0.1
-0.01
Rape
-0.02
Mpml
NWst
FrSt
-0.01
0.005
-0.05
NCpe
Limp
0.01
0.05
0.005
-0.04
ECpe
0.02
-0.05
Arsn
AGBH
Figure 7.5 Two-dimensional CA biplot for the 2007/08 crime contingency table,
constructed from the first two columns of U
1
/
2
1
/
2 . Columns are represented
and V
by axes. Calibrations on axes are in terms of R 1 / 2
( X E ) C 1 / 2 , that is, proportional
by a factor of n 1 / 2 to Pearson residuals. It is shown how to obtain these predictions for
the Western Cape province.
We now illustrate the CA biplot resulting from the SVD R 1 / 2
( X E ) C 1 / 2
=
U V but plotting U
2 . The biplot in Figure 7.7
is obtained by specifying the arguments ca.variant = "PearsonResB" and Pear-
sonRes.scaled.markers = FALSE in our function cabipl .
Although the predictions made from Figure 7.7 are exactly those obtained from
Figure 7.5, the biplot in Figure 7.7 differs in a very obvious way from that in Figure 7.5:
the row points are squeezed towards each other, making graphical prediction difficult.
There is an easy way to address this problem: setting argument lambda = TRUE utilizes
the lambda tool described in Section 2.3 by requiring lambda to satisfy
1
/
2
1
/
and V instead of U
and V
p tr ( VV )
q tr ( U
λ =
2 U ) .
(7.52)
Search WWH ::




Custom Search