[R] A line for SD, colouring simplevioplots, limiting Y axis in barplots, axis in vioplot
Fernando Marmolejo Ramos
fernando.marmolejoramos at adelaide.edu.au
Fri Oct 24 01:55:57 CEST 2008
Dear R people
I have 4 questions about graphics.
-------- QUESTION 1 & 2 -------
I have the following data, which is also set as a dataframe. Then I plot it
using a violin plot
1. How can a draw a line representing 2 SD for each group?
2. How can I colour the simple violin plots in different colours (e.g., Group A
white, Group B grey)?
set.seed(1)
a = rnorm(25, 100, 50)
b = rnorm(25, 300, 50)
# --- dataframing the vectors...
adf <- data.frame(Group = " Group A ", Measure = a)
bdf <- data.frame(Group = " Group B ", Measure = b)
## Combine into a dataframe using rbind
abtruncdata <- rbind(adf, bdf)
attach(abtruncdata)
require(UsingR)
with(abtruncdata, simple.violinplot(Measure~Group))
with(abtruncdata, plotmeans(Measure ~ Group, pch=16, mean.labels=TRUE, col=
c("blue", "red"), digits=2, add=T, p=0.95, barcol=c("blue", "red"), use.t=T,
n.label=F, connect=F, ci.label=F))
-------- QUESTION 3 -------
I have the following bar plot for the same data
3. how can I set the ylim, lets say from 50 to 350 in steps of 25?
tmp = split(Measure, Group)
means = sapply(tmp, mean)
stdev = sapply(tmp, sd)
require(gplots)
plotCI(barcol="blue", pch=16, col="blue", barplot(means, col=c("white",
"gray90"),
ylim=c(0,max(means + stdev)), xlab = "Error bars: ±1SD", ylab="Measure"),
means, stdev, add=TRUE)
-------- QUESTION 4 -------
The same data frame now as violin plot
4. How can I keep a violinplot frame a la simple.violinplot?
## set up frame, without axes
plot(1,1,xlim=c(0.5,2.5),ylim=range(c(a,b)),type="n",
xlab=" ",ylab="Median measure",axes=FALSE)
## bottom axis, with user-specified labels
axis(side=1,at=1:2,labels=c("Group A","Group B"))
axis(side=2)
# but as soon as I call the first violinplot, the upper and right axes appear
again!
In advance, thank you for your answers!
Cheers,
Fer
More information about the R-help
mailing list