mutability<-function (x) { #Description # Calcola gli indici di mutabilità di Gini e di Shannon normalizzati #Usage # mutability(x) #Arguments # x Vettore di reali positivi su cui calcolare gli indici #Value # gini Valore dell'indice di Gini normalizzato # shannon Valore dell'indice di Shannon normalizzato h<-length(x) hh<-sum(x) p<-x/hh som<-sum(p*(1-p)) gini<-h/(h-1)*som shannon<--sum(p*log(p))/log(h) list(gini=gini,shannon=shannon) }