Generate a linear projection decomposition for the model with continuous worker hetergoneity

lin.proja(sdata, y_col = "y", k_col = "k", j_col = "j", do_interacted_reg = 1)