scilab的强大回归

Question

为了实现稳健的线性回归，我想实现一个具有Geman-McLure损失函数的M估计器

M-Estimators类在this document中呈现，Geman-McLure可在第13页找到。 要解决最小化问题，建议使用Iteratively reweighted least squares。如何在scilab中实施此过程？我可以使用optim吗？

Answer 1

您链接的the site of the document还有一些Matlab demo files available in a zip。这个zip中有两个文件我认为很重要：

• utvisToolbox /教程/ lineTut / robustDemo.m
• utvisToolbox /教程/ lineTut / robustIteration.m

在robustDemo.m文件中，有一个鲁棒M估计算法的实现。

回答您的问题如何在SciLab中实现这一点;您可以先使用mfile2sci将这些matlab文件转换为scilab。至少在 sampleRobustObj和robustIteration函数，只使用基本的Matlab内容，应该可以mfile2sci转换。