Дата поступления в стационар: 20.04.2000 г. 4 страница



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

Решение. Создайте рабочую таблицу в Excel. Введите в ячейку А1 слово “Возраст”. Затем в ячейки А1:А12 – соответствующие значения возраста. В ячейку В1 введите “Систолическое давление”. В ячейки               В1:В12 – значения систолического давления. Затем вычисляется значение коэффициента корреляции между выборками. Установите курсор в свободную ячейку (А13). Нажмите кнопку Вставка функции (fx) на панели инструментов. Выберите в диалоговом окне Мастер функций статистические функции, затем выберите функцию КОРРЕЛ.


Появится диалоговое окно КОРРЕЛ. Указателем мыши введите диапазон данных “Возраст” в поле массив 1 (А2:А12). Диапазон данных “Систолическое давление” введите в поле массив 2 (В2:В12).


В ячейке А13 после нажатия кнопки ОК появится значение коэффициента корреляции – 0,61. Если проверить значимость коэффициента корреляции между переменными X и  при уровне значимости a=0,05 (при n=20 tкр.=2,1), то можно сделать вывод, что имеется заметная линейная корреляционная связь между  и X.

 

                        Регрессионный анализ.

 

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

Приведем уравнения регрессии Y на X и X на Y:

M(Y)x=f(x),  M(X)y=j(y),

M(Y)x – условное математическое ожидание величины Y, соответствующее значению x;

M(X)y – условное математическое ожидание величины X, соответствующее значению y.

В результате n независимых опытов получены n пар чисел (x1, y1), (x2, y2), … , (xn, yn).

Найдем по данным наблюдений выборочное уравнение прямой линии регрессии.

Выборочное уравнение линейной регрессии Y на X будем искать в виде

                                                                                       (*)

Угловой коэффициент прямой линии регрессии Y на X называют выборочным коэффициентом регрессии Y на X.

Подберем параметры ryx и b так, чтобы сумма квадратов отклонений ординат всех эмпирических точек от ординат соответствующих точек прямой (*) должна быть минимальной (в этом состоит сущность метода наименьших квадратов).

В результате применения метода наименьших квадратов получим следующие формулы для вычисления ryx и b:

 ;

,

где , ,

, ,

,

.

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

Для получения коэффициентов уравнения регрессии используется процедура Регрессия из пакета анализа. Кроме того, могут быть использованы функция ЛИНЕЙН для получения параметров регрессионного уравнения и функция ТЕНДЕНЦИЯ для получения предсказанных значений Y в требуемых точках.

Для реализации процедуры Регрессия необходимо:

- выполнить команду Сервис, Анализ данных;

- в появившемся диалоговом окне Анализ данных в списке Инструменты анализа выбрать строку Регрессия;

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

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

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

- если необходимо проверить отличие экспериментальных точек от предсказанных по регрессионной модели, следует установить флажок в поле График подбора;

- нажать кнопку ОК.

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

Приводимое значение R – квадрат (коэффициент детерминации) в регрессионной статистике определяет, с какой степенью точности полученное регрессионное уравнение аппроксимирует исходные данные. Если R – квадрат > 0,95, говорят о высокой точности аппроксимации. Если R – квадрат лежит в диапазоне от 0,8 до 0,95, говорят об удовлетворительной аппроксимации. Если R – квадрат < 0,6, то точность аппроксимации недостаточна и модель требует улучшения.

В таблице Дисперсионный анализ оценивается общее качество полученной модели: ее достоверность по уровню значимости критерия Фишера – p, который должен быть меньше, чем 0,05. Значение p определяем в строке Регрессия, в столбце Значимость F.

Значения коэффициентов модели определяются из таблицы в столбце коэффициенты – в строке Y – пересечение приводится свободный член; в строках соответствующих переменных приводятся значения коэффициентов при этих переменных. В столбце p – значение приводится достоверность отличия соответствующих коэффициентов от нуля. В случаях, когда p > 0,5, коэффициент может считаться нулевым. Это означает, что соответствующая независимая переменная практически не влияет на зависимую переменную и коэффициент может быть убран из уравнения.

Пример. Изучали зависимость между объемом Y (мкм3) и диаметром X (мкм) сухого эритроцита у млекопитающих. Результаты наблюдений приведены в таблице:

 

X Y
7,6 87
8,9 81
5,5 50
9,2 112
3,5 18
4,8 37
7,3 71
7,4 69
6,8 54

Необходимо на основании этих данных построить регрессионное уравнение.

Решение.

1. В пункте меню Сервис выберите строку Анализ данных


 и далее укажите курсором мыши на строку Регрессия.

2. В появившемся диалоговом окне задайте Входной интервал Y.

3. Укажите Входной интервал X.

4. Установите флажок в поле График подбора.

