library(tidyverse)
library(GDAtools)
library(ggrepel)
library(showtext)
showtext_auto(TRUE)
日本語版「嗜好データ」を読み込む
load("taste_d_J.rda")
.d_J
## # A tibble: 1,253 × 9
## ID Isup TV Film Art Eat Gender Age Income
## <int> <fct> <fct> <fct> <fct> <fct> <fct> <fct> <fct>
## 1 1 Active TV-メロドラマ 映画-アクション 芸術… 外食… 女性 55-64 £20-29
## 2 2 Active TV-メロドラマ 映画-ホラー 芸術… 外食… 女性 45-54 <£9
## 3 3 Active TV-自然 映画-アクション 芸術… 外食… 女性 55-64 <£9
## 4 4 Active TV-メロドラマ 映画-時代劇 芸術… 外食… 女性 65+ £10-19
## 5 5 Active TV-コメディー 映画-ホラー 芸術… 外食… 女性 35-44 £10-19
## 6 6 Active TV-コメディー 映画-ホラー 芸術… 外食… 女性 18-24 <£9
## 7 7 Active TV-ニュース 映画-アクション 芸術… 外食… 女性 25-34 £10-19
## 8 8 Active TV-ニュース 映画-ドキュメンタ… 芸術… 外食… 男性 65+ £10-19
## 9 9 Active TV-メロドラマ 映画-時代劇 芸術… 外食… 女性 65+ <£9
## 10 10 Active TV-ニュース 映画-アクション 芸術… 外食… 女性 65+ £10-19
## # ℹ 1,243 more rows
Isup == “Active” の1215行を切り出す
.d_J %>% filter(Isup=="Active") -> .d0
speMCAを実行
実行後、1,2軸を反転させて、res.speMCAとする
res.speMCA0 <- speMCA(.d0[,3:6])
res.speMCA0 %>% flip.mca(dim = c(1,2)) -> res.speMCA
追加要素（サプリメンタリ：性別、年齢、収入）を接続する
- 個体座標（1−3軸）と追加要素を接続したデータを、ind_coord_supとする★
res.speMCA$ind$coord[,1:3] -> coord_ind
cbind(coord_ind[,1:3],.d0[,7:9]) -> ind_coord_sup
三つのファイルをrdaで保存
.d0 %>% save(file="d0.rda") # MCAの入力にするdf
res.speMCA %>% save(file="res.speMCA.rda") # .d0[,3:6]にMCAをしたresult
ind_coord_sup %>% save(file="ind_coord_sup.rda") # MCAresultの個体coordと追加変数にする.d0の要素をcbindしたdf