0

私はRプログラミングが初めてです。私は現在、spdep パッケージの mstree() を改善できる研究プロジェクトに取り組んでいます。Polygon()、Polygons()、SpatialPolygons()、SpationPolygonsDataFrame() を機能させる方法について質問があります。以下は私が試したコードです:

 library("spdep")
 library("sp")
 require(maptools)
 Dmat<-matrix(nrow=1,ncol=2)
 rownames(Dmat)<-c("s1")
 Dfram<-data.frame(Dmat)

 #diabetesNum is a dataset that I artificially made
 nrows<-nrow(diabetesNum)
 x1matrix<-data.matrix(diabetesNum)
 x1Frame<-data.frame(diabetesNum)
 x1Frame4<-as(x1Frame,"S4")
 x1X<-prcomp(x1Frame)$x

 #Use the first two principal components for coordinates of the points
 xycoordinates<-data.frame(x1X[,1:2])
 row1<-xycoordinates[1,]
 row1M<-data.frame(row1)
 xycoordinates[nrows+1,]<-c(row1M)
 xycoordinates
            PC1         PC2
  1   285.60120   -6.541733
  2   -78.14211  -30.160212
  3  -263.15002   31.236614
  4    51.00528 -143.520521   
  5  -139.55849   27.842465
  6  -244.09755   31.258001
  7   281.18411   12.072960
  8   580.64105   51.538605
  9  -135.05802  -38.923404
 10 -201.08445   47.932800
 11 -137.34100   17.264425
 12  285.60120   -6.541733

 xycoordinatesF<-data.frame(xycoordinates)
 xycoordinatesM<-data.matrix(xycoordinatesF)
 pl<-Polygon(xycoordinatesM)
 pls<-Polygons(list(pl), ID="1")
 Sr<-SpatialPolygons(list(pls))

 bho.nb<-poly2nb(Sr)
 Error in sp$mbxv[i]:(n * 2) : argument of length 0

...そしてエラー!上記のコードの何が問題になっていますか? Polygon()、Polygons()、SpatialPolygons()、SpatialPolygonsDataFrame()、および poly2nb() を使用する適切な方法を見つけるためにインターネットで検索して検索しましたが、これらの関数を機能させる方法がわかりません! 私はすでに R のマニュアルを読んでいますが、あまり役に立ちません。しかし、興味深いことの 1 つは、Polygons オブジェクトを表示すると、Polygons スロットが空であることです。同様に Sr を表示すると、SpatialPolygons スロットと Polygons スロットは空です。ありがとう、

4

0 に答える 0