なかのヒト

長くなったので”MY NOTE 2″から引っ越しました

多変量解析　R　主成分分析　因子分析　数量化Ⅰ類からⅢ類まで
できるようにがんばります！！

注）全く知識のない私のメモですので他の方には参考になりません・・・^^;

多変量解析概論
な、なんとなく多変量解析の事わかるw　主成分分析のこともわかりやすく説明してある。

Rに関してはRjpWiki
言われた通り、筑波のミラーサイトから。ダウンロードしました。
で？というわけでdemoを使ってみる事に。
使い方から全くわかりません。

Rの簡単な使い方
Rのデモ　demo()と打ってみると，デモ可能な一覧。
例）demo(graphics)
関係ないけど、このヒトの『絵日記でも書いてみようか』にウケた＞＜　眠気がとびました。ありがとうw
で、デモの図がでないんだけど・・・あ、enter押したら見れたw

Rを使って実際に統計解析してみよう
データの入れ方からわかりません。ということでここの例を試してみる。
このデータはどこで保存してあるの？
Excelで作って”test.data”で保存
でもエラーが
Excelで作って.dataで保存してもExcel形式じゃダメだという事で（常識なの？？）フォーマットをCSVに変更。
『スペース区切り』もあるからどっちがいいかいろいろ試してみてねと言われる。。。
意味が分かりませんw試すか、試してみるしかないのか。。。
それでもできない。『マルチバイト文字列が不正』か。
UTFとshift JISの問題？例題をコピペで.dataだとちゃんと動く。
Excelで保存したデータをEditorに持って行くとshift JISになるのが関係あるのでしょうか。眠い。

（CSVについてみて見る・・・いくつかのフィールドをコンマで区切ったテキストデータおよび
テキストファイルExcelで作って『.csv』で保存すればいいのかあ。）

ここにこんなことが書いてありました。
「\ を表すには \\ とするしか無いのですが、今後国際化の過程で \ → 0xA5(UTF-8) とかアレゲな事になるかも
しれませんMS-CP932 , SHIFT_JIS 等からコード変換でもよく過去に話題になりました．（なかまさんより）」
アレゲって何？ああ、なんだか疲れて来てて変な事調べてる。とにかくSHIFT_JISでなんかあったらしい。
なかまさん教えて＞＜

Excelの保存をUnicode テキスト(UTF-16)でしてみる。
マルチバイト文字列が不正です
追加情報: Warning message:
In read.table(“test2.txt”, header = TRUE) :
‘test2.txt’ の readTableHeader で不完全な最終行が見つかりました

Editorに持ってって確認したら漢字が変な文字になってた・・・
やっぱりEditorでUnicodeに変換するしかないのかな。

あああ、って結論で次の行程に進んだら解決策がw　次を読んどくのも大事だよね＞＜

えっとじゃあ、実際のデータ使って主成分分析を・・・で、何からしていいのかw
再び群馬大学の青木先生のところで考えよ。
ソースが打ち込めない・・・改行使っちゃダメなのでしょうか・・・打ち込んじゃうと上の行に戻れないものなんですね。
め、めんどくさいwどっかで作ってペーストが一番いいのかなあ。

xlimの意味から分からない。

もうちょっと調べよう・・・。

お、こっちの方がわかりやすいかな。

データのSpeciesは数字に置き換えないとだめか。また作り直し。
結果らしきものがでたけど、イマイチわからない＞＜
やっぱり１から勉強しないとだめか。（当たり前〜）　

行列の作成

Fisherのirisデータ使ったので一回試してた。コピペで上手くいく（それでいいのかー！）

そして実際のデータ使ってやってみた。結果出たので送った・・・あってるかどきどき・・・結果待ち。

生態学のデータ解析

plotのX,Y軸を０で交差するようにできるのでしょか。
グラフ
rnorm・・・norm 正規分布　頭が拒否するw
　　　　　　rXXX 乱数　rnorm(n) 正規分布に従う乱数をn個生成
axes・・・FALSEだと軸は生成されない
ann=F 軸のラベルを描かない
axis()・・・軸
adjの値はテキスト文字列が整形される仕方。0は左揃い。0.5は中心。1は右揃い。ってaxisのなかにこれ入れるのかなあ^^;
las・・・軸ラベルの形式　las=0 常に軸に平行（規定）
　　　　　　　　　　　　 las=1 常に水平
　　　　　　　　　　　　 las=2 常に軸に垂直
　　　　　　　　　　　　 las=3 常に垂直
pchプロット点の形を指定する

軸はaxisでなんとかなりそう。

PC3の３次元かあpersp()

ここ『統計解析フリーソフト R の備忘録頁』のサイトホントにわかりやすくて助かってます。
基本的な事。
＜オブジェクトについて・・・＞
objects() rm()

＜オブジェクトの表示について＞
作業ディレクトリ指定に print() cat()

＜文字列の表示＞
sprintf()

＜書式付きでオブジェクトの表示＞
この辺でもう意味が・・・また後でわかるかも

＜出力をファイルに送る＞
“sink” 作業ディレクトリの場所に指定した名前をつけたファイルを作り、Rの出力を送ってくれる。

“test_00″というファイルを作った。結果だけが送られている。Rの画面上には表示されない。
sink()で閉じる

＜ベクトル＞
R は基本的データを，要素の個数が一つだけのベクトルとして扱っている。・・・意味わからないw
と、とにかく進もうw
ベクトルを作成する基本的な関数はc()
ベクトルの要素の個数をベクトルの長さと呼ぶ。length()

c()以外の規則的なベクトルを生成する関数　rep() seq()
seq:単調増加（減少）する数値列

seq() sequence()

rep() times=cのcが数列可って数列がよくわかんない

rep() unique() numeric()

＜ベクトルをc()で定義する＞

＜ベクトルの要素＞
ああ！yのベクトルの中身長過ぎてところどころきれてる！！
ベクトルのベクトルの中の数を「要素」と呼ぶ。左から順に 1, 2, ・・・ と番号が振られている。

character・・・文字列（”"や＼を前に置いて表す）

＆と｜

NA・・・欠損値・不定データ
NaN・・・非数（数では表せないもの）
Inf・・・無限大

＜ベクトルの計算＞
◯同じ長さのベクトルの場合（＋，ー，＊，／，％／％整数商，％％剰余，＾べき乗）

◯長さのちがうベクトル
短い方のベクトルの要素が循環して使用される。

◯ベクトル同士の比較

相変わらずバカな事をしてしまう・・・コメントのerrorについてw