2016-04-05 23 views
0

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 
+2

http://sape.inf.usi.ch/quick-reference/ggplot2/shape、それはそれぞれのユニークな形状を見つけることができないので、人々はちょうどFIDの2つの多くのレベルがある – OganM

+2

あなたのオブジェクトをコピー&ペーストすることができますので、dput使用してみてくださいデータを共有します。 – scoa

+0

あなたの答えはありがとうございます。私はグループの数を減らして簡単にしました。 –

答えて

0

私は知らないが、これは、エラーの原因を強調する必要があります。

40レベルのデータまたはグループをプロットしようとしていますが、プロット文字は25個しかありません(pch)。エラーを取り除くためには、グループ数を25以下に制限しなければならないと思います。

require(ggplot2) 

# no error 
xx = 1:25 
ggplot() + geom_point(aes(x = xx,y = xx), col=xx,shape=xx) 

# error 
xx = 1:26 
ggplot() + geom_point(aes(x = xx,y = xx), col=xx,shape=xx) 
+0

あなたの答えはありがとうございます。私はグループの数を減らして簡単にしました。 –

0

実装されていないPCH値 '28':

法的形状値番号0〜25、および番号は32

あなたが40台の形状を有しているようだ127へ(0..39 )?

も参照してください

+0

ありがとうRDW。私は、スクリプトをどのように使って、特にx軸とy軸を定義するものがどのように形に関連するのかを理解していませんでした –

関連する問題