で,
発表分野のうちわけとやらは?
> summary(d$分野1)
01_群落 02_植物個体群 03_植物生理生態
9 4 9
04_植物繁殖 05_植物生活史 08_菌類
3 2 1
10_景観 11_遷移・更新 12_フェノロジー
7 6 1
13_動物と植物の相互関係 14_進化 15_生物多様性
8 12 3
16_数理 17_動物群集 18_動物繁殖
10 10 4
19_動物個体群 20_動物生活史 21_行動
7 3 14
22_動物の社会性 23_分子 24_古生態
4 2 3
25_保全 26_生態系管理 27_外来種
16 11 7
29_物質生産 30_物質循環 31_生態学教育・普及
1 5 2
3 未満は「つぶす」ので,
そのように処理するとこうなる.
> # 発表数 3 未満のグループを廃止
> g <- d$分野1
> table.g <- table(g)
> removed.g1 <- levels.g[table.g < 3]
> people.removed.g1 <- g %in% removed.g1
>
> # 第一希望が通らなかった人たちを第二希望に
> g[people.removed.g1] <- d[people.removed.g1, "分野2"]
> table.g <- table(g) # 再計算
> removed.g <- levels.g[table.g == 0]
> g <- factor(g)
> table.g <- table(g) # 再計算
> table.g
g
01_群落 02_植物個体群 03_植物生理生態
9 5 9
04_植物繁殖 10_景観 11_遷移・更新
4 7 6
13_動物と植物の相互関係 14_進化 15_生物多様性
8 12 3
16_数理 17_動物群集 18_動物繁殖
11 10 4
19_動物個体群 20_動物生活史 21_行動
9 3 14
22_動物の社会性 24_古生態 25_保全
4 3 18
26_生態系管理 27_外来種 30_物質循環
11 7 6