помогите исправить ошибку в скрипте Matlab

ironking77

Всем шалом!
считаю методом максимального правдоподобия оценку для параметра gamma Парето распределния:x^gamma/(1+x^gamma)
вот сам скрипт
 
 
%в эту переменную импортируется выборка
x=sigma1;
>> % Плотность распределения Парето
p=inline('(x^gamma)/(1+x^gamma',... 'x','gamma');
% Функция правдоподобия
L=inline('-prod(feval(p,x,theta',...
'theta','p','x');
% Наиболее правдоподобное значение параметра
m=fminsearch(L,[1],[],p,x)

ошибка следующая
 

Error using ==> inlineeval at 15
Error in inline expression ==> -prod(feval(p,x,theta
Error in inline expression ==> exp(-mu)*mu.^x/factorial(x(1
N must be a matrix of non-negative integers.

Error in ==> inline.subsref at 27
INLINE_OUT_ = inlineeval(INLINE_INPUTS_, INLINE_OBJ_.inputExpr,
INLINE_OBJ_.expr);

Error in ==> fminsearch at 205
fv(:,1) = funfcn(x,varargin{:});

ironking77

ап

Lenochka00884

Перепиши все без всяких inline и будет тебе счастье.
Моё такое мнение.

Vlad128

ну так очевидно же, что ему хочется, чтобы функции вводились извне

ironking77

мне не нужно вводить функции извне...
я взял этот скрипт вот отсюда:http://www.exponenta.ru/educat/systemat/semeriy/lab3/ex3.asp
я не знаю как его переделать под распределение Парето.Если кто сможет это сделать буду очень признателен!
Спасибо за прочтение данных строк.
ПЫСЫ с меня денежное вознаграждение

sonik_23rus

да ты бы полный текст запостил, а так у тебя ошибка

N must be a matrix of non-negative integers.

ironking77

вот весь код ,который нужно переделать под распределение Парето
 
% Плотность распределения Пуассона
p=inline('exp(-mu)*mu.^x/factorial(x(1',...
'x','mu');
% Функция правдоподобия
L=inline('-prod(feval(p,x,theta',...
'theta','p','x');
% Наиболее правдоподобное значение параметра
m=fminsearch(L,[1],[],p,x)
m = 5.6500
% График функции правдоподобия
mu2=m-2:0.1:m+2;

for i=1:length(mu2 L2(i)=L(mu2(ip,x); end
L2=L2/L(m,p,x); plot(mu2,L2)

вот ссылка на выборку данных: ссылка на выборку
данные из екселевского файла я импортирую стандартно т.е. File->Import Data
Оставить комментарий
Имя или ник:
Комментарий: