問題タブ [brms]

brms パッケージで ranef() を使用してランダム効果の結果を管理しています。

「Population2003」の結果を抽出したいと思います。

モデル内の特定の変数の結果を抽出するにはどうすればよいですか?

brms パッケージで順序回帰モデルを構築しました。

次に、結果の条件付き効果をプロットしようとしました (「序数 = TRUE」は廃止されました)。

これらの推定値は、特定のカテゴリで結果が得られる確率であり、特定の条件下ですべてのカテゴリの合計が 1 になるはずです。conditional_effects() からではなく、posterior_epred() から。

同様の状況を検索すると、デフォルトで conditional_effects がサンプルの平均ではなく中央値を返すことがわかりました。

そこで、次のページ（連続予測子を持つ順序モデル）を参照して、もう一度試してみました。 https://cran.r-project.org/web/packages/tidybayes/vignettes/tidy-brms.html#ordinal-model-with-categorical-predictor

しかし、「RDB2000popが見つかりませんでした」というエラーが表示されました。RDB2000pop は説明変数の 1 つです。そして、「dfpop_chenv」にあるはずです。

問題は「add_epred_draws()」に関連しているようです。「add_epre_draws()」を単純に使用しようとすると、同じエラーが表示されます。

各カテゴリの予測確率をプロットするにはどうすればよいですか?

サンプル サイズである標準誤差 (SE) と回帰の重みを組み込んだ加法モデルのモデルを構築しています。私が使用しているデータは文献から得たものであるため、アプローチはメタ分析的です。一部の著者は、生データ、つまり、サンプル サイズ = 1、標準誤差 = 0 の個々の観測値を公開しています。一部の著者は、サンプル サイズ > 1 の平均 +/- SE を報告しています。

問題は、使用している分布のファミリがログ リンクを持つガウス分布であり、0 の SE を入力するとエラーが発生することです。モデルの構造と実行しようとしたときの出力は次のとおりです。

fit1 <- brm(response|weights(Weight) + se(SE) ~ X1 + X2 + X3..., family = gaussian(link="log"), data = Data)

出力：

brms で平均と個々のサンプルを組み合わせる方法があるかどうか、または前進する方法について提案があるかどうか疑問に思っていますか?

お時間をいただきありがとうございます！

brms: 複数のカテゴリ変数を使用してモデルをセットアップするにはどうすればよいですか? すべてのレベルが存在し、一般的な切片に焼き付けられるものはありません。

たとえば、予測子があるとします。

• 性別 (2 レベル)

• educ (教育レベル、3 レベル)

やっている

係数として両方のレベルの性別を存在させることができるようです。ただし、これを次のように変更すると:

個別に 2 つの性別を引き続き持つことにつながりますが、今では 1 つのレベルの教育がインターセプトに組み込まれています。これが望ましいかどうかはわかりません。3 つのレベルすべてを個別に表示し、インターセプトに焼き付けないようにする方法を知りたいです。

やっている

助けにならない。