最適なθを計算するために、ロジスティック回帰では、最適なシータを得るためにfminuncに渡されるcostFunction(最小化される関数)を作成する必要があります。
また、costFunctionの勾配を計算できる場合は、 'GradObj'オプションを 'on'に設定します
options = optimset('GradObj','on');
Xの勾配値gを第2の出力引数として返すようにcostFunctionをコード化する。
それから私は与える
[theta, cost] = fminunc(@(t)(costFunction(t, X, y)), initial_theta, options);
ここで、Xはデータ行列、yは応答です。 Rで上記をどのように実装できますか?