対応分析・多重対応分析
【対応分析】
対応分析(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)