最初に書いた mclustBICとMclustの出力は 本来一緒のはずらしいけど， 今回は異なってしまった理由は不明...なぜなんだろう(´Д` ) あとクラスターの分布が違うのは 使った混合分布の型(球、楕円球)、体積、形と軸の向きが 異なるから． （最初"EEE"，二回目"…

いろいろ調べてたら もう一つやり方見つけたただこれやったら全然結果変わってしまった(´Д` )とりあえずメモ。 > res > plot(res) これよりEEVのクラスター数7が最適 > res > cl > cl[1:10] 1 2 3 4 5 6 7 8 9 10 6 7 6 7 6 6 7 3 6 6 > table(cl) cl 1 2 3 …

モデルに基づいたクラスター分析http://mjin.doshisha.ac.jp/R/29/29.html を参考にやってみたmclustパッケージをダウンロードまず書かれているように plot(EMclust) ってやったけどエラーどうやらmclustパッケージに EMclustは今入っていないらしい最適なモ…

apply関数を使って 行列処理をしてた時の エラーというか困ったこと func list2 for(i in 1:3){ if(list[i]==list[7]) list2[2,i] if(list[i+3]==list[8]) list2[3,i] } if(sum(list2[2,])==1) return(list2[1,list2[2,]==1]) else return(list2[1,list2[3,]…

原因が本当にこれかは不明だけど「はい」「いいえ」みたいな二値データは 基本的には因子分析では扱われなくて、因子分析したいなら カテゴリカル因子分析ってやつを使うらしいhttp://home.hiroshima-u.ac.jp/keiroh/maeda/statsarekore/gedoufa.htmlテトラ…

アンケート結果を因子分析して 因子を抽出しようとしているのだけど、なかなかうまくいかない 使っているのは、 factanal関数 factanal2関数（青木先生） fa関数（psychパッケージ） とかいろいろデータ形式は 「いいものは長く使いたいですか？」：「はい」…

以前に書いたtable()関数を用いて クロス集計が行えるentertainment-lab: Rのtable関数に関する調査ついでにここで扱われているcut関数とは 数値ベクトルを指定した区間に分類するための関数使い方は http://www.okada.jp.org/RWiki/?R%A5%D7%A5%ED%A5%B0%A5…

土曜日は片所さんとMTGそこでやること。１）N社コンテスト内容確認 ・ブランドロイヤル ・バラエティシーカー ・新規購買者 の３つの群に対して・認知 ・意欲 ・購買 の３つの反応を分析した結果はいまいちだけど とりあえず報告データの切り取り方とか もう…

並び替え関数 sort()とorder() 違いが書かれてるサイト。ああああ: R: sort と order の違い

度数分布を表示する関数 table() > table(iris$length) 4.3 4.4 4.5 4.6 4.7 4.8 4.9 1 3 1 4 2 5 6 irisデータのlengthの度数分布 >table(iris[1]) でもいける（lengthが一列目） xtabs() > xtabs(Freq~Inocome,data=dataBox) Income 0-50 50-100 100-150 >…

邪魔な欠損値をどう処理するかベクトルとしてうまく操作して やつらを回避したデータを取り出す dataBoxの二列目がNAじゃないベクトルを取り出す > data データNAを0に置き換える > x > c [1] 0 1 0 1 0 0 1ifelse()は ifelse(条件式,TRUEの場合に返す値,FAL…

sink() sink(file="test.txt") または sink(file = "test.txt", append = TRUE) # 既存のファイルに追記~~~なんかしらの操作~~~sink()で終了

はてなで R用雑記帳開設 なう。