search for: cov1b

Displaying 1 result from an estimated 1 matches for "cov1b".

Did you mean: cov1
2004 Oct 31
3
strange results with dmvnorm
...54,27068,38365,22255,31550,11789,45404,53969 ,13509,36350) Dist <- sqrt(outer(x,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...