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

Методические рекомендации предназначены для врачей СЭС и МСЧ.

Составители: Кемеровская облСЭС (В.А. Зенков), Кемеровский медицинский институт (проф. А.П. Михайлуц).

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

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

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

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

1.1. В соответствии с приказом МЗ РСФСР N 74 от 18.03.88 "О переводе санитарно-эпидемиологической службы на работу в условиях эксперимента" среди главных задач санэпидслужбы в новых условиях названы "изучение, выявление и оценка причинно-следственных связей" между санитарным состоянием условий труда, быта и отдыха, окружающей среды, с одной стороны, и здоровьем населения - с другой стороны.

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

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

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

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

1.5. Следует исходить из того, что в большинстве случаев с повышением степени вредности и опасности условий труда, его тяжести и напряженности возрастает определенность их как причины заболеваний; а наличие профессиональных заболеваний является абсолютным доказательством действия условий труда как причины.

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

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

2.1. При подготовке исходных данных СЭС целесообразно, пользуясь гигиенической классификацией N 4137-86 "Гигиеническая классификация труда (по показателям вредности и опасности факторов производственной среды, тяжести и напряженности трудового процесса)", ранжировать условия труда на 5 классов: оптимальные (1), допустимые (2), вредные и опасные 1 степени (3), вредные и опасные 2 степени (4), вредные и опасные 3 степени (5), а тяжесть и напряженность труда - на 4 общепринятые класса.

2.2. Рекомендуется условия труда вредные и опасные 1 - 3 степени (3 - 5 класса), труд тяжелый (3) и очень тяжелый (4), напряженный (3) и очень напряженный (4) относить к потенциальным причинам заболеваний работающих.

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

2.3. Как следствие, действия причин на заболеваемость работающих целесообразно принимать в настоящее время общепринятые в системе учреждений здравоохранения (МСЧ) показатели:

а) число случаев профессиональных заболеваний на 1000 работающих;

б) суммарные уровни (% болевших лиц, число случаев и число дней на 100 работающих) заболеваемости с временной утратой трудоспособности (ЗВУТ) <*>, а также уровни ЗВУТ по болезням нервной системы и органов чувств (6 кл.), сердечно-сосудистой системы (7 кл.), органов дыхания (8 кл.), органов пищеварения (9 кл.), мочеполовых органов (10 кл.), кожи и подкожной клетчатки (12 кл.), костно-мышечной системы (13 кл.);

--------------------------------

<*> Н.В. Догле, А.Я. Юркевич "Заболеваемость с временной утратой трудоспособности (методы изучения)". М., 1984.

в) суммарные уровни (число случаев на 1000) хронических заболеваний, а также уровни хронических болезней по названным выше классам заболеваний.

2.4. Выбор конкретных производства, цеха, профессии для анализа причинно-следственных связей рекомендуется делать, первоначально ориентируясь на уровни ЗВУТ (табл. 1).

При этом суммарные уровни ЗВУТ коллектива, квалифицируемые как "очень высокие", "высокие", "выше средних", "средние" (табл. 1), являются основанием для установления зависимостей от условий труда, если последние относятся к 3 - 5 классам вредности и опасности или 3 - 4 степени тяжести или 3 - 4 степени напряженности.

2.5. Если трудовой коллектив состоит преимущественно из женщин, то в дополнение к показателям заболеваемости, указанным в п. 2.3, целесообразно получить данные о % беременных-работниц, с различными формами патологии течения беременности и исходов родов.

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

- по возрасту (1 - до 40 лет, 2 - старше 40 лет);

- по полу - (1 - мужчины, 2 - женщины);

- по образу жизни (1 - полезные привычки: физкультура, закаливание; 3 - вредные привычки; курение, алкоголь; 2 - отсутствие как полезных, так в вредных привычек);

- по загрязнению окружающей среды (1 - концентрации вредных веществ в атмосфере или питьевой воде на уровне ПДК; 3 - концентрации вредных веществ в атмосфере или питьевой воде выше ПДК в 5 и более раз; 2 - концентрации вредных веществ в атмосфере или питьевой воде выше ПДК до 5 раз).

2.7. На предприятии (производстве, цехе) необходимо в дополнение к основной изучаемой группе, подбирать 2 - 3 профессиональные группы работающих, включая контрольную с градуированными на классы условиями, тяжестью и напряженностью труда. В качестве контрольной группы рекомендуется брать профессии, имеющие допустимые условия труда, выполняющие труд средней тяжести и умеренно напряженный. Контролем могут быть показатели заболеваемости работающих в целом по отрасли промышленности города, области.

