Электронный учебно-методический комплекс по учебной дисциплине «Прикладная математика» для студентов специальности 6-05-0716-01 «Метрология, стандартизация и контроль качества»

мой. Отклонение отдельной точки от линии регрессии (от предсказанного значения) называется остатком.

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

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

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

Чем меньше разброс значений остатков около линии регрессии по отношению к общему разбросу значений, тем лучше прогноз.

В случае множественной линейной корреляции важную роль играет ве-

личина R2 – коэффициент детерминации (определенности) как показатель качества модели или применимости данного набора предикторов для описания зависимой переменной Y .

Значение R2 является индикатором степени подгонки модели к данным. Коэффициент детерминации непосредственно интерпретируется следующим

образом. Если R2 0,4 , то только 40% от исходной изменчивости ряда Y могут быть объяснены предикторами Xi , а 60% остаются необъясненными.

Таким образом, величина R2 есть доля дисперсии исследуемой переменной Y , объяснённая переменными Xi .

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

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

R2 0,5. Значение R2 , близкое к 1, показывает, что модель объясняет почти всю изменчивость соответствующей переменной.

Интерпретация коэффициента множественной корреляции R. Обычно,

степень зависимости двух или более предикторов Xi с зависимой перемен-

ной Y выражается с помощью коэффициента множественной корреляции

R. Коэффициент множественной корреляции R можно интерпретировать как

61

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

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

Для интерпретации направления связи между переменными смотрят на знаки регрессионных коэффициентов (или B-коэффициентов). Если B- коэффициент положителен, то связь этой переменной с зависимой прямая; если B-коэффициент отрицателен, то связь обратная. Если B-коэффициент равен нулю, связь между переменными отсутствует.

II ПРАКТИЧЕСКИЙ РАЗДЕЛ 2.1. ЛАБОРАТОРНЫЕ РАБОТЫ

Лабораторная работа № 1. Числовые характеристики случайных величин. Вероятностные распределения. Знакомство с пакетом Statistica.

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

Контрольный пример 1.1.

Число α-частиц, достигающих счетчика, является случайной величиной Х, распределенной по закону:

Пользуясь пакетом MS Excel, найти числовые характеристики случайной величины Х: математическое ожидание М(Х), дисперсию D(X), среднее квадратическое отклонение σ(Х), коэффициент асимметрии As(X) и эксцесс Ek(X); указать их размерность; построить многоугольник распределения случайной величины.

Решение.

1. Создадим массивы данных в MS Excel.

Это – массив значений случайной величины (A2:A12) и массив вероятностей (B2:B12), с которыми случайная величина Х принимает свои значения

(рис. 1.1).

62

2. Найдем математическое ожидание

случайной величины X:

В MS Excel для этого можно использовать функцию СУММПРОИЗВ (массив 1, мас-

сив 2), которая позволяет найти сумму произведений диапазонов ячеек:

1) Выделим ячейку B14, в которую будет

Выделим ячейку C2, в которой созадим формулу: =A2-$B$14. После нажатия клавиши «Enter» в ячейке C2 появится число –3,868.

Скопируем формулу протягиванием ячейки C2 указателем мыши за правый нижний угол табличного курсора (при нажатой левой кнопки мыши) до ячейки C12 (см. рис 1.2).

Выделим ячейку B15, в которую будет возвращено значение дисперсии. В выделенную ячейку вводим функцию СУММПРОИЗВ (С2:С12, С2:С12, B2:B12).

63

Выделим ячейку D2, в которой созда-

дим формулу: =C2^3*B2.

После нажатия клавиши «Enter» в ячейке D2 появится число – 1,215.

Скопируем полученную формулу в ячейки D3:D12.

С помощью автосуммирования (кнопка «Σ» на вкладке «Главная») найдем сумму ячеек столбца D.

Результат (3,479) появляется в ячейке

D13 (см. рис. 1.3).

Вычислим коэффициент асимметрии:

As(X ) 1,961 3 3,479 0,462 0,5.

Рис. 1.2 – вычисление M(X) и D(X)

