Расширение "Статистическая связь двух переменных"

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

рассчитать коэффициент корреляции двух переменных;

построить уравнение линейной регрессии;

табулировать значения уравнения линейной регрессии для фактических значений аргумента;

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

получить графическое представление точек фактических данных и уравнения линейной регрессии;

строить тестовые массивы данных, равномерно разбросанных вокруг аналитически заданного тренда;

строить уравнения степенной и экспоненциальной регрессии.

Расширение требует установки файла s2.php (не входит в базовую комплектацию).

Функция s2()

Функция s2(x,y[,p]) возвращает массив, содержащий элементы:

kk - коэффициент корреляции переменных;

a - свободный член уравнения линейной регрессии y на x;

b - коэффициент при независимой переменной уравнения линейной регрессии y на x;

r2 - значение коэффициента R-квадрат;

ulr - текстовое представление уравнения регрессии.

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

Если параметр p=0 (значение по умолчанию), то элемент gxy в массив не включается. Если p>1, то функция не только возвращает массив указанных выше значений, но и выводит график из данного файла на страницу. Вывод графика возможен только в том случае, если в системе установлен файл grafik.php (не входит в базовую комплектацию).

Параметры x и y могут быть строками, содержащими разделённые пробелами фактические значения в формате [ключ:]значение или уже готовыми одномерными массивами.

Разметка

~x=massiv("1 2 3 4 5 6 7")
~y=massiv("2.4 3.7 6.2 7.7 10.5 12.2 13.7")
~u=s2(x,y)
@vms(u)

Вывод

Array
(
    [kk] => 0.99673419461949
    [a] => 0.17142857142857
    [b] => 1.9714285714286
    [r2] => 0.99347905472377
    [ulr] => Y=0.17142857142857+1.9714285714286*X
)

То же самое, но с выводом графика:

Разметка

~u=s2(x,y,2)
@vms(u)

Вывод

Array
(
    [kk] => 0.99673419461949
    [a] => 0.17142857142857
    [b] => 1.9714285714286
    [r2] => 0.99347905472377
    [ulr] => Y=0.17142857142857+1.9714285714286*X
    [gxy] => vf/678f7c8801037.jpeg
)

Полученные данные можно вывести и более изящно, используя возможности системы. Например:

Разметка

~u=s2(x,y,1)
@nt(1)
@nst(1)
@td(u('ulr'),2)
@kst()
@st(1,'Свободный член',u('a'))
@st(1,'Наклон прямой',u('b'))
@st(1,'Коэффициент корреляции',u('kk'))
@st(1,'R-квадрат',u('r2'))
@kt()
@vis(u('gxy'))

Вывод

Y=0.17142857142857+1.9714285714286*X
Свободный член	0,171
Наклон прямой	1,97
Коэффициент корреляции	0,997
R-квадрат	0,993

Функция s2td - формирование тестовых данных

Детальное описание в процессе разработки

s2td(formula,xmin,xmax,po,kz)

Функция s2ur - одновременное построение линейной, степенной и экспоненциальной регрессии

s2ur(vx,vy,p)

Линейный тест

	sko	formula
lr	12,7	1.94353556338+0.35215555617718*x
er	13,0	2.2170345196177exp(0.094382945370789x)
sr	14,6	2.0355081153996*pow(x,0.39306332647977)

sko: Сумма квадратов отклонений расчётных данных от фактических

FZ	Фактические значения
LR	Линейная регрессия
ER	Экспоненциальная регрессия
SR	Степенная регрессия

Экспоненциальный тест

	sko	formula
lr	21,8	1.5748936065932+0.36491721556876*x
er	21,2	1.7980293039081exp(0.11148664730842x)
sr	33,1	2.1569940486014*pow(x,0.30509874633484)

sko: Сумма квадратов отклонений расчётных данных от фактических

FZ	Фактические значения
LR	Линейная регрессия
ER	Экспоненциальная регрессия
SR	Степенная регрессия

Степенной тест

	sko	formula
lr	12,9	1.894300242422+0.15451769573119*x
er	14,2	1.8168904850395exp(0.06775589233912x)
sr	10,4	1.9497789090841*pow(x,0.22176757740004)

sko: Сумма квадратов отклонений расчётных данных от фактических

FZ	Фактические значения
LR	Линейная регрессия
ER	Экспоненциальная регрессия
SR	Степенная регрессия

Оцените публикацию по предложенной шкале

-5 -4 -3 -2 -1 0 +1 +2 +3 +4 +5

Статистика Код