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

KursLektcijTeoriiVeroyat

.pdf
Скачиваний:
24
Добавлен:
16.04.2015
Размер:
444 Кб
Скачать
☆

Большое количество практических задач приводится к задачам однофакторного дисперсионного анализа.

Типичным примером является работа технологической линии в составе которой имеется несколько параллельных рабочих агрегатов.

На выходе имеют место какие-то детали. Эти детали по какому-то параметру можем контролировать.

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

Вопрос: Обусловлены ли эти отличия действием случайных факторов или имеет место влияние конкретного станка агрегата.

В данном случае фактор только один – станок.

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

Имеется m станков, т.о. имеется m совокупностей. Из этих совокупностей мы проводим выборки объёмом n. Так, что значение параметров i-той совокупности i: xi1 , xi2 ,..., xin .

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

i \

j

1

2

.

j

.

n

Ср. выбо-

 

 

 

 

 

 

рочное xi

 

 

 

 

 

 

 

 

1

 

x11

x12

.

x1j

.

x1n

x1

2

 

x21

x22

.

x2j

.

x2n

x2

.

 

.

.

.

.

.

.

.

i

 

xi1

xi2

.

xij

.

xin

xi

.

 

.

.

.

.

.

.

.

m

 

xm1

xm2

.

xmj

.

xmn

xm

Выдвигаем гипотезу Н0 заключающуюся в равенстве средних выборочных.

H0 : x1 = x2 = ... = xm

H1 : x1 ≠ x2 ≠ ... ≠ xm −влияниестанковзначимо

Гипотеза Н0 проверяется сравнением внутригрупповых и межгрупповых дисперсий по F критерию Фишера.

Если расхождение незначительно, то принимается гипотеза Н0, в противном случае гипотеза Н0 отвергается.

 

 

 

n

 

 

 

 

n

 

 

 

 

 

∑x1 j

 

 

 

∑xi j

x

=

j =1

;

x =

j =1

;

 

 

1

 

 

n

 

i

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

m

n

 

1

m

 

 

x =

mn

∑∑

x =

m

∑

x ;

 

 

 

 

 

 

 

 

 

 

 

i =1

j =1

ij

 

 

i =1

i

 

 

 

 

 

 

 

 

 

Далее находят сумму квадратов отклонений от общего среднего:

31

m n

 

m n

 

 

∑∑(xij − x)2 =

∑∑(xij − xi + xi − x)2 =

 

i =1 j =1

 

i =1 j =1

 

 

m n

 

m n

m n

(xi − x)2

= ∑∑(xij − xi

)2 + ∑∑(xi − x)2 + 2∑∑(xij − xi )2

i =1 j =1

 

i =1 j =1

i =1 j =1

 

1442443

1442443

14444244443

Q

 

Q1

0

 

m

Q1 = n∑(xi − x )2

i =1

Ноль потому, что стоит сумма от (xij − xi )(xi − x ) - сумма отклонений переменных одной

14243

0

совокупности от средней арифметической той же совокупности.

m

n

m

m

n

∑∑(xij − x)2 = n∑(xi − x)2 + ∑∑(xij − x)2

i =1

j =1

i =1

i =1

j =1

1442443

14243

1442443

 

Q

Q1

 

Q2

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

Величину Q1 – рассеяние по фактору.

Слагаемое Q2 – представляет собой сумма квадратов разностей между отдельными и средней соответствующей совокупности. Эта сумма называется суммой квадратов отклонений внутри группы.

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

отклонений от общей средней.

Получим оценки дисперсий: S 2 , S 2

, S 2 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

2

 

 

 

 

- дисперсия обусловленная влиянием фактора;

 

 

 

 

2

 

 

Q2

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

S2

=

 

 

 

 

 

= Sост

- остаточная дисперсия – влиянием случайных и других неучтённых

m(n −1)

факторов.

 

 

Q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S 2 =

 

 

 

- полная дисперсия.

 

 

 

 

 

mn −

1

 

S 2

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Далее формируем оценку различия между оценками

иS

:

 

 

 

 

 

 

 

 

 

Q1 (m −1)

 

 

1

2

 

 

S 2

 

 

Sф2

 

 

 

 

 

 

 

 

2

 

 

1

=

 

 

=

 

 

 

 

 

 

