num.lv was not set to zero by default if num.RR was used
bugfix for num.RR + num.lv
added option to turn off the default rotation of LVs to principal direction in ordiplot
Formula interface for constrained model wasn't correctly handling factor predictors due to omitting the intercept. Now, the intercept is first included but the first column of the design matrix is removed (otherwise the categories of a factor are all included as separate predictors), so this should be better.