I have changed theta matrix for It matrix. The formula vcovHAC still confuses me. I am going to ask Kris in the meeting on Thu.