hat function

Return the hat matrix of a tskrr model