search for: cov2b

...x, "-")^2 + outer(y,y, "-")^2) Dist <- Dist/max(Dist) library(spatial) Cov1 <- sphercov(Dist, 0.8, alpha=0, se=sqrt(2)) Cov2 <- sphercov(Dist, 0.6, alpha=0, se=sqrt(0.55)) library(geoR) Cov1b <- cov.spatial(Dist, cov.model= "spherical", cov.pars=c(2, 0.8)) Cov2b <- cov.spatial(Dist, cov.model= "spherical", cov.pars=c(0.55, 0.6)) library(mvtnorm) dmvnorm(rep(0, nrow(Cov)), sigma=Cov1) [1] 37949.22 dmvnorm(rep(0, nrow(Cov)), sigma=Cov1b) [1] 2.920084e-05 dmvnorm(rep(0, nrow(Cov)), sigma=Cov2) [1] 4.722488e+60 dmvnorm(rep(0, nrow(Cov)), sigma=Co...