
Last chance! 50% off unlimited learning
Sale ends in
# make an original plot
plot( 11:20, sample(51:60) )
# add some histograms
subplot( hist(rnorm(100)), 15, 55)
subplot( hist(runif(100),main='',xlab='',ylab=''), 11, 51, hadj=0, vadj=0)
subplot( hist(rexp(100, 1/3)), 20, 60, hadj=1, vadj=1, size=c(0.5,2) )
subplot( hist(rt(100,3)), c(12,16), c(57,59), pars=list(lwd=3,ask=FALSE) )
# augment a map
if( require(TeachingDemos) && require(maptools) ){
plot(state.vbm,fg=NULL)
tmp <- get.Pcent(state.vbm)
for( i in 1:50 ){
tmp2 <- as.matrix(USArrests[i,c(1,4)])
tmp3 <- max(USArrests[,c(1,4)])
subplot( barplot(tmp2, ylim=c(0,tmp3),names=c('',''),yaxt='n'),
x=tmp[i,1], y=tmp[i,2], size=c(.1,.1))
}
}
tmp <- rnorm(25)
qqnorm(tmp)
qqline(tmp)
tmp2 <- subplot( hist(tmp,xlab='',ylab='',main=''),
cnvrt.coords(0.1,0.9,'plt')$usr, vadj=1, hadj=0 )
abline(v=0, col='red') # wrong way to add a reference line to histogram
# right way to add a reference line to histogram
op <- par(no.readonly=TRUE)
par(tmp2)
abline(v=0, col='green')
par(op)
# scatter-plot using images
if(require(rimage)){
data(logo)
x <- runif(10)
y <- runif(10)
plot(x,y,type='n')
for(i in 1:10){
subplot(plot(logo),x[i],y[i],size=c(.3,.3))
}
}
Run the code above in your browser using DataLab