Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Лобоцкая Н.Л. Основы высшей математики учеб. пособие для студентов высш. мед. учеб. заведений

.pdf
Скачиваний:
243
Добавлен:
25.10.2023
Размер:
13.51 Mб
Скачать

Обработку серии измерений следует проводить в следующем порядке.

1)Определить среднее арифметическое (выборочное).

2)Найти среднюю квадратическую отдельного измерения.

3)Определить наибольшую возможную ошибку отдельного

измерения

SH a 6 и убедиться,

что среди результатов

измерений

нет таких,

которые отличались бы от среднего арифметического

более, чем на SH a i l 6 . Если

такие результаты

имеются, их сле­

дует отбросить и начать обработку сначала.

 

 

4) Определить

среднюю

квадратическую

ошибку

среднего

ар ифметического.

 

измерения А х при заданной вероят­

5) Определить

точность

ности а и числе степеней

свободы k.

 

 

Пример 1. Десять измерений диаметра капилляра в стенке легочных альвеол дали следующие результаты: xv мк = 2,83; 2,82; 2,81; 2,85; 2,87; 2,86; 2,83; 2,85; 2,83; 2,84.

Провести обработку серии этих измерений. Р е ш е н и е . 1) По формуле (1)

- _

2,83 + 2,82 +

2,81 +

2,85 +

2,87 +

2,86 + 2,83

+

 

 

 

10

 

 

+

2,85 +

2,8^3 +

2,84

=

^

 

2) По формуле (2) находим среднюю квадратическую ошибку отдельного измерения

5 = | /

(2,83 - 2,84)2 + (2,82 - 2,84)2 + • • • + (2,84 - 2 , 8 4 ) ~ _ Q 0 1 g

3) По формуле (4)

5„аиб = 3-0,019 = 0,057.

4) Средняя квадратическая ошибка среднего арифметического по формуле (3) равна:

S~x

 

=

= - M i = 0,006.

 

8

Vn

]Ло

 

 

 

5) Точность измерения по формуле (5)

АХ = ta,kSxB.

При а = 95% и k = 9 коэффициент Стьюдента t по табл. 11 Приложений t = 2,262 и

Ал: = 2,262 • 0,006 = 0,013. Таким образом, диаметр альвеол

х = (2,84 ± 0,01) мк.

6)

Доверительный интервал, согласно

формуле (6),

заключен

в пределах от 2,83 до

2,85.

Это

значит,

что с

вероятностью

Р = 9 5 %

генеральная

средняя не

выйдет

за пределы

интервала

2,83

и 2,85.

 

 

 

 

 

 

 

При

многократных

непосредственных

измерениях

той

или

иной

физической величины

приборами результаты

могут

полу­

чаться одни и те же. Это значит, что случайные ошибки изме­

рения

меньше или

равны

ошибке

регистрирующего

прибора.

В этом случае критерием точности

измерения является

цена

наименьшего деления

шкалы

прибора.

Тогда

ошибка

измерения

берется равной половине цены наименьшего

деления шкалы

при­

бора.

Например,

при цене

деления

шкалы термометра в

0,1°

ошибка измерения

температуры берется 0,05°.

 

 

§142. ОПРЕДЕЛЕНИЕ ОШИБОК КОСВЕННЫХ ИЗМЕРЕНИЙ

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

непосредственно измеряемые, — ее аргументами ъ

х2, . ..

, хп).

Вид функциональной зависимости

определяется формулой,

свя­

зывающей эти величины. В общем

случае у = Цхъ

х%, ...

, хп).

Определение ошибок косвенных измерений проводится в по­ рядке, указанном выше (§ 141).,

1) Путем подстановки в расчетную формулу средних значе­ ний величин, непосредственно измеренных, вычисляют среднее значение косвенно измеренной величины

У — Ї (ХЪ х2> • • • > ХП)-

2) Вычисляют дисперсию косвенно измеряемой величины по формуле

где

^ — частная производная функции у = / ь х2,

. . . , хп);

Sx;

І

измерен­

— средняя квадратическая ошибка непосредственно

ных

величин.

 

3) Вычисляют среднюю квадратическую ошибку

косвенно

измеренной величины

 

