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 ______________________________________________ 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.