Логарифмическая модель

Исходные значения X и Y

> X: =[1788, 1890, 2064, 2158, 2296, 2330, 2350, 2440, 2510, 2610, 2633, 2688, 2750, 3500, 3620, 3800];

> Y: =[6000, 6000, 6600, 6600, 6800, 7500, 7420, 7580, 7550, 7785, 7855, 7700, 7963, 8500, 8560, 9100];

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

> r: =Correlation(X, Y);

> Sr: =sqrt(1-r^2)/sqrt(Count(X)-2.0);

> t: =r/Sr;

>

> XX: =RandomVariable(StudentT(Count(X)-2.0)):

> 1-CDF(XX, abs(t));

> plot(PDF(XX, x), x=-6..6);

>

Линейная модель

>

f1: =Fit(a*x+b, X, Y, x);

> m1: =Fit(a*x+b, X, Y, x, output=solutionmodule):

> m1: -Results(leastsquaresfunction);

> k: =m1: -Results(parametervector);

> S: =m1: -Results(standarderrors);

> t1: =k[1]/S[1];

> t2: =k[2]/S[2];

> 1-CDF(XX, abs(t2));

>

Квадратичная модель

>

f2: =Fit(a+b*x+c*x^2, X, Y, x);

m: =Fit(a1*x+a2*x^2+b, X, Y, x, output=solutionmodule):

M: -Results(leastsquaresfunction);

k: =m: -Results(parametervector);

S: =m: -Results(standarderrors);

t1: =k[1]/S[1];

t2: =k[2]/S[2];

t3: =k[3]/S[3];

CDF(XX, abs(t1));

CDF(XX, abs(t2));

CDF(XX, abs(t3));

Кубическая модель

> f3: =Fit(a+b*x+c*x^2+d*x^3, X, Y, x);

Warning, model is not of full rank

Экспоненциальная модель

> f4: =ExponentialFit(X, Y, x);

Логарифмическая модель

> f5: =LogarithmicFit(X, Y, x);