Joris,
I found this (http://ceae.colorado.edu/~balajir/r-session-files/) on the web.
It will do exactly what you want. Get the files:
myboxplot-stats.r
myboxplot.r
Leesferry-mon-data.txt <= example data
The usage is:
#Boxplots
#Source the ‘myboxplot’ codes from Balaji’s directory.
source("myboxplot-stats.r")
source("myboxplot.r")
#Define Variable flow3
flow3=as.data.frame(flow2)
#Only one graph per page:
par(mfrow=c(1,1))
#For 12 months all on one graph:
xs=1:12
zz=myboxplot(split(t(flow3),xs),plot=F,cex=1.0)
zz$names=rep(" ",length(zz$names))
z1=bxp(zz,ylim=range(flow3,zmean),xlab="Month",ylab="Monthly Streamflow (cms)",axes=F)
box()
axis(1,at=z1,labels=months)
axis(2)
points(z1,zmean,lty=1,lwd=2, col="red")
title(main=”Monthly Boxplots of Streamflow”)
I hope this helps. I have a complete example for data I am using if you
need something more complete.
Regards,
Tom
Joris Meijerink wrote:
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}}
______________________________________________
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.
--
Thomas E Adams
National Weather Service
Ohio River Forecast Center
1901 South State Route 134
Wilmington, OH 45177
EMAIL: [EMAIL PROTECTED]
VOICE: 937-383-0528
FAX: 937-383-0033
______________________________________________
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.
