Определение значимости коэффициента корреляции

Уровень значимости коэффициента корреляции может быть определён по критерию Стьюдента:

Если , где - уровень значимости =0.95, а - число степеней свободы =n-2, то можно утверждать, что между х и у существует линейная зависимость, в противном случае – линейная зависимость отсутствует.

Значение T табл выбирается по таблице значений критерия Стьюдента.

Отчет выполненной данной работы содержит:

1. Содержательную постановку задачи (выбор значений x и у, полученных в результате эксперимента)

2. Математическую постановку задачи (функции, по которым проводились расчеты)

3. Блок-схему алгоритма решения задач

4. Программу на алгоритмическом языке

5. Результат работы программы

6. Вывод по расчётам

7. Расчет, полученный в Excel

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

Пример выполнения работы

Провести расчет на Qbasic и в Excel.

Содержательная постановка задачи

Проведён эксперимент по определению зависимости коэффициента активности от ионной силы раствора.

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

Ионная сила раствора Х(i) Коэффициент активности Y(i)
0,001 0,002 0,005 0,01 0,04 0,05 0,1 0,2 0,3 0,4 0,5 0,98 0,97 0,95 0,92 0,9 0,84 0,81 0,8 0,81 0,82 0,84

Табличные значения коэффициента Стьюдента = 2.201

Таблица значений критерия Стьюдента

f q
0.80 0.90 0.95 0.98 0.99 0.995 0.998 0.999
3.0770 6.3130 12.7060 31.820 63.656 127.656 318.306 636.619
1.8850 2.9200 4.3020 6.964 9.924 14.089 22.327 31.599
1.6377 2.35340 3.182 4.540 5.840 7.458 10.214 12.924
1.5332 2.13180 2.776 3.746 4.604 5.597 7.173 8.610
1.4759 2.01500 2.570 3.649 4.0321 4.773 5.893 6.863
1.4390 1.943 2.4460 3.1420 3.7070 4.316 5.2070 5.958
1.4149 1.8946 2.3646 2.998 3.4995 4.2293 4.785 5.4079
1.3968 1.8596 2.3060 2.8965 3.3554 3.832 4.5008 5.0413
1.3830 1.8331 2.2622 2.8214 3.2498 3.6897 4.2968 4.780
1.3720 1.8125 2.2281 2.7638 3.1693 3.5814 4.1437 4.5869
1.363 1.795 2.201 2.718 3.105 3.496 4.024 4.437
1.3562 1.7823 2.1788 2.6810 3.0845 3.4284 3.929 4.178
1.3502 1.7709 2.1604 2.6503 3.1123 3.3725 3.852 4.220
1.3450 1.7613 2.1448 2.6245 2.976 3.3257 3.787 4.140
1.3406 1.7530 2.1314 2.6025 2.9467 3.2860 3.732 4.072
1.3360 1.7450 2.1190 2.5830 2.9200 3.2520 3.6860 4.0150
1.3334 1.7396 2.1098 2.5668 2.8982 3.2224 3.6458 3.965
1.3304 1.7341 2.1009 2.5514 2.8784 3.1966 3.6105 3.9216

Значения критерия Стьюдента (t-критерия) для различного уровня значимости q и числа степеней свободы f (f = n-2 или n, n-число опытов).

1. Математическое ожидание: и

2. Дисперсия: и

3. Среднее квадратическое отклонение: и

4 Коэффициент вариации: и

5. Нормированное отклонение: и

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

7. Критерий Стьюдента:

БЛОК-СХЕМА

ПРОГРАММА НА ЯЗЫКЕ QBASIC

CLS

n = 11

INPUT "Введите Коэффициент Стьюдента="; tt

DIM x(1 TO n), y(1 TO n), tx(1 TO n), ty(1 TO n)

DATA 0.001, 0.002, 0.005, 0.01, 0.04, 0.05, 0.1, 0.2, 0.3, 0.4, 0.5

FOR i = 1 TO n

READ x(i)

NEXT i

DATA 0.98, 0.97, 0.95, 0.92, 0.9, 0.84, 0.81, 0.8, 0.8, 0.81, 0.82, 0.84

FOR i = 1 TO n

READ y(i)

NEXT i

mx = 0: my = 0

FOR i = 1 TO n

mx = mx + x(i)

my = my + y(i)

NEXT i

mx = mx / n: my = my / n

dx = 0: dy = 0: S = 0

FOR i = 1 TO n

dx = dx + (x(i) - mx) ^ 2

dy = dy + (y(i) - my) ^ 2

S = S + (x(i) - mx) * (y(i) - my)

NEXT i

dx = dx / (n - 1): dy = dy / (n - 1)

CKOx = SQR(dx): CKOy = SQR(dy)

Vx = CKOx * 100 / mx: Vy = CKOy * 100 / my

FOR i = 1 TO n

tx(i) = (x(i) - mx) / CKOx

ty(i) = (y(i) - my) / CKOy

NEXT i

Kxy = S / ((n - 1) * CKOx * CKOy)

tr = ABS(Kxy) * SQR(n - 2) / SQR(1 - Kxy ^ 2)

PRINT " N", "x", "y", "tx", "ty"

FOR i = 1 TO n

PRINT i, x(i), y(i), tx(i), ty(i)

NEXT i

PRINT "число ответов ="; n

PRINT "табличное значение критерия Стьюдента ="; tt

PRINT "математическое ожидание :"; "по x="; mx, "по y="; my

PRINT "дисперсия :"; "по x="; dx, "по y="; dy

PRINT "среднеквадратическое отклонение :"; "по x="; CKOx, "по y="; CKOy

PRINT "коэффициент вариации :"; "по x="; Vx, "по y="; Vy

PRINT "коэффициент корреляции ="; Kxy

PRINT "рассчитанное значение критерия Стьюдента ="; tr

IF tr >= tt THEN

PRINT "линейная зависимость между x и y Cуществует"

ELSE

PRINT "линейная зависимость между x и y Отсутствует"

END IF

END