diff --git a/test/test.jl b/test/test.jl
index 70988882..f79b01b5 100644
--- a/test/test.jl
+++ b/test/test.jl
@@ -473,7 +473,7 @@ end
     @test anovatable.pval ≈ [0.9176105002855397, 0.6522549061174356, 0.0020109339157131302] atol=1E-4
 end
 
-@testset "  Model: BE RDS 1, 2X2                                     "  begin
+@testset "  Model: BE RDS 1, 2X2 + UN test                           "  begin
     dfrds        = CSV.File(joinpath(path, "csv", "berds2x2", "rds1.csv"), types = Dict(:Var => Float64, :Subject => String, :Period => String, :Sequence => String, :Formulation => String )) |> DataFrame
     dropmissing!(dfrds)
     lmm = Metida.LMM(@formula(log(Var)~Sequence+Period+Formulation), dfrds;
@@ -482,6 +482,12 @@ end
     Metida.fit!(lmm)
     anovatable = Metida.typeiii(lmm)
     @test Metida.m2logreml(lmm)  ≈ -1.0745407333692825 atol=1E-8
+
+    # Unstructured
+    lmm = Metida.LMM(@formula(log(Var)~Sequence+Period+Formulation), dfrds;
+    repeated = Metida.VarEffect(Metida.@covstr(Formulation|Subject), Metida.UN),
+    )
+    Metida.fit!(lmm)
 end