R標準のplot()でカテゴリデータにそれぞれプロットのマーカーを割り当て

Question

前提・実現したいこと

R 3.5.0で主成分分析の結果を2次元平面上にプロットしたのですが、R標準のplot()で3種類のカテゴリデータにそれぞれプロットのマーカーを割り当てたいです。(色分けはできるのですが…)

プロットのマーカーは、
ラベルsetosaには1の○, versicolorには2の△, virginicaには4の×をつけたいです。

該当のソースコード

data(iris)
data <- iris[1:4]
prcomp.obj <- prcomp(data, scale=TRUE) # 主成分分析

pc1 <- prcomp.obj$x[,1] # 第一主成分得点
pc2 <- prcomp.obj$x[,2] # 第二主成分得点

label <- as.factor(iris[,5]) # 分類ラベル
percent <- summary(prcomp.obj)$importance[3,2] * 100 # 累積寄与率

plot(pc1, pc2, col = label, main = paste(percent, "%"))

Answer 1

何度もすみません。訂正です。
長くなったのでコメント欄から回答欄に移動しました。

colはfactor型、pchはnumeric型で認識されているようです。
そのため、factor型を書き直してプロットしてみてはいかがでしょうか。

label2<-as.numeric(label)
label2[label2==3]<-4
plot(pc1, pc2, col = label, pch = label2, main = paste(percent, "%"))

↑最後の1行をこの3行に書きなおすとおそらくご希望のプロットになるかと・・・
（色付け不要なら col = label は削除で。）

Comments

[Mikya]

追加です。こちらの方でプロットの形がカテゴリー別にできました。ありがとうございます。
私の勘違いでコメント欄の方についていた前の方をベストアンサーにしてしまいました…。大変申し訳無いです。高評価とフォローをさせていただきます。

[fermat]

お気になさらず(^^)解決して何よりです～

Answer 2

colの代わりにpchパラメータを使えばよいのでは。
https://stats.biopapyrus.jp/r/graph/scattergram.html#pch

Comments

[Mikya]

ご回答ありがとうございます。
col の代わりにpchをいれると
plot.xy(xy, type, ...) でエラー: プロット記号が不正です
というエラーが出てしまいます。
カテゴリ別にpch パラメータは使えないのでしょうか？

[fermat]

plot(pc1, pc2, col = label, pch = c(1,2,4), main = paste(percent, "%"))
これで通りましたがどうでしょうか。

[Mikya]

試してみたところ解決いたしました。
ありがとうございました。