I unable to run my program.Please help me...
2 次查看(过去 30 天)
显示 更早的评论
y=[repmat(0,1,7840),repmat(1,1,1317),repmat(2,1,239),repmat(3,1,42),repmat(4,1,14),repmat(5,1,4),repmat(6,1,4),repmat(1,1,1)]
n=length(y)
k=tabulate(y)
x=y(1)
syms f j
f =@(x,r,d,j) (factorial(x)/(factorial(x-j)*factorial(j)))*((-1)^j)*((d+r+j+1)/((d+r+j)^2))
s1=symsum(f,j,0,x)
custlogpdf =@(x,d,r) log(factorial(r+x-1)/(factorial(x)*factorial(r-1)))+ log(s1);
phat = mle(x,'logpdf',custlogpdf, 'start',[0.4,5,2,1])
1 个评论
回答(1 个)
José-Luis
2017-9-7
custlogpdf =@(x,d,r) log(factorial(r+x-1)/(factorial(x)*factorial(r-1))+ log(s1));
2 个评论
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Calculus 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!