3.1. Для нахождения количественных зависимостей и их статической достоверности в причинно-следственных отношениях "условия труда - заболеваемость работающих" рекомендуется применять мини, макро ЭВМ, включая персональные компьютеры, с использованием программ стандартизации, вычисления критерия Стьюдента, весовых индексов факторов, дисперсионного, корреляционного и регрессионного анализа. <*>

--------------------------------

<*> Даются в приложениях.

3.2. При установлении количественных зависимостей между ЗВУТ и условиями труда, получив исходные данные (раздел 2), целесообразно последовательно (рис. 1) применить следующие методы статистической обработки: а) нахождение весовых индексов факторов (программа "НИПВИ"); б) стандартизация показателей (программа "СТАНДК"); в) вычисление достоверности различий показателей ЗВУТ на основе критерия Стьюдента (программа "ВУТ"); г) дисперсионный анализ (программа "ГАЛ"); д) коррелляционный - регрессионный анализ (программа "КА").

3.2.1. Определение весовых индексов факторов (Шиган Е.И. "применение современных математических методов и вычислительной техники для изучения здоровья населения". М., 1973) видно на следующем примере:

В цехе работает 280 человек. Из них в 198__ г. имели заболевания органов дыхания с ВУТ 80 человек, у 200 - болезни органов дыхания с ВУТ не зарегистрированы. Предположим, что на заболеваемость с ВУТ могли влиять такие факторы, как условия труда, возраст, пол и загрязнения атмосферного воздуха в местах проживания работающих. Градации факторов и схема определений их весовых индексов представлены в табл. 2, из которой видно, что наибольшее значение (8.1) имеет весовой индекс такого фактора, как класс опасности условий труда. На втором месте по вероятности влияние на возникновение заболеваний органов дыхания стоит фактор - загрязнение атмосферного воздуха в районе проживания работающих (весовой индекс 3.6). Наименьшее влияние оказывает фактор - возраст работающих (весовой индекс 1,8).

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

3.2.2. Выбрав в качестве потенциальной причины повышения заболеваемости условия труда, необходимо устранить влияние на нее возраста и пола. Это делается путем стандартизации, которую из-за трудоемкости целесообразно проводить на персональных компьютерах с использованием программы "СТАНДК". Программа дает стандартизованные по возрасту и полу показатели заболеваемости в изучаемых и контрольной группах.

3.2.3. Для определения статистической достоверности различий в заболеваемости между отдельными профессиональными группами работающих по проценту болевших лиц, числу случаев и дней ЗВУТ рекомендуется использовать программу "ВУТ". Программа запрашивает данные об абсолютных количествах болевших лиц, случаев и дней ЗВУТ в сравниваемых группах и дает ответ о достоверности различий (по критерию Стьюдента) заболеваемости по этим показателям.

3.2.4. Если статистически достоверных различий в заболеваемости работающих в отдельных профессиональных группах в зависимости от выбранных факторов (факторов-причин и факторов-условий реализации причин) не установлено, то процедуру анализа причинно-следственных связей между условиями труда и ЗВУТ стоит прекратить.

3.2.5. При достоверных различиях (по критерию Стьюдента) в показателях заболеваемости сравниваемых групп устанавливается с помощью дисперсионного анализа (на персональном компьютере по программе "ГАЛ") доля влияния в % и статистики как достоверность влияния выбранного фактора на заболеваемость.

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

Предлагаемая программа ГАЛ составлена на основании математического алгоритма (Н.В. Догле, А.Я. Юркевич "Заболеваемость с временной утратой трудоспособности", 1984), запрашивает только абсолютные числа работающих и болевших в группах, а выдает все статистические параметры.

3.2.6. Для установления силы связи и ее статистической достоверности между факторами-причинами и показателями заболеваемости рекомендуется использовать корреляционный регрессионный анализ (программа "КА" для персонального компьютера). При наличии значительной силы и ее статистической достоверности необходимо составить уравнение регрессии y = a_о + a₁x, где y - уровни ЗВУТ (% болевших, число случаев ЗВУТ на 100), x - градация фактора-причины (класс опасности условий труда; степень загрязнения атмосферного воздуха), a_о, a₁ - коэффициенты уравнения регрессии.