= F подчиняется распределению f

 

Фишера.

 

 

 

 

 

 

 

 

 

 

 

S22

 

 

Sост2

 

 

 

Q2

[m(n −1)]

н

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Выбираем уровень значимости α, или доверительной вероятности 1– α = Р и по таблице F- распределения с числом степеней свободы: к1 = m–1; к2 = m(n–1) находим критическое значение Fкр,α Фишера.

P{Fн > Fкр,α }=α P{Fн ≤ Fкр,α }= P =1−α

Сравнивая между собой Fн и Fкр,α мы делаем вывод насколько сильно влияние интересующего нас фактора на исследуемую случайную величину.

В этом и состоит идея дисперсионного анализа.

Однофакторный дисперсионный анализ обычно представляют в виде таблицы.

32

 

Компоненты

 

 

Оценки

 

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

Основной

дисперсии

 

дисперсии

 

свободы

Межгрупповая

S 2 =

 

Q1

 

= S 2

m - 1

фактор

дисперсия

 

1

 

m −1

 

 

 

ф

 

Случайные,

Внутригрупповая

2

 

 

Q2

 

 

 

 

2

 

неучтенные

дисперсия

S2

=

 

 

= Sост

m(n - 1)

m(n −1)

факторы

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Общая дисперсия

 

S 2 =

 

 

Q

 

 

 

 

mn - 1

 

 

 

 

mn −1

 

 

 

 

 

 

 

 

 

 

Основы регрессионного и корреляционного анализа.

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

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

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

X:

M y x = f (x)

Как правило исследуются такие виды статистической связи, при которых значение некоторой случайной переменной зависит в среднем от значений, принимаемых другой случайной переменной:

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

M

x

= f (x)= Y (x)

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

Такое

представление

зависимости

между

X, см

50

 

 

 

переменными X и Y называется полем корреляции.

 

 

 

Можно также построить таблицу корреляции.

 

40

 

 

 

 

Y, см

 

 

 

Проделывая операцию усреднения для всех тех

(Хi, Yi )

30

 

 

 

 

 

 

значений Х,

по которым есть экспериментальный

 

20

 

 

 

 

 

 

 

материал, приходим к тому, что облако исчезает и

 

10

 

 

 

 

 

 

y

получается набор точек,

представляющих

средние

 

30

50

70

значения. Соединяя эти точки,

получаем ломанную,

 

 

 

 

 

называемую эмпирической линией регрессии.

 

Связь между СВ характеризуется формой и теснотой связи.

 

 

 

 

Определение фориы связи и понятие регрессии.

Определить форму связи между СВ – значит выявить механизм получения зависимой случайной величины. При изучении статистических связей, форму связей характеризует функция регрессии:

M (Y X = x)= f (x) - зависимость условного МО Если св Х и Y зависимы, то МО их произведения:

33

y

 

 

x

 

M (xy)= M (x)M

 

= M (y)M

 

 

x

 

 

y

Регрессия св Y относительно Х определяется как:

 

 

 

 

 

 

M (Y

 

 

)=

+∞

 

y

 

 

 

 

 

 

 

X

= x

∫

yf

dy ,

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

−∞

 

 

 

 

где

 

y

 

 

- условная плотность вероятности по формуле Байеса:

f

dy

 

 

 

x

 

 

 

 

 

f (x, y )

 

 

f (x, y )

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

f

 

=

 

 

 

 

=

 

 

 

 

 

 

 

 

 

f (x)

 

 

+∞

 

 

 

 

 

 

 

 

 

x

 

 

 

 

∫ f (x, y )dy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

−∞

 

 

 

 

 

 

=

+∞

xf (x, y)dx - регрессия Х по Y.

 

 

 

 

 

 

M X

 

 

∫

 

 

 

 

 

 

 

 

Y = y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

−∞

 

 

 

 

 

 

 

 

 

 

 

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

Точность прогноза определяется дисперсией условного распределения:

σ 2 (Y X = x)= M {Y X = x − M (Y X = x)}2 = M (Y 2 X = x)− M 2 (Y X = x)

учитывая: σ 2 (x)= M (x − mx )2 = M (x)2 − M 2 (x)

