[R] Boxplot 5% and 95% quantile instead of 25% and 75%
Joris Meijerink
J.Meijerink at witteveenbos.nl
Thu Aug 21 09:37:10 CEST 2008
Hi,
I'm new to the whole R-thing as a replacement for Matlab, not disappointed sofar ;)
I found out how to make nice looking boxplots, but i also would like the make a boxplot with 5% and 95% instead of the standard 25 and 75% quantiles.
My csv input looks something like:
LOCATION FILTER NR DATE VALUE MONTH
Peelhorst01 1 14-Jan-94 23.07 1
Peelhorst01 1 28-Jan-94 23.68 1
Peelhorst01 1 14-Feb-94 23.38 2
Peelhorst01 1 28-Feb-94 23.27 2
Peelhorst01 1 14-Mar-94 23.25 3
Peelhorst01 1 28-Mar-94 23.69 3
Peelhorst01 1 14-Apr-94 23.63 4
Peelhorst01 1 28-Apr-94 23.3 4
Peelhorst01 1 14-May-94 23.14 5
Peelhorst01 1 28-May-94 23.09 5
Peelhorst01 1 14-Jun-94 23.06 6
Peelhorst01 1 28-Jun-94 22.86 6
Peelhorst01 1 14-Jul-94 22.63 7
Peelhorst01 1 28-Jul-94 22.48 7
Peelhorst01 1 14-Aug-94 22.35 8
Peelhorst01 1 28-Aug-94 22.27 8
Peelhorst01 1 14-Sep-94 22.21 9
Peelhorst01 1 28-Sep-94 22.27 9
Peelhorst01 1 14-Oct-94 22.33 10
Peelhorst01 1 28-Oct-94 22.28 10
Peelhorst01 1 14-Nov-94 22.37 11
Peelhorst01 1 28-Nov-94 22.49 11
Peelhorst01 1 14-Dec-94 22.56 12
Peelhorst01 1 28-Dec-94 22.62 12
going on for 13 more years
I used the following to produce a boxplot:
z <- boxplot(VALUE ~ MONTH, data = reeks,
plot = FALSE
)
Then I replace the numbers of the month in jan, feb etc. with
z$names <- c('jan','feb','mrt','apr','mei','jun','jul','aug','sep','okt','nov','dec')
and make the boxplot with the bxp function.
Now I was thinking of using the same sollution by replacing row 2 and 4 in z$stats with the results of the quantile function for 5% and 95% but to be able to calculate that I need the vectors of only 1 month without the other months. How can i do that, or is there even a better/easier sollution to my problem?
kind regards
Joris
------------------------------------------------------------------------------------------
DISCLAIMER:\ This e-mail is strictly confidential and is...{{dropped:16}}
More information about the R-help
mailing list