4) Вычисляют среднюю квадратическую ошибку среднего значения косвенно измеренной величины для п измерений

Sy

5) Точность значения косвенно измеренной величины вычис­ ляют по формуле

Л У = ta,k Sy = ta,k

—7=-

\

V п

6) Доверительный интервал косвенно измеренной величины будет равен:

у — Л у < р <~у 4-А у.

(9)

Пример 2. С помощью колориметра проведено измерение концентрации сх неизвестного раствора путем сравнения с рас­ твором известной концентрации по формуле

Сх — С ° "Т"'

ах

где d0 и dx толщины слоев, одинаково поглощающих моно­ хроматический свет.

В пяти опытах получены результаты: d0 = 5,7 мм, Sd(j = = 0,15 мм, dx = 8,5 мм, Sdx = 0,18 мм. Концентрация раствора с„ = 2% (примем за точное число). Определить среднее значе­ ние концентрации сх и доверительный интервал с доверительной вероятностью а = 0,95.

е ш е н и е . сг = Сп

dx

=

2 • -£-=-

=

1,34 «/о.

 

 

х

 

 

 

8,5

 

 

 

 

Согласно

формуле

(7),

 

 

 

 

 

 

51

 

д 1

do \

2 SI 4-

д '

do V

 

 

 

 

c'kd'k

V 0

/J

 

ddu

 

 

 

 

+ _ddx

 

 

 

 

a*

 

 

 

 

 

 

 

Согласно

формуле

(8)

 

 

 

 

 

51 =

 

 

 

 

 

СЬ

C2

І

c 0^0

 

 

 

 

 

5 5 -

<4 H

dt

 

 

= Yf l/(0,15)2 + |g£-(0,18)2 = 0,04»/o

При

а = 0,95, п = 5, k = 4

по табл.

11

Приложений

найдем

t =

2,776

 

 

 

 

 

 

 

 

 

 

 

Д с, = 2,776 • -Щг-

=

0,05,

сх

= (1,34 ±

0,05) %.

 

 

 

К 5

 

 

 

 

 

 

 

 

Доверительный

интервал

(согласно

выражению

(9)) заключен в

пределах от 1,29 — 1,39%.

 

 

 

 

 

 

 

 

§ 143. ПОНЯТИЕ О СПОСОБЕ НАИМЕНЬШИХ КВАДРАТОВ

Обычные приемы графического

изображения

функциональной

зависимости между

переменными

величинами у =

f(x)

широко

применяют при обработке

результатов

измерений.

 

 

Результаты

измерений

функционально

зависимых

величин

Xi и уі, хг и у2; xs

и уз

и т. д. содержат

ошибки

измерений,

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

Чтобы математически более или менее точно выразить на­

блюдаемую

закономерность, нужно

по возможности устранить

случайные

отклонения

и найти наиболее характерные точки.

Это достигается путем

выравнивания

или сглаживания экспери­

ментальных

данных.

 

 

Известны три способа выравнивания или сглаживания: 1) гра­ фический; 2) арифметический (скользящей средней); 3) алгебраи­

ческий (наименьших квадратов).

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

Пусть производится

опыт, цель которого — исследование

зависимости некоторой

физической величины у от величины X,

например, зависимости температуры кипения воды от давления,

зависимости

веса

людей

от возраста

и др. Предполагается,

что

величины

 

х

и у

связаны

функциональной

зависимостью

у —

=

/ (х),

которую

и

нужно

определить

по

опытным данным.

В

результате

опыта

получены

следующие

данные:

 

 

 

 

 

 

 

 

 

X-t

Х\

Х%

Хз ...

хп

 

 

 

 

Как

подобрать

УІ

Уі

Уг

Уз

• • •

Уп

y f(x),

чтобы

 

математическую

формулу

f(Xi),

f(x2),

f(x3),

• • • ,

/(*:„)'наиболее близко

подходили

к

зна­

чениям уъ

уг,

уз

, . .. ,

уп,

т. е. чтобы

теоретическая

кривая

наилучшим

образом воспроизводила бы зависимость у от х? При

заданном

 

типе зависимости

у = f (х),

полученном из теоретичес­

ких

соображений

и по расположению экспериментальных

точек,

