0

私は R に非常に慣れておらず、〜 500k ポイントを含む形状ファイルがランダムに分散されているかどうかを判断しようとしています。引数に何を入力しても、同じエラーメッセージが表示され続けます。

ここに私が最初から持っているものがあります：

> library(spatstat)

spatstat 1.40-0       (nickname: ‘Do The Maths’) 
For an introduction to spatstat, type ‘beginner’
> as.ppp(area)
marked planar point pattern: 500000 points
Mark variables: 
[1] OBJECTID   Encoded_Ti Time_      Filter     Category   Severity      Action_    Hit_Count  Profile    Encoded_So Source_IP  Source_Por Encoded_De
[14] Dest_IP    Dest_Port  VLAN_Tag   Source_Cou Source_Reg Source_Cit     Source_Lat Source_Lon Dest_Count Dest_Regio Dest_City  Dest_Latit Dest_Longi
window: rectangle = [-159.964, 178.417] x [-46.4, 70.6349] units
Warning message:
some mark values are NA in the point pattern x 
> quadrat.test(area)
Error in UseMethod("quadrat.test") : 
no applicable method for 'quadrat.test' applied to an object of class          "c('SpatialPointsDataFrame', 'SpatialPoints', 'Spatial')"


<bytecode: 0x0000000024e7a660>
<environment: namespace:spatstat>

だから私の次の試みは：

> X <- ppp(x, y, c(-159.964, 178.417), c(-46.4, 70.6349))
Warning message:
In ppp(x, y, c(-159.964, 178.417), c(-46.4, 70.6349)) :
data contain duplicated points
> quadrat.test(X)
Error in rectquadrat.countEngine(X$x, X$y, tess$xgrid, tess$ygrid) : 
xbreaks do not span the actual range of x coordinates in data

そして私の最後の試み：

> quadrat.test(X, nx = 20, ny = 20)
Error in rectquadrat.countEngine(X$x, X$y, tess$xgrid, tess$ygrid) : 
xbreaks do not span the actual range of x coordinates in data
> quadrat.test(X, nx = 20, ny = 20, xbreaks= NULL, ybreaks = NULL)
Error in rectquadrat.countEngine(X$x, X$y, tess$xgrid, tess$ygrid) : 
xbreaks do not span the actual range of x coordinates in data
> data(X)
Warning message:
In data(X) : data set ‘X’ not found
> quadrat.test(X)
Error in rectquadrat.countEngine(X$x, X$y, tess$xgrid, tess$ygrid) : 
xbreaks do not span the actual range of x coordinates in data
> quadrat.test(X, 10)
Error in rectquadrat.countEngine(X$x, X$y, tess$xgrid, tess$ygrid) : 
xbreaks do not span the actual range of x coordinates in data

私が言ったように、私はこれに非常に慣れておらず、アマチュアの python の経験しかありませんが、プロジェクトでこの関数を使用する必要がある大学院生です。どんな助けでも大歓迎です。

乾杯

4

1 に答える 1

0

R は (ほとんど) 入力変数を変更しないことに注意してください。そのため、コマンドに出力を割り当てて、そこから作業する必要があります。具体的には、オブジェクトareaに変換できる変数の場合、結果に名前を付ける必要があります。pppas.ppp

X <- as.ppp(area)

次に、新しく作成されたオブジェクトにquadratcountor関数を適用できます。quadrat.testppp

quadratcount(X, nx=20, ny=20)
quadrat.test(X, nx=20, ny=20)

~500k ポイントの 20 x 20 の四角形のグリッドは、四角形のカウントにとって非常に大きな空間スケールのように見えるかもしれませんが、もちろんそれは特定の設定に依存します。

再現可能な例を提供していないため、これらのコマンドがセットアップで機能するかどうかはわかりませんが、適切なデータで機能します.

于 2015-01-27T22:24:23.730 に答える