Программа запрашивает данные по y и x, и выдает коэффициент корреляции, критерий Стьюдента, коэффициенты регрессии.

3.2.7. Последовательно применяя рекомендуемые методы статистической обработки материалов при изучении причинно-следственных отношений, можно:

- ранжировать все факторы-причины и факторы-условия по влиянию на заболеваемость;

- получить стандартизированные показатели;

- определить статистическую достоверность различий в заболеваемости в зависимости от действующих факторов;

- установить долю влияния (в %) факторов на заболеваемость и ее статистическую достоверность;

- установить силу связи между факторами и заболеваемостью и ее статистическую достоверность;

- получить математические зависимости (уравнения) заболеваемости от действующих факторов.

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

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

Во всех случаях при применении эпидемиологических методов необходимо подбирать группы работающих (основную и контрольную) с одинаковым профессиональным стажем на предприятии (производстве) не менее 8 - 10 лет.

Анализ эпидемиологическими методами может проводиться по двум схемам "воздействие - заболеваемость" и "заболеваемость - воздействие".

3.3.1. При эпидемиологическом анализе по схеме "воздействие - заболеваемость" рассчитывается показатель относительного риска (ОР).

Например, на производстве 400 человек имеют стаж работы 8 - 10 лет. Из них 120 страдает хроническими заболеваниями костно-мышечной системы, 280 человек подобных хронических болезней не имеют. В группе рабочих с хроническими заболеваниями костно-мышечной системы 60 человек (или 50%) выполняют тяжелую работу с нахождением в вынужденном положении. Среди здоровых со стажем также 8 - 10 лет тяжелый физический труд при вынужденном положении тела имеют 56 человек (или 20%) относительный риск равен т.е. риск через 8 - 10 лет работы на данном производстве иметь хроническое заболевание костно-мышечной системы в 2,5 раза больше у лиц, выполняющих тяжелую работу при вынужденном положении тела.

3.3.2. Если эпидемиологический анализ проводится по схеме "заболевание - воздействие", то рассчитываются следующие показатели: показатель абсолютного риска (АР); показатель относительного риска (ОР); разность относительных рисков (РР). Например, в группе работающих (основной) подвергающихся воздействию повышенных концентраций аммиака, через 8 - 10 лет работы распространенность хронических заболеваний верхних дыхательных путей составляет 40 случаев на 100; в контрольной группе, не имеющей воздействия аммиака, при таком же стаже работы на предприятии распространенность хронических заболеваний верхних дыхательных путей равна 10 случаев на 100.

Относительный риск разность абсолютных рисков РР = 40 - 10 = 30 случаев на 100 работающих.

Относительный риск (ОР) характеризует связь между воздействием и заболеваемостью, разность абсолютных рисков (РР) говорит о том, к какому абсолютному повышению заболеваемости приводит профессиональное воздействие на данном производстве.

1 D=SYS(7,0)

10 PRINT "Программа стандартизации косвенной по возрасту и полу"

20 DIM M1(12), G1(12), M2(12), G2(12), M3(12), M4(12), G3(12), G4(12), Y1(12),

25 DIM Y2(12), X1(12), X2(12), C(12), H(12), K(50), A(12), B(12) Z1(12),

27 DIM Z2(12), L1(12), L2(12), L3(12), L4(12)

30 S1=0 \ S2=0 \S3=0 \ S4=0 \ S5=0 \ S6=0 \ S7=0 \ S8=0 \ P1=0 \ P2=0 \ P3=0\ P9=0 \ R1=0\

35 R2=0 \ R3=0 \ R4=0 \ R5=0 \ R6=0 \ R7=0 \ R8=0 \ K1=0 \ K2=0 \ K8=0 \ Q1=0 \ Q2=0

40 PRINT "Введи число возрастных групп"

50 INPUT N

60 FOR I=1 TO N

70 PRINT "Введи число мужчин контрольной гр. в следующем возрасте"

80 INPUT M1(I)

90 S1=S1+M1(I)

100 NEXT I

110 FOR I=1 TO N

120 PRINT "Введи число женщин контрольной гр. в следующем возрасте"

130 INPUT G1(I)

140 S2=S2+G1(I)

150 NEXT I

160 FOR I=1 TO N

170 print "Введи число мужчин изучаемой гр. в следующем возрасте"

180 INPUT M2(I)

190 S<...>=S3+<...>2(I)

200 NEXT I

210 for I=1 TO N

