Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Скачиваний:
15
Добавлен:
27.11.2024
Размер:
902 Кб
Скачать
☆

Стьюдента при n1+n2–2 степенях свободы и вероятности 2α. Если Tнабл < Tкр, то принимается нулевая гипотеза. В противном случае нулевую гипотезу отвергают.

3) При конкурирующей гипотезе Н1: ̅ < ̅ критическую точку

1

2

Tкр(2α, n1 + n2 – 2) находим по таблице критических точек распределения Стьюдента при n1+n2–2 степенях свободы и вероятности 2α. Если Tнабл > –T кр, то принимается нулевая гипотеза. В противном случае нулевую гипотезу отвергают.

ТЕМА 3. РЕГРЕССИОННЫЙ И КОРРЕЛЯЦИОННЫЙ АНАЛИЗ

Лекция 6. Выборочные уравнения регрессии

Зависимость между

переменными величинами, когда каждому

значению одной переменной

может соответствовать множество значений

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

методами корреляционного и регрессионного анализа. Основной задачей

регрессионного анализа является установление формы и изучение

зависимости между переменными, корреляционного анализа – выявление связи между случайными переменными и оценка ее тесноты. В

регрессионном анализе рассматривается зависимость случайного результативного признака y от неслучайных факторных признаков x1, x2, ..., xn. В случае единственного факторного признака x различают следующие виды регрессий: линейную, гиперболическую, показательную, степенную,

логарифмическую, параболическую и т.д. Предположим, что для оценки параметров регрессии взята выборка, содержащая n пар значений (xi, yi), где i

= 1, 2, … , n. Оценкой предложенных выше уравнений регрессии являются выборочные уравнения регрессии:

 

линейное ̂ = 0 + 1 ;

 

 

 

 

гиперболическое ̂ =

+

1

;

 

 

 

0

 

 

 

 

 

 

 

показательное ̂ =

;

 

 

 

 

0

1

 

 

 

 

степенное

̂ =

1;

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

логарифмическое

 

̂ = 0 + 1 ;

 

 

 

 

параболическое

̂ =

0

+ +

2

2

,

 

 

 

 

 

 

 

1

 

 

 

где

параметры

a0,

 

a1,

a2

являются

точечными

оценками

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

Регрессионную модель удобно представлять графически. Для этого на координатной плоскости откладываются точки Pi (xi, yi), (i = 1, 2, … , n) (рис. 4.1). Полученный график называется диаграммой рассеивания.

Y

40

 

 

 

 

 

 

 

35

 

 

 

 

 

 

 

30

 

 

 

 

 

 

 

25

 

 

 

 

 

 

 

20

 

 

 

 

 

 

 

15

 

 

 

 

 

 

 

10

 

 

 

 

 

 

 

5

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

-5

0

5

10

15

20

25

 

 

 

 

 

 

 

X

 

Рис. 4.1. Диаграмма рассеивания

 

Построив диаграмму рассеяния, можно подобрать вид уравнения регрессии. На рис. 4.1 для одних и тех же экспериментальных точек построены линейная и показательная регрессии. Видим, что экспериментальные точки располагаются ближе к линии = 0 1 , чем к прямой. Следовательно, можно сделать вывод, что показательная регрессия более адекватно описывает фактические данные, чем линейная.

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

(значимости) уравнения регрессии на уровне значимости α вычисляют наблюдаемое значение случайной величины

 

 

 

 

 

2( −2)

 

 

 

 

 

 

 

 

 

=

 

 

,

 

 

 

 

(4.3)

 

 

 

2

 

 

 

 

 

набл

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ост

 

 

 

 

 

 

где

остаточная дисперсия

2

 

 

и дисперсия

уравнения регрессии

 

 

 

 

 

 

ост

 

 

 

 

 

2

находятся по формулам

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

=

∑( − ̂ )2

 

,

2

=

∑( ̅− ̂ )2

 

(4.4)

 

 

 

 

ост

 

 

 

−1

 

 

 

−1

 

 

 

 

 

 

 

 

 

 

