Advertisement
Guest User

Untitled

a guest
Aug 26th, 2016
71
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.24 KB | None | 0 0
  1. m = 10^6; lam=5; bet = 1/50; t = numeric(m)
  2. for (i in 1:m) {
  3. n = rpois(1, lam); t[i] = sum(rexp(n+1,bet))
  4. }
  5. mean(t); sd(t); mean(t > 600)
  6. ## 299.7921 # aprx E(T) = 300
  7. ## 165.8158 # aprx SD(T) = 165.83
  8. ## 0.053278 # aprx P(X > 600)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement