問題タブ [mixed-models]
For questions regarding programming in ECMAScript (JavaScript/JS) and its various dialects/implementations (excluding ActionScript). Note JavaScript is NOT the same as Java! Please include all relevant tags on your question; e.g., [node.js], [jquery], [json], [reactjs], [angular], [ember.js], [vue.js], [typescript], [svelte], etc.
r - glmerによるモデル仕様
GLMM をセットアップしようとしていますが、R で階層データ構造を表現するのに問題がありますglmer
。私のデータは次の構造を持っています:
wherevisit
は 内にネストされyear
、これは 内にネストされますsite
(つまり、1 年間に最大 5 回の訪問があり、サイトごとに最大 10 年間のデータがある可能性があります)。
私は on の効果に興味がありますがyear
、ランダムな効果としてy
扱いたいと思ってsite
います (つまり、切片と勾配 [年間] はサイトによって異なります)。
次のコマンドを使用しました。
モデルを正しく表現しましたか?
r - 重み付けされた非線形混合モデルの「適合しない配列」エラー
nlme パッケージの nlme() 関数を使用して非線形混合モデルを正常に適合させましたが、重み引数を含めてそのモデルを改善しようとしても、解釈できないエラーしか得られないため、できませんでした。私の問題の正確な原因を特定するために。
非加重モデルに次のコードを使用しました。
これは加重モデルです。
次のエラーが表示されます。
conLin$Xy * varWeights(object) のエラー: 適合しない配列
これは、行列とベクトルの積に関係していることを理解しています。そのうちの1つは同じ次元ではないため、操作を正常に実行できます。ただし、これらのオブジェクトの寸法を確認する方法や、それらにアクセスして問題を解決しようとする方法さえわかりません。
varFixed() 引数を varIdent に変更すると、メッセージは変わりますが、自分で収集できる限り、基本的に同じことを伝えています。
recalc.varFunc(object[[i]], conLin) のエラー: ディム [製品 3618] がオブジェクトの長さと一致しません [37776] さらに: 警告メッセージ: in conLin$Xy * varWeights(object) : オブジェクトの長さが長くなります短いオブジェクトの長さの倍数ではありません
この問題に関するヘルプは大歓迎です!
これが私のデータのサンプルです:
前もって感謝します。
r - R の部分的にネストされた/ブロックされた実験計画
実験のデザインには 10 人の参加者が含まれます。それらのすべては、治療のために条件 A、B、C、D を通過しますが、参加者 1 ~ 5 は条件 E、F を通過し、参加者 6 ~ 10 は条件 G、H を通過します。
lme 関数を備えた nlme パッケージを使用して、欠落しているデータを処理し、参加者のリストごとの削除を防ぎます。測定変数 = DV、固定効果 = 条件、変量効果 = 参加者)。すべてがちょうど交差したとき、これは私が持っているものです:
最初の部分 (条件 A、B、C、D) が交差し、2 番目の部分 E、F および G、H がネストされている場合に設定される統計は何ですか....ヘルプやガイダンスをいただければ幸いです! ありがとう。
r - lmer の代替最適化アルゴリズム
パッケージ内の関数lmer
は、最適化アルゴリズムとしてパッケージからlme4
デフォルトで使用します。bobyqa
minqa
次の投稿https://stat.ethz.ch/pipermail/r-sig-mixed-models/2013q1/020075.htmlminqa
によると、パッケージ内の他の最適化アルゴリズムも使用できます。
uobyqa
またはnewuoa
の最適化アルゴリズムとしてどのように使用できますlmer
か?