cokrig: Cokriging
In victorvicpal/MGSR: Multivariate Gaussian Subspatial Regression

Description Usage Arguments Value Author(s) References See Also Examples

Subspatial Cokriging estimation.

1	cokrig(RES, xygrid)

`RES`	Resulting lmc fitting
`xygrid`	Grid

Dataframe containing coordinates and their estimated values.

Victor Vicente Palacios

Myers D.E. Matrix Formulation of Co-Kriging. Math Geol. 1982;14(3):249-257.

lmc,GRID_MGSR

data("iris")
Versicolor <- iris[which(iris$Species=='versicolor'),-5]
##Data Standarization
means_vers <- apply(Versicolor,2,mean)
sd_vers <- apply(Versicolor,2,sd)
Versicolor_st <- Versicolor

for (i in 1:length(Versicolor[1,]))
{Versicolor_st[,i] <- (Versicolor[,i]-means_vers[i])/sd_vers[i]}

##PrComp
PC_train <- princomp(Versicolor_st)
## CrossVariogram Calculation
CV_vers <- crossvariogram(as.data.frame(PC_train$scores[,1:2]),as.data.frame(Versicolor_st),11)
## lmc 
RES_vers <- lmc(CV_vers,'Pow',1.6)

## Cokriging (standarized values)
Z_vers_st <- cokrig(RES_vers,xygrid)