Error in `solve.default(tau1^2 * (as.numeric(lambda1) * D + (1 - as.numeric(lambda1)) * 
    I))`: system is computationally singular: reciprocal condition number = 7.20799e-17

test_that("An example for a real data model fitting", {

"An example for a real data model fitting"

An example for a real data model fitting