Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- beta <- 0.5*log(2)
- N <- 1000
- t0 <- rexp(N, exp(beta*(sample(5, N, replace=TRUE)-1)))
- x <- seq(0.1,7,0.1)
- pdf <- rep(0, length(x))
- for (i in 1:length(x)) {
- pdf[i] <- mean(dexp(x[i], exp(beta*(0:4))))
- }
- hist(t0, freq=FALSE)
- lines(pdf~x, lwd=2)
Add Comment
Please, Sign In to add comment