
```{julia;term=true}
using DataFrames, MixedModels, RData
const dat = Dict(Symbol(k)=>v for (k,v) in 
    load(joinpath(dirname(pathof(MixedModels)), "..", "test", "dat.rda")));
```

```{julia;term=true}
mm1 = fit(LinearMixedModel, @formula(Y ~ 1+S+T+U+V+W+X+Z+(1+S+T+U+V+W+X+Z|G)+(1+S+T+U+V+W+X+Z|H)), dat[:kb07])
```

```{julia;term=true}
mm1.optsum
```

```{julia;term=true}
mm1.trms[1].Λ
```

```{julia;term=true}
mm1.trms[2].Λ
```
