|
In applied statistics, regression-kriging (RK) is a spatial prediction technique that combines a regression of the dependent variable on auxiliary variables (such as parameters derived from digital elevation modelling, remote sensing/imagery, and thematic maps) with kriging of the regression residuals. It is mathematically equivalent to the interpolation method variously called ''universal kriging'' and ''kriging with external drift'', where auxiliary predictors are used directly to solve the kriging weights. == BLUP for spatial data == Regression-kriging is an implementation of the best unbiased linear predictor for spatial data, i.e. the best linear interpolator assuming the universal model of spatial variation. Matheron (1969) proposed that a value of a target variable at some location can be modeled as a sum of the deterministic and stochastic components: : which he termed ''universal model of spatial variation''. Both deterministic and stochastic components of spatial variation can be modeled separately. By combining the two approaches, we obtain: : where is the fitted deterministic part, is the interpolated residual, are estimated deterministic model coefficients ( is the estimated intercept), are kriging weights determined by the spatial dependence structure of the residual and where is the residual at location where is the vector of estimated regression coefficients, is the covariance matrix of the residuals, is the vector of measured values of the target variable. The GLS estimation of regression coefficients is, in fact, a special case of the geographically weighted regression. In the case, the weights are determined objectively to account for the spatial auto-correlation between the residuals. Once the deterministic part of variation has been estimated (regression-part), the residual can be interpolated with kriging and added to the estimated trend. The estimation of the residuals is an iterative process: first the deterministic part of variation is estimated using OLS, then the covariance function of the residuals is used to obtain the GLS coefficients. Next, these are used to re-compute the residuals, from which an updated covariance function is computed, and so on. Although this is by many geostatisticians recommended as the proper procedure, Kitanidis (1994) showed that use of the covariance function derived from the OLS residuals (i.e. a single iteration) is often satisfactory, because it is not different enough from the function derived after several iterations; i.e. it does not affect much the final predictions. Minasny and McBratney (2007) report similar results—it seems that using more higher quality data is more important then to use more sophisticated statistical methods. In matrix notation, regression-kriging is commonly written as: : where , is the vector of predictors and where is the sill variation and ) becomes an identity matrix. Likewise, if the target variable shows no correlation with the auxiliary predictors, the regression-kriging model reduces to ordinary kriging model because the deterministic part equals the (global) mean value. Hence, pure kriging and pure regression should be considered as only special cases of regression-kriging (see figure). 抄文引用元・出典: フリー百科事典『 ウィキペディア(Wikipedia)』 ■ウィキペディアで「Regression-kriging」の詳細全文を読む スポンサード リンク
|