analysisk-meansとは非階層クラスタリングの手法の一つです。クラスタの平均を用い、与えられたクラスタ数k個に分類することからk-meansと名付けられています。Rでk-meansによるクラスタリングを行う方法についてまとめます。データの準備クラスタリング結果の検証(validation)1.データの準備あるサンプルについてのデータが行に、あるパラメータについての各サンプルのデータが列に並んだデータフレームまたは行列を用意します。サンプル数をm、パラメータ数をnとすると、m×nのデータフレーム(または行列)となります。このデータフレームにはNAのデータを含めることはできないので予め取り除いておく必要があります。いう関数を使います(パッケージは不要です)。centers=クラスタ数)です。クラスタ数は1からサンプル数m-1までの整数です。centersに数値を指定した場合、クラスタの代表点の初期値はランダムに選ばれますが、ベクトルを与えることで任意に選ぶこともできます。他の引数としては、繰り返しの最大値を与えるiter.max、アルゴリズムを指定するalgorithmなどがあります。それぞれデフォルトではiter.max=10L、algorithm=”Hartigan-Wong”となっています。kmeansの返り値は9つの要素からなるリストです。$clusterでクラスタリングの結果がベクトルで得られます。3.結果の検証(validation)適切なクラスタ数は目的によって異なりますが、少なすぎるとクラスタ内のばらつきが大きくなってしまいますし、多すぎては意味がわからなくなります。分割の適切さの尺度の1つとしてクラスター内平方和があります。$withinssでそれぞれのクラスターについての平方和、$tot.withinssで平方和の合計が得られます。クラスター数kがある程度以上になると平方和はあまり変化しなくなります。Index:SPI)は、様々な時間スケールの干ばつを示す指標です。1.4)は、全球気候モデル(GCMs)のシミュレーションで得られた気候データである。