未来はキミドリイロ

サイコメトリックアイドルを目指す心理学徒の勉強部屋です。勉強・趣味などについて書いています。

R

【R】行列における複数行の繰り返し

R

備忘をかねて,Rでの行列操作,複数行の繰り返しについてまとめておきます。 例えば,手元のデータが次のようになっていることを考えます。 y1 print(y1) [,1] [,2] [,3] [1,] 1 2 3 [2,] 4 5 6 [3,] 7 8 9 このとき,この行列を行方向に3回反復した9行3列の…

mvrnorm関数で「'Sigma' is not positive definite」というエラーが出たときの対処法

R

シミュレーションなどのために,多変量正規分布からのサンプリングを行う関数として{MASS}パッケージのmvrnorm関数があります。 この関数は次のようなプログラム例を通してサンプリングを行います。 library(MASS) #3変量正規分布の例 #平均値として0が3つの…

【R】任意の因子負荷から架空データを発生させる(1因子モデル編)

R

久しぶりの更新です。 更新したいネタはたくさんあるのですが,なかなか筆を執ることができません。そんなどうでもいい個人的事情はおいておきまして,今回は備忘もかねて,統計ソフトRの{psych}パッケージを利用して任意の因子負荷行列から架空データを発生…