fitTD function

Fit single trial mixed model