require(survival)

kma <- survfit( Surv(time, status) ~ rx + adhere, data=colon )

kmplot(kma, mark='', simple=FALSE,
	xaxis.at=c(0,.5,1:9)*365, xaxis.lab=c(0,.5,1:9), # n.risk.at
	lty.surv=c(1,2), lwd.surv=1, col.surv=c(1,1,2,2,4,4), # survival.curves
	col.ci=0, # confidence intervals not plotted
	group.names=c('Obs ','Obs tumor adh','Lev','Lev tumor adh','Lev+5FU ','Lev+5FU tumor adh'),
	group.order=c(5,3,1,6,4,2), # order of appearance in the n.risk.at table and legend.
	extra.left.margin=6, label.n.at.risk=FALSE, draw.lines=TRUE, 
	cex.axis=.8, xlab='Years', ylab='Survival Probability', # labels
	grid=TRUE, lty.grid=1, lwd.grid=1, col.grid=grey(.9), 
	legend=TRUE, loc.legend='bottomleft',
	cex.lab=.8, xaxs='r', bty='L', las=1, tcl=-.2  # other parameters passed to plot()
) 
title(main='Chemotherapy for stage B/C colon cancer', adj=.1, font.main=1, line=0.5, cex.main=1)

set.seed(100)
age <- rnorm(80, rep(c(26,36), c(70,10)), 4)
sex <- factor( sample(c('Female','Male'), 80, TRUE) )
group <- paste('Group ', sample(1:4, 40, prob=c(2,5,4,1), replace=TRUE), sep='')
d <- data.frame(age, sex, group)

tplot( age ~ group, data=d, 
    las=1, cex=1, cex.axis=1, bty='L', 
    show.n=TRUE, dist=.25, jit=.05, type=c('db','db','db','d'),
    group.pch=TRUE, pch=c(15,17,19,8), 
    col=c('darkred','darkblue')[c(sex)],
    boxcol=c('lightsteelblue1','lightyellow1',grey(.9), 0), boxborder=grey(.8),
    boxplot.pars=list(notch=TRUE, boxwex=.5)
    )

 set.seed(274)
 x <- rnorm(100, 20, 5)
 y <- rexp(100)
 levels <- as.factor(sample(c("Male","Female"), 100, TRUE))
 jmplot(x, y, levels, col=levels, las=1, type='db', main='Some scores', jit=.02)

set.seed(6)
x <- round(rnorm(400, 100, 4))
y <- round(rnorm(400, 200, 4))
sex <- sample(c('Female','Male'), 400, T)
dsplot(y ~ x, pch=19, col=1+(sex=='Female'), cex=.6, bkgr=T, 
   xlab='measurement 1', ylab='measurement 2', bty='L')
legend( 'bottomright', pch=19, col=1:2, legend=c('Male','Female'), cex=.8 )

show.colors <- function(){
   par(mfrow=c(1,1))
   #   par(mai=c(.4,.4,.4,.4), oma=c(.2,0,0,.2))
   x <- 22 ; y <- 30
   plot(c(-1,x),c(-1,y), xlab='', ylab='', type='n', xaxt='n', yaxt='n', bty='n')
   for(i in 1:x){for(j in 1:y){
         k <- y*(i-1)+j ; co <- colors()[k]
         rect(i-1, j-1, i, j, col=co, border=grey(.5) )}}

   text(rep(-.5,y),(1:y)-.5, 1:y, cex=1.2-.016*y)
   text((1:x)-.5, rep(-.5,x), y*(0:(x-1)), cex=1.2-.022*x)
   title('col=colors()[n]')
}

> colors()[465]
[1] "mediumorchid3"