2

R Kest 関数 (spatstat) で検索距離 "r" がどのように計算されるかがわかりません。パッケージのドキュメントには次のように記載されています。

オプション。K(r) を評価する引数 r の値のベクトル。この引数は指定しないことをお勧めします。賢明なデフォルトがあります。

計算で使用される「賢明なデフォルト」とは何ですか? どんなドキュメントでも大歓迎です。

4

2 に答える 2

5

Rオープンソースなので、自分でコードを読むことができます。入力Kestすると、コードが一覧表示されます。が割り当てられている場所を見つけるrと、 というサブルーチンによって行われていることがわかりhandle.r.b.argsます。プロンプトでその名前をR入力して手順を読みます。わずか 6 行です。

    rmax <- if (missing(rmaxdefault)) 
        diameter(as.rectangle(window))
    else rmaxdefault
    if (is.null(eps)) {
        if (!is.null(window$xstep)) 
            eps <- window$xstep/4
        else eps <- rmax/512
    }
    breaks <- make.even.breaks(rmax, bstep = eps)

明らかに、領域の直径の推定値は 512 で割られます。

于 2013-01-15T02:53:40.583 に答える
3

@whuber による非常に役立つ回答をフォローアップすると、Kest で r に使用される正確な計算に関する追加の詳細を見つけることができました。

Kest(X) で r 値はどのように決定されますか?

Kest のデフォルトの r 値は次のように計算されます。

r の最大値は、関数 'rmax.rule' によって計算されます。rmax <- rmax.rule("K", W, lambda) ここで、W はデータを含むウィンドウで、lambda は単位面積あたりのポイントの平均密度です。 . 現在、このルールは、Ripley の経験則の最小値を取ります: rmax = 囲んでいる長方形の最小辺の 4 分の 1 大きなサンプル ルール: rmax = sqrt(1000/(pi * lambda)) r 値は 0 から rmax まで等間隔です。ステップ値「eps」。eps が指定されていない場合は、513 個の値または 512 個の間隔があるように、eps = rmax/512 になります。必要に応じて、いつでも「r」値をオーバーライドできます。

出典：spatstat

于 2013-01-15T11:50:25.790 に答える