Научная статья на тему 'Статистический анализ влияния загрязнения среды обитания на заболеваемость населения в Республике Марий Эл'

Статистический анализ влияния загрязнения среды обитания на заболеваемость населения в Республике Марий Эл Текст научной статьи по специальности «Науки о здоровье»

CC BY
221
40
i Надоели баннеры? Вы всегда можете отключить рекламу.
Область наук
Ключевые слова
ЗАБОЛЕВАЕМОСТЬ / ЗАГРЯЗНЕНИЕ СРЕДЫ ОБИТАНИЯ / ПАНЕЛЬНЫЕ ДАННЫЕ / MORBIDITY / ENVIRONMENTAL POLLUTION / PANEL DATA

Аннотация научной статьи по наукам о здоровье, автор научной работы — Коротков П.А., Трубянов А.Б., Авдеева А.А., Гисмиева А.И.

В статье предложен эконометрической подход к анализу связи показателей экологически зависимой заболеваемости населения и показателей загрязнения среды обитания, отличающийся использованием панельных данных.Цель. Установить количественные связи между состоянием среды обитания и здоровьем населения в условиях дифференцированных антропогенных нагрузок, при которых существует значимый риск здоровью населения в Республике Марий Эл.Материалы и методы. Методология исследования основана на подходах к корреляционно-регрессионному анализу панельных данных. Для выявления показателей загрязнения среды обитания, статистически связанных с показателями экологически зависимой заболеваемости, рассчитываются коэффициенты корреляции Пирсона и ранговые коэффициенты корреляции Спирмена. Затем строятся модели регрессии для панельных данных: модель с фиксированными эффектами и модель со случайными эффектами. Источники панельных данных: Маристат, Управление Роспотребнадзора по Республике Марий Эл и Министерство здравоохранения Республики Марий Эл. В исходный набор вошли 6 показателей загрязнения воды и атмосферы и 7 приоритетных показателей заболеваемости населения в разрезе 15 муниципальных образований Республики Марий Эл за период с 2009 по 2017 гг.Результаты. Анализ матриц коэффициентов Пирсона и Спирмена позволил выявить показатели загрязнения среды обитания, наиболее тесно связанные с показателями заболеваемости населения. Указанные показатели вошли в исходные спецификации моделей панельной регрессии. Получено три статистически значимых модели панельной регрессии, описывающих влияние загрязнения питьевой воды из распределительной сети на заболеваемость бронхиальной астмой детей 0-14 лет с диагнозом, установленным впервые в жизни, и выбросов в атмосферу загрязняющих веществ, отходящих от стационарных источников загрязнения на заболеваемость гастритом и дуоденитом подростков 15-17 лет с диагнозом, установленным впервые в жизни.Заключение. Полученные модели обладают биологическим правдоподобием. Этиопатогенетический анализ подтверждает возможность существования выявленных связей. Найденные статистически значимые связи между загрязнениями среды обитания и здоровьем населения - это не доказательство наличия причинно-следственной связи между ними, а лишь статистическое подтверждение гипотезы о возможном ее наличии. Это подтверждение является необходимым этапом работы для перевода гипотезы в разряд твердо установленных фактов. В дальнейшем предлагается использовать дополнительные, более объективные и полные интегральные оценки качества окружающей среды, например, флуктуирующую асимметрию билатеральных признаков биологических объектов.

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по наукам о здоровье , автор научной работы — Коротков П.А., Трубянов А.Б., Авдеева А.А., Гисмиева А.И.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

STATISTICAL ANALYSIS OF ENVIRONMENTAL POLLUTION IMPACT ON POPULATION MORBIDITY IN THE REPUBLIC OF MARI EL

The article considers an econometric approach to the analysis of relation between the population morbidity rate depending on ecology and the environmental pollution index. Panel data are used in this approach.The purpose is to find quantitative relations between the state of the environment and public health under the differentiated man-caused load threatening public health in the Republic of Mari El. Materials and methods. The research methods are based on the approaches to correlation and regression analysis of the panel data. In order to identify the environmental pollution index statistically related to the morbidity rate, Pearson and Spearman's correlation coefficients were calculated. Then the regression models for the panel data were developed: a fixed-effect model and a random-effect model. The sources of the panel data are the following: Regional Statistics Office in the Mari El Republic (Maristat), Office of the Federal Service for Surveillance on Consumer Rights Protection and Human Wellbeing in the Mari El Republic (Rospotrebnadzor) and Ministry of Healthcare of the Mari El Republic. The data include six air and water pollution indexes and seven priority indicators of the population morbidity rate in 15 municipal districts of the Mari El Republic in the period of 2009-2017.Results. The analysis of the Pearson and Spearman's correlation coefficients helped to identify environmental pollution indexes closely related to the population morbidity rate. These indicators were used as input data of the panel regression model. Three statistically significant panel regression models were identified. They describe the impact of pollution of drinking water from the distributed network on bronchial asthma morbidity among 0-14-aged children diagnosed for the first time in their life; and the impact of emission into the atmosphere of pollutants from the point emission sources on gastritis and duodenitis morbidity among 15-17 aged teenagers diagnosed for the first time in their life.Conclusion. The identified models have biological plausibility. The ethiopathogenetic analysis confirms the possibility of existence of the identified relations. The statistically significant relations between environmental pollution and public health do not prove existence of cause-and-effect links between them. It is statistical demonstration of the hypothesis of their possible existence. This demonstration is an essential work stage to make the hypothesis a hard fact. In the future, it is proposed to use additional, more objective and integral evaluation of environmental quality, for example, the fluctuating asymmetry of bilateral features of biological objects.