метод наименьших квадратов дает возможность выбрать число­ вые параметры кривой так, чтобы кривая y = f(x) наилучшим образом отображала экспериментальные данные.

Согласно методу наименьших квадратов, сумма квадратов отклонений экспериментальных точек от сглаживающей кривой

 

 

 

 

 

 

у*ах+Ь

у =

f (х)

(теоретической)

должна

 

 

 

 

 

 

обращаться в минимум, т. е. сум­

 

 

 

 

 

 

 

 

ма

 

квадратов

отклонений наблю­

 

 

 

 

 

 

 

 

даемых

значений

от

f(xt)

должна

 

 

 

 

 

 

 

 

быть наименьшей

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 [ у г - / ( \ ) Р

=

гшп.

 

 

\0C2Xj

 

 

 

Хп

 

 

Пусть выбран общий вид функ­

 

 

 

 

 

ции

f(x,

а, Ь,

с, . . . ) ,

зависящей

 

 

 

 

 

 

 

 

 

 

Р и с .

138

 

 

от

нескольких

числовых

парамет­

параметры

a,

b,

 

с,

 

 

ров

а, Ь, с, . . .

Требуется

выбрать

 

так,

чтобы

выполнялось

условие

 

 

U

=

'2lyl

— f(xl,

 

а,

Ъ,

с,

...)]» ==min.

 

 

(10)

Параметры а,

Ь,

с, ...

можно

найти из

условия

минимума:

 

 

=

0

 

 

 

 

 

 

 

 

 

(И)

 

 

 

да

~

 

' db

 

 

 

 

 

 

 

 

 

Из

системы

уравнений

(11)

для

конкретной

функции у =

= f(x,

а,

Ь,

с,

. . . )

можно

найти

коэффициенты

а,

Ь,

с, ...

Рассмотрим применение способа наименьших квадратов для случая линейной зависимости между двумя величинами. Пусть в опыте зарегистрирована совокупность значений (хг; t/t), при­ веденная в виде экспериментальных точек на рис. 138. Требу­ ется по методу наименьших квадратов подобрать параметры в линейной функции у = ах + Ь, изображающей данную экспери­ ментальную зависимость.

По формуле (10) можно написать ,

V = (У! — axj. — bf + . . . + (yn — axn — bf = min. Значения параметров а и b определим из условия (11).

Ж^Ж

 

a X l

~ 6

) 2

+ ^ 2

— ax2

— b)2

+

---Jr

 

+

п-ахп-Ь)*],]

 

 

 

 

 

 

 

(12)

 

 

 

 

 

 

 

 

 

 

 

 

F =

