R: implementando meu próprio algoritmo de aumento de gradiente

Estou tentando escrever meu próprio algoritmo de aumento de gradiente. Eu entendo que existem pacotes como gbme, xgboost,mas eu queria entender como o algoritmo funciona escrevendo meus próprios. Estou usando o irisconjunto de dados e meu resultado é Sepal.Length(contínuo). Minha função de perda é...