lme4

0熱

1答えて

多項式lmerモデルのlsmeansとポストホックテスト（ペアワイズ比較）を実行するにはどうすればよいですか？

今、私は、時間（10時点：1〜10、ポリ因子）、Age_Group1（大人と子供）、タイプ（2つの意味タイプの単語：反復性と所有性のある粒子）とWord（4つのトーン：T1、T2、T3、T4の4つの単語）とランダムな要素：参加者。子供や大人の作品の種類や言葉、とりわけその傾き（線形傾向）と鮮明度（二次傾向）のピッチカーブにしたい。モデルは次のようである： model <- lmer(Pitch

1熱

1答えて

lmerとglsの値が一致しない

lme4にlmerを使用してHLM分析を行っている縦方向のデータセットがあります。この分析の結果を同じデータの結果と比較したいが、nlmeパッケージのglsを使用している。データセット内の各参加者には複数の尺度があり、参加者の中には1つ以上の時点で欠損値があります。だから私は、（1）どのように二つの質問を持って lmerは、これに伴う問題を持っていないようですが、私はglsを使用して同様の解析

0熱

1答えて

グループ構造指定なしの非線形混合モデル

私は非常に単純な非線形混合モデル（Gompertz曲線）を階層構造なしで（反復測定だけで）フィットしようとしています。まず、固定してからランダムな効果で試してみたい。これはデータである data <- structure(list(CumGDD = c(124.66, 124.66, 124.66, 124.66, 124.66, 124.66, 124.66, 124.66, 124.66,

2熱

1答えて

Rでlmeに区分的混合モデルをコーディングするにはどうすればよいですか？

lmerを使用して区分的混合モデルを実行するためにthis exampleに続き、非常にうまく機能します。しかし、私はモデルをlmeに翻訳するのに問題があります。なぜなら、私は異分散を扱う必要があり、lmerにはその能力がありません。問題を再現するコードはhereです。質問に答える必要があると思われる場合は、コード内の実験デザインの詳細を含めました。ここでは、ブレークポイントなしのモデルです：

0熱

1答えて

モデルのコンバージェンスに関する問題：ans.ret [meth、] < - c（ans $ par、ans $ value、ans $ fevals、ans $ gevals、

）のエラーlme4を使用して、収束の問題は、次のリソースを経由して解決しようとすると、新たなエラーを取得：。。 Error in ans.ret[meth, ] <- c(ans$par, ans$value, ans$fevals, ans$gevals, : number of items to replace is not a multiple of replacement length こ

1熱

2答えて

条件付き平均値（条件付きモード）またはランダム効果のプロットがこのように見えるのはなぜですか？

glmerのランダムエフェクトモデルをlme4からRに設定しています。私の理解では、ランダム効果は平均0、分散1.632（上記参照）の正規分布から来るということです。だから私は、条件付き平均の分布（または条件付きモード、getME(modelfit, 'b')を使用して得られる）がベルカーブを多かれ少なかれ従うべきであると予想していました。私は条件付き手段のヒストグラムをプロットするときは、し

0熱

1答えて

GLMERマルチレベルのシミュレーション "zelig style"

rを使ってロジスティック混合効果回帰を実行します。さて、固定効果と私は希望|： glmer（国）、データ= HDP、家族=二項X2 +（1 Y〜X1 + X2 + X1：）回帰は次のように何とかですZeligと一緒に試してみましたが、これは私がシミュレーションをして確率を予測する最も簡単な方法として学んだものですが、新しいバージョンにはマルチレベルモデルはなく、以前のZelig Multilev