Dear list:

I'm trying to fit a spatio-temporal variogram model in gstat.

I'm reading: http://cran.r-project.org/web/packages/gstat/vignettes/st.pdf

My dataset is: CO, 23 Monitoring stations (irregular), 24 hours.
I create a STFDF object and it looks fine, but when I try to create
a variogramST, I found an error. (I'm attaching the code)

I have read a previous discussion about it but I cann't solve my error.

Please, could you help me?

Thanks in advance,
Gema



##--------------------------------------------------------------------------------------##

> sp = cbind(coords[,1], coords[,2])
> row.names(sp) = paste("point", 1:nrow(sp), sep="")
> sp = SpatialPoints(sp)

> time = as.POSIXct("2008-01-01")+3600*(00:23)
> #select 24 hours
> dia.co=co_stl_23[1:24,]
> 
> mydata = (length(sp)*length(time),dia.co)
> mydata = data.frame(values = signif(mydata,3), ID=IDs)
> stfdf.co = STFDF(sp, time, mydata)
> 
> summary(stfdf.co)

Object of class STFDF
 with Dimensions (s, t, attr): (23, 24, 2)
[[Spatial:]]
Object of class SpatialPoints
Coordinates:
                min     max
coords.x1  436436.6  448850
coords.x2 4470109.5 4481161
Is projected: NA 
proj4string : [NA]
Number of points: 23
[[Temporal:]]
     Index                       timeIndex    
 Min.   :2008-01-01 00:00:00   Min.   : 1.00  
 1st Qu.:2008-01-01 05:45:00   1st Qu.: 6.75  
 Median :2008-01-01 11:30:00   Median :12.50  
 Mean   :2008-01-01 11:30:00   Mean   :12.50  
 3rd Qu.:2008-01-01 17:15:00   3rd Qu.:18.25  
 Max.   :2008-01-01 23:00:00   Max.   :24.00  
[[Data attributes:]]
     values               ID     
 Min.   :-3.12000   ID 1   : 46  
 1st Qu.:-0.75150   ID 10  : 46  
 Median :-0.01895   ID 11  : 46  
 Mean   : 0.03365   ID 12  : 46  
 3rd Qu.: 0.75600   ID 2   : 46  
 Max.   : 2.94000   ID 3   : 46  
                    (Other):276  
> #plot
> stplot(stfdf.co)

> st.var<- variogramST(values ~ 1, data=stfdf.co, tlags=0:5)

Error en if (nrow(x) < nrow(y)) matrix(t(apply(x, 1, spDiN1, y = y, ll = longlat)),  : 
  argumento tiene longitud cero


##--------------------------------------------------------------------------------------##




Gema Fernández-Avilés 

Department of Statistics
Faculty of Law and Social Sciences
University of Castilla-La Mancha at Toledo

http://www.uclm.es/profesorado/gemafaviles/
New book: http://eu.wiley.com/WileyCDA/WileyTitle/productCd-1118413180.html

	[[alternative HTML version deleted]]

