I am looking at the aov function in R. I see that it uses a modified QR factorization routine dqrdc2 based on the Linpack routine dqrdc. Pivoting is done different than the original Linpack function.

My questions:

  *   Why is it necessary to modify the pivoting strategy? Something necessary for aov?
  *   Will the LAPACK function dgeqp3 work well or would the pivoting have o be modified?

Thank you
MRC

	[[alternative HTML version deleted]]

