Colorless Green Ideas:R言語で二変量正規分布を扱う

はじめにR本体だけでは二変量正規分布（二次元正規分布）をはじめとした多変量正規分布（多次元正規分布）を扱うことはできません。このため、自分で定義するか、MASSあるいはmvtnormというライブラリを読み込む必要があります。単に、多次元の正規乱数を作成するだけなら、MASSパッケージを使うのが楽です。確率密度や累積密度を求めたかったら、mvtnormパッケージを使うのがよいでしょう。 以下では、二変量正規分布を中心に扱いますが、三次元以上も同様に処理することができます。 正規分布とは？二変量正規分布を扱う前に、もっとも基本的な正規分布、すなわち一次元の正規分布をRでどう扱うかを確認しましょう。 一変量正規分布（一次元正規分布）は、μ，σ2の2つのパラメータを持ちます。μは平均、σ2は分散に相当します。正規分布は、N(μ, σ2)と表します。例えば、平均が4で、分散が9の正規分布は、N(4,

[dora223]

“library(MASS)”