拟构模型ln(y)=β0+β1lnx+β2(lnx)^2+β3(lnx)^3+e，lny、lnx是时间序列，已经判定是一阶单整。请问接下来怎么做协整检验、格兰杰检验和误差修正模型？我需要判断(lnx)^2和(lnx)^3的平稳性吗？