Несмотря на важность функции регрессии, возможности ее практического использования ограничены, т.к. для ее вычисления необходимо знать аналитический вид двумерной функции {x,y}. Мы же, как правило, имеем выборку ограниченного объема.

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

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

величина Х: Y (x)= f (x) .

Определение кривой регрессии инвариантно закона совместного распределения св Х и Y.

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

Линейная регрессия (ЛР). Метод наименьших квадратов.

Линейная регрессия занимает в технике и теории корреляции особое место. Она обусловлена двумерным нормальным законом распределения СВ Х и Y:

Y (x)= a0 + a1 x , где

а0 и а1 – коэффициенты регрессии, х – независимая случайная величина

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

Строятся квадратичные формы:

n

Q = ∑(xi −ε)2 → min

i=1

xi – измеренное значение переменной,

ε - истинное или теоретическое значение этой величины.

Требуется, чтобы сумма квадратов отклонений измеренных значений относительно истинных была минимальна.

34

В случае линейной регрессии за теоретическое значение принимается значение Y (x), т.е.

ищется такая прямая линия с коэффициентами а0 и а1, чтобы сумма квадратов отклонений от этой линии была минимальна.

n

Q = ∑(yi − a0 − a1 x)2 ,

i=1

уi – измеренное значение переменной Y.

Минимальные квадратичные формы получают, приравнивая к нулю ее производные по а0 и

а1:

 

∂Q

= −2∑(y − a0 − a1 x)= 0

a 0 , a1 = const

 

 

 

 

 

 

 

∂a0

 

 

 

 

 

 

 

 

 

 

 

 

∂Q

 

n

n

n

 

= −2∑(y − a0 − a1 x)x = 0

 

 

 

∑a 0 = na0

∑a1 x = a1

∑x

∂a

 

=

0

=

0

=

0

 

1

 

 

i

i

i

na0 + a1 ∑x = ∑y

a0 ∑x + a1 ∑x2 =∑yx

a0

=

∑y ∑x2 − ∑x∑yx

a1

=

n∑yx -∑x∑y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n∑x2 −

∑x 2

 

 

n∑x2 −

∑x 2

 

 

 

 

 

 

 

 

Нелинейная регрессия (НР).

Форма связи между условными средними определяется уравнениями регрессии. В зависимости от вида уравнений можно говорить о ЛР или НР.

В общем случае эта зависимость может быть представлена в виде полинома степени k:

Y (x)= a0 + a1 x + a1 x2 +... + ak xk

Определение коэффициентов регресии производится по методу наименьших квадратов:

Q = ∑(yi −Y (x))2 → min

∂Q

,

 

∂Q

,

 

∂Q

,...,

 

∂Q

 

 

 

 

 

 

 

 

 

 

∂a

 

 

 

 

 

 

 

 

 

∂a

0

 

 

 

∂a

 

 

∂a

2

 

k

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂Q

= −2

∑

(y − a

 

− a x − a

 

x2

−... − a

 

xk )= 0

 

 

 

 

 

 

0

2

k

∂a0

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂Q

 

= −2

∑

(y − a

 

 

− a x − a

 

 

x2

−... − a

 

 

xk )x = 0

 

 

 

 

 

 

 

 

0

2

k

∂a1

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.......

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂Q

= −2

∑

 

(y − a

 

− a x − a

 

x2

−... − a

 

xk )xk = 0

 

 

 

 

k

 

 

0

2

∂ak

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

35

a0 n + a1 ∑x + a 2 ∑x 2 +... + a k ∑x k = ∑y

 

 

 

 

 

a0 ∑x + a1 ∑x2

+ a2 ∑x3 +... + a k ∑x k+1

= ∑yx

 

 

 

 

 

........

 

 

 

 

 

 

 

 

 

∑xk + a1

∑xk +1 + a2 ∑x k+2 +... + a k ∑x 2k = ∑yxk

a0

 

 

 

 

 

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

Измерение тесноты связи.

Если бы величина Y полностью определялась аргументом Х, все точки лежали бы на линии регрессии. Чем сильнее влияние прочих факторов, тем дальше отстоят точки от линии регрессии. В случае в) связь между Х и Y является более тесной.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

х

 

 

 

х

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а)

 

 

 

в)

 

 

За основу показателя, характеризующего тесноту связи, берется общий показатель

изменчивости дисперсии:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

σ y2 = M

