Check number pairs CV

victorvicpal · Jan 30, 2017 · 7b1560c · 7b1560c
1 parent d9bf3ba
commit 7b1560c
Showing 1 changed file with 14 additions and 0 deletions.
diff --git a/R/n_pairs_opt.R b/R/n_pairs_opt.R
@@ -0,0 +1,14 @@
+n_pairs_opt <- function(coord,values,n_min=12,n_max=20)
+{
+	v <- rep(0,(n_max-n_min+1))
+	n <- 1
+	for (i in n_min:n_max)
+	{
+		CV <- crossvariogram(coord,values,i)
+		aux <- CV[[1]][,3]
+		v[n] <- max(aux)-min(aux)
+		n <- n+1
+	}
+	best_ap <- which.min(v)
+	list(dif_pairs=v,best=best_ap+n_min-1)
+}