[问题] R的规划求解 optim()、constrOptim()

楼主: MrSomebody (MrSomebody)   2016-06-29 21:02:31
[问题类型]:
程式咨询
[软件熟悉度]:
入门
[问题叙述]:
最近在研究SVM,但同样资料用R和matlab跑出来结果不一样,所以想尝试写算法,
在过程中有地方卡住,故上PTT求助各位大大。
http://imgur.com/5c2wQue
解线性kernal的SVM算法,即解出函数max{L_D function}(如图中红框),
x_i为解释变量(已知)、y_i为反应变量(已知),
要求alpha_i带入什么值时使L_D function最大。
可利用optim() 函数求主要L_D最大值,但无法加入条件式(图中蓝框两条件式);
constrOptim() 函数可以加不等式条件,但只限定alpha_i线性组合之不等式,无法设
sum(alpha_i&y_i)=0条件式。
不知R是否有其他函数可以求类似的问题?
谢谢大家耐心看完。
作者: qhair (卡早睡卡有眠)   2016-06-30 00:05:00
sum(\alpha_i*y_i)=0是线性的,可以写成矩阵式
楼主: MrSomebody (MrSomebody)   2016-06-30 07:11:00
谢谢Q大回应,但constrOptim函数已设定不等式为大于等于,但sum条件为等于,解出来得alpha无法满足,悲剧...
作者: qhair (卡早睡卡有眠)   2016-06-30 12:19:00
你的问题用kernlab套件的ipop或quadprog套件应该比较合适
楼主: MrSomebody (MrSomebody)   2016-06-30 19:25:00
谢谢Q大! 我用ipop可以解这个问题!

Links booklink

Contact Us: admin [ a t ] ucptt.com