perc = c(5,10,25,50,75,90,95); P = perc/100; tit <- c("gl",perc)
gl <- c(1,2,3,4,5,10,20,50,100)
ch = function(x) signif(qchisq(x,df=g),2)
write(tit,file="",sep="\t",ncolumns=8); for(g in gl) write(c(g,ch(P)),file="",sep="\t",ncolumns=8)