tests/heuristic.ls.capa.ident.R

In kappalab: Non-Additive Measure and Integral Manipulation Functions

library(kappalab)

## number of criteria
n <- 5

## the number of alternatives
n.a <- 10000

## a randomly generated 5-criteria matrix
C <- matrix(rnorm(n*n.a,10,2),n.a,n)

## the corresponding global scores
g <- numeric(n.a)

## generate a random capacity
x <- runif(2^n-1)
for (i in 2:(2^n-1))
    x[i] <- x[i] + x[i-1]
mu <- normalize(capacity(c(0,x)))
for (i in 1:n.a)
  g[i] <- Choquet.integral(mu,C[i,])

mu.in <- as.capacity(uniform.capacity(n))

## the full solution 
lsc <- heuristic.ls.capa.ident(n,mu.in,C,g)
mu.sol <- lsc$solution

## the difference between mu and mu.sol
mu@data - mu.sol@data

lsc