「統計学関連なんでもあり」の過去ログ--- 044

No.13599　ポリコリック相関係数を用いた因子分析について　　【ＪＪ】　2010/10/12(Tue) 21:30

青木先生

添付のページ参考にさせていただきました。当方も統計，Rともに初心者で，試行錯誤しながら行っていてホームページを見ていてもよくわかからないところがあります。

ポリコリック相関係数を用いた因子分析を行いたいと思っているのですが，自分なりに試行錯誤してみて，作成したのが下のコマンドです。

factanal(polychor(データ),factors=2,rotation=promax,scores="regression")

データは，もともとの生データで17項目の質問紙となります。
ただ，このコマンドを実行してもpolychorは，相関係数を算出する関数ですので，やはり無理です。ここの部分がどうしても上手くいかず，お知恵を拝借したいと思っております。

上記の質問時に，factanal(covmat=cor(iris[1:4]), factors=1)と先生のお答えがありますが，covmat=cor(iris[1:4]）の部分は，iris[1:4]が生のデータで，その他の部分がポリコリック相関係数行列という認識でよろしいのでしょうか？

No.13600　Re: ポリコリック相関係数を用いた因子分析について　　【青木繁伸】　2010/10/12(Tue) 21:51

> このコマンドを実行してもpolychorは，相関係数を算出する関数ですので，やはり無理です。ここの部分がどうしても上手くいかず

どういう風にうまくいかないと？

> 上記の質問時に，factanal(covmat=cor(iris[1:4]), factors=1)と先生のお答えがありますが

スレッド番号は？

> covmat=cor(iris[1:4]

は，単に，相関係数行列を指定してるだけでは？

No.13601　Re: ポリコリック相関係数を用いた因子分析について　　【JJ】　2010/10/12(Tue) 22:49

早速の返信ありがとうございます。
スレッド番号は10969です。

説明不足で申し訳ありません。
非常に端的に言うと，ポリコリック相関係数を用いた因子分析を生データから行うコマンドが知りたいということになります。

↓のコマンドまではなんとか，たどり着けたの
factanal(***(データ),factors=2,rotation=promax,scores="regression")

***に入るコマンドがあるのでしょうか？

No.13602　追加です。　　【JJ】　2010/10/13(Wed) 00:19

http://aoki2.si.gunma-u.ac.jp/R/cfa.html

先生のホームページを参考にさせていただき，

factanal(covmat=polychoric(生データ),factors=2,promax,scores="regression")

としてみました。僕は，"covmat=polychoric(生データ)"で生データのポリコリック相関係数行列を算出しそれをfactanal関数で因子分析しているというように理解し，因子得点も算出できると思っているのですが，
'covmat' は有効な共分散リストではありません　というエラーが出てしまいます。

説明が不十分で申し訳ありません。なんとかご教授いただければありがたく思います。

No.13603　Re: ポリコリック相関係数を用いた因子分析について　　【青木繁伸】　2010/10/13(Wed) 08:22

> ***に入るコマンドがあるのでしょうか？

> "covmat=polychoric(生データ)"で生データのポリコリック相関係数行列を算出しそれをfactanal関数で因子分析しているというように理解し，因子得点も算出できると思っているのですが，
'covmat' は有効な共分散リストではありません　というエラーが出てしまいます。

polychoric 関数って，あなたか誰かが作ったのですか？R には標準ではありませんよね。ないものを指定しても，エラーになるだけでしょう。

なお，ポリコリック相関係数行列を指定できたとしても，因子得点は計算されないでしょう。factanal は，そのようには作られていませんから。

ちなみに，http://aoki2.si.gunma-u.ac.jp/R/cfa.html は確証的因子分析ですから，無関係ですね。

いろいろ試行錯誤をするのはよいのですけど，でたらめにやってはダメです。

No.13604　Re: ポリコリック相関係数を用いた因子分析について　　【JJ】　2010/10/13(Wed) 09:00

ご指導ありがとうございます。
まだまだ根本的なところがわかってないみたいです。

ありがとうございました。

No.13605　Re: ポリコリック相関係数を用いた因子分析について　　【青木繁伸】　2010/10/13(Wed) 21:26

他の所でも同時に質問していたのですね。
なんとか回答を得たいと思う気持ちはわかりますけど，あまり好ましいことではありません。
それと，「結局の所どうなったの？」ということには，最終的に問題が解決したかどうかにかかわらず，答えておいた方がよいですよね。
ネットでの質問における，エチケットです。

No.13607　Re: ポリコリック相関係数を用いた因子分析について　　【JJ】　2010/10/13(Wed) 23:29

焦っており，ご指摘のとおりエチケット不足だったと思います。
結果ですが，質問の回答の意味を調べるのに四苦八苦している状況であり，取り急ぎの期限に間に合いそうもありません。
今回はRを使用してポリコリック相関係数行列を作成することはできました。
しかし，それを使って因子分析を行うことはできなかったので，SPSSに読み込んで行いました。しかし，この方式では因子得点が算出できませんでした。
結果，ポリコリック相関係数を使用することあきらめ，SPSSでパッケージ化された因子分析を行う方向で考えています。
少しでも正確な分析を行おうと思ったのですが，背伸びをしすぎて結果的に無礼な形になってしまったことを反省しています。
ご指導ありがとうございました。