5. Укажите выходной диапазон. Для этого поставьте переключатель в положение Выходной интервал, затем наведите указатель мыши на правое поле ввода Выходной интервал  и, щелкнув левой кнопкой мыши, указатель мыши наведите на левую верхнюю ячейку выходного диапазона (C1). Щелкните левой кнопкой мыши. Нажмите кнопку ОК.


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

 

ВЫВОД ИТОГОВ

             
                 

Регрессионная статистика

             
Множественный R 0,94943              
R-квадрат 0,901418              
Нормированный R-квадрат 0,887335              
Стандартная ошибка 9,446213              
Наблюдения 9              
                 

Дисперсионный анализ

           
df SS MS F

Значимость F

   
Регрессия 1 5711,383 5711,383 64,00676 9,11E-05      
Остаток 7 624,6166 89,23094          
Итого 8 6336            
                 
Коэффициенты Стандартная ошибка t-статистика P-Значение Нижние 95% Верхние 95% Нижние 95,0% Верхние 95,0%
Y-пересечение -32,4752 12,50339 -2,59731 0,035564 -62,041 -2,90936 -62,041 -2,90936
Переменная X 1 14,28322 1,785308 8,000422 9,11E-05 10,06164 18,5048 10,06164 18,5048
                 
                 
                 

ВЫВОД ОСТАТКА

             
                 
Наблюдение Предсказанное Y Остатки            
1 76,07731 10,92269            
2 94,6455 -13,6455            
3 46,08255 3,917447            
4 98,93047 13,06953            
5 17,51611 0,483886            
6 36,0843 0,915701            
7 71,79235 -0,79235            
8 73,22067 -4,22067            
9 64,65074 -10,6507            

 


В таблице Дисперсионный анализ оценивается общее качество полученной модели: ее достоверность по уровню значимости критерия Фишера (строка Регрессия, столбец Значимость F, в примере 0,0000911, то есть p=0,0000911 и модель значима).

Приводимое значение R – квадрат (коэффициент детерминации) в регрессионной статистике определяет степень точности описания моделью процесса. В примере R – квадрат=0,9015. Так как R – квадрат < 0,95, не можем говорить о высокой точности аппроксимации.

Определим значения коэффициентов модели. На пересечении строки Y – пересечение и столбца Коэффициент приводится свободный член. В строке Переменная X1 приводится коэффициент при X1.

Поэтому выражение для определения объема сухого эритроцита у млекопитающих от диаметра будет иметь вид:

 

 

 

Однофакторный дисперсионный анализ.

 

 Для сравнения нескольких средних пользуются дисперсионным анализом. На практике дисперсионный анализ применяют, чтобы установить, оказывает ли существенное влияние некоторый качественный фактор А, который имеет m уровней А1, А2 … А m  на изучаемую величину Х. Например, если требуется выяснить, какая доза рентгеновского излучения наиболее эффективно влияет на темп размножения бактерий, то фактор А – рентгеновское излучение, а его уровни – дозы излучений.

 Основная идея дисперсионного анализа состоит в сравнении факторной дисперсии и остаточной дисперсии. В математической статистике доказывается, что факторная дисперсия характеризует влияние фактора А на величину Х, а остаточная – влияние случайных причин.

Рассмотрим случай, когда число испытаний на различных уровнях различно. Пусть произведено q1 испытаний на уровне А1, q2 испытаний на уровне А2, …, q m   испытаний – на уровне Аm.

Общую сумму квадратов отклонений наблюдаемых зачений от общей средней х находят по формуле :

 

Sобщ = [P1 + P2 + …+ Pm] – (R1 +R2 +… + Rm)2/n,

 где

P1 = - сумма квадратов наблюдавшихся значений признака на уровне А1;

 

P2 = - сумма квадратов наблюдавшихся значений признака на уровне А2;

                   . . .

 

Pm = - сумма квадратов наблюдавшихся значений признака на уровне А m.

 

R1 = , R2 = … Rm =  - суммы наблюдавшихся значений признака соответственно на уровнях А1, А2, … Аm .

 

n= q1 + q2 +… +qm - общее число испытаний ( объем выборки).

 

Факторную сумму квадратов отклонений групповых средних от общей средней, которая характеризует рассеяние " между группами" находят по формуле:

 Sфакт = [ ( R12/q1) + (R22/q2 ) +… + (Rm2/qm )] – [ (R1 + R2 + …+ Rm)2 /n]

 

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

S ост = S общ -S факт

Факторную дисперсию находят по формуле :

 

S2 факт = S факт / (m-1)

 Остаточную дисперсию находят по формуле:

 S2 ост = S ост /(n-m)

 

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

 Если факторная дисперсия окажется меньше остаточной, то фактор не оказывает существенное влияние на величину Х.

Если факторная дисперсия больше остаточной, то применяем критерий Фишера - Снедекора, для чего найдем наблюдаемое значение критерия


Дата добавления: 2019-02-12; просмотров: 136; Мы поможем в написании вашей работы!

Поделиться с друзьями:






Мы поможем в написании ваших работ!