fminsearchの可能性のある使用

Question

私は与えられたデータに対して以下のコードを書いています。おそらくネストされた関数を書くことによって間違いを犯します。誰も私にエラーを示すことができますか？私はどんな種類の答えにも満足しています。

clear; 
    t=[1 3 6 9 12 18]';%time 
    y=[.94 .77 .40 .26 .24 .16]';%probability 
    n=100;%trial 
    x=y.*n;%correct replies 
    p=rand(2,1);%starting parameters 
    f=@(t)p(1,1).*t.^(-p(2,1));%model function 
%%%%%%%%%%%%%%%%%%%%%%% 
    ff=@(p)sum(x.*log(-f)+(n-x(i).*log(-1+f)));%sum of the -loglikelihood 
    ffmin1=fminsearch(ff,p(1,1)) 
    ffmin2=fminsearch(ff,p(2,1)) 

Answer 1

コードには多くの誤りがあります。定義されていないi、ではないインデックスxをすることができます

f=@(t,p) p(1,1).*t.^(-p(2,1));%model function 

あなたの目的関数は、パラメータなしでfを呼び出すことはできません

あなたのモデル関数は、あなたのパラメーターを受け入れなければならない、と：行ずつ行きます。

ff= @(p) sum(x.*log(-f(t,p))+((n-x).*log(-1+f(t,p)))); 

fminsearchは、すべての変数と同時に、いない時に1つの最小化問題を解決します。

ffmin=fminsearch(ff,p) 

これはモデルで発生する可能性のあるエラーは考慮しませんが、少なくともコードは実行する必要があります。