---

[ツリー表示] [留意事項] [ワード検索] [過去ログ]

---

樋口です。以下、バブルプロットに関する補足です。

KH Coderの次のバージョンでは、対応分析の結果をバブルプロットであらわす
ことができるようになります。あるいは、こちらでひっそり公開しているテス
ト版では、すでにサポートされています：
http://khcoder.info/cgi-bin/bbs_khn/khcf.cgi?no=805&mode=allread

バブルプロットの作成には、座標データだけでなく、バブルの大きさのデータ
も必要になります。これも併せて出力するためには、以下のコマンドをお使い
下さい。

#--------------------------------------------------------------------
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 ( exists("b_size") ){                               # バブル対応
    if ( exists("biplot") ){
        if (biplot){
            v_size <- sqrt(n_total);
            if (std_radius){ # 大小差をデフォルメ
                v_size <- v_size / sd(v_size)
                v_size <- v_size - mean(v_size)
                v_size <- v_size * 5 * bubble_var / 100 + 10
                v_size <- neg_to_zero(v_size)
            }
            b_size <- c(b_size, v_size)
        }
    }
    out_data <- cbind(out_data, b_size)
    colnames(out_data) <- c("label","x","y","type","size")
}
if (nchar(savefile)) {                                 # 保存
    write.table(
        out_data,
        savefile,
        sep=",",
        quote=F,
        row.names=F
    )
    print( paste("saved: ",savefile) )
}
#--------------------------------------------------------------------