[(Y − my )]2 = M [(Y +

 

(x)−

 

(x)− my )]2 =

 

 

 

 

 

Y

Y

 

 

 

 

 

 

= M [(Y −

 

(x))]2

+ M [

 

 

(x)− my ]2

+ 2M [(Y −

 

(x))(

 

(x)− my )]

 

Y

Y

Y

Y

 

1442443

 

1442443

144424443

 

 

 

 

 

 

σ 2 y

 

 

 

 

δ

2 y

x

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∑

 

 

 

 

 

 

 

y

 

 

 

y

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(Y −Y (x))p

x

= ∑p

x

 

−Y (x)∑p

x

= 0

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

σ y2

=σ y2

+δ y2

 

 

(*)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

x

 

 

 

 

 

 

 

 

σ y2

x

-

 

 

дисперсия переменной Y относительно теоретической линии

дисперсии,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

определяющей влияние прочих факторов на величину Y.

 

 

δ y2

 

-

 

условная

дисперсия,

характеризует

дисперсию теоретической линии

регрессии

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

относительно условной генеральной средней my. Именно она определяет влияние данного фактора (Х) на величину Y и может быть использована для оценки тесноты связи между величинами Х и

Y.

 

 

 

M (Y (x)− my

)2

 

δ y2

 

 

 

 

 

 

 

 

 

 

 

ηT , y x

=

 

 

 

 

 

=

 

x

 

- теоретическое корреляционное отношение.

 

 

σ 2

 

 

σ

2

 

 

 

 

y

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

что легко доказать, поделив (*) на σу2:

 

 

Изменяется от 0 до 1,

 

 

 

 

 

 

 

 

 

 

 

 

σ 2

 

δ 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 =

x

+

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

σ 2

 

σ

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

σ

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

σ 2

≤ σ 2

 

 

 

 

 

 

 

 

 

 

 

η

=1−

 

x

;

 

 

 

 

 

 

 

 

 

 

 

σ 2

 

 

 

 

 

 

 

 

 

 

 

T , y x

 

 

 

y x

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

1)

Если η

y

=1, то σ y2 = 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T ,

 

x

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Влияние прочих факторов отсутствует. Все распределение будет сконцентрировано на линии регрессии. В этом случае между Х и Y существует простая функциональная зависимость.

36

Sy(x)2

2)Если ηT , y x =0, когда Y (x)= my .

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

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

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

Определение эмпирического корреляционного соотношения.

SY2 = 1n ∑[Y −Yi (x)]2 + 1k ∑[Yi (x)−Y ]2 = S y2(x ) +δ y2(x)

y – измеряемое значение зависимой переменной n – общее количество измерений

Yi (x)- условное среднее (среднее значение зависимой переменной у в i-ом интервале св Х) k – общее количество интервалов

Y- среднее всей совокупности измерений

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

ymax

 

 

 

 

y

 

 

 

Yk(x)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Эмпирическая

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Yср

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

линия

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

регрессии

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ymin

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

х

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

хmin

 

 

 

хma x

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

δy(x)2 – характеризует дисперсию эмпирической линии регрессии относительно среднего всей совокупности, т.е. влияние исследуемого фактора на зависимую переменную Y.

η2 =

δ y2(x )

=1

−

S y2(x )

 

- Эмпирическое корреляционное соотношение

S y2

S y2

 

 

 

 

 

 

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

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

Свойства:

1.0 ≤ η2 ≤ 1

2.если η2 =1, все точки корреляционного поля лежат на линии регрессии – функциональная связь между Х и Y.

3.Если η2 =0 (когда ∑[Yi (x)−Y ]2 = 0 ), отсутствует изменчивость условных средних [Yi (x)],

эмпирическая линия регрессии проходит параллельно оси абсцисс – свзи между Х и Y нет. Эмпирическое корреляционное соотношение η2 завышает тесноту связи между

переменными и случайными величинами, причем тем сильнее, чем меньше число измерений, поэтому η2 рекомендуется использовать для предварительной оценки тесноты связи, а для

окончательной оценки – теоретическое корреляционное соотношение.

Коэфициент корреляции.

Рассмотрим случай вычисления теоретического корреляционного соотношения ηт2 , когда связь между случайными величинами Х и Y является линейной.

