R中的任何功能与SAS中的PROC OPTMODEL相同吗？

Question

您知道是否有任何R包或功能可以帮助解决      ＆＃34; min y = x ^ 3 - 3 * x ^ 2 + 4＆＃34; 像下面的SAS代码一样？

proc optmodel;      
   var x;
   min y=x**3 - 3*x**2 + 4;
 solve;
 print x;
 quit;

PROC OPTMODEL在此处记录： http://support.sas.com/documentation/cdl/en/ormpug/67517/HTML/default/viewer.htm#ormpug_optmodel_toc.htm

我正在寻找最小化y的x的值。 提前谢谢！

Answer 1

您正试图找到多项式的根。执行此操作的R函数是polyroot()。 polyroot()的参数是多项式系数的向量，按递增顺序排列。试试这个：

polyroot(c(4, 0, -3, 1))
[1] -1+0i  2-0i  2+0i

现在，您可以通过设置curve()的绘图范围，使用此信息可视化多项式：

curve(x^3 - 3*x^2 + 4, from=-2, to=3, col="blue")
abline(h=0)