ss <- 1

N <- 10000

for(i in 1:N) {

eta <- rpois(1,4)
xx <- rlnorm(eta,4.26,0.83)
ss[i] <- sum(xx)
rm(xx)

}

mean(ss)

sd(ss)