地理空間加重回帰分析 (GWR) (空間統計ツール)
サマリー
空間的に変化する関係をモデリングするときに使用される局所形の線形回帰である、地理空間加重回帰分析を実行します。
図
使用法
このツールは、空間的に変化する関係をモデリングするときに使用される局所形の線形回帰である、地理空間加重回帰分析を実行します。 GWR ツールは、データセット内のあらゆるフィーチャに回帰方程式をあてはめることで、理解または予測しようとしている変数またはプロセスのローカル モデルを作成します。 GWR ツールでは、各ターゲット フィーチャの近傍内に存在するフィーチャの従属変数と説明変数を組み込んで、このような別々の方程式を作成します。 各近傍の形状と範囲は、近傍タイプ と 近傍選択方法 パラメーターに対する入力に基づきますが、これには 1 つの制約があります。近隣のフィーチャ数が 1000 を超えると、各局所方程式には最近隣の 1000 個だけが組み込まれます。
このツールは、数百のフィーチャを持つデータセットに対して使用します。 このツールは、小規模なデータセットには適していません。 このツールをマルチポイント データに使用することはできません。
入力フィーチャ パラメーターは、モデル化している現象 (従属変数 の値) を表すフィールドと 説明変数 の値を表す 1 つ以上のフィールドとともに使用します。 これらのフィールドは数値で、値の範囲を含む必要があります。 従属変数または説明変数に欠損値があるフィーチャは、解析から除外されます。ただし、欠損値の補完 ツールを使用して、GWR ツールを実行する前にデータセットを完全なものにすることができます。
-
注意:
GWR ツールでは、各種出力を生成できます。 地理空間加重回帰分析モデルのサマリーは、ツールの操作中に ジオプロセシング ウィンドウの下部でメッセージとして利用できます。 このメッセージにアクセスするには、ジオプロセシング ウィンドウで、進行状況バーにポインターを合わせるか、ポップアップ ボタンをクリックするか、メッセージ セクションを展開します。 ジオプロセシング履歴から、以前に実行した GWR ツールのメッセージにアクセスすることもできます。 また、GWR ツールは、出力フィーチャ の値を生成し、ローカルな診断値を報告するフィールドを追加します。 出力フィーチャ の値と関連するチャートは、残余をモデリングできるよう、ホット/コールド レンダリング方式が適用された状態で、コンテンツ ウィンドウに自動的に追加されます。 各出力およびチャートの詳細な説明については、地理空間加重回帰分析の詳細をご参照ください。
The tool accepts points and polygons as input. For polygons, all distances and neighbors are determined using the distance between polygon centroids (points). However—特に大きいポリゴン、細長いポリゴン、またはマルチパートのポリゴンの場合は、単一ポイントは適切にポリゴンを表さないことがあります。 このような場合、ポリゴン間の近傍数と距離が不自然になるか、誤解を招く可能性があります。 たとえば、境界を共有する 2 つのポリゴンは、その重心が離れていると、近傍と見なされない場合があります。 このツールが使用する重心を確認するには、インサイド パラメーターをオフにして フィーチャ → ポイント ツールを使用し、ポリゴンを重心ポイントに変換します。 また、近隣探索を使用して、ポリゴンまたはポイントの重心の近傍を視覚化することもできます。
通常、ラインに対して地理空間加重回帰分析を実行することはお勧めしません。重心がラインを適切に表現することはほとんどないためです。 ただし、ツールでラインを使用するために、フィーチャ → ポイント ツールを使用してラインを重心ポイントに変換し、ツールで重心を使用することができます。 その後、結果を元のラインに結合できます。
指定する モデル タイプ パラメーターの値は、モデル化しているデータによって異なります。 回帰分析の正確な結果を得るには、解析に正しいモデルを使用することが重要です。
投影データを使用することをお勧めします。 これは、距離が分析の要素になっている場合には、近傍タイプ パラメーターに 距離バンド を指定すると、地理空間加重回帰分析にとっても距離が成分となるため、特に重要です。 データの投影法には、(地理座標系ではなく) 投影座標系を使用することをおすすめします。
計算の中には、パフォーマンス向上のために複数の CPU を使用できるものもあり、この場合は自動的に最大 8 つのスレッド/CPU が処理に使用されます。
このツールを使用してデータをローカルに調査する前に、一般化線形回帰 ツールを使用してデータをグローバルに調査する方法が一般的です。
従属変数 および 説明変数 パラメーター値は、さまざまな値を含む数値フィールドである必要があります。 これらの値には、グローバルおよびローカルなばらつきがあります。 このため、地理空間加重回帰分析モデルでは、別の空間様式を表すダミーの説明変数を使用しないでください (都心外の国勢調査地区に 1 という値を割り当て、それ以外にはすべて 0 という値を割り当てるなど)。 GWR ツールでは、説明変数の係数が異なっていてもよいため、このような空間様式の説明変数が不要です。これを含めてしまうと、局所的な多重共線性の問題が生じます。
一般化線形回帰のようなグローバル回帰モデルでは、複数の変数が多重共線性を示すと (複数の変数が冗長であるか、同じ「ストーリー」を語るとき)、結果は信頼できません。 GWR ツールでは、データセットの各フィーチャについて、ローカル回帰方程式が構築されます。 特定の説明変数の値が空間的にクラスタリングしている場合は、ローカル多重共線性が問題となっている可能性が高いといえます。 ローカル多重共線性のために不安定な結果になっているかどうかは、出力フィーチャクラスの条件数フィールド (
COND) で確認できます。 一般的に、条件数が 30 を上回っている場合、NULL に等しい場合、またはシェープ ファイルの場合であれば -1.7976931348623158e+308 に等しい場合の結果には信頼性がありません。 条件数は、モデル内の説明変数の数を補正するようにスケール調整されます。 これによって、異なる数の説明変数を使用するモデル間で条件数を直接比較できるようにします。地理空間加重回帰分析モデルに名義またはカテゴリー データを含める場合には、注意してください。 カテゴリーが空間的にクラスタリングしているようなときには、ローカル多重共線性の問題が生じるリスクがあります。 ローカル共線性が問題となっているかどうかは、地理空間加重回帰分析出力に含まれている条件数からわかります (条件数が 0 未満、30 超、または NULL に設定されている場合)。 ローカル多重共線性が存在する場合の結果は安定していません。
説明変数の係数にみられる地域的なばらつきについて理解を深めるために、GWR ツールによって作成されるオプションのラスター係数サーフェスを確認します。 これらのラスター サーフェスは、その他のオプション の下にある 係数ラスター ワークスペース パラメーターを指定していれば、この中に作成されます。 ポリゴン データの場合は、出力フィーチャ の値の各係数フィールドにグラデーションの色やコールドからホットまでのレンダリングを使用して、分析範囲全体における変化を確認できます。
GWR ツールは、予測位置 の値のフィーチャクラス (多くの場合、このフィーチャクラスは 入力フィーチャ の値と同じです) を指定し、説明変数を照合して、出力予測済みフィーチャ の値を指定することで、予測を行う場合にも使用できます。 入力フィーチャ 値の 照合する説明変数 フィールドが 予測位置からのフィールド フィールドと一致する場合、それらは自動的に設定されます。 一致しない場合は、正しいフィールドを指定します。
重要な説明変数が欠けていると、回帰モデルは誤って指定されます。 回帰残余の空間的自己相関が統計的に有意である場合や、1 つ以上の説明変数の係数に予期せぬ空間的なばらつきがある場合は、モデルが誤って指定されている可能性が考えられます。 これらの欠けている重要な変数を特定し、これらをモデルに含められるよう、(一般化線形回帰残余分析や地理空間加重回帰分析係数差異分析などを通じて) あらゆる努力を行ってください。
説明変数を非定常とすることが理にかなっているかどうかを判断します。 たとえば、ASPECT をはじめとするいくつかの変数を関数として、特定の植物種の密度をモデリングしているとします。 ASPECT 変数の係数が分析範囲全体で変化する場合には、重要な説明変数 (優勢な競合する植物など) が欠けているという事実を示している可能性が高いといえます。 あらゆる努力を行って、回帰モデルにすべての重要な説明変数を含めてください。
計算結果が無限または未定義の場合、シェープファイル以外の結果は NULL に、シェープファイルの結果は -DBL_MAX = -1.7976931348623158e+308 になります。
注意:
シェープファイルは NULL 値を格納できません。 シェープファイルではない入力からシェープファイルを作成するツールやその他の手順では、結果的に NULL 値がゼロまたは非常に小さな負の値 (-DBL_MAX = -1.7976931348623158e+308) として格納されることがあります。 この場合、予期せぬ結果に至ることがあります。 詳細については、ジオプロセシングでのシェープファイル出力の注意事項をご参照ください。
近傍選択方法 パラメーターには、次の 3 つのオプションがあります。 黄金探索 を指定すると、ツールは黄金分割探索方法を使用して、距離バンド または 近傍数 パラメーターに最適な値を見つけます。 手動間隔 オプションは、距離バンドまたは近傍数を少しずつ増やして、指定された距離間の近傍をテストします。 どちらの場合も、使用される近傍サイズは、赤池情報量基準 (AICc) の値を最小化するサイズです。 しかしながら、どちらの方法も、局所的な多重共線性の問題のために、最適な距離バンドまたは近傍数の解を求めることができません。 エラーが表示されたり深刻なモデル設計上の問題が発生した場合は、ユーザー定義 オプションを使用して、特定の距離または近傍数を指定してみてください。 次に、出力フィーチャクラスの条件数を確認し、どのフィーチャが局所的な共線性の問題に関係しているのか判断してください。
深刻なモデル設計上の問題や、局所方程式に十分な近傍数が含まれていないことを示すエラーは、多くの場合、グローバルまたは局所的な多重共線性の問題があることを示しています。 どこに問題があるのかを判断するために 一般化線形回帰 ツールを使用してグローバル モデルを実行し、各説明変数について VIF 値を確認します。 大きな VIF 値がある場合 (たとえば 7.5 以上)、グローバル多重共線性によって地理空間加重回帰分析では解を導けなくなっています。 しかし、もっと可能性が高いのは、ローカル多重共線性の問題です。 各独立変数の主題図を作成してみてください。 主題図で同一値が空間的にクラスターを形成している場合は、モデルからそのような変数を除外するか、別の独立変数と組み合わせて値のばらつきを増やします。 たとえば、住宅の価値をモデル化していて、寝室と風呂に変数がある場合は、両方を組み合わせて値のばらつきを増やすか、それぞれを風呂と寝室の面積として表します。 地理空間加重回帰分析モデルの構築時には、空間様式のダミー変数、空間クラスターとなっているカテゴリーや名義変数、可能性のある値が非常に限られている変数の使用を避けてください。
地理空間加重回帰分析は、一般化線形回帰と同じ要件の適用を受ける線形モデルです。 地理空間加重回帰分析モデルが適切に指定できていることを確認するには、地理空間加重回帰分析の仕組みで説明されている診断を確認してください。 「回帰分析の基礎」の回帰モデルがうまくいかない理由にも、モデルが正確であることを確認するための情報が記載されています。
パラメーター
| ラベル | 説明 | データ タイプ |
|---|---|---|
|
入力フィーチャ |
従属変数と説明変数を格納したフィーチャクラス。 |
Feature Layer |
|
従属変数 |
モデル化される観測値を含む数値フィールド。 |
Field |
|
モデル タイプ |
モデル化されるデータのタイプを指定します。
|
String |
|
説明変数 |
回帰モデルの独立説明変数を表すフィールドのリスト。 |
Field |
|
出力フィーチャ |
従属変数の推定値と残差が含まれている新しいフィーチャクラス。 |
Feature Class |
|
近傍タイプ |
使用される近傍を固定距離として構築するか、フィーチャの密度に応じて空間範囲の変化を許容するかを指定します。
|
String |
|
近傍選択方法 |
近傍サイズの決定方法を指定します。 黄金探索 および 手動間隔 オプションを使用して選択された近傍は、AICc の値を最小化する処理に基づきます。
|
String |
|
最小近傍数 (オプション) |
各フィーチャの計算に含まれる最小の近傍数。 少なくとも 30 の近傍を使用することをお勧めします。 |
Long |
|
最大近傍数 (オプション) |
各フィーチャの計算に含まれる最大の近傍数 (最大 1000)。 |
Long |
|
最小検索距離 (オプション) |
最小近傍検索距離。 各フィーチャに少なくとも 30 の近傍がある距離を使用することをお勧めします。 |
Linear Unit |
|
最大検索距離 (オプション) |
最大近傍検索距離。 フィーチャの近傍が 1000 を超える距離の場合、ツールはターゲット フィーチャの計算に最初の 1000 を使用します。 |
Linear Unit |
|
近傍数の増分 (オプション) |
近傍テストごとに、この近傍数だけ手動間隔が増やされます。 |
Long |
|
検索距離の増分 (オプション) |
近傍テストごとに、この距離だけ手動間隔が増やされます。 |
Linear Unit |
|
増分数 (オプション) |
最小近傍数 または 最小検索距離 パラメーター値から始まる、テストする近傍サイズの数。 |
Long |
|
近傍数 (オプション) |
各フィーチャについて考慮する近隣フィーチャの最も近い数 (最大 1000)。 この数は 2 ～ 1000 の整数の必要があります。 |
Long |
|
距離バンド (オプション) |
近傍の空間範囲。 |
Linear Unit |
|
予測位置 (オプション) |
推定を計算する場合に、位置を表すフィーチャが格納されるフィーチャクラス。 このデータセット内の各フィーチャには、指定されたすべての説明変数の値が含まれています。 これらのフィーチャの従属変数は、入力フィーチャクラス データに対してキャリブレーションされたモデルを使用して推定されます。 予測を行うには、これらのフィーチャの位置が 入力フィーチャ パラメーターの値と同じ分析範囲内か、近く (範囲プラス 15 パーセント内) の必要があります。 |
Feature Layer |
|
照合する説明変数 (オプション) |
予測位置 パラメーターの説明変数を 入力フィーチャ パラメーターの対応する説明変数と照合します。 値テーブル列:
|
Value Table |
|
出力予測済みフィーチャ (オプション) |
各 予測位置 値の従属変数の推定値を受け取る出力フィーチャクラス。 |
Feature Class |
|
ロバスト予測 (オプション) |
予測計算に使用されるフィーチャを指定します。
|
Boolean |
|
ローカル加重方式 (オプション) |
モデルの空間的な加重を提供するために使用されるカーネル タイプを指定します。 カーネルは、各フィーチャと近傍内にある他のフィーチャとの関連性を定義します。
|
String |
|
係数ラスター ワークスペース (オプション) |
係数ラスターが作成されることになるワークスペース。 ワークスペースを設けると、インターセプトとすべての説明変数についてラスターが作成されます。 このパラメーターは Desktop Advanced ライセンスでのみ利用可能です。 |
Workspace |
|
データのスケール (オプション) |
モデルをフィッティングする前に、平均が 0、標準偏差が 1 になるように説明変数と従属変数の値をスケーリングするかどうかを指定します。
|
Boolean |
派生した出力
| ラベル | 説明 | データ タイプ |
|---|---|---|
|
係数ラスター レイヤー |
出力係数ラスター。 |
Raster Layer |
環境
出力座標系, 地理座標系変換, 現在のワークスペース, スクラッチ ワークスペース, セル サイズ, スナップ対象ラスター
特殊なケース
- 出力座標系
-
Feature geometry is projected to the output coordinate system after analysis is complete.
ライセンス情報
- Basic: 限定的
係数ラスターの生成には Advanced ライセンスが必要です。 - Standard: 限定的
係数ラスターの生成には Advanced ライセンスが必要です。 - Advanced: はい