params <- matrix(c(1,1,1,1,-2,1,0,1,0,0,0,0),4,3)
x <- c(-1,0,1)
rec.mat <- recursive.raw(x, params)
#same results with gen.rec.raw
Pij <- array(NA,dim = c(3,2,4))
Pi1 <- irf(params, x)$f
Pij[,2,] <- Pi1
Pij[,1,] <- 1-Pi1
rec.arr <- gen.rec.raw(Pij,x)
Run the code above in your browser using DataLab