220 PRINT "Введи число женщин изучаемой гр. в следующем возрасте"

230 INPUT G2(I)

240 S4=S4+G2(I)

250 NEXT I

260 FOR I=1 TO N

270 E3(I)=M1(I)*100/S1

280 <...>(I)=G1(I)*100/S2

290 <...>(I)=M2(I)*100/S3

300 <...>(I)=G2(I)*100/S4

310 NEXT I

320 S5=S1+S2

330 SG=S3+S4

340 FOR I=1 TO N

350 PRINT "Введи число заболеваний у муж. в контроле в возрасте"

360 INPUT Y1(I)

370 S7=S7+Y1(I)

380 NEXT I

390 FOR I=1 TO N

400 PRINT "Введи число заболеваний у жен. в контроле в возрасте"

410 INPUT Y2(I)

420 S8=S8+Y2(I)

440 NEXT I

450 FOR I=1 TO N

460 PRINT "Введи число заболеваний у муж. в изучаем. гр. в возрасте"

470 INPUT X1(I)

480 P1=P1+X1(I)

490 NEXT I

500 FOR I=1 TO N

510 PRINT "Введи число заболеваний у жен. в изучаемой гр. в возрасте"

520 INPUT X2(I)

530 P2=P2+X2(I)

540 NEXT I

550 P3=S7+S8

560 P4=P1+P2

570 FOR I=1 TO N

580 A(I)=M1(I)+M2(I)

590 B(I)=G1(I)+G2(I)

600 P5=P5+A(I)

610 P6=P6+8(I)

620 NEXT I

630 P7=P5+P6

640 FOR I=1 TO N

650 C(I)=X1(I)+Y1(I)

660 H(I)=Y2(I)+X2(I)

670 P8=P8+C(I)

690 P9=P9+H(I)

700 NEXT I

705 R1=P<...>

710 FOR I=1 <...>

720 Z1(I)=C(I)*100/A(I)

730 Z2(I)= H(I)*100/B(I)

740 NEXT I

750 R2=P8*100/P5

760 R3=P9*100/P6

770 R4=R1*100/P7

780 FOR I=1 TO N

790 L1(I)=M1(I)*Z1(I)/100

800 L2(I)=M2(I)*Z1(I)/100

810 L3(I)=G1(I)*Z2(I)/100

820 L4(I)=G2(I)*Z2(I)/100

830 NEXT I

840 FOR I=1 TO N

850 R5=R5+L1(I)

860 R6=R6+L2(I)

870 R7=R7+L3(I)

880 R8=R8+L4(I)

890 NEXT I

900 K1=R5+R7

910 K2=R6+R8

920 K3=S7*R2/R5

930 PRINT "СТАНДАРТИЗИРОВАННЫЙ ПОКАЗАТЕЛЬ З-ТИ МУЖ. В КОНТРОЛЕ" <...>

940 K4=P1*R2/R6

950 PRINT "СТАНД. ПОКАЗАТЕЛЬ З-ТИ МУЖ. В ИЗУЧАЕМ. ГРУППЕ" <...>

960 K5=S<...>*R3/R7

970 PRINT СТАНД. ПОКАЗАТЕЛЬ З-ТИ ЖЕН. В КОНТРОЛЕ" <...>

980 K6=P2*R3/R8

990 PRINT "СТАНД. ПОКАЗАТЕЛЬ З-ТИ ЖЕН. В ИЗУЧАЕМОЙ ГРУППЕ" <...>

1000 K7=(S7+S8)*R4/K1

1010 PRINT "ОБЩИЙ СТАНДАРТ. ПОКАЗАТЕЛЬ З-ТИ В КОНТРОЛЕ; <...>

1020 K8=(P1+P2)*R4/K2

1030 PRINT "ОБЩИЙ СТАНДАРТ ПОКАЗАТЕЛЬ З-ТИ В ИЗУЧ. ГРУППЕ" <...>

1040 Q1=P3*100/K1

1050 PRINT "СТАНДАРТИЗИРОВАННОЕ ОТНОШЕНИЕ БОЛЕВШИХ В КОНТРОЛЕ" <...>

1060 Q2=P4*100/K2

1070 PRINT "СТАНДАРТИЗОВАННОЕ ОТНОШ-Е К БОЛЕВШЕМУ В ИЗУЧ. ГР." <...>

1080 END

1 U=SYS(7,0)

RUN

