library(ade4) data(aviurba) aviurba$traits$feed.hab plot(aviurba$traits[,1],main="feeding habit",col="blue") plot(aviurba$traits[,2],main="feeding stratum",col="magenta") table(aviurba$traits$feed.hab) table(aviurba$traits$feed.strat) table(aviurba$traits$breeding) table(aviurba$traits$migratory) plot(table(aviurba$traits[,1]),col="blue",main="feeding habit") plot(table(aviurba$traits[,1]),col="blue", main="feeding habit",lwd=50,type="h",lend="butt",ylab="") barplot(table(aviurba$traits[,1]),main="feeding habit",col="blue") barplot(table(aviurba$traits[,2]),main="feeding stratum",col="magenta") barplot(table(aviurba$traits[,3]),main="breeding stratum",col="green") barplot(table(aviurba$traits[,4]),main="migration stragegy",col="red") t<-cbind(table(aviurba$traits[,1]),c(15,15,10)) barplot(t,legend.text=rownames(t),col=c(2:4),main="feeding habit",names=c("sample 1","sample 2"),xlim=c(0,8),width=2) barplot(t,legend.text=rownames(t),col=c(2:4),main="feeding habit",names=c("sample 1","sample 2"),xlim=c(0,20),width=2,beside=TRUE) t<-table(rpois(100,5)) t plot(t,col="blue",main="rpois(100,lambda=5)",lwd=20,type="h",lend="butt",ylab="") barplot(t,col="blue",main="rpois(100,lambda=5)") help(hist) library(MASS) data(crabs) str(crabs) dataB<-crabs[which(crabs$sp=="B"),4] dataO<-crabs[which(crabs$sp=="O"),4] hist(dataB) hist(dataB,plot=FALSE) sort(dataB) hist(dataB,col="grey") hist(dataB,col="blue",breaks=seq(6,22,2)) hist(dataB,col="magenta",breaks=20) hist(dataB,col="green",breaks=seq(6,22,1)) help(nclass.Sturges) nclass.Sturges nclass.scott nclass.FD x<-rnorm(100) NC <- function(x)c(Sturges=nclass.Sturges(x),Scott=nclass.scott(x),FD=nclass.FD(x)) NC(x) par(mfrow=c(3,1)) hist(x,col="grey",main="Sturges") hist(x,col="yellow",breaks="Scott",main="Scott") hist(x,col="magenta",breaks="FD",main="Freedman-Diaconis") hist(dataB) abline(v=quantile(dataB),col="red",lwd=3) hist(dataB,probability=TRUE,ylim=range(0,0.18),main="Frontal lobe size",xlab="FL(mm)",ylab="frequency per mm") lines(density(dataB), col='blue',lwd=3) lines(density(dataB, bw=0.5), col='red',lwd=3) library(Hmisc) options(digits=1) histbackback(dataB,dataO,xlab=c("Blue crabs","Orange crabs"),main="Frontal lobe size") library(UsingR) FL<-crabs[which(crabs$sp=="B"),4] RW<-crabs[which(crabs$sp=="B"),5] scatter.with.hist(FL,RW) library(Hmisc) plot(FL,RW,lwd=3,col="blue") histSpike(FL,3,add=TRUE,lwd=3) histSpike(RW,4,add=TRUE,lwd=3) histSpike(FL,3,type="density",col="red",lwd=3,add=TRUE) histSpike(RW,4,type="density",col="red",lwd=3,add=TRUE) library(gplots) x <- rnorm(2000, sd=4) y <- rnorm(2000, sd=1) h2d <- hist2d(x,y,show=FALSE, same.scale=TRUE, nbins=c(20,30)) persp( h2d$x, h2d$y, h2d$counts,ticktype="detailed", theta=30, phi=30,expand=0.5, shade=0.5, col="cyan",ltheta=-30) library(iplots) ihist(dataB) library(plotrix) l <- list(runif(10)*10,1:10,c(1,1,1,1,4,8)) l multhist(l)