Текст научной работы на тему «Статистический анализ влияния загрязнения среды обитания на заболеваемость населения в Республике Марий Эл»

УДК 314.8

DOI: http://dx.doi.org/10.21686/2500-3925-2020-3-58-66

П.А. Коротков1, А.Б. Трубянов2, А.А. Авдеева3, А.И. Гисмиева4

1 Поволжский государственный технологический университет,

Йошкар-Ола, Россия

2 Марийский государственный университет, Йошкар-Ола, Россия

3 Территориальный орган Федеральной службы государственной

статистики по Республике Марий Эл, Йошкар-Ола, Россия 1 Казанский государственный медицинский университет, Казань, Россия

Статистический анализ влияния загрязнения среды обитания на заболеваемость населения в Республике Марий Эл*

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

Цель. Установить количественные связи между состоянием среды обитания и здоровьем населения в условиях дифференцированных антропогенных нагрузок, при которых существует значимый риск здоровью населения в Республике Марий Эл. Материалы и методы. Методология исследования основана на подходах к корреляционно-регрессионному анализу панельных данных. Для выявления показателей загрязнения среды обитания, статистически связанных с показателями экологически зависимой заболеваемости, рассчитываются коэффициенты корреляции Пирсона и ранговые коэффициенты корреляции Спирмена. Затем строятся модели регрессии для панельных данных: модель с фиксированными эффектами и модель со случайными эффектами. Источники панельных данных: Ма-ристат, Управление Роспотребнадзора по Республике Марий Эл и Министерство здравоохранения Республики Марий Эл. В исходный набор вошли 6 показателей загрязнения воды и атмосферы и 7 приоритетных показателей заболеваемости населения в разрезе 15 муниципальных образований Республики Марий Эл за период с 2009 по 2017 гг.

Результаты. Анализ матриц коэффициентов Пирсона и Спирмена позволил выявить показатели загрязнения среды обитания, наиболее тесно связанные с показателями заболе-

ваемости населения. Указанные показатели вошли в исходные спецификации моделей панельной регрессии. Получено три статистически значимых модели панельной регрессии, описывающих влияние загрязнения питьевой воды из распределительной сети на заболеваемость бронхиальной астмой детей 0-14 лет с диагнозом, установленным впервые в жизни, и выбросов в атмосферу загрязняющих веществ, отходящих от стационарных источников загрязнения на заболеваемость гастритом и дуоденитом подростков 15—17 лет с диагнозом, установленным впервые в жизни.

Заключение. Полученные модели обладают биологическим правдоподобием. Этиопатогенетический анализ подтверждает возможность существования выявленных связей. Найденные статистически значимые связи между загрязнениями среды обитания и здоровьем населения — это не доказательство наличия причинно-следственной связи между ними, а лишь статистическое подтверждение гипотезы о возможном ее наличии. Это подтверждение является необходимым этапом работы для перевода гипотезы в разряд твердо установленных фактов. В дальнейшем предлагается использовать дополнительные, более объективные и полные интегральные оценки качества окружающей среды, например, флуктуирующую асимметрию билатеральных признаков биологических объектов.

Ключевые слова: заболеваемость, загрязнение среды обитания, панельные данные.

Peter A. Korotkov1, Aleksey B. Trubyanov2, Anastasiya A. Avdeeva3, Alina I. Gismieva4

1 Volga State University of Technology, Yoshkar-Ola, Russia 2 Mari State University, Yoshkar-Ola, Russia 3 Maristat, Yoshkar-Ola, Russia 4 Kazan State Medical University, Kazan, Russia

Statistical Analysis of Environmental Pollution Impact on Population Morbidity in the Republic of Mari El

The article considers an econometric approach to the analysis of relation between the population morbidity rate depending on ecology and the environmental pollution index. Panel data are used in this approach.

The purpose is to find quantitative relations between the state of the environment and public health under the differentiated man-caused load threatening public health in the Republic of Mari El.

Materials and methods. The research methods are based on the approaches to correlation and regression analysis of the panel data. In order to identify the environmental pollution index statistically related to the morbidity rate, Pearson and Spearman's correlation coefficients were calculated. Then the regression models for the panel data were developed: a fixed-effect model and a random-effect model. The sources of the panel data are the following: Regional Statistics

'Исследование выполнено при финансовой поддержке РФФИ и Правительства Республики Марий Эл в рамках научного проекта № 19-413-120001.

Office in the Mari El Republic (Maristat), Office of the Federal Service for Surveillance on Consumer Rights Protection and Human Wellbeing in the Mari El Republic (Rospotrebnadzor) and Ministry of Healthcare of the Mari El Republic. The data include six air and water pollution indexes and seven priority indicators of the population morbidity rate in 15 municipal districts of the Mari El Republic in the period of 2009-2017.

Results. The analysis of the Pearson and Spearman's correlation coefficients helped to identify environmental pollution indexes closely related to the population morbidity rate. These indicators were used as input data of the panel regression model. Three statistically significant panel regression models were identified. They describe the impact ofpollution of drinking water from the distributed network on bronchial asthma morbidity among 0-14-aged children diagnosed for the first time in their life; and the impact of emission into the

atmosphere of pollutants from the point emission sources on gastritis and duodenitis morbidity among 15—17 aged teenagers diagnosed for the first time in their life.

