What are the equivalents of MCA variables coordinates and supplementary variables coordinates in mjca?

536 views Asked by At

I would like to use mjca (package 'ca') on my data in order to estimate explained variation more realistically for the dimensions. The problem is that I would like to extract the coordinates of the active and supplementary variables in order to edit them in a data frame. However, the names of the variables and the dimensions are not given in the output of mjca. In MCA (package 'FactoMineR') the output is given as follows:

> mca$var$coord
                       Dim 1        Dim 2        Dim 3       Dim 4       Dim 5
a                0.620468268  0.011534137 -0.542655702  0.47922448  0.15548571
cl               1.231043177  4.591555841 -0.323929172  0.19597918 -0.41446395
np              -0.347646238  0.003735466 -0.006099464  0.02238883  0.16510343
num              0.417635652 -0.351884061 -0.760499677  0.60590774 -0.35647256
pr               0.945109906 -0.227098798  3.411969743  2.70823750 -0.64981046
vp               0.809895398 -0.303805822  0.048900811 -0.50023568 -0.53191069
EMB_no           0.396034450 -0.046768029 -0.058069978  0.05448188  0.06326411
EM_yes          -1.009887848  0.119258474  0.148078445 -0.13892880 -0.16132349
ca              -0.345163332 -0.088791765 -0.222907122  0.16679404 -0.12407031
to               0.375618920  0.096626332  0.242575397 -0.18151117  0.13501769
ART_no          -0.006456155  0.021963298  0.049258256 -0.05919682 -0.07539649
ART_yes          0.044475732 -0.151302718 -0.339334655  0.40780032  0.51939806

> mca$quali.sup$coord
                   Dim 1       Dim 2      Dim 3       Dim 4       Dim 5
ipva            -0.1508708  0.04768873 -0.0233159  0.08795449  0.01645747
isv              0.6731160 -0.21276510  0.1040248 -0.39241234 -0.07342562

Is there a way to extract and paste the names of the variables to the coordinates in mjca? In mjca the output is not easily interpretable:

> mjca$colcoord
             [,1]         [,2]         [,3]        [,4]        [,5]        [,6]        [,7]
 [1,] -1.14877950  -0.03284730   1.85933139  1.71512222  0.59783898  2.33329527 -1.96334559
 [2,] -2.27924173 -13.07598521   1.10989653  0.70140040 -1.59360435 -0.53239520  0.59206669
 [3,]  0.64365721  -0.01063798   0.02089893  0.08012859  0.63481889 -0.18685503  0.16878106
 [4,] -0.77324063   1.00210711   2.60574231  2.16851574 -1.37062878  1.00967151 -7.97223158
 [5,] -1.74984434   0.64673950 -11.69062156  9.69265663 -2.49850629 -0.05490287 -0.25510472
 [6,] -1.49949849   0.86518831  -0.16755157 -1.79032033 -2.04518437 -0.20949668  0.65583755
 [7,] -0.73324662   0.13318755   0.19896839  0.19498813  0.24324906  0.05337135  0.15195814
 [8,]  1.86977887  -0.33962824  -0.50736940 -0.49721974 -0.62028511 -0.13609695 -0.38749326
 [9,]  0.63906018   0.25286414   0.76375906  0.59694816 -0.47704750 -0.12352216 -0.44654719
[10,] -0.69544784  -0.27517568  -0.83114957 -0.64962006  0.51913993  0.13442117  0.48594841
[11,]  0.01195339  -0.06254781  -0.16877630 -0.21186268 -0.28989778 -0.27854165 -0.39280823
[12,] -0.08234556   0.43088491   1.16268119  1.45949845  1.99707361  1.91884250  2.70601223
[13,] -0.76304235   0.48124458   1.23544015  1.02426084  0.57083939  1.79253772 -0.83528819
[14,] -1.64355033  -5.87863606   0.76112230  0.59773026 -0.22298799  0.40707187 -0.23302386
[15,] -0.40920789   0.55803109   1.06737902  0.94054386  0.10598994 -0.60665661 -1.11993976

I think I should first extract the coordinates like this:

coord.mjca<-as.data.frame(mjca$colcoord)
row.names(coord.mjca)<-mjca$levelnames
colnames(coord.mjca)<-c("Dim 1", "Dim 2", "Dim 3", "Dim 4", "Dim 5", "Dim 6", "Dim 7")

Do you think I should do it like this?

Thank you for your help!

1

There are 1 answers

0
Bora On BEST ANSWER

In ca package, you can use cacoord for this.

For example

cacoord(mca, type='rowprincipal', rows=T)
cacoord(mca, type='symmetric', cols=T)

I hope this helps.