10 DIM C(30), Y(25,25), X(25,25), n(25,25), A(25,25), B(25,25), D(25,25), V(25,25)

20 DIM W(25), Z(20,20), P(25), Q(25), S1(25), S2(25)

30 S1=0 \ S2=0 \ E=0 \ G1=0

40 PRINT "Программа вычисления НИП и весовых индексов факторов"

50 PRINT "Введи общее число наблюдений"

60 INPUT N

65 PRINT "Введи общее число заболеваний"

70 INPUT M

80 PRINT "Введи число факторов"

90 INPUT F

100 FOR j=1 TO F

110 PRINT "Введи числа градаций каждого фактора"

120 INPUT G(j)

130 NEXT j

140 G1=M·100/N

145 FOR j=1 TO F

150 K=G(j)

160 FOR i=1 TO K

170 PRINT "Введи число заболеван. по градациям факторов"

180 INPUT X(i,j)

190 PRINT "Введи число здоров. по градациям факторов"

200 INPUT X(i,j)

210 S1(j)=S1(j)+Y(i,j)

220 S2(j)=S2(j)+X(i,j)

230 NEXT i

240 FOR i=1 TO K

250 A(i,j)=Y(i,j)·100/S1(j)

260 B(i,j)=X(i,j)·100/S2(j)

270 R(i,j)=A(i,j)·B(i,j)

280 PRINT "Коэффициент отношения правдоподобия"; "R(i,j)="; R(i,j)

290 INPUT L

300 NEXT i

310 S1(j)=0

320 S2(j)=0

330 P(j)=R(i,j)

340 Q(j)=R(i,j)

345 FOR i=1 TO K

350 IF P(j)<P(i,j) THEN P(j)=R(i,j)

360 IF Q(j)>n(i,j) THEN Q(j)=R(i,j)

370 W(j)=P(j)/Q(j)

380 PRINT "весовой индекс фактора"; "W(j)="; W(j)

390 INPUT L

400 NEXT j

410 FOR j=1 TO P

420 K=G(j)

430 FOR i=1 TO K