37

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(x)= a0 + a1 x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Такая форма связи между Х и Y имеет место в случае, когда случайные величины

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M

[

 

(x)

− my

]2

 

 

 

 

 

∑[

 

 

(x)−

 

 

]2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

η

2

=

 

Y

 

=

 

Y

Y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

т

 

 

 

 

 

σ y2

 

 

 

 

 

 

 

 

 

 

nσ y2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Подставив вместо Y и

 

 

их значения для случая линейной зависимости:

 

 

 

 

 

 

 

 

Y

 

 

 

 

= ∑y =

na0 + a1 ∑x

 

= a

 

 

+

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Y

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Y

(х)=а0 + а1х

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= ∑(a0 + a1 x − a0 − a1

 

 

)2

 

 

= ∑a12 (x −

 

)2

= a2

nσ X2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

η2

 

x

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

т

 

 

 

 

 

 

 

nσY2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

nσY2

 

 

 

1 nσY2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Заменим а1 ее значением, полученным из решения нормальных уравнений:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∑xy

 

∑x ∑y

 

2

 

 

 

 

∑xy

− ∑x ∑y

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

σ

2

n xy −

 

x y

 

 

σ

2

 

 

 

 

 

 

 

 

 

 

 

 

−

 

 

 

 

 

 

 

 

σ

2

 

 

n

n n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

x

 

 

 

∑ ∑ ∑

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

n

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

=

 

 

 

 

 

 

 

 

 

 

 

x

 

=

 

 

 

 

 

 

x

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∑x2

 

 

∑x 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

σ y2

 

n∑x2 − (∑x)2

 

 

 

 

σ y2

 

 

 

 

 

 

 

 

σ y2

 

 

 

 

 

σ x2σYx2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

−

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∑xy

−

∑x∑ y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r = a

 

 

σ x

=

n

 

 

 

 

n

 

 

 

=

∑(x − x)(y − y)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 σY

 

 

σ xσY

 

 

 

 

 

 

 

nσ xσY

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Коэфициент корреляции является частным случаем теоретического корреляционного соотношения ηт2 , когда связь между СВ является линейной. В этом случае r является показателем тесноты связи.

k x, y

=

∑(x − x)(y − y)

- выборочный корреляционный момент

n

 

ρx, y

=

∑(x − x)(y − y)

r = ρ(x, y)ρ(y, x)

nσ x

2

 

 

 

Выборочный коэфициент корреляции обладает свойствами:

1.r=0, если св Х и Y независимы

2.r < 1 - Для любых св Х и Y

3.

 

r

 

= 1 - Для случая линейной зависимости св Х и Y.

1

a1

> 0

 

 

r =

a1

 

 

 

 

 

 

−1

< 0

 

 

 

 

 

Коэфициент корреляции используется для оценки тесноты связи и в случае нелинейной зависимости между случайными величинами.

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

Если модуль коэфициента корреляции r = 0.8 ÷0.9 , то независимо от вида связи можно

считать, что она достаточно тесна, чтобы исследоват ее форму.

Двумерное нормальное распределение.

Его возникновение объясняется центральной предельной теоремой Ляпунова:

38

 

1

 

 

 

 

1

 

 

 

 

x − mx

2

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

− my

 

x − mx y − my

f xy (x, y) =

 

 

exp

−

 

 

2

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

2

−ρ

)

 

σx

 

 

 

 

σy

 

 

− 2ρ

 

σy

 

2πσx σy 1−ρ

 

 

