r语言optim函数给不同的初值为什么结果不同

发布网友发布时间：2022-04-22 14:18

共2个回答

热心网友时间：2022-04-15 22:14

用内置函数optim()optim(par,fun,lower,upper,method)大致用到这5个参数par是初始值，你选离你峰值差不远的xfun是生成你正弦波的函数lower和upper定义域method用"BFGS"牛顿迭代法，或者"L-BFGS-B"升级版牛顿迭代法。以下是得到的结果，我用f(x)=x^2-2*x+1试了以下>optim(3,fun,lower=-5,upper=5,method="BFGS")$par[1]1#x值$value[1]0#y值$countsfunctiongradient#不知道是什么4444$convergence[1]52#算了多少次收敛

热心网友时间：2022-04-15 23:32

-res$value

>##初始似然值

>L0<-sum(s1)*log(0.5)+sum(s2)*log(0.5)

>##拟合度计算

>##结果输出

>##p^2的值

>cat("roh = ",(L0-LL)/L0)

>##修正的p^2值

>cat("rohbar=",(L0-(LL-length(b)))/L0)

>print(res)

>print(tval)
-res$value

>##初始似然值

>L0<-sum(s1)*log(0.5)+sum(s2)*log(0.5)

>##拟合度计算

>##结果输出

>##p^2的值

>cat("roh = ",(L0-LL)/L0)

>##修正的p^2值

>cat("rohbar=",(L0-(LL-length(b)))/L0)

>print(res)

>print(tval)
-res$value

>##初始似然值

>L0<-sum(s1)*log(0.5)+sum(s2)*log(0.5)

>##拟合度计算

>##结果输出

>##p^2的值

>cat("roh = ",(L0-LL)/L0)

>##修正的p^2值

>cat("rohbar=",(L0-(LL-length(b)))/L0)

>print(res)

>print(tval)