tnx for the nice sollutions, I'll look into them because they look much prettier then my current sollution, quick and durty I replaced some lines in boxplot$stats with the following:
#replace quantile for (i in 1:12) { z$stats[1,i] <- max(reeks$VALUE * match(reeks$MONTH, i), na.rm = TRUE) z$stats[2,i] <- quantile(reeks$VALUE * match(reeks$MONTH, i), probs = c(0.95), na.rm = TRUE, names= FALSE) z$stats[4,i] <- quantile(reeks$VALUE * match(reeks$MONTH, i), probs = c(0.05), na.rm = TRUE, names= FALSE) z$stats[5,i] <- min(reeks$VALUE * match(reeks$MONTH, i), na.rm = TRUE) } end #remove outliers z$out <- NA z$group <- NA This did it for me. regards, Joris ------------------------------------------------------------------------------------------ DISCLAIMER:\ This e-mail is strictly confidential and is...{{dropped:14}} ______________________________________________ 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.