2016-04-09 15 views
0

53 obsを持つデータフレームのglmmを修正しようとしています。 17変数のすべての変数は標準化されていますが、正規分布に従わず、欠損値もありません。データフレームのstr()は以下のようなものです。family = gaussian(アイデンティティ・リンク)を使用したglmer()の呼び出しエラー

  • 種:19レベル "SPP1"、 "SPP2"、../Wファクタ:5 18 12 15 19 4 6 14 16 5 ...
  • アソシエーション:ファクタW/4レベル "assoA" 、 "assoB"、..:1 1 2 2 2 3 3 4 4 1 ...
  • site:Factor w/2レベル "site1"、 "site2":1 1 1 1 1 1 1 1 1 1 2。 ...
  • obs.no:int型1 1 1 1 1 1 1 1 1 1 ...
  • trait1:NUM 0.652 0.428 0.535 0.389 0.486 ...
  • trait2:NUM 0.135 0.16 0.134 0.142 0.159
  • (13にtrait3をクリッピング)

Iは、サイト指定された形質の関連クラス間の有意性を確認するために、次のコードを実行します。

model1= glmer(trait1 ~ association+site+ (1 | species),data=df6,family=gaussian) 

以下のエラーが発生しました。 glmerで

(trait1〜関連+サイト+(1 |種)、データ= DF6、:家族=(lmerするためのショートカットとしてガウス(アイデンティティリンク))でglmer()を呼び出す が廃止されました。ご連絡ください。 lmer()を直接

この後、私は、ガウス - エルミート直交してパラメータを推定したい。ガウス - エルミート求積を実行するために、このエラーやコードを修正するために、任意の勧告が非常に高く評価される。

答えて

0

あなたが実際に投稿答え:glmerではなくlmerを使用してください:

model1 = lmer(trait1~association+site+(1|species), data=df6) 
+0

ありがとうございました – Gee

+0

しかし、私は助けが不思議ですか?lmerは、REMLや最尤法を使用して、線形混合効果モデル(LMM)をデータに適合させます。 glmerは一般化線形混合効果モデル(GLMM)に適合します。固定効果とランダム効果の両方は、モデル式で指定します。同じことができますか? – Gee

関連する問題