南京邮电大学多元统计分析实验六-对应分析

使用平台：R、Rstudio

题目描述：选取Degree（学历）和Marital（婚姻状况）两个变量进行对应分析，并对结果进行解释。其中，Degree变量是定序尺度：0-less than high school,1-high school,2-junior college,3-bachelor,4-graduate；Marital变量是定类尺度：1-married,2-widowed,3-divorced,4-separated,5-never married

> library(haven)
> library(factoextra)
> library(FactoMineR)
> library(ca)
> library(ggplot2)
> data=read_sav("对应分析.sav")
> high=matrix(0,nrow=5,ncol=5)
> for (i in 1:nrow(data)) {
+   ii <- data[i,1]
+   jj <- data[i,2]
+   ii <- as.integer(ii) -1
+   jj <- as.integer(jj)
+   high[ii,jj] <- high[ii,jj]+1+ 
}
> rownames(high) <- c("less than high school","high school","junior college","bachelor","graduate")
> colnames(high) <- c("married","widowed","divorced","separated","never married")
> # CA函数
> fit_ca=FactoMineR::CA(high,graph=FALSE)
> print(fit_ca)
> # 方差解释
> fit_ca$eig
> # 行惯量
> fit_ca$row$inertia
> # 行主成分轮廓坐标
> FF=fit_ca$row$coord
> # 列惯量
> fit_ca$col$inertia
> # 列主成分轮廓坐标
> GG=fit_ca$col$coord
> # plot.CA(fit_ca, cex = 0.7)
> plot.CA(fit_ca, axes = c(1, 2), cex = 0.7)