data(welding)
## toy example from Smyth et al.
HLfit(Strength ~ Drying + Material,resid.formula = ~ Material+Preheating ,data=welding)
## toy example of overparameterized model
HLfit(Strength ~ Rods+Thickness*Angle+(1|Rods),resid.formula = ~ Rods+Thickness*Angle ,data=welding)
Run the code above in your browser using DataLab