data <- read.table("InputFile", header=TRUE, row.names=1, sep=" ")
ma <- cor(data, NULL)
ma.pv <- pvclust(ma, method.hclust="single", method.dist="euclidean", nboot=1000)
plot(ma.pv, hang=-1)
InputFile がスペース区切りの CSV ファイルだとすると、
一行目で data というオブジェクトに InputFile のマトリックス情報を代入。
二行目で data の相関行列を計算して、ma に代入。
三行目で 相関行列 ma に対してユークリッド距離での最短距離法を用いて、
クラスタリングを実行し、結果を ma.pv に代入。
四行目で 得られた ma.pv を plot。
0 件のコメント:
コメントを投稿