主成分分析:R言語

分析に使われる場面

ライブラリとデータの読み込み

  • ライブラリの読み込み
    library(mvpart)
  • データの読み込み
    onsen <- read.table("onsen.txt",sep=",",header=T)
    rownames(onsen) <- onsen[,2]

主成分分析

  • prcomp()関数による主成分分析
    onsen.prc <- prcomp(onsen[,3:8])
  • 結果の表示
    onsen.prc
    summary(onsen.prc)
    #onsen.prc$rotation
  • 第一主成分と第二主成分の散布図の作成
    plot(onsen.prc$rotation[,1],onsen.prc$rotation[,2])
    plot(onsen.prc$rotation[,1],onsen.prc$rotation[,2],type="n")
    text(onsen.prc$rotation[,1],onsen.prc$rotation[,2],colnames(onsen[,3:8]))
  • 主成分得点の散布図の作成
    onsen.prc$x
    plot(onsen.prc$x[,1],onsen.prc$x[,2])
    plot(onsen.prc$x[,1],onsen.prc$x[,2], type="n")
    text(onsen.prc$x[,1],onsen.prc$x[,2], rownames(onsen[,3:8]))
  • 散布図の作成
    biplot(onsen.prc)

トップ   編集 凍結解除 差分 バックアップ 添付 複製 名前変更 リロード   新規 一覧 単語検索 最終更新   ヘルプ   最終更新のRSS
Last-modified: 2009-02-03 (火) 19:40:52 (3820d)