Hi, Why this function doesn't work?
function (x) { if (is.factor(x)) { if (!is.ordered(x)) { warning("La mediana non si puo' calcolare!!!") return(NA) } me <- median(unclass(x)) if (me - floor(me) != 0) { warning("Mediana indeterminata") return(NA) } else { levels(x)[me] } } else if(class(x)=="histogram"){ N<-sum(x$counts) cl<-min(which(cumsum(x$counts)>=N/2)) return(x$breaks[cl]+ (N/2-sum(x$counts[1:(cl-1)]))/(x$densit[cl]*N) } else median(x) } thanks Angelo Scozzarella ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.