4

だから私はこれが初めてです。次のデータ マトリックスで PCoA を実行する必要があります。ADE4、labdsv、Ginko、Aabel ソフトウェアを使用して解析を実行できます。気になるのは、散布図でラベルを色分けする方法です。私のマトリックスは、次の順序で存在/不在マトリックスです。

SpecieName Value1 Value2
A1         0      1
A2         1      1
A3         1      1
B1         0      0
B2         0      1
E1         1      0
E2         0      0

私が望むのはA1、 、A2、およびA3を赤でB1B2青で、すべてのEものを黒で表すことです。どんな助けでも大歓迎です。

4

1 に答える 1

6

これらのグループを示す係数をプロット コマンドに渡すだけです。

data = read.table('data.txt', header=T)

data.pca = prcomp(data[,-1])

groups = factor(gsub('(.).', '\\1', data$SpecieName))

plot(data.pca$x, col=groups)

ここに画像の説明を入力

また、特定の色を設定したい場合は、いつでも同じ方法でカスタム リストにインデックスを付けることができます。

cols = c('red', 'blue', 'black')[groups]
plot(data.pca$x, col=cols)
于 2011-12-03T23:25:11.873 に答える