If the kernel type is not one of the six implemented, the function will terminate with an error message.
The spatial two stage least square estimator is based on the matrix of instruments H=[X,WX,W^2X^2].