powered by
MLE_gamma(g)
g_list <- list() for(i in 1:1000){ g_list <- c(g_list, rgamma(1, shape=100, rate=5)) } mle <- MLE_gamma(g_list) mle
Run the code above in your browser using DataLab