2016-11-09 9 views
1

です。私はpcaとggbiplotを使用してpca結果を表示しようとしていましたが、何とか解決できなかったいくつかのエラーが発生しました。コードが他のデータとうまく機能するため、おそらくデータに問題があります。 は、私は、コードと、私はあなたが以下のリンクでシナリオを再作成したい場合に使用するデータファイルを置く: -Rca ggbiplotエラー:置換は36行、データは35

https://drive.google.com/drive/folders/0B2jQ7Vh3S3PaZkt3Y2ZyaV9XaXc

コード:PCA-plot.R データファイル1:DAT1を。 rda(これは正常に動作します) データファイル2:dat2.rda(これには問題があります)

助けをお待ちしています。 私が得たエラーは一番下にあります。

dat2.rdaからあなたのdfall --we

> g <- ggbiplot(tr.pca, obs.scale = 1, var.scale = 1, 
+    groups = Ydfall, 
+    ellipse = TRUE, 
+    circle = TRUE) 
Error in `$<-.data.frame`(`*tmp*`, "groups", value = c(1L, 1L, 1L, 1L, : 
    replacement has 36 rows, data has 35 
> g <- g + scale_color_discrete(name = '') 
Scale for 'colour' is already present. Adding another scale for 'colour',  which will replace the existing scale. 
> g <- g + theme(legend.direction = 'horizontal', 
+    legend.position = 'top') 
> g<- g+ geom_point(size=1, shape=1, color="black", stroke=2) 
> 
> print(g) 
> 

答えて

0

、ありがとうNAを持っている(which(is.na(dfall), arr.ind = T)を試してみてください)、それがあなたの問題が発生します。 prcomp()を使用したときにna.omit()を使用しましたが、Ydfallを作成したときは使用しませんでした。

Ydfall <- na.omit(dfall)[,1] # quick fix 

# but if I were you, I would do first; 
dfall <- na.omit(dfall) 
関連する問題