library(ggplot2)
bn <- read.table("..../Bluenile6.txt",header=T,sep="\t",quote="")
univx<-bn$carat

dev.new(width=1500,height=800)
par(mfcol=c(4,5))

#Basic graphics default
plot(univx)

#Dotplots
stripchart(univx)
stripchart(univx,method="jitter")

#Stem and leaf
stem(univx)

#Histograms
hist(univx)
truehist(univx)
minu<-trunc(min(univx))
maxu<-trunc(max(univx))+1
d<-signif((maxu-minu)/100,1)
hist(univx,breaks=c(seq(minu,maxu+d,d)))

#Barcharts of individual values
barplot(table(univx))

#Distribution function
Fcarat<-ecdf(univx)
plot(Fcarat)
plot(Fcarat,verticals=TRUE,do.points=FALSE)

#Density estimates (and rugplot)
d1<-density(univx)
plot(d1)
rug(univx)
plot(density(univx,d1$bw/2))
plot(density(univx,d1$bw/5))

#Boxplots
boxplot(univx)
boxplot(univx,horizontal=TRUE,notch=TRUE)

#QQPlots
qqnorm(univx)
qqnorm(univx,pch=19)
qqline(univx,pch=19)
qqnorm(sort(univx),typ="l")

#Shorth plots
library(lshorth)
lshorth(univx)

#ggplot2
dev.new()
qplot(univx)
dev.new()
qplot(univx,binwidth=0.01)
dev.new()
p <- ggplot(data=bn,aes(x=carat))
p + geom_histogram()
dev.new()
p + geom_density()
dev.new()
qplot(sample=univx)