Далее находим критическое значение критерия F(α, 1; n – 2) по таблице критических точек распределения Фишера (приложение 4) при k1=1, k2 = n – 2 степенях свободы и уровне значимости α. Если Fнабл > F(α; 1; n

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

Лекция 7. Модель линейной регрессии

Рассмотрим более подробно линейное уравнение регрессии.

Сущность метода наименьших квадратов заключается в нахождении

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

 

 

 

= ∑( − ̂ )2

→ .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для нахождения параметров

a0, a1 линейного уравнения регрессии

̂ = 0 + 1 решается система уравнений

 

 

 

 

 

 

 

 

0

+

∑

= ∑ ;

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

{ ∑

+

∑

 

2 =

 

∑ .

 

(4.1)

 

 

 

0

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

Для параболического

уравнения

регрессии

̂ =

+

+

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

1

2

 

система уравнений для нахождения параметров a0, a1, a2

имеет вид

 

 

 

 

+

∑

+

 

∑ 2 = ∑ ;

 

 

 

 

 

 

0

 

1

 

 

 

2

 

 

 

 

 

 

 

 

 

{

∑

+

∑ 2

+

2

∑

3

= ∑ ;

(4.2)

 

 

0

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

∑ 2

+

∑

3

+

2

∑

4

= ∑ 2

.

 

 

 

 

0

 

 

1

 

 

 

 

 

 

 

 

 

 

 

Все предложенные выше виды нелинейных регрессий (кроме параболической) могут быть сведены к линейной путем какой-либо замены

переменной. Для гиперболической регрессии вводится переменная x′= 1/x,

для логарифмической регрессии x′= lnx , уравнения показательной и степенной регрессии предварительно логарифмируют.

В качестве универсального показателя тесноты связи между величинами x и y используется выборочный линейный коэффициент корреляции

 

 

 

= 1

 

(4.5)

 

 

 

 

Здесь sx и sy – средние квадратические отклонения соответствующих признаков (факторного и результативного).

Линейный коэффициент корреляции изменяется в пределах –1 ≤ r ≤ 1.

Если r > 0, то связь между переменными x и у прямая, если r < 0, то связь между переменными x и у обратная. При r = 0 связь между переменными отсутствует. При |r| = 1 связь между x и у функциональная, т.е.

наблюдаемые значения располагаются точно на прямой.

Пусть вычисленное значение r≠ 0. Проверим гипотезу H0 об отсутствии линейной корреляционной связи между переменными, т.е. H0: ρ= 0 при альтернативной гипотезе H1: ρ ≠ 0. Для проверки этой гипотезы на уровне значимости α вычисляют наблюдаемое значение критерия

набл = √1|−| 2 √ − 2. (4.6)

Критическое значение критерия T(1–α, n–2) находят по таблице критических точек распределения Стьюдента (приложение 2) для числа степеней свободы n – 2 и уровня значимости α. Если Tнабл < T(α, n – 2), то гипотеза H0 принимается, в противном случае гипотеза H0 отвергается, т.е.

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

По уравнению линейной регрессии можно получить точечный и интервальный прогнозы.

Точечный прогноз заключается в получении прогнозного значения уp,

которое определяется путем подстановки в уравнение регрессии ̂ = 0 +

1 соответствующего (прогнозного) значения xp

̂р = 0 + 1 р.

Интервальный прогноз заключается в построении доверительного интервала прогноза, т. е. нижней и верхней границ уpmin, уpmax интервала,

содержащего точную величину для прогнозного значения ̂р (̂р < ̂р <

̂р ). Доверительный интервал всегда

определяется с

заданной

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

(степенью

уверенности),

соответствующей

принятому

значению уровня значимости α. Предварительно вычисляется стандартная

ошибка прогноза

 

 

 

 

 

 

 

 

 

 

̂р

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

(

− ̅)2

 

 

 

 

 

= √1 +

+

 

 

(4.7)

 

 

 

 

2

 

 

̂р

ост

 

 

 

 

 

 

 

 

 

 

 

∑( − ̅)

 

 

 

 

 

 

 

 

 

 

 

 

 

где = √

∑( − ̂)2

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ост

− −1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

̂р и верхняя ̂р границы интервала прогноза

̂р = ̂р − ̂р, ̂р = ̂р + ̂р,

где ̂р = ∙ ̂р, значение t находят по таблице критических точек распределения Стьюдента для числа степеней свободы n–2 и уровня значимости α.

ТЕМА 4. МОДЕЛИРОВАНИЕ ВРЕМЕННЫХ РЯДОВ И

ПРОГНОЗИРОВАНИЕ Лекция 8-9. Моделирование временных рядов

Временным рядом (рядом динамики, динамическим рядом)

называется упорядоченная во времени последовательность численных показателей{(yi,ti), i=1,2,...,n}, характеризующих уровни развития изучаемого явления в последовательные моменты или периоды времени.

Величины yi называются уровнями ряда, а ti – временными метками

(моменты или интервалы наблюдения). Обычно рассматриваются временные ряды с равными интервалами между наблюдениями, в качестве значений ti

берутся порядковые номера наблюдений и временной ряд представляется в

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

, … , , где n – количество наблюдений.

1

2

 

 

 

 

Целью исследования

временного

ряда

является

выявление

закономерностей в изменении уровней ряда и построении его модели в целях прогнозирования и исследования взаимосвязей между явлениями.

Моделирование тенденции временного ряда начинается с проверки

наличия тенденции. Для этого наиболее широко применяются

метод

сравнения средних и метод Фостера-Стюарта.

 

Для получения ряда с меньшим разбросом уровней, что в ряде случаев

позволяет на основе визуального анализа сделать вывод о

наличии

тенденции, применяется сглаживание временного ряда.

Сглаживание временного ряда по методу скользящей средней

заключается

в замене исходных уровней ряда yt сглаженными значениями

y′t, которые получаются

 

как

 

среднее значение

определенного

числа

уровней исходного ряда,

 

симметрично

окружающих

значение yt. В

результате

получается

временной

ряд

y′t,

меньше

подверженный

колебаниям.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для вычисления сглаженных значений y′t

по

методу

простой

скользящей средней используются следующие формулы:

 

 

1) Нечетный интервал сглаживания g = 2p+1 (интервал сглаживания –

количество исходных уровней ряда (yt), используемых для сглаживания):

 

 

∑ +

 

 

 

 

−

+

− +1

+ +

+

 

 

 

′ =

= −

 

=

 

 

 

+ −1

+

(5.5)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 +1

 

 

 

 

 

2 +1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где уt – фактическое значение уровня исходного ряда в момент t; y′t –

значение скользящей средней в момент t; 2р+1- длина интервала сглаживания. Формула (5.5) при интервалах сглаживания g = 3 и g = 5

принимает вид

′ =

 

+ +

; ′ =

 

+

+ +

+

−1

 

+1

−2

−1

 

+1

+2

 

 

3

 

 

 

 

5

 

 

 

 

 

 

 

 

 

 

2) Четный интервал сглаживания g = 2p:

 

1

 

+∑ + −1

 

+

1

 

 

1

 

+

+ +

+

1

 

 

 

 

 

 

 

 

 

′ = 2

−

= − +1

 

2

+

= 2

−

− +1

+ −1

2

+

(5.6)

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Формула (5.6) при интервалах сглаживания g = 2 и g = 4 принимает вид

1

 

+ +

1

 

1

 

+

+ +

+

1

 

′ = 2

−1

 

2

+1

; ′ = 2

−2

−1

 

+1

 

2

+2

 

 

2

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

При использовании скользящей средней с длиной активного участка g = 2p+1 первые и последние р уровней ряда сгладить нельзя, их значения теряются. Для восстановления потерянных значений временного ряда можно использовать следующий прием:

а) Вычисляется средний прирост ∆у на последнем активном участке

