st.prelimFit function

Compute preliminary estimates for a linear model with ST-distributed error term