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

No.10457　Rによる数量化II類分析について　　【水野】　2009/07/23(Thu) 19:35

青木先生のプログラムは，サンプル数，説明変数の数，総カテゴリー数に制約はあるのですか？パソコン（ウィンドウズ版）のメモリー等の制約のみですか？5000サンプル，100変数，総カテゴリー数500（＝100変数＊5カテゴリー）を実行したいと考えております。よろしくご教示をお願いいたします。hironobu_

No.10459　Re: Rによる数量化II類分析について　　【青木繁伸】　2009/07/23(Thu) 21:27

実際にやってみてうまくいかなかったということですか？
まだ実際にやってみていないなら，やってみてうまくいかなかったら，R, JavaScript, Excel どの環境におけるプログラムか（URL で示してください）書いてください。うまくいかなかったということなら，どの環境かと同時に，その他の条件を書いてください。
阿呆タロウではないですが，仮定の質問にはお答えしがたいのです(^_^;)

5カテゴリーを持つ100変数で，4群の判別をする場合，library(MASS) の lda を使えば，分析可能のようですが。
n <- 5000
v <- 100
a <- cbind(
	x=matrix(sample(5, n*v, replace=TRUE), n, v),
	g=(sample(4, n, replace=TRUE)))
a <- data.frame(a)
a <- data.frame(sapply(a, factor))
colnames(a) <- c(paste("X", 1:100, sep=""), "g")
(ans <- lda(g~., data=a))

No.10463　Re: Rによる数量化II類分析について　　【水野】　2009/07/25(Sat) 22:10

実際の調査データが10月頃に作成される予定ですので，その時に実行してみます。
ありがとうございました。