r 在PCAmixdata PC地图中绘制单个观察结果

Posted 2021-05-23

## Taken from tutorial at https://cran.r-project.org/web/packages/PCAmixdata/vignettes/PCAmixdata.html

## Import library
library(PCAmixdata)

## Split qualitative and quantitative variables
split <- splitmix(gironde$housing)
X1 <- split$X.quanti 
X2 <- split$X.quali 

## FAMD
res.pcamix <- PCAmix(X.quanti=X1, X.quali=X2,rename.level=TRUE,
                     graph=FALSE)

## Plot individual observations along first two PCs
plot(res.pcamix,
      choice="ind", 
      coloring.ind=X2$houses,
      label=FALSE,
      posleg="bottomright", 
      main="Observations")