데이터 프레임 열이 제거 된 경우 autoKrige이 다른 결과를 표시합니다.

여기 뭔가 잘못되었을 수도 있지만 관련없는 열을 제거하여 데이터 프레임을 단순화하면 automap 라이브러리의 autoKrige 함수가 다른 결과를 제공합니다. automap 라이브러리의 meuse 데이터로이 문제를 재연 해 보았습니다.데이터 프레임 열이 제거 된 경우 autoKrige이 다른 결과를 표시합니다.

library(automap) 
data(meuse) 
colnames(meuse) 
[1] "x"  "y"  "cadmium" "copper" "lead" "zinc" "elev" 
[8] "dist" "om"  "ffreq" "soil" "lime" "landuse" "dist.m" 

coordinates(meuse) =~ x+y 
data(meuse.grid) 
gridded(meuse.grid) =~ x+y 
kriging_result_01 = autoKrige(zinc~1, meuse) 
plot(kriging_result_01) 

meuse <- NULL 
data(meuse) 
meuse <- meuse[, c(1,2,6)] 

coordinates(meuse) =~ x+y 
data(meuse.grid) 
gridded(meuse.grid) =~ x+y 
kriging_result_02 = autoKrige(zinc~1, meuse) 
plot(kriging_result_02) 

identical(kriging_result_01, kriging_result_02) 
[1] FALSE

플롯의 세부 사항도 다릅니다.

이것은 예상되는 동작입니까?

덕분에, 빌

출처

2014-10-01 William Chivers

필자는 세부적인 이유를 찾지 못했지만 'autoKrige'에 제공되지 않을 때 new_data 객체가 어떻게 계산되는지와 관련이 있다고 생각합니다. 'autoKrige (zinc ~ 1, meuse, meuse.grid)'를 사용하면 동일한 결과를 얻을 수 있습니다. –

'kriging_result_01 $ krige_output @ grid'와'kriging_result_02 $ krige_output @ grid'를 보면 계산 된'new_data' 객체가 약간 다른 것을 볼 수 있습니다. –

문제가되지 제거 된 열이지만, 그리드에 장착 모델을 적용 할 수 있습니다. 귀하의 예제를 사용하면 볼 수 있습니다 : 당신이 볼 수 있듯이

library(automap) 
set.seed(42) 

data(meuse) 
coordinates(meuse) =~ x+y 
data(meuse.grid) 
gridded(meuse.grid) =~ x+y 
kriging_result_01 = autoKrige(zinc~1, meuse) 

meuse <- NULL 
data(meuse) 
meuse <- meuse[, c(1,2,6)] 
coordinates(meuse) =~ x+y 

kriging_result_02 = autoKrige(zinc~1, meuse) 


[email protected] 
          x1   x2 
cellcentre.offset 178635.12844 329744.86186 
cellsize    32.93423  32.93423 
cells.dim    84.00000 118.00000

및

[email protected] 
          x1   x2 
cellcentre.offset 178614.42379 329741.35016 
cellsize    32.93423  32.93423 
cells.dim    85.00000 118.00000

는, 그리드가 약간 다릅니다, 4999 및 4996 격자 점으로 각각. 당신이

kriging_result_01 = autoKrige(zinc~1, meuse, meuse.grid) 
kriging_result_02 = autoKrige(zinc~1, meuse, meuse.grid) 

[email protected] 
         x  y 
cellcentre.offset 178460 329620 
cellsize    40  40 
cells.dim    78 104 

[email protected] 
         x  y 
cellcentre.offset 178460 329620 
cellsize    40  40 
cells.dim    78 104

를 사용하는 경우

은 그리드는 비슷합니다.

출처

2014-10-01 05:51:14

Brilliant, thanks 파스칼! –

@WilliamChivers이 답변으로 질문을 해결 한 경우 점수 아래에있는 회색 눈금을 클릭하여 정답으로 받아 들여야합니다. –

데이터 프레임 열이 제거 된 경우 autoKrige이 다른 결과를 표시합니다.

답변

관련 문제