1. Введение

Внешняя торговля является одной из быстроразвивающихся сфер экономической жизни России. Особенно это касается агропродовольственного сектора в целом и его ключевой составляющей – зерновых. Согласно статистике за 2021 год, объем мирового рынка зерна составил более 180 миллионов тонн, причем 35% этого объема пришлось на пшеницу [1]. По результатам 2022 года Россия заняла семнадцатое место среди основных экспортеров продовольствия в мире. По этой категории доля России составляет 8,8%, а выручка показала рост на 16% [1].  Зерновые являются стратегическим товаром для обеспечения продовольствием и их импортируют более 130 стран мира. Наиболее крупными импортерами являются Китай, Турция, Иран, Египет, Саудовская Аравия [2]. Зерно исторически было основой экспорта России. Согласно рейтинга мировой статистике за 2022 г. [3], Россия вышла на первое место в мире с долей 20% по экспорту только пшеницы. В топ 5 мировых экспортеров также вошли: Канада (13,9%), США (13,7%), Франция (9,3%) и Украина (9%). В структуре же экспорта продукции АПК за 2022 год зерновые занимают 32%. В число крупнейших покупателей нашего продовольствия входят Китай (12,3%) и Турция (12,2%). Что касается зерновых, то это Турция и Египет [4], [5]. В связи со стратегической значимостью зерновых как во внешнеэкономической деятельности, так и для экономики страны в целом, представляет интерес ее анализ и поиск современных подходов к моделированию и прогнозированию. Одним из таких подходов, наряду с классическими статистическими методами [6], [7], выступают методы машинного обучения. В последние годы растет интерес к применению методов машинного обучения  к анализу и прогнозированию финансовых рынков, социально-экономических процессов, технических систем [8], [9]. В некоторых работах [10] так же предприняты попытки применения методов машинного обучения  к внешнеэкономической деятельности. Целью исследования является сравнительный анализ внешней торговли зерновыми и определение возможности использования методов машинного обучения  в поддержке принятия решений агропродовольственного экспорта зерновых.

В качестве методов исследования выступают как классические методы статистического анализа, так и методы машинного обучения [11]. В качестве материалов использовались данные федеральной таможенной службы [5], данные агрегаторов мировой торговой статистики Trade Map [2], и российской статистики Ru-stat [4].

2. Основные результаты

Зерновые являются стратегическим товаром с точки зрения продовольственной безопасности любой страны. В этой связи их производство и налаженный экспорт становятся доминантой развития экономики. Исторически Россия входила в число  ведущих мировых держав по экспорту зерновых. Однако в девяностые годы мы потеряли ключевые позиции в этой сфере, превратившись практически в импортера продовольствия. В начале нового тысячелетия государство обратило пристальное внимание на эту сферу и приступило к реализации аграрной реформы. Благодаря предпринятым мерам внутреннее потребление было практически обеспечено, что позволило в качестве приоритета определить выход на внешние рынки.

Pассмотрим динамику объемов (в млрд. долларов) мирового и российского экспорта зерновых (товарная группа «Злаки») в период 2015-2021 гг (рис. 1).

Уравнение регрессии мирового экспорта зерновых [LATEX_FORMULA]y=6,17∙x+87,88 [/LATEX_FORMULA]   показывает явную тенденцию к росту, что обусловлено возрастающими потребностями мирового сообщества в обеспечении продовольствием. Такая же тенденция наблюдается и в российском экспорте. При этом мировые цены на зерновые постоянно растут. На рисунке 2 представлены средние мировые цены на некоторые зерновые.Проведем сравнительный анализ структуры мирового и российского экспорта. Как видно из приведенной диаграммы (рис. 3), ключевыми здесь являются пшеница, рис и кукуруза, на долю которых, в совокупности приходится около 86% в 2015 году. Эта тенденция сохраняется в динамике и в 2021 году составляет уже почти 90%.

Что касается России, то как следует из приведенной диаграммы (рис. 4), ключевыми здесь являются  пшеница и кукуруза, с преобладанием пшеницы, которая в динамике составляла около 70 % в 2015 году и почти 80% в 2021 году.

На примере 2021 года видим (рис. 5), что крупнейшими экспортерами являются США (24%), Индия (10%), Аргентина (10%), Украина (9%), Австралия (8%) и Россия (7%).

В то время как крупнейшими потребителями (рис. 6) являются Китай (25%), Япония (10%), Мексика (7%), Республика Корея (6%), Египет (6%) и Турция (5%).

