1

MCMC アルゴリズムに反復の配列があります。行は、分布からの抽選を表します。列は、分布のパラメーター (変数) を表します。説明を簡単にするために、2 つの変数、5 回の反復を仮定します。ので、私は持っています：

> draws <- data.frame( iteration = c(1:5),
                       alpha     = rnorm(5,0,1),
                       beta      = rnorm(5,0,1))

  iteration      alpha       beta
1         1 -0.3157940  0.2122465
2         2  1.0087298 -0.2346733
3         3  1.0366165  0.3472915
4         4 -2.4256564  0.9863279
5         5 -0.6089072 -1.1213000

データセットを溶かすと、次のようになります。

> melt(draws)
Using  as id variables
    variable      value
1  iteration  1.0000000
2  iteration  2.0000000
3  iteration  3.0000000
4  iteration  4.0000000
5  iteration  5.0000000
6      alpha -0.1042616
7      alpha  1.0707001
8      alpha  0.2166865
9      alpha  0.0771617
10     alpha -0.8893614
11      beta -0.4846693
12      beta -1.5950729
13      beta -0.7178340
14      beta  1.0149766
15      beta -0.3128256

しかし、（手作業で編集した）と同等のものを得るために、反復を延期したいと思います。

> melt(draws)
Using  as id variables
   iteration  variable      value
1          1     alpha -0.1042616
2          2     alpha  1.0707001
3          3     alpha  0.2166865
4          4     alpha  0.0771617
5          5     alpha -0.8893614
6          1      beta -0.4846693
7          2      beta -1.5950729
8          3      beta -0.7178340
9          4      beta  1.0149766
10         5      beta -0.3128256
4

2 に答える 2

7

変数idを以下に指定しますmelt：

melt(draws, id = "iteration")

与える：

   iteration variable       value
1          1    alpha -0.02765436
2          2    alpha -1.42138702
3          3    alpha  0.83525096
4          4    alpha -1.10677555
5          5    alpha  0.72465936
6          1     beta  0.59269720
7          2     beta -0.32164072
8          3     beta -1.31097204
9          4     beta  0.94993620
10         5     beta  0.20919169
于 2011-07-03T17:51:52.280 に答える
2

ばぁ。私はいつも答えを見つける直前に質問をします...

を読んhelp(melt.array)でいたのですが、data.frame に変換して質問を投稿したところ、最終的に の回答にたどり着きましたhelp(melt.data.frame)。

私が欲しいものを得るために、私は使用します：

myMelt <- melt( draws, id.var = "iteration" );

次に、ファセットプロットを作成できるようにします。

ggplot(myMelt, aes(x = iteration,y = value)) + geom_point() + stat_smooth() + facet_grid(variable ~ ., scales="free")

ここに画像の説明を入力

于 2011-07-03T17:52:54.920 に答える