# distribuzione chi-quadrato: quantili

job <- function(df,p) qchisq(p,df,lower.tail=FALSE)
glib <- c(1:50)
pval <- c(0.95,seq(0.9,0.1,by=-0.1),0.05,0.025,0.01,0.005,0.001)
qntl <- outer(glib,pval,job) # notare l'uso della funzione outer

tbl <- round(qntl,3)
nr <- nrow(tbl)+1
nc <- ncol(tbl)+1
out <- matrix(NA,nr,nc)
out[2:nr,2:nc] <- tbl
out[2:nr,1] <- glib
out[1,2:nc] <- pval
out