6. Аналогичным образом вычислим эксцесс случайной величины (см.

Рис. 1.3 – вычисление коэффициентов асимметрии и эксцесса

64

7. Построим многоугольник распределения случайной величины (рис. 1.4). Для этого выделим ячейки со значениями вероятностей В2:В12. На вкладке «Вставка» в группе «График» (для EXCEL 2016 – «Диаграммы»)

выберем тип графика – график с маркерами.

Щелкнем правой кнопкой мыши на поле графика и выберем «Выбрать данные». В появившемся окне «Выбор источника данных», справа, где «Под-

писи горизонтальной оси (категории)», выбираем «Изменить».

В открывшемся окне «Диапазон подписей оси» добавляем массив значений случайной величины А2:А12 и нажимаем ОK.

На вкладке «Макет» в разделе «Названия осей» выбираем подписи осей: х и p. На той же вкладке выбираем «Название диаграммы» и указываем название графика: «Многоугольник распределения случайной величины Х»

(рис. 1.4).

Нужный макет графика можно выбрать также на вкладке «Конструктор». Редактирование графика производится после щелчка правой кнопкой мыши на нужном элементе диаграммы.

Рис. 1.4 – Многоугольник распределения исходной случайной величины

Контрольный пример 1.2. В пакете MS Excel необходимо:

1. Ввести в таблицу значения аргумента х в диапазоне от –3 до 5 с шагом

0,2.

2.Вычислить значение плотности стандартного нормального распределения, а также плотности нормального распределения с параметрами a = 2, σ

=1; a = 0, σ = 0,5; a = 1, σ = 2.

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

65

Решение. Для вычисления значений плотности и функции нормального распределения в Excel используется встроенная статистическая функция НОРМ.РАСП.

Для построения графика плотности распределения протабулируем функцию f x на отрезке (–3; 5) с шагом 0,2. Для этого вводим в А3 подпись

«х», в ячейки В3:Е3 подписи, соответствующие условию задачи (рис.1.6).

Вводим в A4 значение –3, в А5 – значение –2,8, обводим, выделяя, ячейки A4 и A5 и, захватив нижний правый угол рамки вокруг ячеек A4 и A5, перетягиваем его вниз до ячейки А44, что позволит автоматически занести в столбец значения от –3 до 5 с шагом 0,2 (на рис.1.6 приведена часть значений отрезка).

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

ем категорию «Статистические» и функцию НОРМ.РАСП.

Появляется диалоговое окно функции. Заполняем ее как показано на рис. 1.5.

Рис. 1.5 – Диалоговое окно функции НОРМ.РАСП для N(0,1)

В поле «Интегральная» ставим ноль, нажимаем «ОК». В ячейке B4 появился результат.

Захватив нижний правый угол ячейки B4, автозаполнением растягиваем результат на ячейки B4:B44 (см. рис. 1.6).

Аналогично находим значения функции f x для N(2;1); N(0;0,05);

N(1;2).

Строим график по данным. Выделяем диапазон ячеек A4:E44. Вызываем мастер диаграмм, выбрав пункты меню Вставка/Диаграммы. Выбираем тип диаграммы «Точечная» и вид – «Точечная с гладкими кривыми». Получаем графики плотности нормального распределения (см. рис. 1.6).

66

Делаем вывод: параметр a изменяет положение графика, с увеличением параметра график смещается вправо; параметр σ влияет на ширину графика, с увеличением параметра график растягивается.

Графики функции нормального распределения получим аналогично – при вызове функции НОРМ.РАСП в поле «Интегральная» нужно поставить

1.

Рис. 1.6 – Часть исходных данных и решение контрольного примера 1.2

Контрольный пример 1.3. Исследование основных распределений в па-

кете Statistica.

1.3.1. Основы работы в пакете Statistica. Пакет Statistica состоит из ста-

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

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

Прикладное окно пакета Statistica имеет стандартную для окна приложе-

ния MS Windows структуру (рис. 1.7).

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

67

Набор данных в пакете Statistica – это прямоугольная таблица, столбцам которой соответствуют обрабатываемые переменные (Variables), а строкам отвечают наблюдения (Cases) значений переменных. Для создания нового набора данных нужно, прежде всего, завести файл с трафаретом таблицы нужных размеров, который открывается при запуске пакета (рис. 1.7).

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

Рис. 1.7 – Окно пакета Statistica

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

полняется процедура Probability Calculator:

Statistics – Probability Calculator – Distributions.

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

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

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

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

ласть Long name (label, link or formula with function), которая предназначена для ввода выражений и комментариев. Набор формулы следует производить, начиная со знака равно (=), далее встроенные функции могут быть набраны

68

непосредственно с клавиатуры либо с помощью конструктора (кнопка Function).

1.3.2. Нормальное распределение.

Запустим процедуру Probability Calculator, для чего в меню Statistics выберем модуль Probability Calculator, выделим соответствующую строку и нажмем кнопку ОК. В результате откроется окно калькулятора Probability Distribution Calculator.

Выполним расчеты для нормального распределения со средним значением a 0,5 и 1. Первая задача будет заключаться в поиске квантили для вероятности p 0.8 и построении графиков плотности и функции распределения.

В окне калькулятора в соответствующие поля введем параметры распределения, вероятность p и отметим опции Inverse и Create Graph.

После нажатия кнопки Compute получим результат (рис. 1.8). Значение квантили для заданной вероятности равно 1.341621.

Рис. 1.8 – Графики плотности и функции распределения для нормального распределения

Следующая задача заключается в определении значения функции распределения для заданного значения случайной величины x 1. Введя в поле Х значение, равное 1, и нажав кнопку Compute, в поле p получим значение

функции распределения, равное 0,691462 (рис. 1.9).

69