search

lem-usp / EvolQG

Tools for Evolutionary Quantitative Genetics
http://cran.r-project.org/web/packages/evolqg/
Other
10 stars 8 forks source link

MeanMatrix fails to converge on a simple example #116

Closed diogro closed 7 years ago

diogro commented 7 years ago

From an user email. This hangs, and should probably be easy:

data = array(NA, dim = c(5, 5, 3))
data[,,1] = c(0.8467  ,0.5195  ,0.8089  ,0.2300, -0.0228,
                     0.5195  ,0.4611  ,0.5595  ,0.4020, -0.2681,
                     0.8089  ,0.5595  ,0.9253,  0.5862, -0.2901,
                     0.2300  ,0.4020  ,0.5862  ,1.0472, -0.7404,
                     -0.0228, -0.2681, -0.2901, -0.7404,  2.2454)

data[,,2] = c(0.7395, 0.6315, 0.8940, 0.4118, 0.1810,
                     0.6315, 0.8664, 0.7663, 0.4747, 0.2356,
                     0.8940, 0.7663, 1.1668, 0.7014, 0.2386,
                     0.4118, 0.4747, 0.7014, 0.7541, 0.1393,
                     0.1810, 0.2356, 0.2386, 0.1393, 0.9006)

data[,,3] = c(0.9727 ,0.7276  ,0.9256, 0.2217, -0.0499,
                     0.7276 ,0.8336  ,0.7469, 0.5929,  0.3071,
                     0.9256 ,0.7469  ,1.0325, 0.5313, -0.0542,
                     0.2217 ,0.5929  ,0.5313, 1.0943,  0.4322,
                     -0.0499, 0.3071, -0.0542, 0.4322,  0.9342)
MeanMatrix(data, verbose = TRUE)