この先へ

• 二値表現型の場合の分散分解はまだ改善の余地がありそうなので、それについては調べる必要があるだろう(めぼしい資料が見つからないので、もしかしたら、これから誰かがやらないといけないのかもしれない)
• やるとしたら、混合モデルを引き継ぎつつ、線形ではなく、ロジスティック回帰的にするのがよいだろう
• ロジスティック回帰型の混合モデルはもちろんあり、上記記事にも出てきたlme4パッケージのglmer()関数である
• 一般化線形混合モデルについては、こちらに解説記事、そしてこちらにRでのデータハンドリング例