対応分析・多重対応分析

【対応分析】

対応分析(Correspondence Analysis)を行うためには,

まずcaパッケージをダウンロードして,caライブラリを読み込む。

smokeデータに対して対応分析を行う。

行主成分解でプロットする。

mapを指定しないと,対称解でプロットする。

library(ca)

data(smoke)

ca(smoke)

plot(ca(smoke))

plot(ca(smoke), map = "rowprincipal")

エクセル上のデータをコピーして、それを読み込むこともできる。

対応分析のときは、

行変数名も列変数名も読み込む。

x <- read.table("clipboard")

plot(ca(x))

plot(ca(x), map = "rowprincipal")

【多重対応分析】

多重対応分析(Multiple Correspondence Analysis)を行うためには,

まずcaパッケージをダウンロードして,caライブラリを読み込む。

library(ca)

x <- read.table("file.txt", header=T, sep="\t")

summary(x)

mca <- mjca(x,lambda="indicator")

mca <- mjca(x,lambda="Burt")

plot(mca, what=c("none", "all"))

summary(mca)