Conclusion. The identified models have biological plausibility. The ethiopathogenetic analysis confirms the possibility of existence of the identified relations. The statistically significant relations between environmental pollution and public health do not prove existence of cause-and-effect links between them. It is statistical demonstration of the hypothesis of their possible existence. This demonstration is an essential work stage to make the hypothesis a hard fact. In the future, it is proposed to use additional, more objective and integral evaluation of environmental quality, for example, the fluctuating asymmetry of bilateral features of biological objects.

Keywords: morbidity, environmental pollution, panel data.

Введение

На региональном и муниципальном уровнях проводится недостаточное количество комплексных исследований влияния состояния среды обитания на здоровье населения. Между тем, по данным государственной статистики для субъектов Российской Федерации, влияние антропогенной нагрузки (совокупная оценка выбросов загрязняющих веществ в атмосферу, сбросов сточных вод в водные объекты, образование токсичных отходов) на здоровье населения может достигать 40 % и более.

В Республике Марий Эл сложилась неоднозначная медико-экологическая ситуация. С одной стороны, результаты мониторинга указывают на благоприятную экологическую ситуацию в республике: на ее территории не выявлено высокого уровня загрязнения атмосферного воздуха населённых мест, питьевой воды систем централизованного хозяйственно-питьевого водоснабжения, почвы [1]. С другой стороны, анализ приоритетных заболеваний населения Республики Марий Эл, обусловленных неблагоприятным воздействием факторов среды обитания, показывает, что республика по большинству таких показателей, превышающих средние заболеваемости по РФ в 1,4—1,5 раза, отнесена к территориям риска [1].

Противоречие «благоприятная экологическая ситуация в населенных пунктах республики — превышение показателей заболеваемости, обусловленных неблагоприятным воздействием факторов среды обитания, показателей по РФ в 1,5 раза» свидетельствует об актуальности проведения количественного исследования взаимосвязей между загрязнением отдельных компонентов окружающей среды (атмосферы, воды и почвы) и состоянием здоровья населения.

Одна из основных задач социально-гигиенического мониторинга — выявление корреляционно-регрессионных связей между параметрами здоровья и состояния окружающей среды [2]. Подробный обзор по применению статистических моделей в задачах социально-гигиенического мониторинга представлен в монографии [2]. К типичным ошибкам автор монографии относит ошибки: формального построения моделей; трактовки результатов; некорректное использование временных рядов и др. Несмотря на то, что данный обзор выполнен достаточно давно, его результаты во многом остаются актуальными. Действительно, до сих пор в современных работах по социально-гигиеническому мониторингу при оценке влияния загрязнения компонентов окружающей среды (атмосферного воздуха, воды, почвы) на показатели заболеваемости населения часто при-

водятся лишь коэффициенты корреляции [3—8], которые в отличие от коэффициентов регрессии являются показателями качества, но не силы связи. Вместе с тем регрессионный анализ медико-экологических данных, безусловно, также используется [9—12], но применительно к пространственным совокупностям, либо временным рядам, в то время как панельные данные (пространственно-временные совокупности), обладающие значительными преимуществами не используются.

В связи этим в качестве основного способа количественного описания влияния состояния среды обитания на здоровье населения предлагается использовать метод панельных данных.

Цель работы — установление статистических связей между состоянием среды обитания и здоровьем населения в условиях дифференцированных антропогенных нагрузок, при которых существует значимый риск здоровью населения в Республике Марий Эл.

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

Основные задачи:

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

— выявить статистически значимые статистически значимые корреляции между загрязнения загрязнении среды обитания и заболеваемости;

— построить панельные регрессионные модели, устанавливающие количественные вязи между этими показателями.

Данные и методология

Для сбора исходных данных об экологически зависимой заболеваемости [13] и загрязнении среды обитания (атмосферы, воды, почвы) в разрезе 15 муниципальных образований за период с 2009 по 2017 гг. использовались официальные источники статистики: Ма-ристат [14] и Управление Ро-спотребнадзора по Республике Марий Эл [1]. В результате прямого обращения в Минздрав Республики Марий Эл данные по приоритетным заболеваниям, обусловленным неблагоприятным воздействием факторов среды обитания были уточнены.

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

Зависимые переменные:

Y1 — младенческая смертность по муниципальным образованиям Республики Марий Эл на 1000 родившихся живыми;

Y1 — заболеваемость анемией детей 0—14 лет с диагнозом, установленным впервые в жизни (на 100 тыс. соответствующего населения);

Y3 — заболеваемость бронхиальной астмой детей 0—14 лет с диагнозом, установленным впервые в жизни (на 100

тыс. соответствующего населения);

Y4 — заболеваемость анемией подростков 15—17 лет с диагнозом, установленным впервые в жизни (на 100 тыс. соответствующего населения);

Y5 — заболеваемость бронхиальной астмой подростков 15—17 лет с диагнозом, установленным впервые в жизни (на 100 тыс. соответствующего населения);

Y6í — заболеваемость гастритом и дуоденитом подростков 15—17 лет с диагнозом, установленным впервые в жизни (на 100 тыс. соответствующего населения);

Y7 — заболеваемость язвой желудка и 12-перстной кишки взрослых от 18 лет и старше с диагнозом, установленным впервые в жизни (на 100 тыс. соответствующего населения).

Независимые переменные:

Х1 — удельный вес нестандартных проб воды из поверхностных водоемов по микробиологическим показателям, %;

Х2 — удельный вес нестандартных проб воды из поверхностных водоемов по санитар-но-химическим показателям, %;

