R/coregionalization_model_gauss.R

In tleblevecIMP/TransPGS: Study the indicator transiograms obtained from multi-gaussian integrals

#' this function computes the cross correlations for a distance h
#' using a gaussian covariance


coregionalization_model_gauss<-function(h,r1,r2,rho,shift,rho1a){

  rho1h=gaussian_cov(r1,h)
  rho2h = ((rho^2)/(rho1a^2)) * rho1h + (1-((rho^2)/(rho1a^2))) * gaussian_cov(r2,h)
  rho12h =  rho*gaussian_cov(r1,h+shift) / rho1a
  rho21h =  rho*gaussian_cov(r1,h-shift) / rho1a

  return (c(rho1h,rho2h,rho12h,rho21h))
}