#bivariate gamma distribution
data(flood)
est_res1 <- MLEce(flood, "BiGam")
summary(est_res1)
#bivariate Weibull distribution
datt = rBiWei(n=50, c(2,3,3,4,0.4))
est_res2 <-MLEce(datt, "BiWei")
summary(est_res2)
#Dirichilet distribution
data(fossil_pollen)
fossil_data <- fossil_pollen/rowSums(fossil_pollen)
eps <- 1e-10
fossil_data <- (fossil_data +eps)/(1+2*eps)
est_res3 <- MLEce(fossil_data, "Dirichlet")
summary(est_res3)
Run the code above in your browser using DataLab