Х3 — удельный вес проб питьевой воды из распределительной сети, не соответствующих гигиеническим требованиям по санитарно-хи-мическим показателям, %;

Х4 — удельный вес проб питьевой воды из распределительной сети, не соответствующих гигиеническим требованиям по микробиологическим показателям, %;

Х5 — удельный вес проб воды из нецентрализованных источников водоснабжения, не отвечающих гигиеническим требованиям по микробиологическим показателям, %;

Х6 — выбросы в атмосферу загрязняющих веществ, отходящих от стационарных источников загрязнения, т/км2.

Для характеристики загрязнения атмосферы используется единственный доступный показатель «Выбросы в атмос-

феру загрязняющих веществ, отходящих от стационарных источников загрязнения, т/км2». Количество выбросов нормируется на площадь муниципальных образований [15]. Правомерность использования такого показателя обусловлена небольшим размером территории муниципальных образований Республики Марий Эл.

Анализ распределения переменных выявил наличие «скошенности» вправо. Для уменьшения асимметрии распределения в дальнейшем анализе использовались логарифмы переменных: у1, ..., у7, х1, ... х6. В связи с тем, что многие показатели загрязнения воды имеют нулевые значения, что при их логарифмировании приводит к потере значений переменных, дополнительно использовались первые разности: dY1, ..., dY7, йХ1, ..., йХ6.

Традиционно в статистическом анализе данных о состоянии окружающей среды и заболеваемости населения различают временные ряды и пространственные совокупности [2]. Регрессионный анализ позволяет получать достоверные результаты при использовании длинных временных рядов и предметном (медицинском) обосновании возможной причинно-следственной связи между загрязнением и заболеваемостью [16]. При коротких рядах, когда связь конкретного типа загрязнения с уровнем конкретной патологии явно не просматривается, можно получить ложные корреляции. При анализе пространственных совокупностей вероятность получения ложной корреляции значительно меньше, поскольку условие независимости наблюдений в разных точках, как правило, соблюдается.

Методология исследования основана на подходах к корреляционно-регрессионному анализу панельных данных [17—22]. Панельные данные представляют собой двумерные — пространственно-вре-

менные — массивы. Панельные данные позволяет получать более эффективные оценки [22], поскольку характеризуются большим количеством наблюдений, что увеличивает число степеней свободы и уменьшает мультиколлинеарность факторов за счет учета индивидуальных различий.

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

На первом этапе для выявления показателей загрязнения среды обитания, оказывающих влияние на показатели экологически зависимой заболеваемости, рассчитываются коэффициенты корреляции Пирсона (г) и ранговые коэффициенты корреляции Спир-мена (г5).

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

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

Модель регрессии с фиксированными эффектами имеет вид:

у1, =а + р1+81г + РХ, + £,.,, (1)

