53 個の obs を持つデータフレームの glmm を修正しようとしています。17 変数の すべての変数は標準化されていますが、正規分布に従っておらず、欠損値はありません。データフレームの str() は以下のようなものです。
- 種 : 19 レベルの因子 "spp1","spp2",..: 5 18 12 15 19 4 6 14 16 5 ...
- association : 4 水準の因数 "assoA","assoB",..: 1 1 2 2 2 3 3 4 4 1 ...
- site : 2 水準の因数 "site1","site2": 1 1 1 1 1 1 1 1 1 2 ...
- obs.no : int 1 1 1 1 1 1 1 1 1 ...
- 特性 1: 数値 0.652 0.428 0.535 0.389 0.486 ...
- 特性 2 : 数値 0.135 0.16 0.134 0.142 0.159
- (trait3 から 13 を切り取った)
次のコードを実行して、指定された特性の関連クラスとサイト間の重要性を確認しました。
model1= glmer(trait1 ~ association+site+ (1 | species),data=df6,family=gaussian)
以下のエラーを受け取りました。
glmer(trait1 ~ association+site+ (1 | species),data=df6, : lmer() へのショートカットとして family=gaussian (identity link) を指定して glmer() を呼び出すことは非推奨です; lmer() を直接呼び出してください
この後、ガウス・エルミート求積法でパラメータを推定したいと思います。このエラーを修正するための推奨事項と、ガウス・エルミート求積法を実行するためのコードは非常に高く評価されています。