This version produces correct output as per figures 5.3, 5.4(a) & 5.4(b)
of Jim's book.
In this code I add the intercept to cest as WLS does not appear to
return expected intercept despite my inserting a row of 1s