Constrained ordination analysis with flexible response functions