R
備忘をかねて,Rでの行列操作,複数行の繰り返しについてまとめておきます。 例えば,手元のデータが次のようになっていることを考えます。 y1 print(y1) [,1] [,2] [,3] [1,] 1 2 3 [2,] 4 5 6 [3,] 7 8 9 このとき,この行列を行方向に3回反復した9行3列の…
シミュレーションなどのために,多変量正規分布からのサンプリングを行う関数として{MASS}パッケージのmvrnorm関数があります。 この関数は次のようなプログラム例を通してサンプリングを行います。 library(MASS) #3変量正規分布の例 #平均値として0が3つの…
久しぶりの更新です。 更新したいネタはたくさんあるのですが,なかなか筆を執ることができません。そんなどうでもいい個人的事情はおいておきまして,今回は備忘もかねて,統計ソフトRの{psych}パッケージを利用して任意の因子負荷行列から架空データを発生…