問題タブ [lmer]

にmer呼び出されて作成されたオブジェクトがありますlmer()。

でランダム効果を取得できますが、ランダム効果ranef()ごとに対応する数の観測値が必要です-それを行う簡単な方法はありますか？

追加情報：

私は上記で自分自身を完全に明確にしなかったかもしれません。たとえば、患者が病院内にクラスター化され、病院のランダム切片が含まれる単純な2レベルのモデルがある場合、各病院内の患者数ranef() とともに、各病院のランダム効果を抽出したいと思います。現時点では、

これは私に次のようなものを与えます：

私は次のようなものを手に入れたいです：

これを行うために、私は使用しています

そして、cbindこれをからの結果にranef(fullmodel)[[1]]

ただし、これは洗練されておらず、エラーが発生しやすいようです。

Raudenbush/Bryk(1992/2002)、Hierarchical Linear Models、p178 - 179 に基づいて、データを使用して区分的な成長回帰を実行したいと考えています。

これを行うには、マルチレベル モデルを適用する必要があります。

I は、使用する必要がある lme4 パッケージです。

私のデータには、ID ごとの経時的な売上高が含まれています。sales>10 と sales>25 の 2 つのブレークポイントがあります。ブレークポイント前後の売上の伸びを分析したいと考えています。

ID ごとに 2 つの変数があります。 0)

私のマルチレベルモデルの形式は次のとおりです。

Y= b0i + b1i*a1ti + b2i*a2ti + eti

(申し訳ありませんが、画像と数式を投稿することは許可されていません :( )

a1ti と a2ti は、ピースを指定するコード化された変数 (var1 と var2) です。

レベル 1: 個々のレベルごとの反復観察
レベル 2: 被験者間レベル

lmer() コードでレベル 1 とレベル 2 を指定する必要がある場所と、2 つのピースの 2 つの勾配を取得する方法を誰かが説明してくれませんか? インターセプトのみを取得するために、ブラケットの間に仕様を O + ... または 1 + ... として追加する必要がありますか?

私のコードは次のようになります。

またはこれ：

ありがとうございました！

おめでとうアレクサ

G'dayみんな、

'lme4'の'lmer'関数を実行すると、多くのスロットを持つS4オブジェクトが生成されます。これらのスロットの1つ、つまりmodel @ Xと、この「X」スロット出力がどのように生成されるかに興味があります。この出力を、この「X」出力を自動的に生成しない別のモデル関数（glmmPQL）で再現してみたい（FYI「lmer」はクラス「mer」のオブジェクトを生成し、スロット「X」はモデルです）固定効果の行列）。以下のコードは、私が話していることを示しています。

私が理解したいのは、この「X」データがどのように生成されたかです。ターミナルに「（）」を付けずに「lmer」のコードを記述して調べましたが、何も見つかりませんでした。showMethod（'lmer'）も試しましたが、関数'lmer'：と表示されます。

「X」スロットが特に実行していること（またはS4オブジェクトの任意のスロット）のソースコードを取得する方法があるかどうか疑問に思っていますか？または、これを再現する方法を知っている人はいますか？あなたの助けと時間をありがとう。

ここに私の問題があります。混合効果モデルを実行する必要がある R のデータ セットがあります。コードは次のとおりです。

スクリプトを実行しようとすると、次のエラーが返されます。

「follows」パラメーターまでエラーを追跡しました。これは、stress_limit と word_position を使用するだけで正常に動作するためです。「follows」のデータは、n または l、子音、母音の 3 つの文字列のみです。スペースを _ に置き換えてみましたが、成功しませんでした。この場合、「フォロー」の使用を妨げている lmer() 関数の内部動作について何かありますか? どんな助けでも素晴らしいでしょう！

詳細情報: intdiff には数値が含まれ、stress_limit は文字列 (Stressed または Unstressed) であり、単語の位置も文字列 (Word Medial または Word Initial) です。

編集: エラーを再現するデータ サンプルを次に示します。

lme() 関数も試しましたが、次のエラーが返されました。

元の投稿のコードは、私が使用している正確なコードから library(lme4) 呼び出しを除いたものであるため、考えられる情報を省略していません。

私のRバージョンは2.15.2です

2 つのカテゴリ変数によって説明される MAM のグラフィカル表現を作成しようとしています。2つの間の相互作用は重要ではないため、各説明変数の効果を示すために2つの箱ひげ図を作成する必要があることに気付きました。手段などを抽出できないようで、lmerを使用しているためかどうかわからないため、ボックスプロットをグラフ化するにはどうすればよいですか。

これまでの私のコードは次のとおりです....

2 つのグラフをプロットするには、午前または午後には m4 を使用し、天気には m5 を使用する必要があると思います。

私はこれまでのところこれを持っていますが、行き詰まり、「$」を使用すると理解できないエラーが発生するため、役立つ答えがどこにも見つからないようです。

ありがとう

lmer 関数を使用して、3 つの異なる条件 (条件 = 0、1、2) とターゲットの存在 (ターゲット = False または True) の間の反応時間 (RT) に相互作用の影響があるかどうかを調査しようとしています。患者（患者）。

次の式を書きました。

私の問題は、この関数のデフォルトの切片が cond = 0 および target = False であるのに対し、切片を cond= 0 および target=True にしたいことです (cond0*target= の間に大きな違いがあるかどうかを確認するため)。 True および cond1*target=True)。

よろしくお願いします。

ここに私が持っている出力があります

私のデータに基づいて、選択されている切片はcond0:targetTRUEで、出力の他のレベルはcond1:targetFALSEとcond2:targetFALSEです。

lme4 ライブラリを使用して anova モデルを実行しようとしています。実際のデータがないので、MASS から Cars93 データをヤンクしました。lmer(...) を機能させることはできましたが、anova テーブルを印刷できないようです。

何が起こったかは次のとおりです。

どんな助けや提案も大歓迎です。

次の機能を備えた R パッケージはありますか。

(1) 交互作用変数のさまざまな値をシミュレートし、(2) 交互作用の項のさまざまな値に対する Y に対する交互作用の影響を示すグラフをプロットし、(3) lmer( ) lme4 パッケージの機能は?

arm、ez、coefplot2、および fanovaGraph パ​​ッケージを調べましたが、探しているものが見つかりませんでした。

次のコードを使用して、データベースから特定の行を削除しました。

これはダミー変数であり、その特定のダミー変数の値が 1 である行は正常に削除されました。（前後に「dim」機能を使用してデータベースのサイズを確認しましたが、すべて正常に見えました。）

しかし、今回は新しいデータ セット「df2」を使用して回帰モデルを実行したところ、自由度が急激に低下したことがわかりました。これは、削除された行の数をはるかに超えていました!

どうしてこうなったのだろうと思いました。その後、新しいデータセットには、NA のみを持つ行が多数あることに気付きました。確率変数に欠損値がある各行で、R は NA 値の完全な行を作成しました。

上記のコードが行を削除する最良の方法ではないことに気付いた後、次のことを試しました。

私はもう同じ問題を抱えていないので、うまくいったようです。しかし、上記の新しいコードには、私が今のところ気付いていない (類似またはその他の) 問題がある可能性があると思いますか?