Re: 対応分析の座標データの取り出し (HIGUCHI Koichi) KH Coder 旧掲示板
[ツリー表示] [留意事項] [ワード検索] [過去ログ]

  [No.825] Re: 対応分析の座標データの取り出し 投稿者:HIGUCHI Koichi  投稿日:2011/08/06(Sat) 18:34:46

はじめまして、樋口です。書き込みありがとうございます。

対応分析にはRを使っていますので、Rからデータを取り出す形となります。手
順は以下の通りです。

[1] KH Coderで対応分析を実行します。

[2] 「保存」ボタンをクリックして、「R Source」形式で対応分析の結果を保
存します。

[3] kh_coder.exeと同じ場所にあるRgui.batをダブルクリックして、KH Coder
に付属のRを起動します。

[4] Rのコンソール画面に[2]で保存したファイルをドラッグします。※これに
よってR上で、同様の対応分析が実行されます。

[5] 以下のコマンドを実行する(Rのコンソール画面に貼り付けてエンター・
キーを押す)ことで、対応分析の入力となったデータ表がC:\temp0.csvに保存
されます。

#--------------------------------------------------------------------
savefile <- "C:\\temp0.csv"
out_data <- cbind(rownames(d), d)
if (nchar(savefile)) {
    write.table(
        out_data,
        savefile,
        sep=",",
        quote=F,
        row.names=F
    )
    print( paste("saved: ",savefile) )
}
#--------------------------------------------------------------------

[6] 以下のコマンドを実行する(Rのコンソール画面に貼り付けてエンター・
キーを押す)ことで、プロット作成に必要な座標データ等がC:\temp1.csvに保
存されます。

#--------------------------------------------------------------------
savefile <- "C:\\temp1.csv"
out_data <- cbind(c$cscore[,d_x], c$cscore[,d_y], 1 )  # 語
if ( exists("biplot") ){                               # 変数・見出し
    if (biplot){
        out_data <- rbind(
            out_data, 
            cbind(c$rscore[,d_x], c$rscore[,d_y], 2)
        )
    }
}
out_data <- cbind(rownames(out_data), out_data)        # 仕上げ
colnames(out_data) <- c("label","x","y","type")
if (nchar(savefile)) {                                 # 保存
    write.table(
        out_data,
        savefile,
        sep=",",
        quote=F,
        row.names=F
    )
    print( paste("saved: ",savefile) )
}
#--------------------------------------------------------------------


なお、[6]についても[7]についても、最初の1行「savefile <- "C:\\temp1.cs
v"」を次の内容に変更すると、GUIで保存ファイル名を決めることができます。

#--------------------------------------------------------------------
savefile <- NULL
require(tcltk)
savefile <- tclvalue(
    tkgetSaveFile(
        filetypes = "{{CSV Files} {.csv}}",
        defaultextension=".txt"
    )
)
#--------------------------------------------------------------------


- 関連一覧ツリー (■ をクリックするとツリー全体を一括表示します)