R/shapley_mfoc.R

In Inventorymodel: Inventory Models

shapley_mfoc <-
function(n=NA,a=NA,d=NA,K=NA){
  if (is.na(a)==T|sum(is.na(d)==T)==length(d)|sum(is.na(K)==T)==length(K)){ 
    cat("Values for a, d and K are necessary. Please, check them.", sep="\n")
  } else {
    cat("Shapley-Value", sep="\n")
    dk<-order(d/K)
    d<-d[dk];K<-K[dk]
    cind<-as.vector(mfoc(n,a,d,K,cooperation=0))
    shapley<-c();shapley[1]<-cind[1]/n
    for (i in 2:n){
      aux<-0
      for (j in 2:i){aux<-aux+(cind[j]-cind[j-1])/(n-j+1)}
      shapley[i]<-shapley[1]+aux
    }
    return(shapley)
  }
}