Problemas com o kriging comum

8

Eu estava seguindo este artigo relacionado à krigagem comum

insira a descrição da imagem aqui

Agora minha matriz de covariância se parece com isso, para 4 variáveis

1   0.740818220681718   0.548811636094027   0.406569659740599
0.740818220681718   1   0.740818220681718   0.548811636094027
0.548811636094027   0.740818220681718   1   0.740818220681718
0.406569659740599   0.548811636094027   0.740818220681718   1

Bem, a relação entre semvariograma e variograma é dada por

γ(h)/(C0)=1C(h)/C(0)

Então, calculei o também. Agora, quando tento calcular os pesos comoγ(h)

A = 1.0000    0.7408    0.5488    1.0000
    0.7408    1.0000    0.7408    1.0000
    0.5488    0.7408    1.0000    1.0000
    1.0000    1.0000    1.0000         0 

B =  0.4066
    0.5488
    0.7408
    1.0000

Estou tomando a quarta variável como ausente

 [W;mu] = inv(A)*B =  0.1148
                      0.0297
                      0.8555
                     -0.1997

O acima foi usando covariância. Agora usando semi-variância eu tinha

A = 0         0.2592    0.4512    1.0000
    0.2592         0    0.2592    1.0000
    0.4512    0.2592         0    1.0000
    1.0000    1.0000    1.0000         0

B = 0.5934
    0.4512
    0.2592
    1.0000


inv(A)*B =  0.1148
            0.0297
            0.8555
            0.1997

Como você pode ver, os últimos termos não são iguais. Quando de acordo com a derivação, eles são equacionados ou considerados iguais. Algum esclarecimento?

user34790
fonte
Qualquer um caras. Isso está me matando. O que estou fazendo errado?
user34790
Não é uma solução (eu não sabia como postar isso na seção de comentários em um bom formato legível), mas você notou a estrutura do inverso de A nos dois casos diferentes? > A = matriz (c (1.0000,0.7408,0.5488,1.0000, + 0,7408,1.0000,0.7408,1.0000, + 0,5488,0.7408,1.0000,1.0000, + 1.0000,1.0000,1.0000,0), nrow = 4)>> (A) [, 1] [, 2] [, 3] [, 4] [1,] 1.9619812 -1.7076503 -0.2543309 0.4426230 [2,] -1.7076503 3.4153005 -1.7076503 0.1147541 [3,] -0.2543309 -1.7076503 1.9619812 0.4426230 4,] 0.4426230 0.1147541 0.4426230 -0.7705443>>> A = matriz (c (0,0.2592,0.4512,1.0000, + 0.2592,0,0.2592
2
Não há nada na derivação que diga deve ser o mesmo nas formulações de covariância e semivariância. μ
whuber

Respostas:

2

γxx0

μwnγ ΓnxnΓ= [ γ ( x i ,

[Γ110][wμ]=[γ1]
Γn×nγγ= [ γ( x ,Γ=[γ(xi,xj)]i,jγx1nγ=[γ(x,xi)]i envolvendo o novo local e é um vetor daqueles com comprimento .x1n

Consulte (alterações de até notações) Statistics for Spatial Data por N. Cressie p. 121 na edição revisada.

## using the covariance 
Acov <-  matrix(c(1.0000, 0.7408, 0.5488, 1.0000,
                  0.7408, 1.0000, 0.7408, 1.0000,
                  0.5488, 0.7408, 1.0000, 1.0000,
                  1.0000, 1.0000, 1.0000, 0.0000),
                nrow=4) 
Bcov <- c(0.4066, 0.5488, 0.7408, 1.0000)
## using the variogram 
Avario <- matrix(-1, nrow = 4, ncol = 4)
Avario[1:3, 1:3] <- 1 - Acov[1:3, 1:3]
Avario[4, 4] <- 0
Bvario <- 1 - Bcov
Bvario[4] <- -1
## compare
cbind(cov = solve(Acov, Bcov), vario = solve(Avario, Bvario))
Yves
fonte