VPF::[MatLab] - Форум программистов

Metalex

Дата 18.2.2011, 15:55 (ссылка)

(нет голосов)

Загрузка ...

Опытный

Профиль
Группа: Участник
Сообщений: 635
Регистрация: 22.10.2008
Где: Украина-ZPсity

Репутация: нет
Всего: нет

Код


function [p, Q]= chi2test(x)
    
    % Check the arguments.
    if(nargin ~= 1),            error('One and only one argument required!');                    end
    if(ndims(x) ~= 2),            error('The argument (x) must be a 2d matrix!');                    end
    if(any(size(x) == 1)),        error('The argument (x) must be a 2d matrix!');                    end
    if(any(~isreal(x))),        error('All values of the argument (x) must be real values!');    end
    
    % Calculate Q = sum( (a-np*)^2/(np*(1-p*)) )
    s=        size(x, 1);
    r=        size(x, 2);
    np=        sum(x, 2)/sum(sum(x)) * sum(x);        % p=sum(x, 2)/sum(sum(x)) and n=sum(x)
    Q=        sum(sum((x-np).^2./(np)));
    
    % Calculate cdf of chi-squared to Q. Degrees of freedom, v, is (r-1)*(s-1).
    p=        1 - gammainc(Q/2, (r-1)*(s-1)/2);
end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

clear all      % Ochistit pamyat
clc
sf='LR2.txt';  % Imya faila dannux
x=load(sf);    % Vvodim vuxodnue dannue
x=sort(x(:));  % Pereformatiruem stolbets i razsortiruem
H2=chi2test(x);
fprintf('Критерій Пірсона H2=%d\n',H2)

Ошибка:

Цитата

One and only one argument required!

В чем дело?

Файл LR2.txt состоит из 900 действительных чисел в столбик

--------------------

Don't let the system get you down.

0 Пользователей читают эту тему (0 Гостей и 0 Скрытых Пользователей)
0 Пользователей:
« Предыдущая тема \| Математические пакеты \| Следующая тема »