time format in r axis

x=c(paste(“81m”, 1:12, sep=”"),paste(“82m”, 1:12, sep=”"))
ind=1:24
y=cumsum(rnorm(24,0))
par(lwd=2)
plot(y~ind,type=’b',xaxt=’n',col=4,ylab=”,xlab=”);grid()

axis(1,seq(1,24,2),labels=F,tcl=-0.3)
#down axis
axis(1,seq(2,24,2),labels=F,tcl=-0.5,lwd.ticks=2)
mtext(x[seq(1,24,2)],1,at=seq(1,24,2),padj=3)
#down axis names
mtext(x[seq(2,24,2)],1,at=seq(2,24,2),padj=1)

R barplot

 x1 <- c(23.2,34.5,76.3,65.8,12.6)
  x2 <- c(15.6,12.4,21.8,20,5.2)
  A <- gl(5,1,5,labels=c("a1","a2","a3","a4","a5"))
  data <- cbind(x1,x2)
  rownames(data) <- levels(A)
  barplot(x1,names.arg=levels(A))
  barplot(t(data),beside=T,ylim=c(0,100),legend.text=colnames(data),
    col=c("grey50","grey80"),ylab="Fréquence")

different y axis unit

Stata
twoway connected y t ||connected var2 t, yaxis(2)

R
par(mar=c(5,4,2,4),bg=’snow2′)
plot(unemp~tse,type=’b',pch=16,
col=4,ylab=”,xaxt=’n',xlab=’Time’)
#minor axis
axis(1,at=stse,labels=stse,
col.axis=1, las=1,tck=-0.009)
#major axis
axis(1,at=ltse,labels=ltse,tck=0.009)
abline(v=1958,col=’purple’,lwd=2)
mtext(‘Unemployment rate’,2,padj=-3.5)
#y2 axis
par(new=T)
plot(GNP~tse,xaxt=’n',ann=F,yaxt=’n',col=2,
xlab=”,ylab=’GNP’,type=’b',pch=15)
axis(4,at=seq(80,120,7),
labels=seq(80,120,7))
mtext(‘GNP’,4,padj=3.5)

/*SAS*/

%LET PATH=C:\DATA;
DATA LONG;
INFILE “&PATH\LONGLY.CSV” DLM=’,';
INPUT GNPD    GNP    Unemp ArmedF Population    Year Employed;
RUN;

goptions reset=all cback=WHITE;
/*V=DOT,STAR,CIRCLE,SQUARE*/

SYMBOL1 C=BLUE V=DOT L=1 W=1 I=JOIN;
SYMBOL2 C=RED V=STAR L=2 W=1 I=JOIN;

/* Create axis definitions  F:TIME FORMAT(‘Helvetica’),
H:HEIGHT*/

AXIS1  LABEL=(F=times H=4 PCT
ANGLE=90 ‘GNP’) OFFSET=(0,0)   ORDER=(200 TO 600 BY 50);
AXIS2  LABEL=(F=times H=4 PCT ANGLE=90 ‘UNEMP’)
OFFSET=(0,0) ORDER=(100 TO 500 BY 50);
AXIS3 ORDER=(1947 TO 1962 BY 2);

/*cREATE LEGEND CONTENT*/
LEGEND1 LABEL=none VALUE=(‘GNP’ ‘unemp’)
position=(top left inside)offset=(3,-1);
LEGEND2 LABEL=none  VALUE=(‘UNEMP’)
position=(top left inside) offset=(3,-2.5);

TITLE F=times H=4 PCT ‘SAS GPLOT’;
PROC GPLOT DATA=LONG;
PLOT GNP*YEAR/OVERLAY VAXIS=AXIS1 HAXIS=AXIS3 LEGEND=LEGEND1 HREF = 1958;
PLOT2 UNEMP*YEAR/OVERLAY VAXIS=AXIS2 LEGEND=legend2;
RUN;
QUIT;