Детальный анализ ежегодных объемов экспорта (табл. 1) показывает явную тенденцию к росту (рис. 7). Уравнения регрессии [LATEX_FORMULA]y=0,5624∙x+4,6649[/LATEX_FORMULA]  как общего объема , так и его ключевой составляющей (пшеницы и меслина)  [LATEX_FORMULA]y=0,5335∙x+3,2092[/LATEX_FORMULA] подтверждают эту динамику.

Зерновые, входящие в структуру российского экспорта, как показано в таблице 1, можно разделить на три группы: первая это пшеница и меслин, вторая – кукуруза, ячмень и овес, и, третья – рожь, сорго, гречиха и прочие злаки.

Наиболее интересным представляется анализ структуры экспорта в разрезе стран и его динамики, а так же анализ значимости российского экспорта для страны в которую он осуществляется. Такой анализ необходим для выявления тенденций выбора партнеров для определения устойчивости страновой диверсификации экспорта.

Проанализировав данные за период 2013-2020 годы, нами была получена следующая таблица (таблица 2), характеризующая динамику доли двадцати ключевых стран (ежегодно), в которые производился экспорт зерновых. Как видно, таких стран получилось – 34. Причем, какие-то страны весь исследуемый период попадали в двадцатку ключевых импортеров российских зерновых (например, Турция и Египет), часть стран имела смешанную тенденцию (например, Азербайджан и Южная Африка), а часть стран появилась только в конце исследуемого периода (например, Беларусь, Казахстан, Оман).

Визуальный анализ показывает, что возможно выделение трех категорий импортеров российских зерновых: Ключевые (Турция, Египет, Саудовская Аравия, Иран, Азербайджан), Стабильные (Грузия, Израиль, Иордания, Йемен, Кения, Ливан, Ливия, Нигерия, ОАЭ, Судан, Танзания, Южная Африка, Южная Корея) и с неопределенной стабильностью (Бангладеш, Беларусь, Вьетнам, Греция, Индонезия, Испания, Казахстан, Латвия, Марокко, Мексика, Мозамбик, Нидерланды, Оман, Пакистан, Перу, Филиппины).

Однако, традиционные методы анализа, не всегда позволяют выявить существующие тенденции, так как, число стран в которые производится экспорт российских зерновых гораздо больше рассмотренный нами выше тридцати четырех, а во вторых – среди стран не входящих в число экспортеров могут быть потенциальные экспортеры. Поэтому задача сводится к классификации стран импортеров зерновых в разрезе импорта зерновых из России и предсказанию, в какую экспортную российскую категорию попадет страна. При этом нас в первую очередь интересуют страны возможные импортеры, то есть те, которые потенциально могут войти в число крупнейших импортеров российских зерновых.

2.1. Машинное обучение

Описание исходных данных для задачи машинного обучения: источником для получения этих данных послужили сайты «https://russian-trade.com/» и данные агрегаторов мировой торговой статистики Trade Map. Датасет имеет название «Trade-statistics», содержит в себе информацию о торговых операциях с различными товарами  в разных странах в период с 1995 год по 2022 год. Он состоит из 8226871 строк и 10 столбцов, которые включают в себя следующие характеристики (рисунок 8 слева):

· country_or_area: название страны или региона, с которым происходила торговля; Тип данных: object

· year: год, в котором была совершена торговая операция; Тип данных: int64

· comm_code: код товара по классификации торговых товаров; Тип данных: object

· commodity: название товара; Тип данных: object

· flow: тип торговой операции (импорт или экспорт); Тип данных: object

· trade_usd: объем торговых операций в долларах США; Тип данных: int64

· weight_kg: вес товара в килограммах; Тип данных: float64

· quantity_name: единица измерения количества товара; Тип данных: object

· quantity: количество товара; Тип данных: float64

·  category: категория товара; Тип данных: object

Для дальнейшей работы  методами машинного обучения, датасет был предобработан, были добавлены столбцы, определены столбцы, которые будут кодироваться и масштабироваться. Закодировано было название страны (country_or_area). Так как мы рассматривали только зерновые, то столбец название товара (commodity) был исключен. Были рассмотрены – объем импорта из России, объем импорта всего по зерновым (weight) и объем торговых операций в долларах США (trade_usd) так же из России и всего по зерновым.

2.2. Построение моделей и алгоритмов машинного обучения

Для начала, мы разобьем датасет на обучающую и тестовую выборку. Параметр X послужит для нас матрицей признаков, а Y будет вектором целевой переменной. С помощью библиотеки мы разобьем данные так, чтобы они были перемешаны на тестовой и обучающей выборке, чтобы на тестовых данных модель сталкивалась с ранее неизвестными данными, и давала более точные результаты работы, а так же демонстрировала свою устойчивость.  Было принято решение к построению нескольких моделей машинного обучения единовременно (RandomForestClassifier, AdaBoostClassifier, GradientBoostingClassifier, BaggingClassifier, StackingClassifier, VotingClassifier) с использованием метода кросс-валидации из библиотеки, а также определением среднего значения оценки производительности в нескольких метриках [12]. Данный код позволяет сравнить производительность различных моделей машинного обучения на реальных данных агроэкспорта и выбрать наиболее эффективный алгоритм. Фрагмент кода приведен на рисунке 8 (справа).

