В чем разница между z~1 и z~LON+LAT при построении вариограммы в R?

library(sp)
data(meuse)
coordinates(meuse)<-~x+y
v<-variogram(log(zinc)~1,meuse)
v
         np       dist     gamma dir.hor dir.ver   id
    1   57   79.29244 0.1234479       0       0 var1
    2  299  163.97367 0.2162185       0       0 var1
    3  419  267.36483 0.3027859       0       0 var1
    4  457  372.73542 0.4121448       0       0 var1
    .
    .

v1<-variogram(log(zinc)~x+y,meuse)
v1
            np       dist     gamma dir.hor dir.ver   id
        1   57   79.29244 0.1060834       0       0 var1
        2  299  163.97367 0.1829983       0       0 var1
        3  419  267.36483 0.2264256       0       0 var1
        4  457  372.73542 0.2847192       0       0 var1
        .
        .

Из приведенного выше кода и вывода я могу обнаружить, что log(zinc)~1 и log(zinc)~x+y возвращают разные значения гаммы для v и v1. В чем принципиальная разница между этими двумя операциями?

help(variogram) сказал мне, что "формула, определяющая вектор отклика и (возможные) регрессоры, в случае отсутствия регрессоров, использует, например, z~1;". Но я не совсем понял это предложение! может ли кто-нибудь подробно рассказать мне, когда я должен использовать z~1 или когда я должен использовать z~LON+LAT?


r sp spatial gstat
person Orpheus    schedule 28.07.2015    source источник

Ответы (1)


arrow_upward
1
arrow_downward

Первый аргумент variogram определяет модель средней структуры. в случае постоянного среднего значения модель содержит только точку пересечения, следовательно, log(zinc)~1. Если среднее значение смоделировано как модель линейной регрессии в x и y, используйте log(zinc)~x+y. В этом случае в качестве основы для значений вариограммы вычисляются обычные остатки методом наименьших квадратов, а не измеренные log(zinc) значения.

person Edzer Pebesma    schedule 13.08.2015