m1 <- matrix(
c(3685,18344,
4204,18517,
4723,18171,
5242,18517,
5761,17825,
6280,17479,
6799,17652,
7318,17306,
7836,16961,
8355,16096,
8874,15231,
9393,16096,
9912,15923,
10431,15404,
10950,12982,
11467,15404,
11988,17134,
12505,17306,
13024,17825,
13543,18517,
3685,18517,
4204,18517,
4723,18517,
5242,18517,
5761,18344,
6280,18344,
6799,17825,
7318,17998,
7836,16788,
8355,17134,
8874,15923,
9393,15231,
9912,14366,
10431,16096,
10950,16096,
11467,17306,
11988,16788,
12505,18171,
13024,18171,
13543,18517),2,)

Count <- t(matrix(round(60*(18517-m1[2,])/(18727-8351)),,2))
colnames(Count) <- NULL

stop()

postscript(file="fig1.eps",family="NimbusSan",width=8,height=8,
horiz=F,onefile=F,pointsize=15,paper="special")
par(mar=c(5,4,4,5))
barplot(Count[,-1],beside=T,col=c(2,4),xlab="",ylab="",mar=c(5,4,4,5))
legend(4,30,legend=c("ISR","USA"),title="NAMECNTRY",fill=c(2,4))
axis(4,at=(0:12)*3.169494,labels=(0:12)*.01)
mtext("Proportion per bar",side=4,line=2.5)
mtext(1:11,side=1,line=0,at=6*0:10-1.2)
mtext("GRPLNDAVG7",side=1,line=2)
mtext("Count",side=2,line=2.5)
dev.off()