lm.diallel function

Fitting diallel linear models