440 Z((i,j)=y(i,j)·100/((X(i,j)-Y(i,j))

450 V((i,j)=Z(i,j)/01

460 PRINT "НИП"; "V(i,j)="; V(i,j)

470 IMPUT L

480 NEXT i

490 NEXT j

500 FOR i=1 TO F

510 F=E+W(j)

520 NEXT j

530 FOR j=1 TO F

540 K=G(j)

550 FOR i=1 TO K

560 D(i,j)=V(i,j)·W(j)·01/B

570 PRINT "НИП "весовой индекс"; "D(i,j)="; D(i,j)

580 INPUT L

590 NEXT i

600 NEXT j

610 PRINT "Общий % заболеваний"; "01=", 01

620 END

1 U=SYS(7,0)

RUN

10 N1=0 \ N2=0 \ N3=0 \ D4=0 \ I1=0 \ I2=0 \ I3=0 \ D3=0 \ X1=0 \ X2=0 \ M1=0

20 M3=0 \ M4=0 \ T2=0 \ T3=0 \ C=0 \ D1=0 \ D2=0

30 PRINT "Программа оценки достоверности различий показателей ВУТ"

40 PRINT "Введи кол-во работающих в контрольной группе"

50 INPUT N1

60 PRINT "Введи количество болевших в контроле"

70 INPUT N2

80 PRINT "Введи количество случаев заболеваний в контроле"

90 INPUT N3

100 PRINT "Введи кол-во дней нетрудоспособности в контроле"

110 INPUT D4

120 PRINT "Введи количество работающих в изучаемой группе"

130 INPUT i1

140 PRINT "Введи количество болевших в изучаемой группе"

150 INPUT i2

160 PRINT "Введи количество случаев заболеваний в изучаемой группе"

170 INPUT i3

180 PRINT "Введи количество дней нетрудоспособности в изучаемой группе"

190 INPUT D3

200 X1=N2·100/N1

210 X2=N3·100/N1

220 M1= SQR((X1·(100-X1))/N1

230 Y1=i2·100/I1

240 Y2=i3·100/I1

250 M2=SQR((Y1·(100-Y1))/I1)

260 T1=(X1-Y1)/(SQR(M12+M22))

270 M3=0·(SQR(X2/I))

280 M4=10·(SQR(Y2/I1))

290 T2=(X2-Y2) / (SQR(M32+M42))

300 D1=D4·100/N1

310 D2=D3·100/I1

320 T3=(D1-D2) / (SQR(((D1·(36500-D1)) / N1) ((D2·(36500-D2)) / I1)))

330 O=((N2·(i1-i2))-((N1-N2)·i2)) / (SQR(N1·I1+(N2+i2)·((N1-N2) + (i1-I2))

340 PRINT "Процент болевших в контр. и изуч. группе"; "X1="; X1, "Y1="; Y1

350 PRINT "Число случаев на 100 в контр. и изуч. группе"; "X2="; X2, "Y2="; Y2

360 PRINT "Число дней на 100 в контр. и изуч. гр"; "D1="; D1, "D2="; D2

370 PRINT "Достоверность различий по % болевших", "T1="; T1

300 PRINT "Достоверность разл-й по числу сл. заб."; "T2="; T2

390 PRINT "Достоверность разл-й по числу дней"; "T3="; T3

400 PRINT "Коэффициент сопряженности"; "C="; C

410 END

1 U=SYS(7,0)

RUN

5 DIM N(100), M(100), L(200), H(200), P(200)

10 PRINT "Программа дисперсионного анализа"

15 PRINT "Введи число групп"

20 INPUT R

30 FOR i=1 TO R

40 PRINT "Введи число работающих в группе"

45 INPUT N(i)

50 PRINT "Введи число болевших в группе"

55 INPUT M(1)

60 NEXT i

70 FOR i=1 TO R

80 S1=S1+N(i)

90 S2=S2+M(i)

100 

110 S3=S3+L(i)

120 

130 P(i)=M(i)/N(i)

140 S4=S4+H(i)

150 NEXT i

160 S5=S3/B1

170 C1=S4-S5

100 C2=S2-S4

190 C3=S2-S5

200 B1=R-1

210 B2=S1-R

220 B3=S3-1

230 D1=C1/B1

240 D2=C2/B2

250 F=D1/D2

253 PRINT "Факториальная, случайная и общая дисперсии"

260 PRINT "C1="; C1, "C2="; C2, "C3="; C3

270 PRINT "Доля влияния организованного фактора"; "D1="; D1

275 PRINT "Доля влияния случайных факторов"; "D2="; D2

280 PRINT "Критерий Фишера"; "F="; F

230 END

1 U=SYS(7,0)

RUN

10 DIM N (250), X(250), Y(250), A(250), B(250), C(250)

20 S1=0 \ S2=0 \ S3=0 \ S4=0 \ S5=0

30 PRINT "Программа корреляционного регрессионного анализа"

40 PRINT "Введи число наблюдений"

50 INPUT N

60 FOR i=1 TO N

70 PRINT "Введи игрек"

80 INPUT Y(i)

90 PRINT "Введи икс"

100 INPUT X(i)

110 NEXT i

120 FOR i=1 TO N

130 S1=S1+X(i)

140 

150 S2=S2+A(i)

160 NEXT i

170 FOR i=1 TO N

180 S3=S3+Y(i)

190 

200 S4=S4+B(i)

210 NEXT i

220 FOR i=1 TO N

230 C(i)=X(i)·Y(i)

240 S5=S5+C(i)

250 NEXT i

260 

270 

280 PRINT "Коэффициент корреляции"; "R="; R

290 PRINT "Критерий Стьюдента";, "T=";, T

300 

310 K0=(S3/N)-((K1·S1)/N)

320 PRINT "КоэфФициенты регрессии";, "KI =";, KI

330 PRINT "Продолжать? Да - введи 1, нет - введи 5".

340 INPUT j

345 iF=5 THEN 380

350 FOR i=1 TO N

360 PRINT 

370 NEXT i

380 PRINT "Введи масштаб"

390 INPUT H

400 FOR i=1 TO N

410 PRINT TAB (ABS(KO+K1·X(1))·H); "*"

420 NEXT i

430 OPEN "LP" FOR OUTPUT AS FILE # 1

440 PRINT "Выдать на печать? Да - введи 1, нет - введи 5".

450 INPUT D

460 iF D=5 THEN 530

470 FOR i=1 TO N

480 PRINT # 1, "*"

490 NEXT i

500 PRINT # 1, "Коэффициент корреляции"; "R="; R

510 PRINT # 1, "Критерий Стьюдента";, "T=";, T

515 PRINT # 1, "Коэффициенты регрессии"; "K1=";, K1

520 CLOSE # 1

530 END