@Chute

n <- 200 lambda_that <- 3.2 x_call <- rpois(n, lambda_that) T <- sim(x_call) lambda_MLE <- T/n #mean(x_call) lambda_MLE loglik_pois <- function(lambda, x) { if (lambda <= 0) return(-Inf) sum(dpois(x, lambda = lambda, log = TRUE))

#VD : (X1, X2,...Xn) ~ N (miu, sigma^2) #(Xbar - mu)/(S/sqrt(n)) ~ Student(n - 1) N = 10^4 Sample = rep(0,N) for (i in 1:N){ X = rnorm(10, 0, 1) Sample[i] = mean(X)/(sd(X)/sqrt(10)) }