Alisson,
Se você espera que a dependência da variável resposta para sua variável quantitativa independente seja linear¹ (ou seja proporcional aos valores intermediários nos intervalos dos valores que você tem nos seus pontos coletados) então o melhor seria tratá-la como variável intervalar e tecnicamente você passou a fazer uma ANCOVA.
Algumas vantagens do ponto de vista "matemático" nessa abordagem são a redução de "perda" de graus de liberdade nos resíduos e maior potência do teste (no fundo você está testando se as funções afim são iguais (hipótese nula) ou diferentes (hipótese alternativa).
Quando esse relacionamento funcional não é possível de fazer, aí um compromisso pode ser feito (tratando os diversos valores quantitativos em fatores) mas lembre-se da questão dos graus de liberdade e que a informação de "distância" do intervalo entre os fatores será perdida, ele considerará o nível do fator para fert=50 apenas diferente de fert=150 e não que eles distam 100 unidades na escala que mede fertilidade.
HTH
--
Cesar Rabak
[1] ou possam ser linearizados por alguma transformação em função do conhecimento do problema, p. ex., como me parece ser o caso na função que você tem fert + fert²