speMCA: Performs a 'specific' MCA

In GDAtools: A Toolbox for the Analysis of Categorical Data in Social Sciences, and Especially Geometric Data Analysis

Description

Performs a 'specific' Multiple Correspondence Analysis, i.e. a variant of MCA that allows to treat undesirable categories as passive categories.

Usage

speMCA(data, excl = NULL, ncp = 5, row.w = rep(1, times = nrow(data)))

Arguments

data	data frame with n rows (individuals) and p columns (categorical variables)
excl	numeric vector indicating the indexes of the "junk" categories (default is NULL). See "getindexcat" to identify these indexes.
ncp	number of dimensions kept in the results (default is 5)
row.w	an optional numeric vector of row weights (by default, a vector of 1 for uniform row weights)

Details

Undesirable categories may be of several kinds: infrequent categories (say, <5 percents), heterogeneous categories (e.g. 'others') or uninterpretable categories (e.g. 'not available'). In these cases, 'specific' MCA may be useful to ignore these categories for the determination of distances between individuals (see Le Roux and Rouanet, 2004 and 2010).

Value

Returns an object of class 'speMCA', i.e. a list including:

eig	a list of vectors containing all the eigenvalues, the percentage of variance, the cumulative percentage of variance, the modified rates and the cumulative modified rates
call	a list with informations about input data
ind	a list of matrices containing the results for the individuals (coordinates, contributions)
var	a list of matrices containing all the results for the categories and variables (weights, coordinates, square cosine, categories contributions to axes and cloud, test values (v.test), square correlation ratio (eta2), variable contributions to axes and cloud

Author(s)

Nicolas Robette

References

Le Roux B. and Rouanet H., Multiple Correspondence Analysis, SAGE, Series: Quantitative Applications in the Social Sciences, Volume 163, CA:Thousand Oaks (2010).

Le Roux B. and Rouanet H., Geometric Data Analysis: From Correspondence Analysis to Stuctured Data Analysis, Kluwer Academic Publishers, Dordrecht (June 2004).

See Also

getindexcat, plot.speMCA, varsup, contrib, modif.rate, dimdesc.MCA, MCA, csMCA

Examples

## Performs a specific MCA on 'Music' example data set
## ignoring every 'NA' (i.e. 'not available') categories.
data(Music)
getindexcat(Music[,1:5])
mca <- speMCA(Music[,1:5],excl=c(3,6,9,12,15))
str(mca)

Example output

Loading required package: FactoMineR
Loading required package: nleqslv
Loading required package: nnet
 [1] "FrenchPop.No"  "FrenchPop.Yes" "FrenchPop.NA"  "Rap.No"       
 [5] "Rap.Yes"       "Rap.NA"        "Rock.No"       "Rock.Yes"     
 [9] "Rock.NA"       "Jazz.No"       "Jazz.Yes"      "Jazz.NA"      
[13] "Classical.No"  "Classical.Yes" "Classical.NA" 
List of 5
 $ eig :List of 5
  ..$ eigen    : num [1:10] 0.287 0.235 0.188 0.165 0.124 ...
  ..$ rate     : num [1:10] 28.4 23.3 18.6 16.4 12.2 ...
  ..$ cum.rate : num [1:10] 28.4 51.7 70.3 86.7 98.9 ...
  ..$ mrate    : num [1:2] 86 14
  ..$ cum.mrate: num [1:2] 86 100
 $ call:List of 7
  ..$ X        :'data.frame':	500 obs. of  5 variables:
  .. ..$ FrenchPop: Factor w/ 3 levels "No","Yes","NA": 1 2 2 2 2 1 2 2 2 2 ...
  .. ..$ Rap      : Factor w/ 3 levels "No","Yes","NA": 1 1 1 1 1 1 1 1 1 1 ...
  .. ..$ Rock     : Factor w/ 3 levels "No","Yes","NA": 1 3 1 1 2 1 2 1 2 1 ...
  .. ..$ Jazz     : Factor w/ 3 levels "No","Yes","NA": 2 1 1 1 1 1 1 1 2 1 ...
  .. ..$ Classical: Factor w/ 3 levels "No","Yes","NA": 2 1 1 1 1 1 1 1 2 1 ...
  ..$ marge.col: Named num [1:10] 0.0692 0.1304 0.178 0.0188 0.1508 ...
  .. ..- attr(*, "names")= chr [1:10] "FrenchPop.No" "FrenchPop.Yes" "Rap.No" "Rap.Yes" ...
  ..$ marge.row: Named num [1:500] 4e-04 4e-04 4e-04 4e-04 4e-04 4e-04 4e-04 4e-04 4e-04 4e-04 ...
  .. ..- attr(*, "names")= chr [1:500] "1" "2" "3" "4" ...
  ..$ ncp      : num 5
  ..$ quali    : int [1:5] 1 2 3 4 5
  ..$ excl     : num [1:5] 3 6 9 12 15
  ..$ row.w    : num [1:500] 1 1 1 1 1 1 1 1 1 1 ...
 $ ind :List of 2
  ..$ coord  : num [1:500, 1:5] 0.847 -0.25 -0.359 -0.359 0.125 ...
  .. ..- attr(*, "dimnames")=List of 2
  .. .. ..$ : chr [1:500] "4791" "582" "4198" "3704" ...
  .. .. ..$ : chr [1:5] "dim.1" "dim.2" "dim.3" "dim.4" ...
  ..$ contrib: num [1:500, 1:5] 0.499 0.0435 0.0899 0.0899 0.0108 ...
  .. ..- attr(*, "dimnames")=List of 2
  .. .. ..$ : chr [1:500] "4791" "582" "4198" "3704" ...
  .. .. ..$ : chr [1:5] "dim.1" "dim.2" "dim.3" "dim.4" ...
 $ var :List of 9
  ..$ weight    : Named num [1:10] 173 326 445 47 377 113 417 80 370 125
  .. ..- attr(*, "names")= chr [1:10] "FrenchPop.No" "FrenchPop.Yes" "Rap.No" "Rap.Yes" ...
  ..$ coord     : num [1:10, 1:5] -0.271 0.138 -0.088 0.812 -0.294 ...
  .. ..- attr(*, "dimnames")=List of 2
  .. .. ..$ : chr [1:10] "FrenchPop.No" "FrenchPop.Yes" "Rap.No" "Rap.Yes" ...
  .. .. ..$ : chr [1:5] "dim.1" "dim.2" "dim.3" "dim.4" ...
  ..$ contrib   : num [1:10, 1:5] 1.762 0.866 0.479 4.313 4.526 ...
  .. ..- attr(*, "dimnames")=List of 2
  .. .. ..$ : chr [1:10] "FrenchPop.No" "FrenchPop.Yes" "Rap.No" "Rap.Yes" ...
  .. .. ..$ : chr [1:5] "dim.1" "dim.2" "dim.3" "dim.4" ...
  ..$ ctr.cloud :'data.frame':	10 obs. of  1 variable:
  .. ..$ ctr.cloud: num [1:10] 13.08 6.96 2.2 18.12 4.92 ...
  ..$ cos2      : num [1:10, 1:5] 0.0387 0.0357 0.0626 0.0684 0.2644 ...
  .. ..- attr(*, "dimnames")=List of 2
  .. .. ..$ : chr [1:10] "FrenchPop.No" "FrenchPop.Yes" "Rap.No" "Rap.Yes" ...
  .. .. ..$ : chr [1:5] "dim.1" "dim.2" "dim.3" "dim.4" ...
  ..$ v.test    : num [1:10, 1:5] -4.4 4.22 -5.59 5.84 -11.49 ...
  .. ..- attr(*, "dimnames")=List of 2
  .. .. ..$ : chr [1:10] "FrenchPop.No" "FrenchPop.Yes" "Rap.No" "Rap.Yes" ...
  .. .. ..$ : chr [1:5] "dim.1" "dim.2" "dim.3" "dim.4" ...
  ..$ eta2      : num [1:5, 1:5] 0.044 0.0691 0.2944 0.6387 0.3988 ...
  .. ..- attr(*, "dimnames")=List of 2
  .. .. ..$ : chr [1:5] "FrenchPop" "Rap" "Rock" "Jazz" ...
  .. .. ..$ : chr [1:5] "dim.1" "dim.2" "dim.3" "dim.4" ...
  ..$ v.contrib :'data.frame':	5 obs. of  5 variables:
  .. ..$ dim.1: num [1:5] 2.63 4.79 20.39 44.44 27.75
  .. ..$ dim.2: num [1:5] 21.6378 40.8602 16.7148 0.0185 20.7686
  .. ..$ dim.3: num [1:5] 73.31 5.72 6.88 5.29 8.79
  .. ..$ dim.4: num [1:5] 1.206 46.1 47.671 0.198 4.825
  .. ..$ dim.5: num [1:5] 1.21 2.53 8.36 50.03 37.87
  ..$ vctr.cloud:'data.frame':	5 obs. of  1 variable:
  .. ..$ vctr.cloud: num [1:5] 20 20.3 20.4 20.1 20.2
 $ svd :List of 3
  ..$ vs: num [1:10] 0.536 0.485 0.434 0.407 0.352 ...
  ..$ U : num [1:500, 1:10] 0.0706 -0.0208 -0.03 -0.03 0.0104 ...
  ..$ V : num [1:10, 1:10] -0.1327 0.093 -0.0692 0.2077 -0.2127 ...
 - attr(*, "class")= chr [1:2] "speMCA" "list"

GDAtools documentation built on May 29, 2017, 9:48 p.m.