Rに別の `` rnorm`の引数として `` rnorm`を入れるのはどういう意味ですか？

Question

rnormが別の引数の1つとして使用されている場合、意味が分かりません。rnorm？ （詳細は以下で説明します）

たとえば、私のRコードの最初の行では、私はrnorm()を使用し、これをrnorm()：muと呼びます。

muは、xからなる。

ここで、mu自体をmeanという引数に「distribution」という新しいrnorm()という名前を付けてみましょう。

私自身の質問は、meanというこの新しいrnorm()の引数として、xが10,000となっていますか？

P.S .: mean任意のnormal distributionの引数は単一の数値でもよく、単一の単一の平均であっても、単一の完全な法線を持ちます。さて、どのように来て、10,000のmu値を使用しても、結果は単一の正常値になりますか？

mu <- rnorm(1e4 , 178 , 20)   ; plot(density(mu)) 
distribution <- rnorm(1e4 , mu , 1) ; plot(density(distribution)) 

Answer 1

あなたdistributionは条件付き密度です。 plot(density(distribution))で描画する濃度は、限界濃度です。

統計的に言えば、最初に通常のランダム変数mu ~ N(178, 20)、次に別のランダム変数y | mu ~ N(mu, 1)があります。あなたが生産するプロットは、限界密度がyです。

P(y)は、muを統合している数学的には、ジョイント分布P(y | mu) * p(mu)の積分です。

@李哲源ZheyuanLi、ああ！だから、我々がvnをrnormの平均引数またはsd引数として使うとき、最終的な単一のプロットは積分の結果です、そうですか？

つまり、周囲の分布からサンプリングしています。密度推定値は、サンプルからのモンテカルロ積分に近似する。

---

この種のものは、多くの場合、ベイズ計算で見られます。 Toy R code on Bayesian inference for mean of a normal distribution [data of snowfall amount]は完全な例ですが、積分は数値積分によって計算されます。