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))
}

lambda_MLE_so <- optimize(function(lam) - loglik_pois(lam, x_call))
lambda_MLE_so



                          
#bai 3: mtcars
#Uoc luong trung bing mpg cua cac xe trong boj du lieu bang mo hinh chuan
#Tinh khoang tin cay cho trung binh

mpg_data <- mtcars$mpg
n
mean_mpg <- mean(mpg_data)