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@index) == nrow(object@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@metu.edu.tr

 


	[[alternative HTML version deleted]]

