我看到了R MASS软件包的源代码以了解其过程。
(GitHub URL https://github.com/cran/MASS/blob/master/R/mca.R)
对于此代码,我有一个问题。 在矩阵SVD之后,源代码将序列设置如下(第36行)
sec <- 1 + (1L:nf)
此序列用于提取主要成分:
rs <- X %*% X.svd$v[, sec]/p
或
structure(list(rs=rs, cs=cs, fs=fs, d=X.svd$d[sec], p=p, call=Call),
class="mca")
由于序列sec
是开始2
,因此上述提取表示MASS程序包中的MCA丢弃了第一个主成分。
此过程是否正确?
我找不到有关此过程的论文或文章。
如果您知道原因,请告诉我吗?