I unable to run my program.Please help me...
1 回表示 (過去 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!