Greetings all,
Ran into a strange problem with the krige function from geoR. The
problem that I am having is that while the krige function seems to
work well, the resulting predicted values are all NAs. Given the size
of the datasets I am working with can't attach it, but I can provide
snippets of the datasets.
> casedata
station year month day obs mpe bias type lat lon
140 319147 2011 8 28 0.00 0.000000000 0.00000000 COOP 35.48667 -82.96833
141 319354 2011 8 28 0.02 0.001305799 -0.01869420 COOP 34.25722 -78.68722
142 319357 2011 8 28 0.00 0.045194085 0.04519409 COOP 34.40944 -78.79139
143 319440 2011 8 28 6.92 5.530313974 -1.38968603 COOP 35.85278 -77.03056
144 319461 2011 8 28 0.98 0.949224103 -0.03077590 COOP 34.14250 -77.87861
145 319467 2011 8 28 0.30 1.492852601 1.19285260 COOP 34.32083 -77.92056
146 319476 2011 8 28 2.95 3.284362345 0.33436235 COOP 35.69389 -77.94556
147 319555 2011 8 28 0.00 0.000000000 0.00000000 COOP 36.13083 -81.22750
148 319667 2011 8 28 0.00 0.014472943 0.01447294 COOP 35.85583 -80.36083
149 319675 2011 8 28 0.00 0.000000000 0.00000000 COOP 36.13056 -80.54806
> v.fit = vgm(phi=1.2914,model="Mat",range=14.1611,nugget=0.23,kappa=83.1047)
> krige.pred.out =
> krige(formula=bias~1,locations=~lon+lat,data=casedata,newdata=gridmeta,model=v.fit)
[using ordinary kriging]
> krige.pred.out
lon lat var1.pred var1.var
1 -86.28771 33.91335 NaN NaN
2 -86.27388 33.94722 NaN NaN
3 -86.26003 33.98111 NaN NaN
4 -86.24615 34.01500 NaN NaN
5 -86.23225 34.04890 NaN NaN
6 -86.21835 34.08280 NaN NaN
7 -86.20441 34.11672 NaN NaN
8 -86.19045 34.15064 NaN NaN
9 -86.17648 34.18457 NaN NaN
10 -86.16249 34.21851 NaN NaN
> summary(krige.pred.out)
lon lat var1.pred var1.var
Min. :-86.29 Min. :30.51 Min. : NA Min. : NA
1st Qu.:-82.38 1st Qu.:33.73 1st Qu.: NA 1st Qu.: NA
Median :-79.88 Median :35.22 Median : NA Median : NA
Mean :-79.91 Mean :35.26 Mean :NaN Mean :NaN
3rd Qu.:-77.48 3rd Qu.:36.76 3rd Qu.: NA 3rd Qu.: NA
Max. :-73.16 Max. :40.16 Max. : NA Max. : NA
NA's :44408 NA's :44408
Any ideas are greatly appreciated, I'm not getting any warnings or
errors from krige, so this is perplexing to me.
--
Adrienne Wootten
Graduate Research Assistant
State Climate Office of North Carolina
Department of Marine, Earth and Atmospheric Sciences
North Carolina State University
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
