問題タブ [lme4]

私はlmerをやっていて、RにNAセルを無視させたいと思っています。私はこのコマンドを入れました：

ただし、の要約ではlmer、元のExcelファイルの行数と同じ数の行（obsの数）が得られるためNA、式で行が除外されていないようです。これに関する推奨事項はありますか？

多変量線形混合モデルを実行しようとしていますが、計算時間を短縮するためにリモート ワークステーションを使用する必要があります。

パーソナル コンピューター (R バージョン 2.15.1、lme4 バージョン 0.999999-0、64 ビット Unix) で lme4 から glmer() を実行すると、モデルが適切に実行されます。

削除ワークステーション (R バージョン 3.0.2、lme4 バージョン 1.0-6、64 ビット Linux) に切り替えると、モデルが実行されず、エラーと警告メッセージが表示されます。

興味深いことに、従属変数のいずれかで lmer() を実行すると、モデルはエラーや警告なしで実行されます。

サンプルコードは次のとおりです。

ここに私のパソコンからの出力があります：

以下は、新しいバージョンの lme4 を使用したリモート ワークステーションからの出力です。

モデルが正しく形成されていないのか、それとも lme4 の 1 つのバージョンに問題があるのか​​疑問に思っています。ご意見やご提案は大歓迎です。ありがとう！！

ポアソンモデルを使用した glmer と difflsmeans の間の結果の不一致を理解するのに苦労しています。どちらの関数も lmerTest パッケージからのものです。基本的に、glmer は 2 つの係数が p < 0.05 で有意であると教えてくれますが、difflsmeans を使用すると異なる結果が得られます。私はポアソン リンクとモデル カウント データへのオフセットを使用しており、実験のバッチには処理の固定効果と変量効果を含めています。

ガウス リンクを使用する他の分析では、lmer と difflsmeans の両方で同じ結果が得られます。

ポアソン リンクを持つ混合モデルに difflsmeans を使用することは有効ですか?

係数の有意性を確認する別の方法はありますか? （これは以前に尋ねられたことは知っていますが、この分析の文脈で意味します）

前もって感謝します

GLMER の結果

difflsmeans の結果

以下のいくつかの係数のサイズが大きいことに疑問を持っている場合に備えて、私は金融/経済データを扱っています... 私の一般的な質問は、R の線形ランダム効果モデルから出力されるパラメーター係数のシミュレーションに関するものです。 R の同じモデルからモデル係数と分散共分散 (VCOV) 行列を使用して、ベータ係数のランダム サンプルを生成しようとしています。 mvtnorm{} パッケージの rmvnorm() 関数? この警告/問題に対処するにはどうすればよいですか?

警告メッセージ: In sqrt(ev$values) : NaN が生成されました

次の検索文字列を Google で検索すると、「rmvnorm、「警告メッセージ: sqrt(ev$values) 内: NaNs が生成されました」と表示されました: http://www.nickfieller.staff.shef.ac.uk/sheff-only /mvatasksols6-9.pdf 4 ページで、このエラーは「負の固有値」を示していますが、概念的または実際的に、負の固有値が何であるか、またはこの場合に生成される理由がわかりません。

2 番目の検索結果: [ http://www.r-tutor.com/r-introduction/basic-data-types/complex 2 -1 の平方根を取得しようとしたためにこのエラーが発生したことを示します。 「複雑な値ではありません」(-1 の平方根を取ることはできません)。

ここで、ベータのランダム生成で何が起こっているのか、そしてこれをどのように修正できるのかという疑問が残ります。

sessionInfo() R バージョン 3.0.2 (2013-09-25) プラットフォーム: x86_64-apple-darwin10.8.0 (64 ビット)

以下のパッケージ/バージョンを使用 mvtnorm_0.9-9994、lme4_1.1-5、Rcpp_0.10.3、Matrix_1.1-2-2、lattice_0.20-23

そのため、同様のテーマの投稿を多数見てきましたが、どれも私が必要としているものではないようです。または、提供されたソリューションを本当に理解していないだけです...

チンパンジーのデータを調べるために、lme4 で混合効果モデルを実行しました。私には、扶養者 (摂食時間) に影響を与える 2 つの要因 (攻撃率、交尾率) があります。

各予測変数と結果変数の間の関係を示す 2 つの散布図を作成したいと思いますが、モデル推定値から導出された線を引きたいと思います ((lm(y ~ x の abline ではありません) )) 完全な LMM に基づくものではなく、単純な回帰直線のみを与えるタイプ)。

これはggplot2でのみ可能であると感じていますが、実際にこれを行う方法を理解できていません。1 日のほとんどを本やフォーラムを調べて過ごしてきたので、彼らが何をしているのかを知っていれば、これがかなり簡単な答えになることを期待していました。

事前にヒントをありがとう！

アレックス