Bugfix concerning autoregressive models. Thanks EC.

famuvie · Jun 2, 2014 · c6674ee · c6674ee
1 parent a8e8896
commit c6674ee
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 2 deletions.
diff --git a/DESCRIPTION b/DESCRIPTION
@@ -1,7 +1,7 @@
 Package: breedR
 Type: Package
 Title: Statistical methods for forest genetic resources analysts
-Version: 0.6-2
+Version: 0.6-3
 Date: 2013-11-12
 Encoding: UTF-8
 Authors@R: c(person("Facundo", "Muñoz",

diff --git a/R/ar.R b/R/ar.R
@@ -65,7 +65,8 @@ build.ar.model <- function (coord, rho, autofill) {
   # from the precision matrix?
   B <- sparseMatrix(i = 1:length(data.ordering),
                     j = data.ordering,
-                    x = 1)
+                    x = 1,
+                    dims = c(length(data.ordering), dimUinv[1]))
 
   Uinv <- Uinv * gmean(Matrix::diag(B %*% solve(Uinv) %*% Matrix::t(B)))