Prediction with krigeST
Dear Jun, that is hard to tell without knowing the structure of your "data" and "newdata" objects. I was not able to reproduce your issue with our test cases and would suspect some NA values to cause it. Please note that the latest version of gstat (available from r-forge), does as well support STSDF and STIDF allowing for non regular data sets. In case the problem persists in the latest version of gstat (1.0-18) using a STSDF, please provide me (offlist) with some additional information on your R set-up (sessionInfo()), the structure of your data and a minimal example reproducing the issue. Best wishes, Ben
On 06.11.2013 10:18, Jun HU wrote:
Dear all, I'm new in R and Spatiotemporal kriging, I've a naive question, I've constructed a STFDF object with the first 4 hours historical data, and want to forecast the fifth hour data. By using krigeST in R, I've constructed vgm model. But after launching: pred = krigeST(formula=values~1, data=mystfdf, newdata = My_pred, modelList=myVgm) everything seems working fine, but the pred object only contains historical data, without prediction data. Any idea? Thank you in advance. Sincerely, Jun [[alternative HTML version deleted]]
_______________________________________________ R-sig-Geo mailing list R-sig-Geo at r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo
Benedikt Gr?ler ifgi - Institute for Geoinformatics University of Muenster http://ifgi.uni-muenster.de/graeler Phone: +49 251 83-33082 Mail: ben.graeler at uni-muenster.de -------------- next part -------------- A non-text attachment was scrubbed... Name: signature.asc Type: application/pgp-signature Size: 560 bytes Desc: OpenPGP digital signature URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20131106/baf6abc2/attachment.bin>