编辑:为了帮助阐明我的问题,我希望与以下数据相符:
我可以使用cftool函数进行拟合,但是使用最小二乘法对我的二进制数据没有意义。只是为了说明...
所以,我的目标是使用fmincon函数拟合这些数据。
原始帖子:
我有一个运动控制实验的数据,其中参与者在执行任务时被计时,并根据他们的表现给出分数(失败或成功)。如您所料,我们假设参与者将有更多时间执行任务,因此犯错的次数会减少。
我正在尝试使用fmincon将函数拟合到此数据,但出现错误“使用fmincon出错(第609行) 提供的目标函数必须返回标量值。” 我不明白a)这意味着什么,或b)我该如何解决。
我在下面提供了一些示例数据和代码。任何帮助表示赞赏。
%Example Data:
time = [12.16 11.81 12.32 11.87 12.37 12.51 12.63 12.09 11.25
7.73 8.18 9.49 10.29 8.88 9.46 10.12 9.76 9.99 10.08
7.48 7.88 7.81 6.7 7.68 8.05 8.23 7.84 8.52 7.7
6.26 6.12 6.19 6.49 6.25 6.51 6 6.79 5.89 5.93 3.97 4.91 4.78 4.43
3.82 4.72 4.72 4.31 4.81 4.32 3.62 3.71 4.29 3.46 3.9 3.73 4.15
3.92 3.8 3.4 3.7 2.91 2.84 2.7 2.83 2.46 3.19 3.44 2.67 3.49 2.71
3.17 2.97 2.76 2.71 2.88 2.52 2.86 2.83 2.64 2.02 2.37 2.38
2.53 3.03 2.61 2.59 2.59 2.44 2.73 ]
error = [0 0 0 0 0 0 0 0 0 1 0 0 0 0 1 1 1 0 0 0 0 1 1 1 1 1 1 0 0 0 0 1 1 1 0 1 0 1 0 1 1 0 1 1 0 1 1 1 1 1 1 1 1 1 1 1 1 0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1];
%Code:
% initial parameters - a corresponds to params(1), b corresponds to params(2)
a = 3.0;
b = -0.01;
LL = @(params) 1/1+params(1)*(log(time).^params(2));
LL([a b]);
pOpt = fmincon(LL,[a b],[],[]);
答案 0 :(得分:1)
错误来自函数LL,该函数返回等于time
长度的多个值。
要正确使用fmincon
,您需要有一个仅返回一个值的函数。
答案 1 :(得分:0)
我相信逻辑回归将很好地适合您的数据和目的。在那种情况下,为什么不简单地使用Matlab的内置函数进行多项逻辑回归呢?
B = mnrfit(time,error)
关于函数LL
,确定要正确输入函数并且没有遗漏括号吗?
LL = @(params) 1/(1+params(1)*(log(time).^params(2)));
不带括号,您的功能等同于1 + a * log(x)^ b