この分析には青木先生のサイト http://aoki2.si.gunma-u.ac.jp/R/ のプログラムと，選挙の総合サイトhttp://go2senkyo.com/blog/archives/6265 のデータにお世話になっています。
まず青木先生のサイトから，双対尺度法のコードとダミー変数にするコードをとってくる。
source("http://aoki2.si.gunma-u.ac.jp/R/src/dual.R", encoding = "euc-jp")
source("http://aoki2.si.gunma-u.ac.jp/R/src/summary.dual.R", encoding = "euc-jp")
source("http://aoki2.si.gunma-u.ac.jp/R/src/plot.dual.R", encoding = "euc-jp")
source("http://aoki2.si.gunma-u.ac.jp/R/src/make.dummy.R", encoding = "euc-jp")
次にデータを作る。
manifest <- data.frame(nuclear = c(1, 3, 3, 3, 2, 3, 3, 3, 3, 3, 1), salesTax = c(1,
1, 1, 1, 1, 2, 2, 2, 2, 2, 2), TPP = c(1, 1, 2, 2, 2, 3, 1, 3, 3, 3, 1))
manifest.dummy <- make.dummy(manifest)
colnames(manifest.dummy) <- c("nuclear.○", "nuclear.△", "nuclear.×", "salesTax.○",
"salesTax.×", "TPP.○", "TPP.△", "TPP.×")
rownames(manifest.dummy) <- c("民主党", "国民新党", "自民党", "公明党",
"日本維新の会", "共産党", "みんなの党", "社民党", "新党大地",
"新党日本", "幸福実現党")
双対尺度法の実行とプロット。
ans <- dual(manifest.dummy)
plot(ans)