где у, ( — зависимая переменная, характеризующая экологически зависимую заболеваемость в /-м муниципальном образовании Республики Марий Эл в период /; а — свободный член; X1 — Х-мерный вектор объясняющих переменных, характеризующих загрязнение среды

обитания в /-м муниципальном образовании Республики Марий Эл в период /; в — вектор соответствующих коэффициентов панельной регрессии;

— постоянные во времени индивидуальные эффекты для каждого в /-го муниципального образовании Республики Марий Эл, улавливающие влияние неучтенных переменных; 8^ — индивидуальные тренды, которые при необходимости могут быть включены в модель; еу — остатки панельной регрессии, некоррелированные с объясняющими переменными и одинаково распределенные по муниципальным образованиям и времени с нулевым средним и дисперсией ае2; / = 1, 2, ..., N — номер муниципального образования, N = 14; / = 1, 2, ..., Т — номер периода, Т — число наблюдений для /-го муниципального образования, Т = 9.

Для оценки параметров модели с фиксированными эффектами (1) используется панельный метод наименьших квадратов (ПМНК).

В модели с фиксированными эффектами (1) оценка коэффициентов в основана на информации об изменении экологически зависимой заболеваемости во времени внутри одного муниципального образования в зависимости от изменения независимых переменных, характеризующих загрязнение среды обитания. При этом не используется информация о различии в среднем уровне экологически зависимой заболеваемости между муниципальными образованиями — оно полностью объясняется постоянными индивидуальными эффектами ц.

Модель регрессии со случайными эффектами имеет вид:

Уц =Р+Рх[, + и , (2)

где ^ — свободный член; и{ — индивидуальные эффекты, которые рассматриваются как случайные величины, не за-

висящие от времени и некоррелированные с остатками ец, имеющие нулевое условное (по объясняющим переменным) математическое ожидание и дисперсиию а2

Параметры модели случайных эффектов (2) оцениваются обобщенным методом наименьших квадратов (ОМНК), а для дисперсии случайных индивидуальных эффектов и дисперсии остатков используются оценки Swamy—Arora [23].

Оценки модели случайных эффектов (2) являются более эффективными, чем оценки модели (1), поскольку учитывают как изменение экологически зависимой заболеваемости внутри одного муниципального образования, так и его различие между муниципальными образованиями. Оценки модели случайных эффектов являются состоятельными, если верно предположение о независимости случайных эффектов иI от объясняющих переменных и остатков модели ец.

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

Расчеты выполняются в статистическом пакете EViews 11.

Результаты и обсуждение

Результаты корреляционного анализа связей показателей загрязнения среды обитания и экологически зависимой заболеваемости представлены в таблицах 1—4.

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

При построении моделей панельной регрессии использовались логарифмированные значения переменных, по-

Таблица 1

Матрица парных коэффициентов корреляции Пирсона

Х1 Х2 Х3 Х4 Х5 Х6 У1 У2 У3 У4 У5 У6 У7

Х1 1,000

х2 0,404 1,000

х3 0,256 -0,144 1,000

Х4 0,152 -0,147 0,132 1,000

Х5 0,233 -0,123 0,367 0,148 1,000

Х6 -0,039 -0,041 -0,150 -0,110 0,168 1,000

У1 0,109 -0,055 -0,005 0,044 -0,144 -0,350 1,000

У2 -0,051 0,040 0,167 0,288 0,054 -0,170 0,039 1,000

У3 -0,039 0,161 -0,258 0,210 0,102 0,203 0,000 -0,058 1,000

У4 -0,173 -0,225 0,032 0,068 0,015 -0,267 0,087 0,657 -0,087 1,000

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

У5 -0,404 0,369 -0,286 -0,041 0,035 -0,077 0,019 -0,266 0,021 -0,130 1,000

У6 -0,020 0,062 0,111 -0,081 0,138 0,275 -0,136 -0,140 0,041 -0,094 0,249 1,000

У7 -0,322 -0,089 -0,399 0,027 -0,326 0,066 -0,062 0,166 0,227 0,138 0,036 -0,023 1,000

Примечание: логарифмированные переменные.

Таблица 2

Матрица парных коэффициентов корреляции Пирсона

ёХ1 ёХ2 ёХ3 йХ4 ёХ5 ёХ6 ёУ1 ёУ2 ёУ3 ёУ4 ёУ5 ёУ6 ёУ7

ёХ1 1,000

ёХ2 0,261 1,000

ёХ3 -0,117 -0,130 1,000

ёХ4 0,103 -0,149 -0,187 1,000

ёХ5 -0,028 -0,174 0,139 0,081 1,000

ёХ6 0,040 0,001 0,045 0,060 0,075 1,000

ёУ1 0,041 0,036 -0,113 -0,029 0,041 -0,010 1,000

ёУ2 -0,028 -0,157 0,062 0,044 0,055 0,011 0,134 1,000

ёУ3 0,101 0,016 0,017 0,169 -0,031 -0,075 0,134 0,101 1,000

ёУ4 -0,118 -0,155 0,092 -0,082 -0,074 -0,060 0,019 0,431 0,121 1,000

ёУ5 -0,064 0,010 -0,064 0,060 0,096 -0,001 -0,082 0,036 -0,161 0,027 1,000

ёУ6 0,022 -0,012 0,136 0,076 0,023 0,018 -0,028 0,094 -0,195 0,229 0,038 1,000

ёУ7 -0,059 -0,100 0,010 0,014 -0,058 -0,026 0,014 0,292 0,016 0,319 0,049 0,100 1,000

Примечание: первые разности.

Таблица 3

Матрица коэффициентов корреляции Спирмена

Х1 Х2 Х3 Х4 Х5 х6 ^2 ^4 ^7

Х1 1,000

Х2 0,115 1,000

Х3 0,236 0,167 1,000

Х4 -0,048 -0,022 -0,002 1,000

Х5 0,207 0,060 0,481 0,128 1,000

Х6 0,292 0,016 0,258 0,194 -0,026 1,000

-0,135 -0,028 0,082 -0,048 0,081 -0,176 1,000

-0,078 0,023 -0,142 0,060 0,006 -0,046 0,021 1,000

0,167 0,037 0,038 0,089 0,128 0,208 -0,006 -0,105 1,000

-0,157 -0,088 -0,147 -0,076 0,002 -0,158 0,042 0,695 -0,107 1,000

0,009 -0,044 0,026 0,116 -0,059 0,201 -0,077 -0,400 0,141 -0,341 1,000

-0,144 0,243 0,230 0,100 -0,032 0,302 -0,088 -0,095 0,051 0,004 0,308 1,000

^7 -0,061 -0,116 -0,262 0,081 -0,273 0,039 -0,097 0,176 0,127 0,189 -0,094 -0,024 1,000

Примечание: исходные переменные.

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

В результате эконометри-ческого моделирования были получены окончательные панельные регрессии.

Модель с фиксированными эффектами:

Уз,,, = 4,71 + 0,24х4,, + е,, (3)

Модель случайных эффектов:

Уз ,,, = 4,71 + 0,22х4,, + е,, (4)

Модель случайных эффектов:

К,., = 7,5 + 0,07 х6,, (5)

Все модели статистически значимы, а знаки при коэффициентах регрессии имеют ожидаемые знаки. Параметры моделей панельной регрессии (3) и (4), оцененные разными методами, имеют близкие значения.

Тесты на отсутствие постоянных индивидуальных эффектов подтверждают корректность выбранной спецификации модели (3). Тест Ха-усмана не отвергает нулевую гипотезу о некоррелированности объясняющих переменных и случайных эффектов для моделей (4—5).

На заболеваемость бронхиальной астмой детей 0—14 лет с диагнозом, установленным впервые в жизни (на 100 тыс. соответствующего населения) (уц ) влияет удельный вес проб питьевой воды из распределительной сети, не соответствующих гигиеническим требованиям по микробиологическим показателям, % (х4г^).

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

проявление бронхиальной астмы) [24-25].

На заболеваемость гастритом и дуоденитом подростков 15-17 лет с диагнозом, установленным впервые в жизни (на 100 тыс. соответствующего населения) (Уб^) влияют выбросы в атмосферу загрязняющих веществ, отходящих от стационарных источников загрязнения, т/км2 (х6г^).

Исходя из этиологии и патогенеза заболеваний желудочно-кишечного тракта можно сделать вывод, что при вдыхании загрязненного воздуха через рот (например, при физической активности), взвеси химических веществ вместе со слюной попадают в желудочно-кишечный тракт. Тем самым химические вещества повреждают защитный барьер желудка, вследствие чего под действием соляной кислоты повреждаются клетки слизистой оболочки, что, в свою очередь, приводит к воспалительному процессу, т.е. гастриту. Аналогичное влияние данные вещества оказывают и на двенадцатиперстную кишку (она является следующим после желудка отделом желудочно-кишечного тракта), провоцируя дуоденит [24, 26].

Заключение

В результате параметрического и непараметрического корреляционного анализа панельных данных для 15 муниципальных образований Республики Марий Эл за период с 2009 по 2017 гг. выявлены показатели загрязнения среды обитания, наиболее тесно связанные с показателями экологически зависимой заболеваемости в муниципальных образованиях Республики Марий Эл.

В соответствии с принятым подходом выявленные показатели включались в исходные спецификации моделей панельной регрессии.

Были получены три статистически значимые модели

панельной регрессии, характеризующие связь заболеваемости бронхиальной астмой детей 0-14 лет с диагнозом, установленным впервые в жизни (на 100 тыс. соответствующего населения) и удельным весом проб питьевой воды из распределительной сети, не соответствующих гигиеническим требованиям по микробиологическим показателям, %; между заболеваемостью гастритом и дуоденитом подростков 15-17 лет с диагнозом, установленным впервые в жизни (на 100 тыс. соответствующего населения) и выбросами в атмосферу загрязняющих веществ, отходящих от стационарных источников загрязнения, т/км2.

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

Найденные статистически значимые связи между загрязнениями среды обитания и здоровьем населения - это не доказательство наличия причинно-следственной связи между ними, а лишь статистическое подтверждение гипотезы о возможном ее наличии [2]. Это подтверждение является необходимым (но недостаточным) этапом работы для перевода гипотезы в разряд твердо установленных фактов.

Для полного доказательства нужны еще подтверждения. Так, планируется использовать дополнительные, более объективные и полные интегральные оценки качества окружающей среды, например, по состоянию живых существ. Одним из способов такой оценки является флуктуирующая асимметрия билатеральных признаков биологических объектов. Особенность такого интегрального показателя заключается в том, что исследуется не типы загрязнений (почвы, воды, воздуха и др.), а реакция живых организмов на условия среды.

Литература

1. О состоянии санитарно-эпидемиологического благополучия населения Республики Марий Эл в 2018 году: Доклад. Йошкар-Ола: Управление Федеральной службы по надзору в сфере защиты прав потребителей и благополучия человека по Республике Марий Эл, 2019. 227 с.

2. Вараксин А.Н. Статистические модели регрессионного типа в экологии и медицине. Екатеринбург: Гощицкий, 2006. 255 с.

3. Ситникова О.О. Оценка экологических рисков здоровью населения в медико-демографических процессах при антропогенном загрязнении атмосферного воздуха (на примере Белгородской области) // Материалы VI Международной студенческой научной конференции «Студенческий научный форум». 2014.

4. Веремчук Л.В., Черпак Н.А., Гвозденко Т.А., Волкова М.В. Методология оценки влияния загрязнения атмосферного воздуха на формирование уровней общей заболеваемости бронхиальной астмой // Гигиена и санитария. 2015. Т. 94. № 3. С. 119-122.

5. Куркатов С.В., Тихонова И.В., Иванова О.Ю. Оценка риска воздействия атмосферных загрязнений на здоровье населения г. Норильска // Гигиена и санитария. 2015. Т. 94. № 2. С. 28-31.

6. Канатникова Н.В., Егорова Н.А. Влияние жесткости питьевой воды на заболеваемость населения г. Орла // Гигиена и санитария. 2017. Т. 96. № 3. С. 235-240.

7. Лучанинова В.Н., Цветкова М.М., Веремчук Л.В., Крукович Е.В., Мостовая И.Д. Состояние здоровья детей и подростков и факторы, влияющие на его формирование // Гигиена и санитария. Т. 96. № 6. С. 561-568.

8. Малыш Н.Г., Доан С.И. Использование факторного анализа при исследовании эпидемического процесса острых кишечных инфекций // Гигиена и санитария. 2017. Т. 96. № 6. С. 519-523.

9. Кику П. Ф., Измайлова О.А. Горборуко-ва Т.В., Ананьев В.Ю. Влияние эколого-гиги-енических факторов среды обитания на распространение болезней органов дыхания у населения Приморского края // Гигиена и санитария. 2012. № 5.

10. Трифонова Т.А., Марцев А.А. Оценка влияния загрязнения атмосферного воздуха на заболеваемость населения Владимирской области // Гигиена и санитария. 2015. Т. 94. № 4.

11. Зайкова З.А., Бурдуковская А.В., Белых А.И. Определение приоритетных неблагоприятных факторов окружающей среды // Гигиена и санитария. 2016. Т. 95. № 12.

12. Измайлова О.А., Кику П.Ф., Ярыги-на М.В., Морева В.Г., Ананьев В.Ю., Косола-пов А.Б. Гигиенические аспекты распространения экологозависимых заболеваний детей и подростков Приморского края // Гигиена и санитария. 2016. Т. 95. № 11.

13. Щербо А.П. Окружающая среда и здоровье: подходы к оценке риска. СПб.: СПбМАПО, 2002. 376 с.

14. Городские округа и муниципальные районы Республики Марий Эл: статистический сборник. Йошкар-Ола: Маристат, 2017. 266 с.

15. Республика Марий Эл в цифрах: Краткий статистический сборник. Йошкар-Ола: Ма-ристат. 2019. 385 с.

16. Кацнельсон Б.А., Кошелева А.А., Привалова Л.И., Кузьмин С.В., Малых О.Л., Халь-фин Р.А., Никонов Б.И., Озкайнак Х., Шу Дж. Влияние кратковременных повышений загрязнения атмосферного воздуха на смертность населения // Гигиена и санитария. 2000. № 1. С. 15—18.

17. Коротков П.А., Трубянов А.Б., Загай-нова Е.А., Никоноров К.Н. Сопоставительный анализ моделей оценки экологической эффективности крупных городов // Современные проблемы науки и образования. 2015. № 2. С. 328.

18. Коротков П.А., Загайнова Е.А. Взаимосвязь уровня распространенности самоубийств и продолжительности рабочего времени // Статистика и экономика. 2017. Т. 14. № 4. С. 41—53.

19. Шиманович Г. Внешний долг Беларуси: опыт постсоциалистических стран. Рабочий материал Исследовательского центра ИПМ WP/09/01, 2009. 24 с.

20. Копнова Е.Д., Розенталь О.М. Экономе-трический анализ экологического менеджмента рыбных ресурсов // Прикладная эконометрика. 2010. № 2 (18). С. 90-100.

21. Ларин А.В., Тарунина Е.Н. Предпринимательская активность и уровень экономического развития: форма зависимости // Прикладная эконометрика. 2015. № 37 (1). С. 3-26.

22. Вербик М. Модели, основанные на панельных данных // Прикладная эконометрика. 2006. № 1. С. 94-135.

23. Swamy P.A. V.B. and Arora S.S. The exact finite sample properties of the estimators of coefficients in the error components regression models // Econometrica. 1972. № 40. С. 261-275.

24. Серов В.В., Пальцев М.А. Патологическая анатомия: курс лекций. Учебное пособие. М.: 1998. 640 с.

25. Бронхиальная астма: клинические рекомендации. M.: Министерство здравоохранения Российской Федерации, 2016. 76 с.

26. Патофизиология: учебник. М.: ГЭОТАР-Медиа, 2007. 496 с.

References

1. O sostoyanii sanitarno-epidemMogicheskogo blagopoluchiya naseleniya Respubliki Many El v 2018 godu: Doklad. Yoshkar-Ola: Upravleniye Federal'noy sluzhby po nadzoru v sfere zashchity prav potrebiteley i blagopoluchiya cheloveka po Respublike Many El = On the state of the sanitary-epidemiological well-being of the population of the Republic of Mari El in 2018: Report. Yoshkar-Ola: Office of the Federal Service for Supervision of Consumer Rights Protection and Human Well-Being in the Republic of Mari El, 2019. 227 p. (In Russ.)

2. Varaksin A.N. Statisticheskiye modeli regressionnogo tipa v ekologii i meditsine = Statistic models of regression type in ecology and medicine. Yekaterinburg: Goshchitsky; 2006. 255 p. (In Russ.)

3. Sitnikova O.O. Assessment of environmental risks to public health in medical and demographic processes during anthropogenic air pollution (on the example of the Belgorod region). Materialy VI Mezhdunarodnoy studencheskoy nauchnoy konferentsii «Studencheskiy nauchnyy forum» = Materials of the VI International Student Scientific Conference "Student Scientific Forum". 2014. (In Russ.)

4. Veremchuk L.V., Cherpak N.A., Gvozden-ko T.A., Volkova M.V. Methodology for assessing the effect of air pollution on the formation of levels of the general incidence of bronchial asthma. Gigiyena i sanitariya = Hygiene and sanitation. 2015; 94; 3: 119-122. (In Russ.)

5. Kurkatov S.V., Tikhonova I.V., Ivanova O.Yu. Risk assessment of the effects of atmospheric pollution on the health of the population of Norilsk. Gigiyena i sanitariya = Hygiene and sanitation. 2015; 94; 2: 28-31. (In Russ.)

6. Kanatnikova N.V., Yegorova N.A. The effect of hardness of drinking water on the incidence of the population of the city of Orel. Gigiyena i sanitariya = Hygiene and sanitation. 2017; 96; 3: 235-240. (In Russ.)

7. Luchaninova V.N., Tsvetkova M.M., Veremchuk L.V., Krukovich Ye.V., Mostovaya I.D. The health status of children and adolescents and factors affecting its formation. Gigiyena i sanitariya = Hygiene and sanitation. 96; 6: 561-568. (In Russ.)

8. Malysh N.G., Doan S.I. The use of factor analysis in the study of the epidemic process of acute intestinal infections. Gigiyena i sanitariya = Hygiene and sanitation. 2017; 96; 6: 519-523. (In Russ.)

9. Kiku P. F., Izmaylova O.A. Gorborukova T.V., Anan'yev V.Yu. The influence of environmental and hygienic environmental factors on the spread of respiratory diseases in the population of Primorsky Krai. Gigiyena i sanitariya = Hygiene and sanitation. 2012; 5. (In Russ.)

10. Trifonova T.A., Martsev A.A. Evaluation of the effect of air pollution on the incidence of the population of the Vladimir region. Gigiyena i sanitariya = Hygiene and sanitation. 2015; 94; 4. (In Russ.)

11. Zaykova Z.A., Burdukovskaya A.V., Belykh A.I. Determination of priority adverse environmental factors. Gigiyena i sanitariya = Hygiene and sanitation. 2016; 95; 12. (In Russ.)

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

12. Izmaylova O.A., Kiku P.F., Yarygina M.V., Moreva V.G., Anan'yev V.Yu., Kosolapov A.B. Hygienic aspects of the spread of environmentally dependent diseases in children and adolescents of the Primorsky Territory. Gigiyena i sanitariya = Hygiene and sanitation. 2016; 95; 11. (In Russ.)

13. Shcherbo A.P. Okruzhayushchaya sreda i zdorov'ye: podkhody k otsenke riska = Environment and health: risk assessment approaches. Saint Petersburg: SPbMAPO; 2002. 376 p. (In Russ.)

14. Gorodskiye okruga i munitsipal'nyye rayony Respubliki Mariy El: statisticheskiy sbornik = . Urban districts and municipal districts of the Republic of Mari El: statistical compilation. Yoshkar-Ola: Maristat; 2017. 266 p. (In Russ.)

15. Respublika Mariy El v tsifrakh: Kratkiy statisticheskiy sbornik = The Republic of Mari El in numbers: A brief statistical compilation. Yoshkar-Ola: Maristat. 2019; 385 p. (In Russ.)

16. Katsnel'son B.A., Kosheleva A.A., Privalova L.I., Kuz'min S.V., Malykh O.L., Khal'fin R.A., Nikonov B.I., Ozkaynak KH., Shu Dzh. The effect of short-term increases in air pollution on population mortality. Gigiyena i sanitariya = Hygiene and sanitation. 2000; 1: 15-18. (In Russ.)

17. Korotkov P.A., Trubyanov A.B., Zagayno-va Ye.A., Nikonorov K.N. Comparative analysis of models for assessing the environmental performance of large cities. Sovremennyye problemy nauki i obrazovaniya = Modern problems of science and education. 2015; 2: 328. (In Russ.)

18. Korotkov P.A., Zagaynova Ye.A. The relationship between the level of suicide prevalence and working hours. Statistika i ekonomika = Statistics and Economics. 2017; 14; 4: 41—53. (In Russ.)

19. Shimanovich G. Vneshniy dolg Belarusi: opyt postsotsialisticheskikh stran. Rabochiy material Issledovatel'skogo tsentra IPM WP/09/01 = External debt of Belarus: the experience of post-socialist countries. Work material of the IPM Research Center WP / 09/01. 2009; 24 p. (In Russ.)

20. Kopnova Ye.D., Rozental' O.M. Econometric analysis of ecological management of fish resources. Prikladnaya ekonometrika = Applied Econometrics. 2010; 2 (18): 90-100. (In Russ.)

21. Larin A.V., Tarunina Ye.N. Entrepreneurial activity and the level of economic development: a form of dependence. Prikladnaya ekonometrika = Applied Econometrics. 2015; 37 (1): 3-26. (In Russ.)

22. Verbik M. Models based on panel data. Prikladnaya ekonometrika = Applied Econometrics. 2006; 1: 94-135. (In Russ.)

23. Swamy P.A. V.B. and Arora S.S. The exact finite sample properties of the estimators of coefficients in the error components regression model. Econometrica. 1972; 40: 261-275.

24. Serov V.V., Pal'tsev M.A. Patologicheskaya anatomiya: kurs lektsiy. Uchebnoye posobiye = Pathological anatomy: course of lectures. Tutorial. Moscow: 1998. 640 p. (In Russ.)

25. Bronkhial'naya astma: klinicheskiye rekomendatsii = Bronchial asthma: clinical

Сведения об авторах

Петр Анатольевич Коротков

Поволжский государственный технологический университет, Йошкар-Ола, Россия Эл. почта: korotkovpa@volgatech.net

Алексей Борисович Трубянов

Марийский государственный университет, Йошкар-Ола, Россия Эл. почта: true47@mail.ru

Анастасия Анатольевна Авдеева

Территориальный орган Федеральной службы государственной статистики по Республике Марий Эл, Йошкар-Ола, Россия Эл. почта: nastya_87_02@mail.ru

Алина Ильдаровна Гисмиева

Казанский государственный медицинский

университет, Казань, Россия

Эл. почта: alina.gismieva@icloud.com

recommendations. Moscow: Ministry of Health of the Russian Federation; 2016. 76 p. (In Russ.)

26. Patofiziologiya: uchebnik = Pathophysiology: a textbook. Moscow: GEOTAR-Media; 2007. 496 p. (In Russ.)

Information about the authors

Peter A. Korotkov

Volga State University of Technology,

Yoshkar-Ola, Russia

E-mail: korotkovpa@volgatech.net

Aleksey B. Trubyanov

Mari State University, Yoshkar-Ola, Russia E-mail: true47@mail.ru

Anastasiya A. Avdeeva

Maristat,

Yoshkar-Ola, Russia E-mail: nastya_87_02@mail.ru

Alina I. Gismieva

Kazan State Medical University, Kazan, Russia

E-mail: alina.gismieva@icloud.com

i Надоели баннеры? Вы всегда можете отключить рекламу.