ggplot2とちょっと混乱します。私は自分のデータをプロットし、そのプロットにいくつかの欠けているデータを持っていますが、伝説ではありません。私は実際に40のグループ(FID)がプロットされ、いくつかは欠落しています。ggplot2を使ったプロットのデータがありません
Data48a <- read.table("Data48a.eigenvec", header = T, sep = "")
PlotData48a <- ggplot(Data48a, aes(x= PC1,y= PC2, group= FID, shape= FID, color= FID))
PlotData48a + geom_point(aes(color=FID)) +
geom_hline(aes(yintercept=0)) + geom_vline(aes(xintercept=0))+
scale_shape_manual(values=1:nlevels(Data48a$FID)) +
scale_color_manual(values=1:nlevels(Data48a$FID)) +
coord_equal(ratio = 1)
、ここでエラーです:私のコードがある
In grid.Call.graphics(L_points, x$x, x$y, x$pch, x$size) : unimplemented pch value '28'
私のデータは以下のようである:これはあなたの全体の問題を解決する場合
FID IID PC1 PC2 PC3
Ami AX-AM-013607-1-01 0.0302889 -0.0236855 0.0294569
Ami AX-AM-013608-1-01 0.0346382 -0.0232841 0.0175366
Ami AX-AM-013609-1-01 0.0396438 -0.0157996 0.0163948
Ami AX-AM-013610-1-01 0.0386903 -0.017313 0.0146025
Ami AX-AM-013611-1-01 0.0382764 -0.0151733 0.0222798
Ami AX-AM-013613-1-01 0.0416853 -0.00826403 0.0245685
Ami AX-AM-013614-1-01 0.0388888 -0.0113679 0.0281065
Ami AX-AM-013615-1-01 0.0438035 -0.0102871 0.0156296
Ami AX-AM-013616-1-01 0.0364263 -0.0211585 0.0246802
BatakKaro ID-KR-000001-1-01 0.0321929 0.00520535 0.00480784
BatakKaro ID-KR-000002-1-01 0.0345291 -0.00101008 -0.00264379
BatakKaro ID-KR-000012-1-01 0.0352078 -0.000998926 -0.00847158
BatakKaro ID-KR-000014-1-01 0.0306035 -0.00150602 -0.00365927
http://sape.inf.usi.ch/quick-reference/ggplot2/shape、それはそれぞれのユニークな形状を見つけることができないので、人々はちょうどFIDの2つの多くのレベルがある – OganM
あなたのオブジェクトをコピー&ペーストすることができますので、dput使用してみてくださいデータを共有します。 – scoa
あなたの答えはありがとうございます。私はグループの数を減らして簡単にしました。 –