2(1

 

 

 

 

 

 

 

 

 

 

σx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ρ – коэффициент корреляции. Х и У по отдельности распределены нормально (mx,σx) и

 

(my,σy).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В частном случае независимых СВ Х и У ρ=0:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

x

− m

 

 

2

y

− my

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

x

 

 

 

 

 

 

f xy (x, y) =

 

 

 

 

exp −

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2πσx σy

2

 

σx

 

 

 

 

σy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Исходные плотности одномерных нормальных распределений Х и У:

+∞

1

 

−(x − mx )

2

 

f x (x) = ∫−∞

− f xy (x, y)dx =

 

 

 

2πσx

exp

2

 

 

 

 

 

2σx

 

 

 

 

 

 

 

 

 

 

 

 

 

σy

 

 

 

 

 

 

 

 

 

 

 

− my

 

 

(x − mx )

2

f ( y / x) =

f xy (x, y)

=

1

 

 

 

1

y

−ρ σx

 

f x (x)

2πσy 1 −ρ2

exp −

2

 

 

σy 1−ρ2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Условное распределение – нормальное с условиями:

 

 

 

 

 

 

 

 

M ( y / x) = my + ρ

σy

(x − mx ) и σy / x = σy 1−ρ2 .

 

 

 

 

 

 

 

 

σx

 

 

 

 

 

 

 

 

Первое условие является уравнением функции регрессии.

 

 

 

 

 

 

M (x / y) = mx + ρ

σx ( y − my ) и σx / y = σx 1−ρ2 .

 

 

 

 

 

σy

 

 

 

 

 

 

 

 

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

Z=fxy(x,y) – трехмерная поверхность, сечения которой плоскостями XZ и YZ представляют собой графики плотности одномерных распределений.

Коэффициент множественной корреляции

 

 

R =

β r

yx1

+β

2

r

yx2

+... +β

n

r

yxn

=

* /

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

r12

...

 

r1n

 

 

 

 

 

= (−1)n+1 r21

1

...

 

r2n

 

 

 

 

 

 

 

 

 

 

 

... ... ... ...

 

 

 

 

 

 

 

 

 

 

 

rn1

rn2

...

 

1

 

* – это с добавочными верхней строкой и правым столбцом, состоящих из свободных

членов уравнений.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример: Вычислить КМК:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

0,57 = 0,675

 

 

− 0,875

 

− 0,69

 

 

 

0

 

=

*

=

1

 

 

 

 

 

0,57

 

− 0,69

= 0,5 R = * / = 0,86

 

0,57

1

 

 

 

0,57

 

 

 

1

 

− 0,825

 

 

 

 

 

 

 

 

 

 

 

 

Коэффициент корреляции рангов (объединенные ранги)

Анализ информации неподдающейся количественной оценке.

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

39

способностей и ранжируют. Корреляция между рангами значительно точнее отражает взаимосвязь.

Есть n учащихся и ранги по некоторому фактору А: X1…Xn и по фактору B: Y1…Yn. Xi, Yi – перестановки n первых натуральных чисел.

Xk-Yk=dk – мера тесноты связи A и B. Если все dk=0, то A и B полностью соответствуют.

∑X k =

n(n +1)

,

 

 

=

n +1

 

X

n

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

n +1

 

 

 

n +1

 

 

 

 

 

 

 

 

X

 

−

 

 

 

 

 

Y

−

 

 

 

 

 

 

 

∑(X − X )(Y −Y )

 

 

 

 

 

 

 

 

 

 

ρ =

∑

 

 

k

 

 

 

2

 

k

 

2

 

 

 

nσX σY

=

 

 

 

n

+1 2

 

 

 

 

n +1

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X

 

−

 

 

 

 

 

Y

 

−

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

2

 

 

 

 

 

 

 

∑

k

 

 

 

 

 

∑

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ρ =1−

6∑d 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n3 − n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

∑xy

=

12∑xy

=1

−

6∑d 2

∑x2

∑y2

n

3 − n

n3 − n

 

 

 

 

Последнее выражение – коэффициент корреляции рангов Спирмена.

Существуют и другие показатели тесноты связи:

ККР Кендела: удобен для углубленных исследований, когда невозможно установит ранговые различия. Строятся объединенные усредненные ранги и

 

n3 − n

−(Tx +Ty )− ∑d 2

 

ρ =

6

 

 

 

 

 

 

 

 

n3 − n

 

 

n3

− n

 

 

 

 

 

− 2T

 

 

− 2T

 

 

 

 

 

6

 

x

6

 

y

 

 

 

 

 

 

 

 

 

l

ti3 −ti

 

 

 

Tx = Ty =

∑

12

 

 

 

 

 

 

 

i=1

ti – число объединенных рангов.

Метод ранговой корреляции

Позволяет анализировать множество факторов и выделять доминирующие.

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

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

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

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

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

12∑d 2

W = m2 (n3 − n)

m – число специалистов n – число факторов.

Чем больше W, тем больше степень согласованности. Если W=0, то согласованность отсутствует. При W=1 – полная согласованность.

40