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

★ R で主成分分析 ★

　220　R で主成分分析　　青木繁伸　　2002/01/10 (木) 19:08
　　222　Re: R で主成分分析　　へむへむ　　2002/01/10 (木) 20:49
　　　223　Re^2: R で主成分分析　　青木繁伸　　2002/01/10 (木) 21:03

220.　R で主成分分析　　青木繁伸　　2002/01/10 (木) 19:08

my.pca <- function(d) { nc <- ncol(d) result <- eigen(cor(d)) eval <- result$values evec <- result$vectors temp3 <- matrix(eval, nc, nc, byrow=TRUE) sqrt(temp3)*evec } d <- matrix(c(1,2,3,4,5, 3,2,1,2,5, 3,2,1,4,6),5,3) my.pca(d)
みたいになります。

でも，これってほとんどプログラムじゃないですね(^_^;)
簡単すぎます!!!(^_^)

さらには，きっとどこかに（まだ捜していない）pca なんていう関数が作られているに違いない。

すると，主成分分析は
pca(d)
で終わってしまうわけです。

　　　　　[このページのトップへ]

222.　Re: R で主成分分析　　へむへむ　　2002/01/10 (木) 20:49

> さらには，きっとどこかに（まだ捜していない）pca なんていう関数が作られているに違いない。

mva packageの中のprincomp の事でしょうか?

　　　　　[このページのトップへ]

223.　Re^2: R で主成分分析　　青木繁伸　　2002/01/10 (木) 21:03

> mva packageの中のprincomp の事でしょうか?

それは，候補に入っていません(^_^)
普通の人が使うには難しいです。。。って，普通の人ってどういう人か(^_^)

この一連の操作を関数としてひとまとめにすれば使いやすいのでしょう。
ただ，それは，ブラックボックスを作ることになりますね。

なかなか，難しいですね。

　　　　　[このページのトップへ]