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