res2<-inla     Poisson model fit 

csi<-res2$marginals.linear.predictor
theta<-lapply(csi, function(x) inla.emarginal(exp,x))
prob.csi<-lapply(csi, function(x) {1-inla.pmarginal(0,x)})