fileread [name='gen.dt'; imethod=read] fgroup=no model [dist=binomial] Kyp; nbinomial=1 fit s(Age;3) + s(Start;3) + s(Number;3) stop