initialization for j = 0...Nobs do xj = (1,...,1)'