medie<-function (x, r) { #Description # Calcola le medie potenziate di un vettore di dati #Usage # medie(x,r) #Arguments # x Vettore su cui calcolare le medie # r Indice della media potenziata da calcolare #Value # Valore della media potenziata if (r == 0) { prod(x)^(1/length(x)) } else { pd<-r-trunc(r/2) if(pd==0){ mean(x^r)^(1/r) } else{ segno<-sign(sum(x^r)) valore<-abs(mean(x^r))^(1/r) valore*segno } } }