diff --git a/docs/src/validation.md b/docs/src/validation.md
index 7051573d..f358a28b 100644
--- a/docs/src/validation.md
+++ b/docs/src/validation.md
@@ -5,7 +5,7 @@
 | Model  |  REML  Metida | REML SPSS|
 |--------|--------|-------|
 | 1 | 1729.4925602367025 | 1729.492560 |
-| 2 | 1904.3265170722132 | 1662.172084 |
+| 2 | 1904.3265170722132 | 1904.327 |
 | 3 | 1772.0953251997046 | 1772.095 |
 | 4 | 1730.1895427398322 | 1730.189543 |
 
diff --git a/test/lme4.jl b/test/lme4.jl
index 56159497..0fb325a3 100644
--- a/test/lme4.jl
+++ b/test/lme4.jl
@@ -23,7 +23,7 @@ transform!(df, :Days => categorical, renamecols=false)
 end
 
 @testset "  CS sleepstudy.csv                                        " begin
-    #REML 1903.327
+    #REML SPSS 1904.327
     # 1662.172084
     # 296.693108
     #Model 2
@@ -58,7 +58,7 @@ end
     lmm = Metida.LMM(@formula(Reaction~1), df;
     random = Metida.VarEffect(Metida.@covstr(Days), Metida.ARH, subj = :Subject)
     )
-    Metida.fit!(lmm)
+    Metida.fit!(lmm; init = [37.9896, 41.1392, 34.1041, 48.1435, 52.2191, 72.4237, 83.3405, 76.7782, 90.2571, 102.617, 0.900038, 6.83327])
     @test lmm.result.reml ≈ 1730.1895427398322 atol=1E-6
 end