Would be great if hts returned **fitted** and **residual** values when using **comb nonnegative** method