[R-sig-Geo] space-time kriging

Edzer Pebesma edzer.pebesma at uni-muenster.de
Fri Jul 1 13:50:26 CEST 2011


Hi Pinar,

as you say, some data are missing, but the way you create index suggests
none are missing. Does your table "prec" have information for each
observation on which location and at which time it has been taken? Could
you maybe show how prec[1:10,] looks like?

Wbr,

On 07/01/2011 01:12 PM, Pınar Aslantas Bostan wrote:
> Dear all,
> 
>  
> 
> I am working about space-time (ST) kriging. I want to make ST ordinary
> kriging with sparse data set. Dependent data is annual precipitation
> measured from 257 meteorological station. Observation period starts from
> 1970 and finishes at 2008, so temporal entity of my study contains 39
> year-value. Precipitation data comprised from 5365 value. It is not equal to
> 257 (location) *39 (time), because some stations has no value for some
> years. Therefore my data set fits to STSDF. An index table is required for
> STSDF object. But I could not create correct index object, my code does not
> work. 
> 
>  
> 
> I tried to summarize my code:
> 
>  
> 
> # creating STSDF object
> 
>  
> 
> # sp is spatial entity (257 spatial location)
> 
> # time is temporal entity (39 year period)
> 
> # prec is data attribute (5365 number of annual precipitation value)
> 
>  
> 
> index <- as.matrix(expand.grid(1:length(sp),1:nrow(time)))
> 
> w= STSDF(sp, time, data.frame(values = as.vector(t(prec))),index)
> 
>  
> 
> Error: nrow(object at index) == nrow(object at data) is not TRUE
> 
>  
> 
> Can you please help me about this error, how can I create correct index
> object?
> 
>  
> 
> Thank you all in advance!
> 
> With my best regards,
> 
> Pinar
> 
>  
> 
>  
> 
> Pinar Aslantas Bostan
> 
> Research Assistant
> 
> Department of Geodetic and
> 
> Geographic Information Technologies (GGIT) Middle East Technical University
> 
> 06531 Ankara/TURKEY
> 
> aslantas at metu.edu.tr
> 
>  
> 
> 
> 	[[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

-- 
Edzer Pebesma
Institute for Geoinformatics (ifgi), University of Münster
Weseler Straße 253, 48151 Münster, Germany. Phone: +49 251
8333081, Fax: +49 251 8339763  http://ifgi.uni-muenster.de
http://www.52north.org/geostatistics      e.pebesma at wwu.de



More information about the R-sig-Geo mailing list