n = 500
phi<-c(.2,.7)
thet <- 0.2
sd <- .2
x <- matrix(arima.sim(n=n,list(order=c(2,0,1),ar=phi,ma=thet,sd=sd)),ncol=1)
y <- x[7:n]
ym1 <- x[6:(n-1)]
ym2 <- x[5:(n-2)]
H <- cbind(x[4:(n-3)],x[3:(n-4)],x[2:(n-5)],x[1:(n-6)])
g <- y~ym1+ym2
x <- H
res <- gmm(g,x)
# The test is Thets(ym1)=.3 and Theta(ym2)=.6
R <- cbind(c(0,0),diag(2))
c <- c(.3,.6)
lintest(res,R,c)
Run the code above in your browser using DataLab