[(уі~~

a x

 

b¥

+ (У2 — ах2

— Ь)2

+

. . . -t

 

"+I

(Уп axn

ibf].

 

 

 

 

 

 

 

Продифференцируем

выражения (12)

 

 

 

2(уі

ах1 — Ь)(—х1)

+

2(у2

ах2

— Ь)(—х2)

+ ...

+

+ 2(yn-axn-b)(-xn)

 

 

 

=

0,

 

 

 

 

2(yl-ax1-b)(-

 

 

\)

+

2(y2-ax2-b)(-

 

1) + . . .

+ .

+ 2(уп-ахп-Ь)

( - 1 ) = 0.

Раскрыв скобки и произведя суммирование, получим

2 * 4

& аїх?—

bixt

= 0, Іуі

aZXi bn^O.

/=1

j= l

і'=1

і=1

г'=1

Из системы уравнений (13) находим значения а и Ь:

л2 xi 2 г/;

1=1

 

t=i <=1

,

6

і=1

г=1

г=1 1=1

П

 

•<2*г)*

л

2 ж?

 

п 2

xf

 

 

 

(=1

'

(=1

 

 

 

г=1

 

 

 

 

 

 

(13)

(14)

Подставив значения а и b в уравнение прямой, получим урав­ нение г/ = ах + &, график которого и будет расположен наи­ более близко к экспериментальным данным.

Пример 3. В «Основах

химии» Д. И. Менделеева приводятся

данные о растворимости азотнокислого

натрия

NaN03

в зави­

симости от температуры воды [32].

 

(у),

 

 

 

 

Предполагая, что количество NaN03

 

 

которое

растворя­

ется в 100 частях воды,

зависит линейно от

температуры (л:)

раствора, найти параметры а и b в формуле

 

у = ах + b по ме­

тоду наименьших квадратов. Результаты приведены ниже.

№ п. п.

xt

 

У і

 

 

 

 

 

Х;УІ

1

0

 

66,7

 

0

 

 

 

0

2

4

 

71,0

 

16

 

 

284,0

3

10

 

76,3

 

100

 

 

763,0

. 4

15

 

80,6

 

225

 

 

 

1209,0

5

21

 

85,7

 

441

 

 

 

1799,7

6

29

 

92,9

 

841

 

 

 

2694,1

7

36

 

39,4

 

1296

 

 

3578,4

8

51

 

113,6

 

2601

 

 

 

5793,6

9

68

 

125,1

 

4624

 

 

8506,8

2

234

 

811,3

10144

 

 

 

24628,6

Р е ш е н и е .

Согласно

формулам (14)

 

 

 

 

а =

9 • 24628,6 — 234 -811,3

= 0,87,

 

9 • 10144 —(234)2

 

 

Ь =

10144 • 811,3 — 234 • 24628,6

=

67,5.

 

9 • 10144 — (234)2

 

 

 

 

 

 

 

 

 

 

 

Следовательно, зависимость у от х имеет вид: у = 0,87л; + 67,5.

§ 144. НУЛЕВАЯ ГИПОТЕЗА. СРАВНЕНИЕ ВЫБОРОЧНЫХ СРЕДНИХ

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

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

Так, нулевая гипотеза предполагает, что полученная.в опыте разность случайна. Исследование заключается в испытании ну­ левой гипотезы в отношении вероятности того, что она верна.

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

Все статистические показатели, а следовательно, и различия между ними, характеризуются определенными уровнями значи­ мости, поэтому отбрасывание нулевой гипотезы должно быть связано с принятием определенного уровня значимости (§ 137). Пусть нулевая гипотеза отбрасывается, например, при уровне значимости Р = 0,01 (доверительная вероятность 0,99). Если вероятность разности между статистическими показателями ниже 0,99 (например 0,97; 0,91) или Р > 0,01, то нулевая гипотеза принимается, различия являются случайными. Ее надо считать правильной по крайней мере до тех пор, пока новые данные не

дадут

возможности

ее

опровергнуть,

доказав, что существую­

щие различия не являются чисто случайными.

 

 

 

 

Однако

и в

том

случае,

когда

нулевая

гипотеза

считается

опровергнутой,

остается

какая-то возможность, что

она

в

дей­

ствительности верна. При уровне значимости Р =0,01

есть

один

случай

из

100, при котором отбрасывание нулевой

гипотезы

было бы ошибкой. Если

достигнут

уровень значимости

0,001,

то уверенность

в

том,

что

нулевая

гипотеза

действительно

правильно

отвергнута, резко возрастает (лишь один

случай на

тысячу, что она верна).

 

 

 

 

 

 

 

 

 

Нулевая

гипотеза

в

биологических

исследованиях

принима­

ется,

если

уровень значимости полученных

данных

Р >

0,05.

Если

Р < 0 , 0 1 , то нулевую

гипотезу

можно

отбрасывать.

В случае, когда уровень значимости

исследований лежит в

пределах

Между 0,05

и 0,01,

необходимы

дополнительные

опы­

ты, чтобы решить,

следует

ли отбрасывать

нулевую

гипо­

тезу.

 

 

 

 

 

 

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

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

ны

на непосредственном

сравнении

варьирующих

величин

или

их

порядковых

значений.

 

 

 

 

 

 

 

 

 

 

Параметрическим критерием является критерий t

(тест),

зна­

чения которого

оцениваются

по таблицам (табл.

8

Приложений

для

объема

выборки п >

30

и табл. 9 Приложений

для объема

выборки га< 30).

 

 

 

 

 

t.

 

 

 

 

 

 

Рассмотрим

случаи

применения

критерия

 

 

 

 

 

 

1. Известно

среднее

значение

генеральной

совокупности р..

Определить,

имеются

ли существенные отклонения

 

между

сред­

ним

выборочным значением хв из объема выборки

п и

средним

значением генеральной

совокупности

р. Нулевая

гипотеза

соот­

ветствует условию хв = р. Оценка

разности

между

хв и р. мо­

жет

быть осуществлена

при помоши теста t

по формуле

 

 

 

 

 

 

t

 

Хя

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример

4. Вес тела

измерен

с

высокой

степенью

точности

на аналитических весах

и равен

[і = 100 г.

При

взвешивании

его

на обыкновенных

весах

п =

10 раз получили

значение

 

 

 

 

хв =

100,26 г,

S-B = 0,078 г.

 

 

 

 

 

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

 

~

0

 

,

хв

=

100,26 — 100

 

 

Р е ш е н и е . Значение

теста t = —^=

 

§~$T&

=

 

 

 

 

 

 

 

 

 

=

3,33. Пользуясь

табл.

10 Приложений

для k — п— 1 = 9 и

^ =

3,33,

находим

уровень

значимости

Р < 0 , 0 1 ,

следовательно,

нулевая гипотеза о достаточной точности обычных весов отвер­ гается.

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

2

х в 1

х

в 2

d

~

У Si

+

S 2

~ sd'

 

XBl

 

ХВ2

 

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

При сравнении двух групп с малыми объемами п, в особен­ ности с неодинаковыми п, ошибка разницы определяется по формуле

I ~ПХ ~ Щ ~

О

І /

t=\

 

i ^ l

 

« 1

+ »

2

d

f

 

( r t ! - l )

+

( n 2

- 1)

*

щпг

 

или с использованием

тождества

(14) гл. XX

 

 

 

/

Ъх*а

-

«і (*В 1 )2 +

SX? -

пг в 2 )2

 

 

 

г > " ~ ^

 

 

(«і — 1) Н - ( « 2 — 1)

 

И і П 2

^

Уровень значимости в этом случае определяется по табл. 10 Приложений при числе степеней свободы

(«1 — 1) + («2 — 1) = «1 + п2 — 2.

Пример 5. Для сравнения активностей двух аналогичных препаратов были проведены две серии испытаний, в каждой из которых использовалось семь животных [20]. Результаты при­ ведены ниже.

х п

35

83

53

60

71

62

39

хв1

= 57,6

х?

1225

6889

2809

3600

5041

3844

1521

2

=

24929

h

 

 

 

 

 

 

 

 

 

 

Хц

60

63

99

95

78

85

72

Xqo

= 78,9

 

3600

3969

9801

9025

6084

7225

5184

2

=

44888

Следует ли

отбрасывать

нулевую

гипотезу?

 

 

 

Р е ш е н и е .

в1)2

= (57,6)2 = 3317,8; (хв 2 )2

= (78,9)2 = 6225,2.

Согласно формуле (15),

24929 — 7 • 3317,8 +

44888 — 7 • 6225,2

7 + 7 _ о 4 f i

fiXfi '

'

7 • 7

-

 

Критерий

j .

хві

хв2

78,9 — 57,6

о с о

1

~~

Sd

_

8,46

_

Число степеней

свободы

k = 7 +

7 — 2 — 12. По табл. 10 При­

ложений для k = 12 и t — 2,52 Р > 0,05, но Р < 0,01. Отсюда можно заключить, что вероятность того, что срав­

ниваемые препараты имеют различную активность, превышает 95% (но меньше чем 99%). Нужны дополнительные опыты, чтобы решить, следует ли отбрасывать нулевую гипотезу.

§ 145. ТОЧНОСТЬ РЕЗУЛЬТАТА И ЧИСЛО ИСПЫТАНИИ. КОЭФФИЦИЕНТ ВАРИАЦИИ

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

Желаемая точность (обозначим ее через Л) это возможное при принятой вероятности отклонение среднего выборочного хв от среднего значения генеральной совокупности а.

Лв

у п

откуда

д2

Для определения объема выборки по формуле (16) при при­ нятой вероятности значение t определяем из табл. 8 Приложе­ ний. Значение А берется заранее в зависимости от условий задачи. Стандарт отклонения S берется из известных исследо­ ваний или выбирается его наибольшее возможное значение.

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

Соседние файлы в папке книги из ГПНТБ