StatisticsandGraphicsinR.pptVIP

Statistics and Graphics in R Barplot y - as.data.frame(matrix(runif(36), ncol=3, dimnames=list(month.abb, LETTERS[1:3]))) barplot(as.matrix(y[1:4,]), ylim=c(0,max(y[1:4,])+0.1), beside=T) text(x=seq(1.5, 13, by=1)+sort(rep(c(0,1,2), 4)), y=as.vector(as.matrix(y[1:4,]))+0.02,labels=round(as.vector(as.matrix(y[1:4,])),2)) ysub - as.matrix(y[1:4,]); myN - length(ysub[,1]); mycol1 - gray(1:(myN+1)/(myN+1))[-(myN+1)]; mycol2 - sample(colors(),myN); barplot(ysub, beside=T, ylim=c(0,max(ysub)*1.2), col=mycol2, main=Bar Plot, sub=data: example); legend(topright, legend=row.names(ysub), cex=1.3, bty=n, pch=15, pt.cex=1.8, col=mycol2, ncol=myN) par(mar=c(10.1, 4.1, 4.1, 2.1)); par(xpd=TRUE); barplot(ysub, beside=T, ylim=c(0,max(ysub)*1.2), col=mycol2, main=Bar Plot); legend(x=4.5, y=-0.6, legend=row.names(ysub), cex=1.3, bty=n, pch=15, pt.cex=1.8, col=mycol2, ncol=myN) Barplot Barplot with Confidence Interval stdev=data.frame(Sep.Len=tapply(iris[,1], iris$Species,sd), Sep.Wid =tapply(iris[,2], iris$Species,sd), Pet.Len=tapply(iris[,3], iris$Species,sd), Pet.Wid=tapply(iris[,4], iris$Species,sd) ) bar - barplot(ysub, beside=T, ylim=c(0,max(ysub)*1.2), col=mycol2, main=Bar Plot) arrows(as.vector(bar), as.vector(ysub), as.vector(bar), as.vector(ysub)+ as.vector(as.matrix(stdev)), length=0.15, angle = 90); arrows(as.vector(bar), as.vector(ysub), as.vector(bar), as.vector(ysub)-as.vector(as.matrix(stdev)), length=0.15, angle = 90) legend(topright, legend=row.names(ysub), cex=1.3, bty=n, pch=15, pt.cex=1.8, col=mycol2, ncol=myN) require(gplots) mybarcol - “gray20”; ci.l - as.matrix(ysub-stdev); ci.u - as.matrix(ysub+stdev) mp - barplot2(ysub, beside = TRUE,col = c(lightblue, mistyrose,lightcyan), legend = rownames(ysub), ylim = c(0, 10), main = IRIS DataSet, font.main = 4, col.sub = mybarcol, cex.names = 1.5, plot.ci = TRUE, ci.l = ci.l, ci.u = ci.u,plot.grid = TRUE)