RRによるやさしい統計学を参考に,乱数を使って任意のデータセットを作る関数を作ってみた。任意の相関係数を持つ二変数を生成する関数。return(result)}使い方は簡単で,相関係数と作るデータセットのサイズを引数で渡すだけ。cor2val(size=500,rho=0.3)summary(obs)cor(obs)次は任意の相関行列をもつデータセットを作る関数。return(result)}引数に生成するデータサイズ,変数の数,目標となる相関行列を渡す。corMatVal(500,5,target.Mat)ここでちょっと一工夫した。相関行列は対称で対角が1に決まってるんだから,全部入れるのめんどくさい。3変数の相関行列なら,r_12,r_13,r_23の情報だけでいいじゃない。いうことで,そういう横着をしたい人は,目標行列のところに額各要素のベクトルを渡して,corvecオプションをTRUEにすればいいようにした。corMatVal(500,3,c(0.2,0.3,0.4),corvec=TRUE)最後に,任意の因子負荷行列を持つデータセットの生成関数。return(result)}こんな風にして使う。FloadVal(size=50000,load.mat)備え付けの因子分析関数で確認。factanal(obs4,2)もっとも,因子間相関がない仮定。実際に斜交回転で分析すると相関が出て負荷量もおかしくなっちゃう。これについてはちょっと今度手を入れないと行けないですなー。AMデータセットを発生させる関数を利用したいと思いましたが一か所エラーが出て進めませんご教示願えませんでしょうか#引数に生成するデータサイズ,変数の数,目標となる相関行列を渡す。AMこれは失礼しました。変数の数の指定ができてないので,行列を作れてないのです。コメントをキャンセルメールアドレスが公開されることはありません。コメントデータの処理方法の詳細はこちらをご覧ください。2016で発表してきました文系学生に対する心理統計教育,回帰分析編で発表してきました。キャンセル投稿を送信できませんでした。メールアドレスを確認してください。メール送信チェックに失敗しました。このブログではメールでの投稿共有はできません。