Как сохранить данные для кривых отклика, связанных с анализом дерева ускоренной регрессии, с использованием "dismo"?

Question

Я моделирую среду обитания потенциальных видов, используя Boosted Regression Tree на "dismo". Код, используемый мной -

lox.BRT0 <- gbm.step(data = training.data, gbm.x = 2:8, gbm.y = 1, 
                     family = "bernoulli", tree.complexity = 7, 
                     learning.rate = 0.005, bag.fraction = 0.75)

Я построил кривые отклика, используя следующее уравнение -

gbm.plot(lox.BRT0, n.plots = 7)

Тем не менее, я хотел дать несколько прогонов, а затем сделать усреднение модели для конечного результата. В то же время я хотел сохранить данные, связанные с переменным откликом для каждого прогона, а затем, наконец, получить кривые отклика как среднее значение для всех прогонов.

Было бы очень хорошо, если бы кто-то мог дать предложения по этому аспекту.