OUTPUT
 solutions='esthobu1_sol';

DATA

  datfile='hobudat1.dat'
    format='(t19,3f6.0,t1,2f6.0,6x)'
    header=0
    dep = tk ry ky
    indep = loom sugu;

  pedfile='hobuped1.ped'
    format='(3i6.0)'
    link=loom
    header=0;

MODEL
 tk = INT + sugu + loom;

COVARIANCE
 loom;

SYSTEM
 non_zero=500000
 total=5000000
 method='AG'
 missing_value=0

END