R程序求救!

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

请问用nlm()函数求最优解时出现如下语句是什么错误啊？该往哪方面改呢？麻烦解答！谢谢了！

“错误于solve.default(cov, ...) : Lapack例行程序dgesv: 系统正好是奇异的”

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：R程序 Default Lapack dgesv fault 程序 求救

迭代过程中矩阵出现不可逆的情形（行列式为零）。

如果目标函数简单的话，建议自己编写优化程序，呵呵。

谢谢楼上的回答!

可是我的目标函数不简单，自己编程恐怕很难，难道没有现成的R函数或方法可用了吗？拜托再帮帮忙吧！谢谢了！

可以在R的贡献包里找其他的非线性优化包试试，

本人的体会是，都不大好用。

或者多选一些初始点试试，good luck！

做非线性问题跟迭代方法和初始值的选取有很大的关系的，需要不断尝试才行！

那可不可以推荐一些比较好的介绍非线性问题的书呢?

另外,请问如果最优解与迭代方法和初始值的选取密切相关，那岂不是可能会有很多解？那怎么求得完呢？怎么确定哪个最优呢？

谢谢！

现在非线性的书很多,对于非线性,迭代方法一般很难找到所有的最优解.一般尝试多种方法,多个初始值,从中选取相对最优.

请问在R语言里如果不能用牛顿法求解非线性最优化问题，还有其他什么方法吗?谢谢!

R基本安装里就有的优化函数：
optim()
constrOptim()

建议你 自己写出一阶条件和二阶条件，可以Taylor展开到一阶，稍微简化一下，看看效果如何。如果结果可以，再展到2阶。如果是似然函数，一般连续可微都是满足的。
具体可以用R的D()函数进行符号求解，简单用法见下（或者matlab的符号计算，直接把结果贴到R script里)。

初始值可以用格点遍历的方法，跑遍可能的初值组合，记得考虑正负域的不同！

> f=expression(x^2+x^5 +x*log(x))
> D(f,'x')
2 * x + 5 * x^4 + (log(x) + x * (1/x))
>

可是要是目标函数很复杂，很难求一阶或二阶条件，用optim()这类函数又解不出来（老是出现“系统正好是奇异”的错误），到底要怎么最优化呢？