Computing Monte Carlo Standard Error for an arbitrary function

オープン
#57 コメント 10 件 リアクション 0 件 担当者 1 名 GitHub で見る

@MansMeg がすでに取り組んでいます。

2020年1月3日 から。

評価

この issue はまだ評価されていません。

説明

feature

Hi!

Me, @paul-buerkner and @avehtari discussed the need for computing the MCSE for different functions of interest, such as the marginal variances, skewness and kurtosis of individual parameters. I think this functionality is probably best to implement directly into the posterior package. Aki suggested to use the batch means and subsampling approaches of https://projecteuclid.org/download/pdfview_1/euclid.ejs/1417615759 and I could implement this into posterior since I need to implement it anyway.

The most important part for me to know is how you would like the API to look. This is my first suggestion, feel free to comment or come up with a better suggestion.

monte_carlo_standard_error(x, FUN = mean)

where x is a draws object and FUN is an arbitrary function that operates on the individual parameters (rv?) of the draws object. It returns an object similar to summary() but with column with the names:
[parameter_name]_[function_name] and [parameter_name]_[function_name]_se

Any thoughts or comments before I get started?

主要言語
R
スター
172
フォーク
27
平均マージ
5日 12時間
マージ済み PR(30日)
5

環境構築

はじめの一歩

  1. issue を最後まで読み、次にプロジェクトのコントリビューションガイドを読みます。
  2. 着手することを issue にコメントします — 二人が同じ作業をするのを防げます。
  3. リポジトリをフォークし、ブランチを切って変更します。
  4. issue 番号を参照したプルリクエストを送ります。

stan-dev/posterior のほかの issue

stan-dev/posterior の issue をすべて見る

似ている issue

R の issue をもっと見る

新しい issue をメールで受け取る

初心者向けの GitHub issue を短くまとめたダイジェスト。