Update polyreg.R