特定の生殖行動に費やされた時間 (秒、したがって整数カウント) に対する男性グループ (治療、コントロール) の効果をテストするために、ポアソン誤差分布とログ リンクを使用してglmer モデルをフィッティングしています。同じ女性が治療群と対照群の両方に結合されているため、私の確率因子は女性のアイデンティティによって表されます。
glmer<-glmer(behaviour~group + (1|ID_female), family=poisson, data=behaviour)
モデルは過分散であるため、自由度 12 で 3364.749 (比率: 280.396) であるため、2 番目の観測レベルの確率因子を導入します。
glmer<-glmer(behaviour~group + (1|ID_female) + (1|observation),family=poisson, data=behaviour)
このようにして、過分散を修正します (残差偏差: 11 自由度で 3.08、比率: 0.28) が、特異適合の警告が表示されます。
境界 (特異) フィット: 参照?isSingular
どうすればこれを別の方法で処理できますか? どうもありがとうございました。