Rick:

A couple more 'issues' arose as I was trying to duplicate this code. First, the third argument to qntl() that you addressed in the previous comment should not contain a comma so that it represents a row vector rather than a column vector.

Finally, I did not find a function var() for a matrix needed to calculate the standard error of the mean for each airline's arrival times. So I manufactured my own:

dev = x - meanDelay;
varx = dev[##,]/(nrow(x)-1);