vec <- c(1,12,10,8,5,4,7,8,5,9,6,4,5,4,5,8,1) # задаю вектор
N <- sum(vec)
N_1 <- factorial(N)
N_l <- log(N_1)
n_1 <- NULL
for(i in 1:length(vec)) {
n <- factorial(vec[i])
n_1[i] <- log(n)
}
n_l <- sum(n_1)
HB <- (N_l - n_l) / N
S <- length(vec)
NS <- N %/% S
r <- N - S * NS
HB_max <- (1 / N) * log((factorial(N)) / ((factorial(NS) ^ (S - r)) * ((factorial(NS + 1)) ^ r)))
E <- -(HB / HB_max)
N <- sum(vec)
N_1 <- factorial(N)
N_l <- log(N_1)
n_1 <- NULL
for(i in 1:length(vec)) {
n <- factorial(vec[i])
n_1[i] <- log(n)
}
n_l <- sum(n_1)
HB <- (N_l - n_l) / N
S <- length(vec)
NS <- N %/% S
r <- N - S * NS
HB_max <- (1 / N) * log((factorial(N)) / ((factorial(NS) ^ (S - r)) * ((factorial(NS + 1)) ^ r)))
E <- -(HB / HB_max)