(

, … , )

 

 

 

−

 

 

 

 

 

 

−

 

 

∆ =

 

−

,

 

− 1

 

 

 

 

где g – длина активного участка.

 

 

 

б) Определяются значения последних р=(g–1)/2 уровней сглаженного

временного ряда с помощью последовательного прибавления среднего абсолютного прироста ∆у к последнему сглаженному значению y′n–p

′ − +1 = ′ − + ∆ , ′ − +2 = ′ − +1 + ∆ , … , ′ = ′ −1+∆ .

Аналогичная процедура применяется для восстановления первых р уровней временного ряда.

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

нахождение аналитической функции ŷ = f(t), характеризующей основную

тенденцию изменения уровней ряда с течением времени. Сама функция f(t)

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

Чаще всего в качестве кривой роста применяются следующие

функции:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

линейная

=

 

+ ;

 

 

 

 

(5.7 )

 

 

 

 

 

0

 

 

1

 

 

 

 

 

парабола второго и более высоких порядков

 

 

 

=

0

+

+

2

2

+ +

;

(5.9)

 

 

 

1

 

 

 

 

 

 

 

 

 

гиперболическая

 

 

 

=

0

+

1

;

(5.10)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

экспонента

= 0+ 1 ;

 

(5.11)

 

 

 

 

 

 

 

 

 

показательная

 

=

×

;

(5.12)

 

 

 

 

0

1

 

 

 

степенная

=

0

× 1 .

 

(5.13)

 

 

 

 

 

 

 

Построение таких функций аналогично построению уравнений парной

регрессии (линейной или нелинейной) с учетом того, что в качестве зависимой переменной используются фактические уровни временного ряда yt

, а в качестве независимой переменной моменты времени t = 1,2, ..., n.