Выбор метрик для оценки результатов модели является критически важным, поскольку это позволяет определить, насколько хорошо модель работает для конкретной задачи машинного обучения. В данном случае были выбраны следующие четыре метрики: Accuracy (точность), Precision (точность), Recall (полнота), F1 (мера F1) — это гармоническое среднее между точностью Precision и полнотой.  Выбор данных метрик позволяет оценить как точность модели, так и ее способность правильно идентифицировать классы, что может быть критически важным в задаче оценки российского агроэскпорта зерновых. Исходя из основной цели исследования и простоты интерпретации результатов мы остановимся на бинарной классификации, при которой  имеем вероятности принадлежности объектов к классам. Исходя из проведенного ранее анализа экспорта зерновых, мы в качестве порога разделения классов возьмем долю (%) импортеров российских зерновых в общем экспорте России по этой категории и определим ее равной 1%. Сводные характеристики приведены в таблице 3.

Из анализа было выявлено, что GradientBoostingClassifier и ансамблевая модель  StackingClassifier демонстрируют наибольшую производительность по всем метрикам, что говорит о том, что они правильно классифицируют большинство примеров. Однако StackingClassifier также показывает более высокие значения Precision и Recall, чем остальные модели, однако требовала значительно больше времени для обучения. Это означает, что она имеет более высокую точность и полноту по сравнению с другими моделями.

Оценим визуально кривые обучения, чтобы сравнить эти две модели не только по численному признаку, выраженному в метриках, но и по визуальному. Красным будет помечаться обучающая переменная, а зеленым – тестовая (рисунок 9).

Таким образом, в ходе исследования были определены две модели: базовый классификатор GradientBoostingClassifier и ансамблевая модель StackingClassifier, которые показали наилучшую производительность с точки зрения метрик и кривых обучения. В то же самое время, ансамблевая модель достаточно сильно уступает базовому классификатору в плане времени обучения, практически в пять с половиной раз.

Исходная задача сводилась к классификации стран импортеров зерновых в разрезе импорта зерновых из России и предсказанию, в какую экспортную российскую категорию попадет страна: страны-импортеры, которые потенциально входят в число крупнейших импортеров российских зерновых и  прочие страны.

Имеющиеся в наличии данные позволяют оценить работу модели для периода 2013–2021, так как данные 2022 и 2023 года еще не доступны. На рисунке 10 представлена диаграмма числа стран, вошедших в число крупнейших импортеров российских зерновых фактические за указанные периоды и предсказанные моделью.

Фактически в 2021 году четыре страны Камерун, Тунис, Алжир и Китай входят в число крупнейших импортеров российских зерновых. Моделью же так же были выделены еще три страны Сенегал, Демократическая республика Конго и Италия, что показывает достаточно хорошую точность. Рассмотрим подробно  рост стоимости экспорта в период 2019–2020 гг. и 2020–2021 гг. в процентах по странам, фактически вошедшим в двадцатку крупнейших импортеров российских зерновых и странам добавленным моделью,  а так же объем экспорта в 2021 г., тыс. долларов США (по данным trade map) (таблица 4).

Как видно из представленного анализа, все страны, за исключением Алжира показали значительный  процент роста стоимости экспорта. Более того, такие страны как Демократическая республика Конго и Италия,  классифицированные моделью показывают значительный рост стоимости в период 2020-2021 гг. Объем же в стоимостном выражении отличается незначительно, что говорит о потенциале этих стран, особенно если их сравнивать, например с Китаем.

Если поставить вопрос классификации несколько шире и задать в качестве порога разделения классов меньшую долю (%), например 0,5%, то мы очевидно сможем получить большее число стран потенциально выгодных импортеров наших зерновых. Помимо этого, мы так же можем выбрать другой параметр или параметры классификации.

Нами продемонстрированы только некоторые возможности машинного обучения для анализа внешней торговли зерновыми и использования моделей машинного обучения для поддержки принятия решений.

3. Заключение

Проведенный сравнительный анализ внешней торговли зерновыми показывает устойчивый рост российского экспорта зерновых, что обусловлено возрастающими потребностями мирового сообщества в обеспечении продовольствием.

The additional file for this article can be found as follows: