series l_YD = log(YD) series l_CE = log(CE) series y = l_YD series c = 0 u = normal() c[1] = l_CE[1] series c = c(-1) + y + y(-1) + u # print y c ols c 0 c(-1) y y(-1) matrix vmat = $vcv series u2 = $uhat